Transcriptome Analysis of the Signalling Networks in Coronatine-Induced Secondary Laticifer Differentiation from Vascular Cambia in Rubber Trees

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

Transcriptome analysis of Taraxacum kok-saghyz reveals the role of exogenous methyl jasmonate in regulating rubber biosynthesis and drought tolerance

Taraxacum kok-saghyz has been identified as one of the most promising alternative rubber crops, with laticifer cells that produce high-quality rubber. To uncover the underlying molecular mechanisms regulating natural rubber biosynthesis under MeJA induction, a reference transcriptome was constructed from nine samples of T. kok-saghyz. MeJA treatment was applied for 0 h (control), 6 h, and 24 h. A total of 7452 differentially expressed genes (DEGs) were identified in response to MeJA stress, relative to the control. Functional enrichment showed that these DEGs were primarily related to hormone signaling, defensive responses, and secondary metabolism. Combined analysis of the DEGs induced by MeJA and high-expression genes in laticifer cells further identified seven DEGs related to natural rubber biosynthesis that were upregulated in latex tissue, suggesting that these candidate genes could prove valuable in studying the mechanism of MeJA-mediated natural rubber biosynthesis. In addition, 415 MeJA-responsive DEGs were from several transcription factor families associated with drought resistance. This study helps to elucidate the mechanism of natural rubber biosynthesis in T. kok-saghyz in response to MeJA stress and identifies key candidate MeJA-induced DEGs in laticifer tissue, as well as a candidate drought-response target gene, whose knowledge will promote the breeding of T. kok-saghyz in the aspect of rubber yields and quality, and drought tolerance.

Genome-Wide Identification and Expression Analysis of MYC Transcription Factor Family Genes in Rubber Tree (Hevea brasiliensis Muell. Arg.)

Myelocytomatosis (MYC) transcription factors play a core regulator in the jasmonic acid signaling pathway, which regulates the secondary laticifer differentiation and rubber biosynthesis in rubber tree (Hevea brasiliensis). However, there are currently no reports on the MYC gene family in rubber trees, an important industrial raw material crop worldwide. In the present study, 32 HblMYCs were isolated and identified. The diversity in gene structure and presence of various cis-regulatory elements in promotors suggest that HblMYCs participate in various biological processes. Based on the expression patterns in the cambium region and laticifer in, respectively, response to coronatine (COR) and tapping, and the phylogenetic relationship with the MYCs that have been functionally identified in other plants, the HblMYC24 and HblMYC30 may be related to laticifer differentiation while the HblMYC6, HblMYC11 and HblMYC15, as well as HblMYC16 and HblMYC21, may positively regulate rubber biosynthesis. The results provide a foundation for understanding the molecular mechanism of jasmonate signaling in regulating laticifer differentiation and rubber biosynthesis in rubber tree.

Regulatory Potential of bHLH-Type Transcription Factors on the Road to Rubber Biosynthesis in Hevea brasiliensis

Natural rubber is the main component of latex obtained from laticifer cells of Hevea brasiliensis. For improving rubber yield, it is essential to understand the genetic molecular mechanisms responsible for laticifer differentiation and rubber biosynthesis. Jasmonate enhances both secondary laticifer differentiation and rubber biosynthesis. Here, we carried out time-course RNA-seq analysis in suspension-cultured cells treated with methyljasmonic acid (MeJA) to characterize the gene expression profile. Gene Ontology (GO) analysis showed that the term "cell differentiation" was enriched in upregulated genes at 24 hours after treatment, but inversely, the term was enriched in downregulated genes at 5 days, indicating that MeJA could induce cell differentiation at an early stage of the response. Jasmonate signaling is activated by MYC2, a basic helix-loop-helix (bHLH)-type transcription factor (TF). The aim of this work was to find any links between transcriptomic changes after MeJA application and regulation by TFs. Using an in vitro binding assay, we traced candidate genes throughout the whole genome that were targeted by four bHLH TFs: Hb_MYC2-1, Hb_MYC2-2, Hb_bHLH1, and Hb_bHLH2. The latter two are highly expressed in laticifer cells. Their physical binding sites were found in the promoter regions of a variety of other TF genes, which are differentially expressed upon MeJA exposure, and rubber biogenesis-related genes including SRPP1 and REF3. These studies suggest the possibilities that Hb_MYC2-1 and Hb_MYC2-2 regulate cell differentiation and that Hb_bHLH1 and Hb_bHLH2 promote rubber biosynthesis. We expect that our findings will help to increase natural rubber yield through genetic control in the future.

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.

