Molecular Identification and Characterization of a Gene Associated with the Onset of Tapping Panel Dryness (TPD) Syndrome in Rubber Tree (Hevea brasiliensis Muell.) by mRNA Differential Display

Integrated physiology, transcriptome and proteome analyses highlight the potential roles of multiple hormone-mediated signaling pathways involved in tapping panel dryness in rubber tree

Currently, one of the most serious threats to rubber tree is the tapping panel dryness (TPD) that greatly restricts natural rubber production. Over-tapping or excessive ethephon stimulation is regarded as the main cause of TPD occurrence. Although extensive studies have been carried out, the molecular mechanism underlying TPD remains puzzled. An attempt was made to compare the levels of endogenous hormones and the profiles of transcriptome and proteome between healthy and TPD trees. Results showed that most of endogenous hormones such as jasmonic acid (JA), 1-aminocyclopropanecarboxylic acid (ACC), indole-3-acetic acid (IAA), trans-zeatin (tZ) and salicylic acid (SA) in the barks were significantly altered in TPD-affected rubber trees. Accordingly, multiple hormone-mediated signaling pathways were changed. In total, 731 differentially expressed genes (DEGs) and 671 differentially expressed proteins (DEPs) were identified, of which 80 DEGs were identified as putative transcription factors (TFs). Further analysis revealed that 12 DEGs and five DEPs regulated plant hormone synthesis, and that 16 DEGs and six DEPs were involved in plant hormone signal transduction pathway. Nine DEGs and four DEPs participated in rubber biosynthesis and most DEGs and all the four DEPs were repressed in TPD trees. All these results highlight the potential roles of endogenous hormones, signaling pathways mediated by these hormones and rubber biosynthesis pathway in the defense response of rubber trees to TPD. The present study extends our understanding of the nature and mechanism underlying TPD and provides some candidate genes and proteins related to TPD for further research in the future.

Bark transcriptome analyses reveals molecular mechanisms involved in tapping panel dryness occurrence and development in rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) has become the mostimportant limiting factor for increasing natural rubber yield, whereas illuminating the molecular mechanisms underlying TPD is the prerequisite for solving the problem of TPD. However, molecular mechanisms underlying TPD are largely unknown. In this study, healthy and different stages of TPD-affected rubber trees were utilized to analyze TPD for the first time. We found that the changing tendencies of key latex physiological parameters were closely related to TPD occurrence and development. To reveal the molecular mechanisms underlying TPD, we sequenced and compared bark transcriptomes among healthy rubber tree, and TPD-affected ones at initial and advanced stages. In total, 8607 genes were identified as TPD-related genes in contrast to healthy rubber tree. According to gene expression profiles, the five samples were divided into three groups including healthy rubber tree, and TPD-affected rubber tree in the initial and advanced stages, which was consistent with the stages of TPD occurrence and development. Interestingly, only asmall proportionof the TPD-related genes were constantly down- or up-regulated with TPD occurrence and development. The TPD-related genes in KEGG pathways significantly enriched were closely associated with protein metabolism, cell division and differentiation, PCD, stress responses, terpene biosynthesis, and various metabolism processes. Moreover, overexpression of HbAPX2 identified as a TPD-related gene enhanced oxidative stress tolerance in S. cerevisiae. The typical symptoms of TPD, partial or complete dry zone (no latex flow) on tapping panel, might attribute to lower IPP available for rubber biosynthesis, and downregulation of the genes in post-IPP steps of rubber biosynthesis and the genes involved in latex flow. Our results not only provide new insights into molecular mechanisms underlying TPD occurrence and development but also contribute to developing effective measures to control TPD in rubber trees.

Downregulation of HbFPS1 affects rubber biosynthesis of Hevea brasiliensis suffering from tapping panel dryness

