Summary: The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high‐quality, chromosome‐level genome sequence of the wild rubber tree, the first report on selection signatures and a genome‐wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15‐year field trial, GWAS identified 155 marker–trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic‐assisted breeding by offering valuable resources in this new domesticated crop.

The rubber tree (Hevea brasiliensis Muell. Arg) is a tropical, perennial, woody plant that is susceptible to cold stress. In China, cold stress has been found to severely damage rubber plants in plantations in past decades. Although several Hevea clones that are resistant to cold have been developed, their cold hardiness mechanism has yet to be elucidated. For the study reported herein, we subjected the cold-resistant clone CATAS93-114 and the cold-sensitive clone Reken501 to chilling stress, and characterized their transcriptomes at 0, 2, 8 and 24 h after the start of chilling. We found that 7870 genes were differentially expressed in the transcriptomes of the two clones. In CATAS93-114, a greater number of genes were found to be up- or downregulated between 2 h and 8 h than in Reken501, which indicated a more rapid and intensive response by CATAS93-114 than by Reken501. The differentially expressed genes were grouped into seven major clusters, according to their Gene Ontology terms. The expression profiles for genes involved in abscisic acid metabolism and signaling, in an abscisic acid-independent pathway, and in early signal perception were found to have distinct expression patterns for the transcriptomes of the two clones. The differential expression of 22 genes that appeared to have central roles in response to chilling was confirmed by quantitative real-time PCR.

Rubber trees (Hevea brasiliensis) are susceptible to low temperature and therefore are only planted in the tropical regions. In the past few decades, although rubber trees have been successfully planted in the northern margin of tropical area in China, they suffered from cold injury during the winter. To understand the physiological response under cold stress, we isolated a C-repeat binding factor 1 (CBF1) gene from the rubber tree. This gene (HbCBF1) was found to respond to cold stress but not drought or ABA stress. The corresponding HbCBF1 protein showed CRT/DRE binding activity in gel shift experiment. To further characterize its molecular function, the HbCBF1 gene was overexpressed in Arabidopsis. The HbCBF1 over expression (OE) line showed enhanced cold resistance and relatively slow dehydration, and the expression of Arabidopsis CBF pathway downstream target genes, e.g. AtCOR15a and AtRD29a, were significantly activated under non-acclimation condition. These data suggest HbCBF1 gene is a functional member of the CBF gene family, and may play important regulation function in rubber tree.

Rubber tree (Hevea brasiliensis) is the primary source of natural rubber and economically important. We present the haplotype-resolved, telomere-to-telomere, gap-free genome assembly of the cultivar CATAS 7-33–97, with both haplotypes containing complete telomeric and centromeric regions. Structural variations, including a 32.71 Mb inversion on chromosome 8, are identified. The fully assembled 36 chromosomes enable comprehensive identification of rubber biosynthesis genes and their allele-specific expression. By integrating transcriptomic and metabolomic data, we reconstruct the rubber biosynthesis pathway and confirm the mevalonate (MVA) pathway as the major carbon source for rapid latex regeneration during tapping. Jasmonic acid (JA) plays a key role in promoting rubber yield by enhancing biosynthetic activity in response to mechanical wounding. We propose a model where JA-induced myelocytomatosis proteins 2 activate mevalonate kinase 1 expression, boosting MVA synthesis and rubber production. These findings provide insights into rubber tree genomics and its molecular response to tapping.

Rubber trees are economically important tropical tree species and the major source of natural rubber, which is an essential industrial material. This tropical perennial tree is susceptible to cold stress and other abiotic stresses, especially in the marginal northern tropics. Recent years, the genome sequencing and RNA-seq projects produced huge amount of sequence data, which greatly facilitated the functional genomics study. However, the characterization of individual functional gene is in urgent demands, especially for those involved in stress resistance. Here we identified and characterized the rubber tree gene ErbB-3 binding protein 1, which undergoes changes in expression in response to cold, drought stress and ABA treatment. HbEBP1 overexpression in Arabidopsis increased organ size, facilitated root growth and increased adult leaf number by delaying the vegetative-to-reproductive transition. In addition, HbEBP1 overexpression enhanced the resistance of the Arabidopsis plants to freezing and drought stress, demonstrating that this gene participates in the regulation of abiotic stress resistance. RD29a, RD22 and CYCD3;1 expression was also greatly enhanced by HbEBP1 overexpression, which explains its regulatory roles in organ size and stress resistance. The regulation of drought stress resistance is a novel function identified in plant EBP1 genes, which expands our understanding of the roles of EBP1 gene in response to the environment. Our results provide information that may lead to the use of HbEBP1 in genetically engineered crops to increase both biomass and abiotic stress resistance.