Opposite physiological effects upon jasmonic acid and brassinosteroid treatment on laticifer proliferation and co-occurrence of differential expression of genes involved in vascular development in rubber tree

During growth of woody plant-trunk, the secondary meristem functions in giving rise the xylem and phloem. Rubber tree (Hevea brasiliensis Muell. Arg.), in addition, contains laticifers (latex producing vessels) in the vicinity of phloem. Insights into regulatory mechanisms of gene networks underlying laticifer proliferation in rubber tree has remained very limited. The candidate vascular development-related genes were selected to investigate for expression profile in phloem and xylem tissues of high latex yield- and high wood yield-clones of rubber tree. The differential gene expression between the mature branch-xylem and -phloem tissues was clearly observed. The cis-regulatory motif analysis revealed the existent of putative jasmonic acid (JA)- and brassinosteroid (BR)-responsive regulatory motifs in promoter regions of these genes, and consequently the effect of exogenous application of JA, BR or their respective signaling inhibitors, on the formation of laticifers in rubber tree was demonstrated. Interestingly, the laticifer numbers were significantly increased in JA-treatment, correlated with up-regulation of phloem development-related genes in both rubber tree clones. On the contrary, the laticifers were decreased in BR-treatment accompanying by up-regulation of xylem development-related genes, especially in high wood yield-rubber tree clone. BR-inhibitor treatment also enhanced laticifer numbers, while JA-inhibitor suppressed laticifer differentiation. Taken together, this study unveils the molecular interplay between JA/BR on vascular development in rubber tree and how this impacts the appearance of laticifers in this plant. This process is vital for a better understanding on laticifer differentiation and its impact in the manipulation of wood and latex yield in rubber tree improvement program.

Genome-Wide Identification and Characterization of the JAZ Gene Family in Rubber Tree (Hevea brasiliensis)

Jasmonate signaling plays a vital role in the regulation of secondary laticifer differentiation and natural rubber biosynthesis in Hevea brasiliensis. Jasmonate ZIM-domain (JAZ) proteins are the master regulators of jasmonate signaling. Although several JAZs have been reported in the laticifer cells of H. brasiliensis, the genome-wide screening of HbJAZ members has not yet been explored. In the present study, 18 HbJAZs were identified based on the recent H. brasiliensis genome. Phylogenetic construction revealed that the HbJAZs were clustered into five subgroups and that members within the same subgroup shared highly conserved gene structures and protein motifs. Cis-element analysis of HbJAZ promoters suggested the presence of hormone, stress and development-related cis-elements. HbJAZ1.0, HbJAZ2.0, and HbJAZ5.0 interacted with CORONATINE INSENSITIVE1 (COI1) in the presence of coronatine (COR, a JA mimic). HbJAZ1.0, HbJAZ2.0, HbJAZ5.0, and HbJAZ12.0 could also interact with each other. Of the 18 HbJAZs, transcripts of 15 HbJAZs were present in the vascular cambium region except for that of HbJAZ7.0, HbJAZ8.0d, and HbJAZ13.0. Fourteen of the 15 HbJAZs were significantly up-regulated upon COR treatment. The transcripts of three genes that were absent from vascular cambium region were also absent from the latex. Among the 15 HbJAZs in the latex, the expression patterns of 13 HbJAZs were different between the tapping and ethrel treatments. Eight of the 14 COR-up-regulated HbJAZs in the vascular cambium region were also activated by tapping in latex. Of the eight tapping-activated HbJAZs, 5 HbJAZs were repressed by ethrel application. Based on the computational analyses and gene expression patterns described in this study, the HbJAZ5.0 and HbJAZ10.0b may be associated with laticifer differentiation while the HbJAZ8.0b is a negative regulator for natural rubber biosynthesis in H. brasiliensis.

