Profiling Ethylene-Responsive Genes Expressed in the Latex of the Mature Virgin Rubber Trees Using cDNA Microarray.

Ethylene is commonly used as a latex stimulant of Hevea brasiliensis by application of ethephon (chloro-2-ethylphosphonic acid); however, the molecular mechanism by which ethylene increases latex production is not clear. To better understand the effects of ethylene stimulation on the laticiferous cells of rubber trees, a latex expressed sequence tag (EST)-based complementary DNA microarray containing 2,973 unique genes (probes) was first developed and used to analyze the gene expression changes in the latex of the mature virgin rubber trees after ethephon treatment at three different time-points: 8, 24 and 48 h. Transcript levels of 163 genes were significantly altered with fold-change values ≥ 2 or ≤ -2 (q-value < 0.05) in ethephon-treated rubber trees compared with control trees. Of the 163 genes, 92 were up-regulated and 71 down-regulated. The microarray results were further confirmed using real-time quantitative reverse transcript-PCR for 20 selected genes. The 163 ethylene-responsive genes were involved in several biological processes including organic substance metabolism, cellular metabolism, primary metabolism, biosynthetic process, cellular response to stimulus and stress. The presented data suggest that the laticifer water circulation, production and scavenging of reactive oxygen species, sugar metabolism, and assembly and depolymerization of the latex actin cytoskeleton might play important roles in ethylene-induced increase of latex production. The results may provide useful insights into understanding the molecular mechanism underlying the effect of ethylene on latex metabolism of H. brasiliensis.

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.

In-depth proteome analysis of the rubber particle of Hevea brasiliensis (para rubber tree).

The rubber particle is a special organelle in which natural rubber is synthesised and stored in the laticifers of Hevea brasiliensis. To better understand the biological functions of rubber particles and to identify the candidate rubber biosynthesis-related proteins, a comprehensive proteome analysis was performed on H. brasiliensis rubber particles using shotgun tandem mass spectrometry profiling approaches-resulting in a thorough report on the rubber particle proteins. A total of 186 rubber particle proteins were identified, with a range in relative molecular mass of 3.9-194.2 kDa and in isoelectric point values of 4.0-11.2. The rubber particle proteins were analysed for gene ontology and could be categorised into eight major groups according to their functions: including rubber biosynthesis, stress- or defence-related responses, protein processing and folding, signal transduction and cellular transport. In addition to well-known rubber biosynthesis-related proteins such as rubber elongation factor (REF), small rubber particle protein (SRPP) and cis-prenyl transferase (CPT), many proteins were firstly identified to be on the rubber particles, including cyclophilin, phospholipase D, cytochrome P450, small GTP-binding protein, clathrin, eukaryotic translation initiation factor, annexin, ABC transporter, translationally controlled tumour protein, ubiquitin-conjugating enzymes, and several homologues of REF, SRPP and CPT. A procedure of multiple reaction monitoring was established for further protein validation. This comprehensive proteome data of rubber particles would facilitate investigation into molecular mechanisms of biogenesis, self-homeostasis and rubber biosynthesis of the rubber particle, and might serve as valuable biomarkers in molecular breeding studies of H. brasiliensis and other alternative rubber-producing species.

Identification of the Hevea brasiliensis AP2/ERF superfamily by RNA sequencing.

BACKGROUND: Rubber tree (Hevea brasiliensis) laticifers are the source of natural rubber. Rubber production depends on endogenous and exogenous ethylene (ethephon). AP2/ERF transcription factors, and especially Ethylene-Response Factors, play a crucial role in plant development and response to biotic and abiotic stresses. This study set out to sequence transcript expressed in various tissues using next-generation sequencing and to identify AP2/ERF superfamily in the rubber tree. RESULTS: The 454 sequencing technique was used to produce five tissue-type transcript libraries (leaf, bark, latex, embryogenic tissues and root). Reads from all libraries were pooled and reassembled to improve mRNA lengths and produce a global library. One hundred and seventy-three AP2/ERF contigs were identified by in silico analysis based on the amino acid sequence of the conserved AP2 domain from the global library. The 142 contigs with the full AP2 domain were classified into three main families (20 AP2 members, 115 ERF members divided into 11 groups, and 4 RAV members) and 3 soloist members. Fifty-nine AP2/ERF transcripts were found in latex. Alongside the microRNA172 already described in plants, eleven additional microRNAs were predicted to inhibit Hevea AP2/ERF transcripts. CONCLUSIONS: Hevea has a similar number of AP2/ERF genes to that of other dicot species. We adapted the alignment and classification methods to data from next-generation sequencing techniques to provide reliable information. We observed several specific features for the ERF family. Three HbSoloist members form a group in Hevea. Several AP2/ERF genes highly expressed in latex suggest they have a specific function in Hevea. The analysis of AP2/ERF transcripts in Hevea presented here provides the basis for studying the molecular regulation of latex production in response to abiotic stresses and latex cell differentiation.

Proteome analysis of the large and the small rubber particles of Hevea brasiliensis using 2D-DIGE.