Tapping panel dryness (TPD) is a century-old problem that has plagued the natural rubber production of Hevea brasiliensis. TPD may result from self-protective mechanisms of H. brasiliensis in response to stresses such as excessive hormone stimulation and mechanical wounding (bark tapping). It has been hypothesized that TPD impairs rubber biosynthesis; however, the underlying mechanisms remain poorly understood. In the present study, we firstly verified that TPD-affected rubber trees exhibited lower rubber biosynthesis activity and greater rubber molecular weight compared to healthy rubber trees. We then demonstrated that HbFPS1, a key gene of rubber biosynthesis, and its expression products were downregulated in the latex of TPD-affected rubber trees, as revealed by transcriptome sequencing and iTRAQ-based proteome analysis. We further discovered that the farnesyl diphosphate synthase HbFPS1 could be recruited to small rubber particles by HbSRPP1 through protein-protein interactions to catalyze farnesyl diphosphate (FPP) synthesis and facilitate rubber biosynthesis initiation. FPP content in the latex of TPD-affected rubber trees was significantly decreased with the downregulation of HbFPS1, ultimately resulting in abnormal development of rubber particles, decreased rubber biosynthesis activity, and increased rubber molecular weight. Upstream regulator assays indicated that a novel regulator, MYB2-like, may be an important regulator of downregulation of HbFPS1 in the latex of TPD-affected rubber trees. Our findings not only provide new directions for studying the molecular events involved in rubber biosynthesis and TPD syndrome and contribute to rubber management strategies, but also broaden our knowledge of plant isoprenoid metabolism and its regulatory networks.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis

Background

Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea.

Methods

Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20–22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators.

Results

Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Characterization of the rubber tree metallothionein family reveals a role in mitigating the effects of reactive oxygen species associated with physiological stress

Metallothioneins (MTs) as reactive oxygen species (ROS) scavengers play important roles in stress response and heavy metal homeostasis. In Hevea brasiliensis (the para rubber tree that is the source of commercial natural rubber) and in other trees, the functions of MTs are not well understood. Latex exudes when the rubber tree is tapped. The flow of latex and its regeneration can be enhanced by tapping, wounding and ethylene treatment, all of which produce ROS as a by-product. Here, we show the presence of four MT genes in H. brasiliensis, comprising three Type 2 (HbMT2, -2a and -2b) and one Type 3 (HbMT3L) isoforms, representing one of the smallest MT gene families among angiosperms. The four HbMTs exhibited distinct tissue expression patterns: HbMT2 and HbMT3L mainly in leaves, HbMT2a specifically in flowers and HbMT2b in diverse tissues. The expression of HbMT2b, an isoform present in latex, decreased significantly in the latex following the stress-inducing treatments of tapping, wounding and ethephon (an ethylene generator). The expressions of the leaf-abundant isoforms, HbMT2 and -3L were up-regulated following pathogenic fungus infection and high-temperature stress, but down-regulated by low-temperature stress. These reactions were consistent with multiple defense- and hormone-responsive cis-acting elements in the HbMT promoters. Nine transcription factors were shown to implicate in the high-temperature responsiveness of HbMT2 and -3L in leaves. Overexpression of HbMT2 in Escherichia coli enhanced the bacterium's tolerance to heavy metals and ROS, consistent with its predicted role as an ROS scavenger. Taken together, our results, along with other relevant studies, suggest an important role of HbMTs in latex regeneration as well as species adaptation via the regulation of ROS homeostasis.

Over-expression of 3-hydroxy-3- methylglutaryl-coenzyme A reductase 1 (hmgr1) gene under super-promoter for enhanced latex biosynthesis in rubber tree (Hevea brasiliensis Muell. Arg.)

