Dejun Li

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.

Transcriptome analyses reveal molecular mechanism underlying tapping panel dryness of rubber tree ( Hevea brasiliensis )

Tapping panel dryness (TPD) is a serious threat to natural rubber yields from rubber trees, but the molecular mechanisms underlying TPD remain poorly understood. To identify TPD-related genes and reveal these molecular mechanisms, we sequenced and compared the transcriptomes of bark between healthy and TPD trees. In total, 57,760 assembled genes were obtained and analyzed in details. In contrast to healthy rubber trees, 5652 and 2485 genes were up- or downregulated, respectively, in TPD trees. The TPD-related genes were significantly enriched in eight GO terms and five KEGG pathways and were closely associated with ROS metabolism, programmed cell death and rubber biosynthesis. Our results suggest that rubber tree TPD is a complex process involving many genes. The observed lower rubber yield from TPD trees might result from lower isopentenyl diphosphate (IPP) available for rubber biosynthesis and from downregulation of the genes in post-IPP steps of rubber biosynthesis pathway. Our results not only extend our understanding of the complex molecular events involved in TPD but also will be useful for developing effective measures to control TPD of rubber trees.

Gene expression analysis and SNP/InDel discovery to investigate yield heterosis of two rubber tree F1 hybrids.

As an important industrial material, natural rubber is mainly harvested from the rubber tree. Rubber tree breeding is inefficient, expensive and time-consuming, whereas marker-assisted selection is a feasible method for early selection of high-yield hybrids. We thus sequenced and analyzed the transcriptomes of two parent rubber trees (RRIM 600 and PR 107) and their most productive hybrids (RY 7-33-97 and RY 7-20-59) to understand their gene expression patterns and genetic variations including single nucleotide polymorphisms (SNPs) and small insertions/deletions (InDels). We discovered >31,000 genetic variations in 112,702 assembled unigenes. Our results showed that the higher yield in F1 hybrids was positively associated with their higher genome heterozygosity, which was further confirmed by genotyping 10 SNPs in 20 other varieties. We also showed that RY 7-33-97 and RY 7-20-59 were genetically closer to RRIM 600 and PR 107, respectively, in agreement with both their phenotypic similarities and gene expression profiles. After identifying ethylene- and jasmonic acid–responsive genes at the transcription level, we compared and analyzed the genetic variations underlying rubber biosynthesis and the jasmonic acid and ethylene pathways in detail. Our results suggest that genome-wide genetic variations play a substantive role in maintaining rubber tree heterosis.

Molecular cloning, expression profiles, and characterization of a novel polyphenol oxidase (PPO) gene in Hevea brasiliensis

The polyphenol oxidase (PPO) is involved in undesirable browning in many plant foods. Although the PPOs have been studied by several researchers, the isolation and expression profiles of PPO gene were not reported in rubber tree. In this study, a new PPO gene, HbPPO, was isolated from Hevea brasiliensis. The sequence alignment showed that HbPPO indicated high identities to plant PPOs and belonged to dicot branch. The cis-acting regulatory elements related to stress/hormone responses were predicted in the promoter region of HbPPO. Real-time RT-PCR analyses showed that HbPPO expression varied widely depending on different tissues and developmental stages of leaves. Besides being associated with tapping panel dryness, the HbPPO transcripts were regulated by ethrel, wounding, H2O2, and methyl jasmonate treatments. Moreover, the correlation between latex coagulation rate and PPO activity was further confirmed in this study. Our results lay the foundation for further analyzing the function of HbPPO in rubber tree. Graphical Abstract The WRKY71OS and transcription initiation site were highlighted with square frame and bold letter, respectively. The other cis-acting regulatory elements were underlined.

Expression Profiles, Characterization and Function of HbTCTP in Rubber Tree (Hevea brasiliensis)

As a highly conserved protein, the translationally controlled tumor protein (TCTP) carries out vital roles in various life processes. In rubber tree, two TCTP genes, HbTCTP and HbTCTP1, were cloned, but only HbTCTP1 was studied in details. In this study, cis-acting regulatory elements, expression patterns, subcellular localization, interacting proteins, and antioxidant activity of HbTCTP were systematically analyzed. Besides the common cis-acting regulatory elements, HbTCTP promoter also harbored various known cis-elements that respond to hormone/stresses. Being consistent with the aforementioned results, HbTCTP was regulated by drought, low temperature, high salt, ethylene (ET), wounding, H2O2, and methyl jasmonate (MeJA) treatments. HbTCTP was expressed throughout different tissues and developmental stages of leaves. In addition, HbTCTP was associated with tapping panel dryness (TPD). HbTCTP was localized in the membrane, cytoplasm and the nucleus, and interacted with four proteins rubber elongation factor (REF), 17.5 kDa heat shock family protein, annexin, and REF-like stress related protein 1. Being similar to HbTCTP1, HbTCTP also indicated antioxidant activity in metal-catalyzed oxidation (MCO) system. Our results are useful for further understanding the molecular characterization and expression profiles of HbTCTP, but also lay a solid foundation for elucidating the function of HbTCTP in rubber tree.

Next-generation sequencing, assembly, and comparative analyses of the latex transcriptomes from two elite Hevea brasiliensis varieties

