The Para rubber tree (Hevea brasiliensis) is an economically important tropical tree species that produces natural rubber, an essential industrial raw material. Here we present a high-quality genome assembly of this species (1.37 Gb, scaffold N50 = 1.28 Mb) that covers 93.8% of the genome (1.47 Gb) and harbours 43,792 predicted protein-coding genes. A striking expansion of the REF/SRPP (rubber elongation factor/small rubber particle protein) gene family and its divergence into several laticifer-specific isoforms seem crucial for rubber biosynthesis. The REF/SRPP family has isoforms with sizes similar to or larger than SRPP1 (204 amino acids) in 17 other plants examined, but no isoforms with similar sizes to REF1 (138 amino acids), the predominant molecular variant. A pivotal point in Hevea evolution was the emergence of REF1, which is located on the surface of large rubber particles that account for 93% of rubber in the latex (despite constituting only 6% of total rubber particles, large and small). The stringent control of ethylene synthesis under active ethylene signalling and response in laticifers resolves a longstanding mystery of ethylene stimulation in rubber production. Our study, which includes the re-sequencing of five other Hevea cultivars and extensive RNA-seq data, provides a valuable resource for functional genomics and tools for breeding elite Hevea cultivars.

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.

A homolog of MOTHER OF FT AND TFL1 (MFT) was isolated from Hevea brasiliensis and its biological function was investigated. Protein multiple sequence alignment and phylogenetic analysis revealed that HbMFT1 conserved critical amino acid residues to distinguish MFT, FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1)-like proteins and showed a closer genetic relationship to the MFT-like group. The accumulation of HbMFT1 was generally detected in various tissues except pericarps, with the highest expression in embryos and relatively higher expression in roots and stems of seedlings, flowering inflorescences, and male and female flowers. HbMFT1 putative promoter analysis showed that tissue-specific, environmental change responsive and hormone-signaling responsive elements were generally present. HbMFT1 was strongly induced under a short-day condition at 28 °C, with the highest expression after the onset of a day. Overexpression of HbMFT1 inhibited seed germination, seedling growth, and flowering in transgenic Arabidopsis. The qRT-PCR further confirmed that APETALA1 (AP1) and FRUITFULL (FUL) were drastically down-regulated in 35S::HbMFT1 plants. A histochemical β-glucuronidase (GUS) assay showed that HbMFT1::GUS activity was mainly detected in stamens and mature seeds coinciding with its original expression and notably induced in rosette leaves and seedlings of transgenic Arabidopsis by exogenous abscisic acid (ABA) due to the presence of ABA cis-elements in HbMFT1 promoter. These results suggested that HbMFT1 was mainly involved in maintenance of seed maturation and stamen development, but negatively controlled germination, growth and development of seedlings and flowering. In addition, the HbMFT1 promoter can be utilized in controlling transgene expression in stamens and seeds of rubber tree or other plant species.

The recurring growth of bacterium in newly developed resistant cells and a minimal level of bacterial infection rate are the main limiting factors of Agrobacterium-mediated transformation experiments in Hevea brasiliensis. The current study aimed to optimize crucial factors of the transformation protocol in order to obtain an efficient transformation experimental model for Hevea using cotyledonary somatic embryos as explants. Transformation conditions such as antibiotic concentration, preculture duration, Agrobacterium concentration, sonication and cocultivation conditions were analyzed using the binary vector pCAMBIA2301. Transient transformation was confirmed by GUS histochemical staining. The best transformation efficiency was observed when the explants were not cultured on a preculture medium that contained acetosyringone at a level of 100 μM. The best results were obtained using a bacterial density of 0.45 at OD 600 nm, 50 s of sonication of explants in a bacterial liquid culture and a total incubation time of 18 min in the same bacterial suspension. Transmission electron microscopical analysis confirmed the impacts of sonication on bacterial infection efficiency. Cocultivation conditions of 22 °C and 84 h of darkness were optimal for the transfer of T-DNA. Agrobacterium was eliminated with 500 mg/L of timentin, and the selection of transformants was performed using 100 mg/L of kanamycin in the selection medium. The presence of transgene was confirmed in the resistant embryos by polymerase chain reaction (PCR). The improved method of genetic transformation established in the present study will be useful for the introduction of foreign genes of interest into the Hevea genome for the breeding of this economically important plant species in the future.

Background: The processabilities and mechanical properties of natural rubber depend greatly on its molecular weight (MW) and molecular weight distribution (MWD). However, the mechanisms underlying the regulation of molecular weight during rubber biosynthesis remain unclear. Results: In the present study, we determined the MW and particle size of latex from 1-year-old virgin trees and 30-year-old regularly tapped trees of the Hevea clones Reyan7-33–97 and RRIM600. The results showed that both the MW and the particle size of latex varied between these two clones and increased with tree age. Latex from RRIM600 trees had a smaller average particle size than that from Reyan7-33–97 trees of the same age. In 1-year-old trees, the Reyan7-33–97 latex displayed a slightly higher MW than that of RRIM600, whereas in 30-year-old trees, the RRIM600 latex had a significantly higher MW than the Reyan7-33–97 latex. Comparative analysis of the transcriptome profiles indicated that the average rubber particle size is negatively correlated with the expression levels of rubber particle associated proteins, and that the high-MW traits of latex are closely correlated with the enhanced expression of isopentenyl pyrophosphate (IPP) monomer-generating pathway genes and downstream allylic diphosphate (APP) initiator-consuming non-rubber pathways. By bioinformatics analysis, we further identified a group of transcription factors that potentially regulate the biosynthesis of IPP. Conclusions: Altogether, our results revealed the potential regulatory mechanisms involving gene expression variations in IPP-generating pathways and the non-rubber isoprenoid pathways, which affect the ratios and contents of IPP and APP initiators, resulting in significant rubber MW variations among same-aged trees of the Hevea clones Reyan7-33–97 and RRIM600. Our findings provide a better understanding of rubber biosynthesis and lay the foundation for genetic improvement of rubber quality in H. brasiliensis.