Comparative Transcriptome Analysis of Latex Reveals Molecular Mechanisms Underlying Increased Rubber Yield in Hevea brasiliensis Self-Rooting Juvenile Clones

Reprogramming of the Hevea brasiliensis Epigenome and Transcriptome in Response to Cold Stress

Low temperature is a key factor limiting the rubber plantation extending to high latitude area. Previous work has shown that cold-induced DNA demethylation was coordinated with the expression of cold-responsive (COR) genes in Hevea brasiliensis. In this work, reduced representation bisulphite sequencing analysis of H. brasiliensis showed that cold treatment induced global genomic DNA demethylation and altered the sequence contexts of methylated cytosines, but the levels of mCG methylation in transposable elements were slightly enhanced by cold treatment. Integrated analysis of the DNA methylome and transcriptome revealed 400 genes whose expression correlated with altered DNA methylation. DNA demethylation in the upstream region of gene seems to correlate with higher gene expression, whereas demethylation in the gene body has less association. Our results suggest that cold treatment globally change the genomic DNA methylation status of the rubber tree, which might coordinate reprogramming of the transcriptome.

Differential Expression of lncRNAs and miRNAs Between Self-Rooting Juvenile and Donor Clones Unveils Novel Insight Into the Molecular Regulation of Rubber Biosynthesis in Hevea brasiliensis

The rubber tree (Hevea brasiliensis Muell. Arg.) is a tropical tree species that produce natural rubber. Self-rooted juvenile clones (SRJCs) are novel rubber tree planting materials developed through primary somatic embryogenesis. SRJCs have a higher rubber yield compared with donor clones (DCs). The molecular basis underlying increased rubber yield in SRJCs remains largely unknown. Here, the latex from SRJCs and DCs were collected for strand-specific and small RNA-seq methods. A total of 196 differentially expressed long noncoding RNAs (DELs), and 11 differentially expressed microRNAs were identified in latex between SRJCs and DCs. Targeted genes of DELs were markedly enriched for various biological pathways related to plant hormone signal transduction, photosynthesis, glutathione metabolism, and amino acids biosynthesis. DELs probably acted as cis-acting regulation was calculated, and these DELs relevant to potentially regulate rubber biosynthesis, reactive oxygen species metabolism, and epigenetic modification. Furthermore, the DELs acting as microRNA targets were studied. The interaction of microRNA and DELs might involve in the regulation of natural rubber biosynthesis.

Identification and characterization of jasmonic acid- and linolenic acid-mediated transcriptional regulation of secondary laticifer differentiation in Hevea brasiliensis

Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.

Identification and characterization of the MADS-box genes highly expressed in the laticifer cells of Hevea brasiliensis

MADS-box transcription factors possess many functions in plant reproduction and development. However, few MADS-box genes related to secondary metabolites regulation have been identified. In Hevea brasiliensis, natural rubber is a representative cis-polyisoprenoids in secondary metabolism which occurs in the rubber laticifer cells, the molecular regulation basis of natural rubber biosynthesis is not clear. Here, a total of 24 MADS-box genes including 4 type I MADS-box genes and 20 type II MADS-box genes were identified in the transcriptome of rubber tree latex. The phylogenetic analysis was performed to clarify the evolutionary relationships of all the 24 rubber tree MADS-box proteins with MADS-box transcription factors from Arabidopsis thaliana and Oryza sativa. Four type I MADS-box genes were subdivided into Mα (3 genes) and Mβ (1 gene). Twenty type II MADS-box genes were subclassified into MIKC* (8 genes) and MIKC^c (12 genes). Eight MADS-box genes (HblMADS3, 5, 6, 7, 9, 13, 23, 24) were predominant expression in laticifers. ABA up-regulated the expression of HblMADS9, and the expression of HblMADS3, HblMADS5, HblMADS24 were up-regulated by MeJA. The function of HblMADS24 was elucidated. HblMADS24 bound HbFPS1 promoter in yeast and HblMADS24 activated HbFPS1 promoter in tobacco plants. Moreover, we proposed that HblMADS24 is a transcription activator of HbFPS1 which taking part in natural rubber biosynthesis.

Identification of histone methylation modifiers and their expression patterns during somatic embryogenesis in Hevea brasiliensis

Histone methylation plays a crucial role in various biological processes, from heterochromatin formation to transcriptional regulation. Currently, no information is available regarding histone methylation modifiers in the important rubber-producing plant Hevea brasiliensis. Here, we identified 47 histone methyltransferase (HMT) genes and 25 histone demethylase (HDM) genes as possible members of the histone methylation modifiers in the rubber tree genome. According to the structural features of HMT and HDM, the HbHMTs were classified into two groups (HbPRMs and HbSDGs), the HbHDMs have two groups (HbLSDs and HbJMJs). Expression patterns were analyzed in five different tissues and at different phases of somatic embryogenesis. HbSDG10, 21, 25, 33, HbJMJ2, 18, 20 were with high expression at different phases of somatic embryogenesis. HbSDG10,14, 20, 21, 33 and HbPRMT4 were expressed highly in anther, HbSDG14, 20, 21, 22, 23, 33, 35 and HbPRMT1 HbJMJ7 and HbLSD1, 2, 3, 4 showed high expression levels in callus. HbSDG1, 7, 10, 13, 14, 18, 19, 21, 22, 23, 35, HbPRMT1, 8, HbJMJ5, 7, 11, 16, 20 and HbLSD2, 3, 4 were expressed highly in somatic embryo. HbSDG10, 21, 25, 33, HbLSD2, 3 were expressed highly in bud of regenerated plant. The analyses reveal that HbHMTs and HbHDMs exhibit different expression patterns at different phases during somatic embryogenesis, implying that some HbHMTs and HbHDMs play important roles during somatic embryogenesis. This study provide fundamental information for further studies on histone methylation in Hevea brasiliensis.

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

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