Papain-like cysteine protease encoding genes in rubber (Hevea brasiliensis): comparative genomics, phylogenetic, and transcriptional profiling analysis

43 HbPLCPs representing nine subfamilies or 20 orthologous groups were found in rubber, where paralogs were resulted from the recent WGD and local duplication. Several senescence-associated genes were also identified. Papain-like cysteine proteases (PLCPs) comprise a large family of proteolytic enzymes involved in plant growth and development, seed germination, organ senescence, immunity, and stress response. Despite their importance and the extensive research in the model plant Arabidopsis thaliana, little information is available on rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study performed a genome-wide identification of PLCP family genes in rubber, resulting in a relatively high number of 43 members. The phylogenetic analysis assigned these genes into nine subfamilies, i.e., RD21 (6), CEP (4), XCP (4), XBCP3 (2), THI (1), SAG12 (18), RD19 (4), ALP (2), and CTB (2). Most of them were shown to have orthologs in Arabidopsis; however, several members in SAG12, CEP and XBCP3 subfamilies form new groups as observed in other core eudicots such as Manihot esculenta, Ricinus communis, Populus trichocarpa, and Vitis vinifera. Based on an expert sequence comparison, 20 orthologous groups (OGs) were proposed for core eudicots, and rubber paralogs were shown to be resulted from the recent whole-genome duplication (WGD) as well as local duplication. Transcriptional profiling showed distinct expression pattern of different members across various tissues, e.g., root, leaf, bark, laticifer, flower, and seed. By using the senescence-specific HbSAG12H1 as the indicator, the transcriptome of senescent rubber leaves was deeply sequenced and several senescence-associated PLCP genes were identified. Results obtained from this study provide valuable information for future functional analysis and utilization of PLCP genes in Hevea and other species.

Transcriptional analysis of genes involved in competitive nodulation in Bradyrhizobium diazoefficiens at the presence of soybean root exudates

Nodulation competition is a key factor that limits symbiotic nitrogen fixation between rhizobia and their host legumes. Soybean root exudates (SREs) are thought to act as signals that influence Bradyrhizobium ability to colonize roots and to survive in the rhizosphere, and thus they act as a key determinant of nodulation competitiveness. In order to find the competitiveness-related genes in B. diazoefficiens, the transcriptome of two SREs treated B. diazoefficiens with completely different nodulation abilities (B. diazoefficiens 4534 and B. diazoefficiens 4222) were sequenced and compared. In SREs treated strain 4534 (SREs-4534), 253 unigenes were up-regulated and 204 unigenes were down-regulated. In SREs treated strain 4534 (SREs-4222), the numbers of up- and down-regulated unigenes were 108 and 185, respectively. There were considerable differences between the SREs-4534 and SREs-4222 gene expression profiles. Some differentially expressed genes are associated with a two-component system (i.g., nodW, phyR-σEcfG), bacterial chemotaxis (i.g., cheA, unigene04832), ABC transport proteins (i.g., unigene02212), IAA (indole-3-acetic acid) metabolism (i.g., nthA, nthB), and metabolic fitness (i.g., put.), which may explain the higher nodulation competitiveness of B. diazoefficiens in the rhizosphere. Our results provide a comprehensive transcriptomic resource for SREs treated B. diazoefficiens and will facilitate further studies on competitiveness-related genes in B. diazoefficiens.