The great progress has been made in rubber tree breeding, but the molecular mechanisms underlying high yield are not well understood. Here, we reported the sequencing, assembly, and comparative analyses of latex transcriptome from two rubber tree varieties. In total, 33,852 unigenes were generated with de novo assembly. The blastx results indicated that 27,886 and 15,704 unigenes showed significant similarities to known proteins from NCBI nr and Swissprot databases, respectively. Among these annotated unigenes, 21,841 and 9010 ones were separately assigned to Gene Ontology (GO) functional categories and Clusters of Orthologous Groups (COGs). Of 126 KEGG pathways, metabolic pathway was the biggest one, suggesting that active metabolic processes happen in rubber tree latex. In contrast to RRIM 600, 2513 and 1391 genes were separately up- and downregulated in RY 7-20-59. The expression profiles of 25 unigenes were further confirmed by real-time RT-PCR, suggesting that the differently expressed genes (DEGs) identified by RNA-seq were accurate and reliable in this study. The DEGs between RRIM 600 and RY 7-20-59 were significantly enriched in plant-pathogen interactions, phenylpropanoid biosynthesis, phenylalanine metabolism, ubiquinone and other terpenoid-quinone biosynthesis, biosynthesis of secondary metabolites, and photosynthesis. Interestingly, the genes involved in rubber biosynthesis pathway, such as CPT, GPPS, HMGR, HMGS, FPPS and DXS, were differently expressed between RRIM 600 and RY 7-20-59. It was the first time that the latex transcriptomes of two rubber tree varieties have been compared and analyzed on a transcriptome-wide scale. Our results not only enrich the transcriptome data of rubber tree but also provide new insights into understanding latex transcriptome and molecular mechanisms underlying high yielding in rubber tree.

Genome-wide identification and expression analysis of the metacaspase gene family in Hevea brasiliensis

Metacaspases, a family of cysteine proteases, have been suggested to play important roles in programmed cell death (PCD) during plant development and stress responses. To date, no systematic characterization of this gene family has been reported in rubber tree (Hevea brasiliensis). In the present study, nine metacaspase genes, designated as HbMC1 to HbMC9, were identified from whole-genome sequence of rubber tree. Multiple sequence alignment and phylogenetic analyses suggested that these genes were divided into two types: type I (HbMC1-HBMC7) and type II (HbMC8 and HbMC9). Gene structure analysis demonstrated that type I and type II HbMCs separately contained four and two introns, indicating the conserved exon-intron organization of HbMCs. Quantitative real-time PCR analysis revealed that HbMCs showed distinct expression patterns in different tissues, suggesting the functional diversity of HbMCs in various tissues during development. Most of the HbMCs were regulated by drought, cold, and salt stress, implying their possible functions in regulating abiotic stress-induced cell death. Of the nine HbMCs, HbMC1, HbMC2, HbMC5, and HbMC8 displayed a significantly higher relative transcript accumulation in barks of tapping panel dryness (TPD) trees compared with healthy trees. In addition, the four genes were up-regulated by ethephon (ET) and methyl jasmonate (MeJA), indicating their potential involvement in TPD resulting from ET- or JA-induced PCD. In summary, this work provides valuable information for further functional characterization of HbMC genes in rubber tree.

Molecular cloning, expression profiles and characterization of a glutathione reductase in Hevea brasiliensis

Glutathione reductase (GR; EC 1.8.1.7) is an important oxidoreductase that can protect organisms against various oxidative stresses. In this study, a new GR gene, named as HbGR2, was isolated from Hevea brasiliensis. The HbGR2 cDNA contained a 1674-bp open reading frame encoding 557 amino acids and the deduced HbGR2 protein showed high identities to the chloroplastic GRs from other plant species. HbGR2 was localized in the chloroplasts of tobacco mesophyll protoplasts. The cis-acting regulatory elements related to stress or hormone responses were predicted in the promoter region of HbGR2. The results from real-time RT-PCR analyses showed that HbGR2 was expressed throughout different tissues and developmental stages of leaves. Besides being related to tapping panel dryness (TPD), HbGR2 was regulated by several treatments including ethephon (ET), methyl jasmonate (MeJA), drought, low temperature, high salt, wounding and hydrogen peroxide (H2O2). The Escherichia coli (E. coli) cells overexpressing HbGR2 markedly increased their tolerance and survival at high concentrations of H2O2, suggesting that HbGR2 might play an important role in oxidative stress response in Hevea brasiliensis.

Transcriptomic analyses reveal molecular mechanisms underlying growth heterosis and weakness of rubber tree seedlings

BackgroundBreeding rubber tree seedling with growth heterosis is vital for natural rubber production. It is the prerequisites for effectively utilizing growth heterosis to elucidate its molecular mechanisms, but the molecular mechanisms remain poorly understood in rubber tree. To elucidate seedling growth heterosis, we conducted comparative transcriptomic analyses between the two hybrids and their parents.ResultsBy identifying and comparing differently expressed genes (DEGs), we found that the hybrids (BT 3410 and WC 11) show significantly differential expression profiles from their parents (PR 107 and RRIM 600). In BT 3410-parent triad, 1092 (49.95%) and 1094 (50.05%) DEGs indicated clear underdominance or overdominance, respectively. Whereas in WC 11-parent triad, most DEGs (78.2%, 721) showed low- or high-parent dominance; 160 (17.35%) exhibited expression patterns that are not statistically distinguishable from additivity, and 8 (0.87%) and 33 (3.58%) DEGs exhibited underdominance and overdominance, respectively. Furthermore, some biological processes are differentially regulated between two hybrids. Interestingly, the pathway in response to stimulus is significantly downregulated and metabolic pathways are more highly regulated in BT 3410.ConclusionsTaken together, the genotypes, transcriptomes and biological pathways (especially, carbohydrate metabolism) are highly divergent between two hybrids, which may be associated with growth heterosis and weakness. Analyzing gene action models in hybrid-parent triads, we propose that overdominance may play important roles on growth heterosis, whereas dominance on hybrid weakness in rubber tree seedlings. These findings bring new insights into our understanding of growth heterosis of rubber tree seedling.