Chaorong Tang

The rubber tree genome reveals new insights into rubber production and species adaptation.

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.

Structure and expression profile of the sucrose synthase gene family in the rubber tree: indicative of roles in stress response and sucrose utilization in the laticifers

Sucrose synthase (Sus, EC 2.4.1.13) is widely recognized as a key enzyme in sucrose metabolism in plants. However, nothing is known about this gene family in Hevea brasiliensis (para rubber tree). Here, we identified six Sus genes in H. brasiliensis that comprise the entire Sus family in this species. Analysis of the gene structure and phylogeny of the Sus genes demonstrates evolutionary conservation in the Sus families across Hevea and other plant species. The expression of Sus genes was investigated via Solexa sequencing and quantitative PCR in various tissues, at various phases of leaf development, and under abiotic stresses and ethylene treatment. The Sus genes exhibited distinct but partially redundant expression profiles. Each tissue has one abundant Sus isoform, with HbSus3, 4 and 5 being the predominant isoforms in latex (cytoplasm of rubber‐producing laticifers), bark and root, respectively. HbSus1 and 6 were barely expressed in any tissue examined. In mature leaves (source), all HbSus genes were expressed at low levels, but HbSus3 and 4 were abundantly expressed in immature leaves (sink). Low temperature and drought treatments conspicuously induced HbSus5 expression in root and leaf, suggesting a role in stress responses. HbSus2 and 3 transcripts were decreased by ethylene treatment, consistent with the reduced sucrose‐synthesizing activity of Sus enzymes in the latex in response to ethylene stimulation. Our results are beneficial to further determination of functions for the Sus genes in Hevea trees, especially roles in regulating latex regeneration.

Comparative analysis of latex transcriptome reveals putative molecular mechanisms underlying super productivity of Hevea brasiliensis.

Increasing demand for natural rubber prompts studies into the mechanisms governing the productivity of rubber tree ( Hevea brasiliensis ). It is very interesting to notice that a rubber tree of clone PR107 in Yunnan, China is reported to yield more than 20 times higher than the average rubber tree. This super-high-yielding (SHY) rubber tree (designated as SY107), produced 4.12 kg of latex (cytoplasm of rubber producing laticifers, containing about 30% of rubber) per tapping, more than 7-fold higher than that of the control. This rubber tree is therefore a good material to study how the rubber production is regulated at a molecular aspect. A comprehensive cDNA-AFLP transcript profiling was performed on the latex of SY107 and its average counterparts by using the 384 selective primer pairs for two restriction enzyme combinations (ApoI/MseI and TaqI/MseI). A total of 746 differentially expressed (DE) transcript-derived fragments (TDFs) were identified, of which the expression patterns of 453 TDFs were further confirmed by RT-PCR. These RT-PCR confirmed TDFs represented 352 non-redundant genes, of which 215 had known or partially known functions and were grouped into 10 functional categories. The top three largest categories were transcription and protein synthesis (representing 24.7% of the total genes), defense and stress (15.3%), and primary and secondary metabolism (14.0%). Detailed analysis of the DE-genes suggests notable characteristics of SHY phenotype in improved sucrose loading capability, rubber biosynthesis-preferred sugar utilization, enhanced general metabolism and timely stress alleviation. However, the SHY phenotype has little correlation with rubber-biosynthesis pathway genes.

Transcriptional and post-transcriptional regulation of the jasmonate signalling pathway in response to abiotic and harvesting stress in Hevea brasiliensis