Natural rubber (cis-1, 4-polyisoprene) is being produced from bark laticifer cells of Hevea brasiliensis and the popular high latex yielding Indian rubber clones are easily prone to onset of tapping panel dryness syndrome (TPD) which is considered as a physiological syndrome affecting latex production either partially or completely. This report describes an efficient protocol for development of transgenic rubber plants by over-expression of 3-hydroxy 3-methylglutaryl Co-enzyme A reductase 1 (hmgr1) gene which is considered as rate limiting factor for latex biosynthesis via Agrobacterium-mediated transformation. The pBIB plasmid vector containing hmgr1 gene cloned under the control of a super-promoter was used for genetic transformation using embryogenic callus. Putatively transgenic cell lines were obtained on selection medium and produced plantlets with 44% regeneration efficiency. Transgene integration was confirmed by PCR amplification of 1.8 kb hmgr1 and 0.6 kb hpt genes from all putatively transformed callus lines as well as transgenic plants. Southern blot analysis showed the stable integration and presence of transgene in the transgenic plants. Over expression of hmgr1 transgene was determined by Northern blot hybridization, semi-quantitative PCR and real-time PCR (qRT-PCR) analysis. Accumulation of hmgr1 mRNA transcripts was more abundant in transgenic plants than control. Increased level of photosynthetic pigments, protein contents and HMGR enzyme activity was also noticed in transgenic plants over control. Interestingly, the latex yield was significantly enhanced in all transgenic plants compared to the control. The qRT-PCR results exhibit that the hmgr1 mRNA transcript levels was 160-fold more abundance in transgenic plants over untransformed control. These results altogether suggest that there is a positive correlation between latex yield and accumulation of mRNA transcripts level as well as HMGR enzyme activity in transgenic rubber plants. It is presumed that there is a possibility for enhanced level of latex biosynthesis in transgenic plants as the level of mRNA transcripts and HMGR enzyme activity is directly correlated with latex yield in rubber tree. Further, the present results clearly suggest that the quantification of HMGR enzyme activity in young seedlings will be highly beneficial for early selection of high latex yielding plants in rubber breeding programs.

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.

Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis

Background

Rubber tree (Hevea brasiliensis) is an important industrial crop cultivated in tropical areas for natural rubber production. Treatment of the bark of rubber trees with ehephon (an ethylene releaser) has been a routine measure to increase latex yield, but the molecular mechanism behind the stimulation of rubber production by ethylene still remains a puzzle. Deciphering the enigma is of great importance for improvement of rubber tree for high yield.

Results

De novo sequencing and assembly of the bark transciptomes of Hevea brasiliensis induced with ethephon for 8 h (E8) and 24 h (E24) were performed. 51,965,770, 52,303,714 and 53,177,976 high-quality clean reads from E8, E24 and C (control) samples were assembled into 81,335, 80,048 and 80,800 unigenes respectively, with a total of 84,425 unigenes and an average length of 1,101 bp generated. 10,216 and 9,374 differentially expressed genes (DEGs) in E8 and E24 compared with C were respectively detected. The expression of several enzymes in crucial points of regulation in glycolysis were up-regulated and DEGs were not significantly enriched in isopentenyl diphosphate (IPP) biosynthesis pathway. In addition, up-regulated genes of great regulatory importance in carbon fixation (Calvin cycle) were identified.

Conclusions

The rapid acceleration of glycolytic pathway supplying precursors for the biosynthesis of IPP and natural rubber, instead of rubber biosynthesis per se, may be responsible for ethylene stimulation of latex yield in rubber tree. The elevated rate of flux throughout the Calvin cycle may account for some durability of ethylene-induced stimulation. Our finding lays the foundations for molecular diagnostic and genetic engineering for high-yielding improvement of rubber tree.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2587-4) contains supplementary material, which is available to authorized users.

Identification and characterization of MAGO and Y14 genes in Hevea brasiliensis

Mago nashi (MAGO) and Y14 proteins are highly conserved among eukaryotes. In this study, we identified two MAGO (designated as HbMAGO1 andHbMAGO2) and two Y14 (designated as HbY14aand HbY14b) genes in the rubber tree (Hevea brasiliensis) genome annotation. Multiple amino acid sequence alignments predicted that HbMAGO and HbY14 proteins are structurally similar to homologous proteins from other species. Tissue-specific expression profiles showed that HbMAGO and HbY14 genes were expressed in at least one of the tissues (bark, flower, latex, leaf and root) examined. HbMAGOs and HbY14s were predominately located in the nucleus and were found to interact in yeast two-hybrid analysis (YTH) and bimolecular fluorescence complementation (BiFC) assays. HbMAGOs and HbY14s showed the highest transcription in latex and were regulated by ethylene and jasmonate. Interaction between HbMAGO2 and gp91^phox (a large subunit of nicotinamide adenine dinucleotide phosphate) was identified using YTH and BiFC assays. These findings suggested that HbMAGO may be involved in the aggregation of rubber particles in H. brasiliensis.

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.

Identification and Characterization of the Glucose-6-Phosphate Dehydrogenase Gene Family in the Para Rubber Tree, Hevea brasiliensis