The rubber particle is a specialized organelle in which natural rubber is synthesised and stored in the laticifers of Hevea brasiliensis (para rubber tree). It has been demonstrated that the small rubber particles (SRPs) has higher rubber biosynthesis ratio than the large rubber particles (LRPs), but the underlying molecular mechanism still remains unknown. In this study, LRPs and SRPs were firstly separated from the fresh latex using differential centrifugation, and two-dimensional difference in-gel electrophoresis (2D-DIGE) combined with MALDI-TOF/TOF was then applied to investigate the proteomic alterations associated with the changed rubber biosynthesis capacity between LRPs and SRPs. A total of 53 spots corresponding to 22 gene products, were significantly altered with the |ratio|≥2.0 and T value ≤0.05, among which 15 proteins were up-regulated and 7 were down-regulated in the SRPs compared with the LRPs. The 15 up-regulated proteins in the SRPs included small rubber particle protein (SRPP), 3-hydroxy-3-methylglutaryl-CoA synthase (HMGCS), phospholipase D alpha (PLD α), ethylene response factor 2, eukaryotic translation initiation factor 5A isoform IV (eIF 5A-4), 70-kDa heat shock cognate protein (HSC 70), several unknown proteins, etc., whereas the 7 up-regulated proteins in the LRPs were rubber elongation factor (REF, 19.6kDa), ASR-like protein 1, REF-like stress-related protein 1, a putative phosphoglyceride transfer family protein, ß-1,3-glucanase, a putative retroelement, and a hypothetical protein. Since several proteins related to rubber biosynthesis were differentially expressed between LRPs and SRPs, the comparative proteome data may provide useful insights into understanding the mechanism involved in rubber biosynthesis and latex coagulation in H. brasiliensis.

Identification of novel microRNAs in Hevea brasiliensis and computational prediction of their targets.

BACKGROUND: Plants respond to external stimuli through fine regulation of gene expression partially ensured by small RNAs. Of these, microRNAs (miRNAs) play a crucial role. They negatively regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs (mRNAs). In Hevea brasiliensis, environmental and harvesting stresses are known to affect natural rubber production. This study set out to identify abiotic stress-related miRNAs in Hevea using next-generation sequencing and bioinformatic analysis. RESULTS: Deep sequencing of small RNAs was carried out on plantlets subjected to severe abiotic stress using the Solexa technique. By combining the LeARN pipeline, data from the Plant microRNA database (PMRD) and Hevea EST sequences, we identified 48 conserved miRNA families already characterized in other plant species, and 10 putatively novel miRNA families. The results showed the most abundant size for miRNAs to be 24 nucleotides, except for seven families. Several MIR genes produced both 20-22 nucleotides and 23-27 nucleotides. The two miRNA class sizes were detected for both conserved and putative novel miRNA families, suggesting their functional duality. The EST databases were scanned with conserved and novel miRNA sequences. MiRNA targets were computationally predicted and analysed. The predicted targets involved in "responses to stimuli" and to "antioxidant" and "transcription activities" are presented. CONCLUSIONS: Deep sequencing of small RNAs combined with transcriptomic data is a powerful tool for identifying conserved and novel miRNAs when the complete genome is not yet available. Our study provided additional information for evolutionary studies and revealed potentially specific regulation of the control of redox status in Hevea.

Differential gene expression in different types of Hevea brasiliensis roots.

Three types of roots (taproots, first order laterals and second order laterals) were functionally characterized on 7-month-old in vitro plantlets regenerated by somatic embryogenesis in Hevea brasiliensis. A histological analysis revealed different levels of differentiation depending on root diameter. A primary structure was found in first and second order lateral roots, while taproots displayed a secondary structure. The expression of 48 genes linked to some of the regulatory pathways acting in roots was compared in leaves, stems and the different types of roots by real-time RT-PCR. Thirteen genes were differentially expressed in the different organs studied in plants grown under control conditions. Nine additional other genes were differentially regulated between organs under water deficit conditions. In addition, 10 genes were significantly regulated in response to water deficit, including 8 regulated mainly in lateral roots types. Our results suggest that the regulation of gene expression in lateral roots is different than that in taproots, which have a main role in nutrient uptake and transport, respectively.

Gene expression pattern in response to wounding, methyl jasmonate and ethylene in the bark of Hevea brasiliensis.

Natural rubber production in Hevea brasiliensis is determined by both tapping and ethephon frequencies. It is affected by a complex physiological disorder called tapping panel dryness. This syndrome is likely to be induced by environmental and latex harvesting stresses. Defence responses, including rubber biosynthesis, are dramatically mediated by wounding, jasmonate and ethylene (ET), among other factors. Using real-time RT-PCR, the effects of wounding, methyl jasmonate (MeJA) and ET on the relative transcript abundance of a set of 25 genes involved in their signalling and metabolic pathways were studied in the bark of 3-month-old epicormic shoots. Temporal regulation was found for 9 out of 25 genes. Wounding treatment regulated the transcript abundance of 10 genes. Wounding-specific regulation was noted for the HbMAPK, HbBTF3b, HbCAS1, HbLTPP and HbPLD genes. MeJA treatment regulated the transcript abundance of nine genes. Of these, the HbMYB, HbCAS2, HbCIPK and HbChi genes were shown to be specifically MeJA inducible. ET response was accompanied by regulation of the transcript abundance of eight genes, and six genes, HbETR2, HbEIN2, HbEIN3, HbCaM, HbPIP1 and HbQM, were specifically regulated by ET treatment. Additionally, the transcript level of the HbGP and HbACR genes was enhanced by all three treatments simultaneously. Overall, a large number of genes were found to be regulated 4 h after the treatments were applied. This study nevertheless revealed some jasmonic acid-independent wound signalling pathways in H. brasiliensis, provided a general characterization of signalling pathways and will serve as a new base from which to launch advanced studies of the network of pathways operating in H. brasiliensis.