BackgroundLatex harvesting in Hevea brasiliensis amounts to strong abiotic stress that can cause a halt in production in the most susceptible clones. Although the role of jasmonic acid has been suggested in laticifer differentiation, its role in latex production and in the response to harvesting stress has received very little attention. Only a few key genes acting in the COI-JAZ-MYC module have been isolated and studied at transcriptional level.ResultsUse of a reference transcriptome obtained on rubber clone PB 260 covering a large number of tissues under different environmental conditions enabled us to identify 24 contigs implicated in the jasmonate signalling pathway in the rubber tree. An analysis of their expression profile by qPCR, combined with hierarchical clustering, suggested that the jasmonate signalling pathway is highly activated in laticifer cells and, more particularly, in the response to harvesting stress. By comparison with their genomic sequences, the existence of regulation by alternative splicing was discovered for JAZ transcripts in response to harvesting stress. Lastly, positive transcriptional regulation of the HbJAZ_1405 gene by MYC was demonstrated.ConclusionThis study led to the identification of all actors of jasmonate signalling pathway and revealed a specific gene expression pattern in latex cells. In-depth analysis of this regulation showed alternative splicing that has been previously shown in Arabidopsis. Interestingly, genotypic variation was observed in Hevea clones with contrasting latex metabolism. This result suggests an involvement of jasmonate signalling pathway in latex production. The data suggest that specific variability of the JA pathway may have some major consequences for resistance to stress. The data support the hypothesis that a better understanding of transcriptional regulations of jasmonate pathway during harvesting stress, along with the use of genotypic diversity in response to such stress, can be used to improve resistance to stress and rubber production in Hevea.

The small RNA profile in latex from Hevea brasiliensis trees is affected by tapping panel dryness.

Natural rubber is harvested by tapping Hevea brasiliensis (Willd. ex A. Juss.) Müll. Arg. Harvesting stress can lead to tapping panel dryness (TPD). MicroRNAs (miRNAs) are induced by abiotic stress and regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs. This study set out to sequence miRNAs expressed in latex cells and to identify TPD-related putative targets. Deep sequencing of small RNAs was carried out on latex from trees affected by TPD using Solexa technology. The most abundant small RNA class size was 21 nucleotides for TPD trees compared with 24 nucleotides in healthy trees. By combining the LeARN pipeline, data from the Plant MicroRNA database and Hevea EST sequences, we identified 19 additional conserved and four putative species-specific miRNA families not found in previous studies on rubber. The relative transcript abundance of the Hbpre-MIR159b gene increased with TPD. This study revealed a small RNA-specific signature of TPD-affected trees. Both RNA degradation and a shift in miRNA biogenesis are suggested to explain the general decline in small RNAs and, particularly, in miRNAs.

Regulation of MIR Genes in Response to Abiotic Stress in Hevea brasiliensis

Increasing demand for natural rubber (NR) calls for an increase in latex yield and also an extension of rubber plantations in marginal zones. Both harvesting and abiotic stresses lead to tapping panel dryness through the production of reactive oxygen species. Many microRNAs regulated during abiotic stress modulate growth and development. The objective of this paper was to study the regulation of microRNAs in response to different types of abiotic stress and hormone treatments in Hevea. Regulation of MIR genes differs depending on the tissue and abiotic stress applied. A negative co-regulation between HbMIR398b with its chloroplastic HbCuZnSOD target messenger is observed in response to salinity. The involvement of MIR gene regulation during latex harvesting and tapping panel dryness (TPD) occurrence is further discussed.

Ethylene response factors are controlled by multiple harvesting stresses in hevea brasiliensis

Tolerance of recurrent mechanical wounding and exogenous ethylene is a feature of the rubber tree. Latex harvesting involves tapping of the tree bark and ethephon is applied to increase latex flow. Ethylene is an essential element in controlling latex production. The ethylene signalling pathway leads to the activation of Ethylene Response Factor (ERF) transcription factors. This family has been identified in Hevea brasiliensis. This study set out to understand the regulation of ERF genes during latex harvesting in relation to abiotic stress and hormonal treatments. Analyses of the relative transcript abundance were carried out for 35 HbERF genes in latex, in bark from mature trees and in leaves from juvenile plants under multiple abiotic stresses. Twenty-one HbERF genes were regulated by harvesting stress in laticifers, revealing an overrepresentation of genes in group IX. Transcripts of three HbERF-IX genes from HbERF-IXc4, HbERF-IXc5 and HbERF-IXc6 were dramatically accumulated by combining wounding, methyl jasmonate and ethylene treatments. When an ethylene inhibitor was used, the transcript accumulation for these three genes was halted, showing ethylene-dependent induction. Subcellular localization and transactivation experiments confirmed that several members of HbERF-IX are activator-type transcription factors. This study suggested that latex harvesting induces mechanisms developed for the response to abiotic stress. These mechanisms probably depend on various hormonal signalling pathways. Several members of HbERF-IX could be essential integrators of complex hormonal signalling pathways in Hevea.