As a key enzyme in the pentose phosphate pathway (PPP), glucose-6-phosphate dehydrogenase (G6PDH) provides nicotinamide adenine dinucleotide phosphate (NADPH) and intermediary metabolites for rubber biosynthesis, and plays an important role in plant development and stress responses. In this study, four Hevea brasiliensis (Para rubber tree) G6PDH genes (HbG6PDH1 to 4) were identified and cloned using a genome-wide scanning approach. All four HbG6PDH genes encode functional G6PDH enzymes as shown by heterologous expression in E. coli. Phylogeny analysis and subcellular localization prediction show that HbG6PDH3 is a cytosolic isoform, while the other three genes (HbG6PDH1, 2 and 4) are plastidic isoforms. The subcellular locations of HbG6PDH3 and 4, two latex-abundant isoforms were further verified by transient expression in rice protoplasts. Enzyme activity assay and expression analysis showed HbG6PDH3 and 4 were implicated in PPP during latex regeneration, and to influence rubber production positively in rubber tree. The cytosolic HbG6PDH3 is a predominant isoform in latex, implying a principal role for this isoform in controlling carbon flow and NADPH production in the PPP during latex regeneration. The expression pattern of plastidic HbG6PDH4 correlates well with the degree of tapping panel dryness, a physiological disorder that stops the flow of latex from affected rubber trees. In addition, the four HbG6PDHs responded to temperature and drought stresses in root, bark, and leaves, implicating their roles in maintaining redox balance and defending against oxidative stress.

Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD)

Background

Tapping panel dryness (TPD) involves in the partial or complete cessation of latex flow thus seriously affect latex production in rubber tree (Hevea brasiliensis). Numerous studies have been conducted to define its origin and nature, but the molecular nature and mechanism of TPD occurrence remains unknown. This study is committed to de novo sequencing and comparative analysis of the transcriptomes of healthy (H) and TPD-affected (T) rubber trees to identify the genes and pathways related to the TPD.

Results

Total raw reads of 34,632,012 and 35,913,020 bp were obtained from H and T library, respectively using Illumina Hiseq 2000 sequencing technology. De novo assemblies yielded 141,456 and 169,285 contigs, and 96,070 and 112,243 unigenes from H and T library, respectively.

Among 73597 genes, 22577 genes were identified as differential expressed genes between H and T library via comparative transcript profiling. A majority of genes involved in natural rubber biosynthesis and jasmonate synthesis with most potential relevance in TPD occurrence were found to be differentially expressed.

Conclusions

In TPD-affected trees, the expression of most genes related to the latex biosynthesis and jasmonate synthesis was severely inhibited and is probably the direct cause of the TPD. These new de novo transcriptome data sets provide a significant resource for the discovery of genes related to TPD and improve our understanding of the occurrence and maintainace of TPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1562-9) contains supplementary material, which is available to authorized users.

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.

Molecular cloning and expression of a novel MYB transcription factor gene in rubber tree

MYB family proteins regulate a variety of cellular processes in plants. Tapping panel dryness (TPD) in rubber tree (Hevea brasiliensis Muell. Arg.) affects latex biosynthesis and causes serious losses to rubber producers. In this study, a novel SANT/MYB transcription factor gene down-regulated in TPD rubber tree, named as HbSM1, was isolated from rubber tree. The complete HbSM1 open reading frame (ORF) was 948 bp in length. The deduced HbSM1 protein is 315 amino acids. HbSM1 belonged to 1RMYB subfamily with a single SANT domain. Sequence alignment revealed that HbSM1 had high homology with MYB members from Ricinus communis and Manihot esculenta, with 72 and 78 % identity, respectively. Moreover, HbSM1 shared 56 % identity with Glycine max GmMYB176. Phylogenetic analysis revealed that HbSM1, GmMYB176, rice OsMYBS2, and OsMYBS3 fell into the same cluster with 93 % bootstrap support value. Comparing expression among different tissues demonstrated that HbSM1 was ubiquitously expressed in all tissues, but it appeared to be preferentially expressed in leaf and latex. Furthermore, HbSM1 transcripts were significantly induced by various phytohormones (including gibberellic acid, ethephon, methyl jasmonate, salicylic acid, and abscisic acid) and wounding treatments. These results suggested that HbSM1 might play multiple roles in plant development via different phytohormones signaling pathways.