BACKGROUND

Breeding 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.

RESULTS

By 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.

CONCLUSIONS

Taken 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.

Identification of the Sex-Biased Gene Expression and Putative Sex-Associated Genes in Eucommia ulmoides Oliver Using Comparative Transcriptome Analyses

Eucommia ulmoides is a model representative of the dioecious plants with sex differentiation at initiation. Nevertheless, the genetic mechanisms of sexual dimorphism and sex determination in E. ulmoides remain poorly understood. In this study de novo transcriptome sequencing on Illumina platform generated >45 billion high-quality bases from fresh leaves of six male and female individuals of E. ulmoides. A total of 148,595 unigenes with an average length of 801 base-pairs (bp) were assembled. Through comparative transcriptome analyses, 116 differentially expressed genes (DEGs) between the males and the females were detected, including 73 male-biased genes and 43 female-biased genes. Of these DEGs, three female-biased genes were annotated to be related with the sexually dimorphic gutta content in E. ulmoides. One male-biased DEG was identified as putative MADS box gene APETALA3, a B class floral organ identity gene in the flowering plants. SNPs calling analyses further confirmed that the APETALA3-like gene was probably involved in the sex determination in E. ulmoides. Four other male-biased DEGs were potential sex-associated genes as well with segregated SNPs in accord with sex type. In addition, the SNPs density was 1.02 per kilobase (kb) in the expressed genes of E. ulmoides, implying a relatively high genetic diversity.

Transcriptome-Wide Identification and Characterization of MYB Transcription Factor Genes in the Laticifer Cells of Hevea brasiliensis

MYB transcription factors hold vital roles in the regulation of plant secondary metabolic pathways. Laticifers in rubber trees (Hevea brasiliensis) are of primary importance in natural rubber production because natural rubber is formed and stored within these structures. To understand the role of MYB transcription factors in the specialized cells, we identified 44 MYB genes (named HblMYB1 to HblMYB44) by using our previously obtained transcriptome database of rubber tree laticifer cells and the public rubber tree genome database. Expression profiles showed that five MYB genes were highly expressed in the laticifers. HblMYB19 and HblMYB44 were selected for further study. HblMYB19 and HblMYB44 bound the promoters of HbFDPS1, HbSRPP, and HRT1 in yeast. Furthermore, the transient overexpression of HblMYB19 and HblMYB44 in tobacco plants significantly increased the activity of the promoters of HbFDPS1, HbSRPP, and HRT1. Basing on this information, we proposed that HblMYB19 and HblMYB44 are the regulators of HbFDPS1, HbSRPP, and HRT1, which are involved in the biosynthesis pathway of natural rubber.

Genome-Wide Identification of MicroRNAs and Their Targets in the Leaves and Fruits of Eucommia ulmoides Using High-Throughput Sequencing

MicroRNAs (miRNAs), a group of endogenous small non-coding RNAs, play important roles in plant growth, development, and stress response processes. Eucommia ulmoides Oliver (hardy rubber tree) is one of the few woody plants capable of producing trans-1, 4-polyisoprene (TPI), also known as Eu-rubber, which has been utilized as an industrial raw material and is extensively cultivated in China. However, the mechanism of TPI biosynthesis has not been identified in E. ulmoides. To characterize small RNAs and their targets with potential biological roles involved in the TPI biosynthesis in E. ulmoides, in the present study, eight small RNA libraries were constructed and sequenced from young and mature leaves and fruits of E. ulmoides. Further analysis identified 34 conserved miRNAs belonging to 20 families (two unclassified families), and 115 novel miRNAs seemed to be specific to E. ulmoides. Among these miRNAs, fourteen conserved miRNAs and 49 novel miRNAs were significantly differentially expressed and identified as Eu-rubber accumulation related miRNAs. Based on the E. ulmoides genomic data, 202 and 306 potential target genes were predicted for 33 conserved and 92 novel miRNAs, respectively; the predicted targets are mostly transcription factors and functional genes, which were enriched in metabolic pathways and biosynthesis of secondary metabolites. Noticeably, based on the expression patterns of miRNAs and their target genes in combination with the Eu-rubber accumulation, the negative correlation of expression of six miRNAs (Eu-miR14, Eu-miR91, miR162a, miR166a, miR172c, and miR396a) and their predicted targets serving as potential regulators in Eu-rubber accumulation. This study is the first to detect conserved and novel miRNAs and their potential targets in E. ulmoides and identify several candidate genes potentially controlling rubber accumulation, and thus provide molecular evidence for understanding the roles of miRNAs in regulating the TPI biosynthesis in E. ulmoides.