De novo Transcriptome Analysis Reveals Distinct Defense Mechanisms by Young and Mature Leaves of Hevea brasiliensis (Para Rubber Tree)

Along with changes in morphology in the course of maturation, leaves of Hevea brasiliensis become more resistant to leaf diseases, including the South American Leaf Blight (SALB), a devastating fungal disease of this economically important tree species. To understand the underlying mechanisms of this defense, and to identify the candidate genes involved, we sequenced the Hevea leaf transcriptome at four developmental stages (I to IV) by Illumina sequencing. A total of 62.6 million high-quality reads were generated, and assembled into 98,796 unique transcripts. We identified 3,905 differentially expressed genes implicated in leaf development, 67.8% (2,651) of which were during the transition to leaf maturation. The genes involved in cyanogenic metabolism, lignin and anthocyanin biosynthesis were noteworthy for their distinct patterns of expression between developing leaves (stages I to III) and mature leaves (stage IV), and the correlation with the change in resistance to SALB and the Oidium/Colletotrichum leaf fall. The results provide a first profile of the molecular events that relate to the dynamics of leaf morphology and defense strategies during Hevea leaf development. This dataset is beneficial to devising strategies to engineer resistance to leaf diseases as well as other in-depth studies in Hevea tree.

Molecular identification and characterization of the pyruvate decarboxylase gene family associated with latex regeneration and stress response in rubber tree.

In plants, ethanolic fermentation occurs not only under anaerobic conditions but also under aerobic conditions, and involves carbohydrate and energy metabolism. Pyruvate decarboxylase (PDC) is the first and the key enzyme of ethanolic fermentation, which branches off the main glycolytic pathway at pyruvate. Here, four PDC genes were isolated and identified in a rubber tree, and the protein sequences they encode are very similar. The expression patterns of HbPDC4 correlated well with tapping-simulated rubber productivity in virgin rubber trees, indicating it plays an important role in regulating glycometabolism during latex regeneration. HbPDC1, HbPDC2 and HbPDC3 had striking expressional responses in leaves and bark to drought, low temperature and high temperature stresses, indicating that the HbPDC genes are involve in self-protection and defense in response to various abiotic and biotic stresses during rubber tree growth and development. To understand ethanolic fermentation in rubber trees, it will be necessary to perform an in-depth study of the regulatory pathways controlling the HbPDCs in the future.

Validation of reference genes for quantitative real-time PCR during latex regeneration in rubber tree.

In rubber tree, latex regeneration is one of the decisive factors influencing the rubber yield, although its molecular regulation is not well known. Quantitative real-time PCR (qPCR) is a popular and powerful tool used to understand the molecular mechanisms of latex regeneration. However, the suitable reference genes required for qPCR are not available to investigate the expressions of target genes during latex regeneration. In this study, 20 candidate reference genes were selected and evaluated for their expression stability across the samples during the process of latex regeneration. All reference genes showed a relatively wide range of the threshold cycle values, and their stability was validated by four different algorithms (comparative delta Ct method, Bestkeeper, NormFinder and GeNorm). Three softwares (comparative delta Ct method, NormFinder and GeNorm) exported similar results that identify UBC4, ADF, UBC2a, eIF2 and ADF4 as the top five suitable references, and 18S as the least suitable one. The application of the screened references would improve accuracy and reliability of gene expression analysis in latex regeneration experiments.