Differential expression of microRNAs and their targets reveals a possible dual role in physiological bark disorder in rubber tree

Trunk phloem necrosis (TPN), a physiological bark disorder of the rubber tree (Hevea brasiliensis), is a serious problem that affects the yield of natural rubber. The resultant bark dryness occurs in up to half of a plantation's trees in almost every rubber tree plantation region, causing a great annual loss of dry rubber for natural rubber production. Different types of injury and physical damage caused by mechanical activation as well as environmental stresses cause physiological bark disorder in tree. Due to the essential role of miR166, miR393 and miR167 in vascular development and abiotic stress response in diverse plant species, it was interesting to investigate the role of these miRNAs in rubber trees, particularly during development of a physiological bark disorder. In this study, the expression pattern of miR166, miR393 and miR167; and their target genes, HD-ZIP III; TIR1 and ARF8, respectively; was demonstrated in healthy tree and different TPN trees. Their existence and function in vivo was validated using RNA ligase-mediated 5' rapid amplification of cDNA ends. Taken together, the results suggest a possible dual role of these three miRNAs in maintaining normal bark regeneration in healthy trees, coping with overtapping by affecting the wound healing system leading to abnormal bark regeneration in overtapped-TPN trees, and act as additional forces that enhance the attenuation of vascular development resulting in bark necrosis and cell death in the natural-TPN tree. This is the first study to address the molecular events of miRNAs involved in the physiological bark disorder TPN in rubber tree. Further study will open the possibility to better understanding of physiological and molecular perspectives during TPN development, and lead to improvement of monitoring the exploitation of rubber tree plantations.

Sequence and expression analyses of ethylene response factors highly expressed in latex cells from Hevea brasiliensis

The AP2/ERF superfamily encodes transcription factors that play a key role in plant development and responses to abiotic and biotic stress. In Hevea brasiliensis, ERF genes have been identified by RNA sequencing. This study set out to validate the number of HbERF genes, and identify ERF genes involved in the regulation of latex cell metabolism. A comprehensive Hevea transcriptome was improved using additional RNA reads from reproductive tissues. Newly assembled contigs were annotated in the Gene Ontology database and were assigned to 3 main categories. The AP2/ERF superfamily is the third most represented compared with other transcription factor families. A comparison with genomic scaffolds led to an estimation of 114 AP2/ERF genes and 1 soloist in Hevea brasiliensis. Based on a phylogenetic analysis, functions were predicted for 26 HbERF genes. A relative transcript abundance analysis was performed by real-time RT-PCR in various tissues. Transcripts of ERFs from group I and VIII were very abundant in all tissues while those of group VII were highly accumulated in latex cells. Seven of the thirty-five ERF expression marker genes were highly expressed in latex. Subcellular localization and transactivation analyses suggested that HbERF-VII candidate genes encoded functional transcription factors.

Impact of the timing and duration of weed control on the establishment of a rubber tree plantation

Rubber tree production is reduced by weeds that compete for environmental resources; therefore, the timing and duration of weed control influences weed interference. The objectives of this study were to evaluate the growth of rubber tree (Hevea brasiliensis) plants, to determine the critical period for weed control, and to evaluate the growth recovery of rubber trees that coexisted with weeds for different periods of time after planting. Two groups of treatments were established under field conditions in the first year of the investigation: one group contained crescent periods of weed infestation, while the other contained crescent periods of weed control, also including a weed-free check and a total weedy check. In the second year of the investigation, the weeds were totally controlled. Urochloa decumbens was the dominant weed (over 90% groundcover). Crop growth was greatly reduced due to the weed interference. Plant height decreased more rapidly than did any other characteristic. Plant height, leaf dry mass, and leaf area decreased by 99%, 97% and 96%, respectively, and were the most reduced characteristics. Plant height also recovered more rapidly than did any characteristic when the period of weed control was lengthened. However, stem dry mass increased by 750%, making it the most recovered characteristic. The critical period for weed control was between 4 and 9½ months after planting in the first year; however, the rubber trees showed an expressive growth recovery when the weeds were controlled throughout the second year.

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.

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.

Over-expression of a cytosolic isoform of the HbCuZnSOD gene in Hevea brasiliensis changes its response to a water deficit

Hevea brasiliensis is the main commercial source of natural rubber. Reactive oxygen species (ROS) scavenging systems are involved in various biotic and abiotic stresses. Genetic engineering was undertaken to study the strengthening of plant defences by antioxidants. To that end, Hevea transgenic plant lines over-expressing a Hevea brasiliensis cytosolic HbCuZnSOD gene were successfully established and regenerated. Over-expression of the HbCuZnSOD gene was not clearly related to an increase in SOD activity in plant leaves. The impact of HbCuZnSOD gene over-expression in somatic embryogenesis and in plant development are presented and discussed. The water deficit tolerance of two HbCuZnSOD over-expressing lines was evaluated. The physiological parameters of transgenic plantlets subjected to a water deficit suggested that plants from line TS4T8An displayed lower stomatal conductance and a higher proline content. Over-expression of the HbCuZnSOD gene and activation of all ROS-scavenging enzymes also suggested that protection against ROS was more efficient in the TS4T8An transgenic line.

Development and characterization of intron-flanking EST-PCR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

With a long-term goal of constructing a linkage map enriched with gene-specific markers in rubber tree (Hevea brasiliensis Muell. Arg.), we utilized rubber tree ESTs associated with tapping panel dryness (TPD) to develop intron-flanking PCR markers. After downloading and assembling the rubber tree ESTs associated with TPD, we predicted the exon/exon junction sites (E/E) by aligning rubber tree transcripts with the genomic sequences of castor bean (Ricinus communis L.). Based on the predicted E/E, the primers flanking intron(s) and no intron were designed. Compared with the markers designed by conventional method, the PCR success rate of the markers designed with the predicted E/E increased 28-30%, whereas the polymorphism rate of intron-flanking EST-PCR markers was approximately 3.43-fold increase. Therefore, the intron-flanking marker was more polymorphism-generating efficient than the markers designed by conventional methods. In addition, analyzing the polymorphic information content (PIC) among Hevea germplasm showed that the polymorphism of wild rubber tree accessions was higher than one of cultivated rubber tree clones and Hevea species. This study enriches the categories and numbers of molecular markers in rubber tree, and the markers developed in this research will have a wide application in DNA fingerprinting, marker-assisted selection and genetic mapping in rubber tree. This research also indicates that it is possible to develop intron-flanking EST-PCR markers of rubber tree with castor bean genome as reference sequences, which provides new insights into developing intron-flanking EST-PCR markers for rubber tree or other plant species without genomic information.

De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

Background

In rubber tree, bark is one of important agricultural and biological organs. However, the molecular mechanism involved in the bark formation and development in rubber tree remains largely unknown, which is at least partially due to lack of bark transcriptomic and genomic information. Therefore, it is necessary to carried out high-throughput transcriptome sequencing of rubber tree bark to generate enormous transcript sequences for the functional characterization and molecular marker development.

Results

In this study, more than 30 million sequencing reads were generated using Illumina paired-end sequencing technology. In total, 22,756 unigenes with an average length of 485 bp were obtained with de novo assembly. The similarity search indicated that 16,520 and 12,558 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 6,867 and 5,559 unigenes were separately assigned to Gene Ontology (GO) and Clusters of Orthologous Group (COG). When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 12,097 unigenes were assigned to 5 main categories including 123 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (9,043, 74.75%), suggesting the active metabolic processes in rubber tree bark. In addition, a total of 39,257 EST-SSRs were identified from 22,756 unigenes, and the characterizations of EST-SSRs were further analyzed in rubber tree. 110 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among 13 Hevea germplasms, PCR success rate and polymorphism rate of 110 markers were separately 96.36% and 55.45% in this study.

Conclusion

By assembling and analyzing de novo transcriptome sequencing data, we reported the comprehensive functional characterization of rubber tree bark. This research generated a substantial fraction of rubber tree transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation, and microarrays development in rubber tree. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding in rubber tree. Moreover, this study also supported that transcriptome analysis based on Illumina paired-end sequencing is a powerful tool for transcriptome characterization and molecular marker development in non-model species, especially those with large and complex genomes.

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.

HbMyb1, a Myb transcription factor from Hevea brasiliensis, suppresses stress induced cell death in transgenic tobacco

Tapping panel dryness (TPD) is a complex physiological syndrome found widely in rubber tree (Hevea brasiliensis) plantations that causes severe yield loss in natural rubber-producing countries. In an earlier study, we confirmed that there is a negative correlation between HbMyb1 expression and TPD severity. To further investigate the function of HbMyb1 in TPD, HbMyb1 was over-expressed in tobacco controlled by a CaMV 35S promoter. In transgenic plants expressing HbMyb1, cell death induced by UV-B irradiation, paraquat and the hypersensitive reaction to necrotrophic fungal infection (Botrytis cinerea) was suppressed with a close correlation between HbMyb1 protein levels and the extent of suppression. In addition the nuclear condensation and degradation were observed in laticifer cells of TPD trees, while the nucleus of laticifer cells of healthy trees was morphologically normal. On the basis of the results described above, we propose that HbMyb1 maybe suppress stress induced cell death in rubber trees.

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.

HbMT2, an ethephon-induced metallothionein gene from Hevea brasiliensis responds to H(2)O(2) stress

Metallothioneins (MTs) are the cysteine-rich proteins with low molecular weight, which play important roles in maintaining intracellular ion homeostasis, detoxification of heavy metal ions and protecting against intracellular oxidative damages. In this study a novel ethephon-induced metallothionein gene, designated as HbMT2, was isolated and characterized from Hevea brasiliensis. The HbMT2 cDNA contained a 237 bp open reading frame encoding 78 amino acids and the deduced protein showed high similarity to the type 2 MTs from other plant species. Expression analysis revealed more significant accumulation of HbMT2 transcripts in leaves and latex than in roots and barks. The transcription of HbMT2 in latex was strongly induced by ethephon and hydrogen peroxide (H(2)O(2)) stress. Overproduction of recombinant HbMT2 protein gave the Escherichia coli cells more tolerance on Cu(2+) and Zn(2+), and the recombinant HbMT2 could scavenge the reactive oxidant species (ROS) in vitro. All these results indicated that HbMT2 could respond to ethephon stimulation and H(2)O(2) stress as a ROS scavenger in H. brasiliensis. It is also suggested that HbMT2 function in improving the tolerance of rubber trees to heavy metal ions, and repressing the ethephon-induced senilism and tapping panel dryness (TPD) development by ROS scavenge system in H. brasiliensis.

Identification and characterization of genes associated with tapping panel dryness from Hevea brasiliensis latex using suppression subtractive hybridization

BACKGROUND

Tapping panel dryness (TPD) is one of the most serious threats to natural rubber production. Although a great deal of effort has been made to study TPD in rubber tree, the molecular mechanisms underlying TPD remain poorly understood. Identification and systematical analyses of the genes associated with TPD are the prerequisites for elucidating the molecular mechanisms involved in TPD. The present study is undertaken to generate information about the genes related to TPD in rubber tree.

RESULTS

To identify the genes related to TPD in rubber tree, forward and reverse cDNA libraries from the latex of healthy and TPD trees were constructed using suppression subtractive hybridization (SSH) method. Among the 1106 clones obtained from the two cDNA libraries, 822 clones showed differential expression in two libraries by reverse Northern blot analyses. Sequence analyses indicated that the 822 clones represented 237 unique genes; and most of them have not been reported to be associated with TPD in rubber tree. The expression patterns of 20 differentially expressed genes were further investigated to validate the SSH data by reverse transcription PCR (RT-PCR) and real-time PCR analysis. According to the Gene Ontology convention, 237 unique genes were classified into 10 functional groups, such as stress/defense response, protein metabolism, transcription and post-transcription, rubber biosynthesis, etc. Among the genes with known function, the genes preferentially expressed were associated with stress/defense response in the reverse library, whereas metabolism and energy in the forward one.

CONCLUSIONS

The genes associated with TPD were identified by SSH method in this research. Systematic analyses of the genes related to TPD suggest that the production and scavenging of reactive oxygen species (ROS), ubiquitin proteasome pathway, programmed cell death and rubber biosynthesis might play important roles in TPD. Therefore, our results not only enrich information about the genes related to TPD, but also provide new insights into understanding the TPD process in rubber tree.