HbSnRK2.6 Functions in ABA-Regulated Cold Stress Response by Promoting HbICE2 Transcriptional Activity in Hevea brasiliensis.

Low temperature remarkably limits rubber tree (Hevea brasiliensis Muell. Arg.) growth, latex production, and geographical distribution, but the underlying mechanisms of Hevea brasiliensis cold stress response remain elusive. Here, we identified HbSnRK2.6 as a key component in ABA signaling functions in phytohormone abscisic acid (ABA)-regulated cold stress response in Hevea brasiliensis. Exogenous application of ABA enhances Hevea brasiliensis cold tolerance. Cold-regulated (COR) genes in the CBF pathway are upregulated by ABA. Transcript levels of all five HbSnRK2.6 members are significantly induced by cold, while HbSnRK2.6A, HbSnRK2.6B, and HbSnRK2.6C can be further activated by ABA under cold conditions. Additionally, HbSnRK2.6s are localized in the cytoplasm and nucleus, and can physically interact with HbICE2, a crucial positive regulator in the cold signaling pathway. Overexpression of HbSnRK2.6A or HbSnRK2.6B in Arabidopsis extensively enhances plant responses to ABA and expression of COR genes, leading to increased cold stress tolerance. Furthermore, HbSnRK2.6A and HbSnRK2.6B can promote transcriptional activity of HbICE2, thus, increasing the expression of HbCBF1. Taken together, we demonstrate that HbSnRK2.6s are involved in ABA-regulated cold stress response in Hevea brasiliensis by regulating transcriptional activity of HbICE2.

HbWRKY40 plays an important role in the regulation of pathogen resistance in Hevea brasiliensis.

KEY MESSAGE: Overexpression of HbWRKY40 induces ROS burst in tobacco and increases disease resistance in Arabidopsis; RNA-seq and ChIP assays revealed the regulatory network of HbWRKY40 in plant defense. WRKY, a family of plant transcription factors, are involved in the regulation of numerous biological processes. In rubber tree Hevea brasiliensis, the roles of WRKYs remain poorly understood. In the present study, a total of 111 genes encoding putative HbWRKY proteins were identified in the H. brasiliensis genome. Among these genes, HbWRKY40 transcripts were significantly induced by Colletotrichum gloeosporioides and salicylic acid. To assess its roles in plant defense, HbWRKY40 was over-expressed in Nicotiana benthamiana and Arabidopsis thaliana. The results showed that HbWRKY40 significantly induced reactive oxygen species burst in N. benthamiana and increased resistance of Arabidopsis against Botrytis cinerea. Transient expression in mesophyll cell protoplasts of H. brasiliensis showed that HbWRKY40 localizes at nuclei. In addition, transcripts of 145 genes were significantly up-regulated and 6 genes were down-regulated in the protoplasts over-expressing HbWRKY40 based on the RNA-seq analysis. Among these potential downstream targets, 12 genes contain potential WRKY-binding sites at the promoter regions. Further analysis through chromatin immunoprecipitation revealed that 10 of these 12 genes were the downstream targets of HbWRKY40. Taken together, our findings indicate that HbWRKY40 plays an important role in the disease resistance by regulating defense-associated genes in H. brasiliensis.

The ICE-like transcription factor HbICE2 is involved in jasmonate-regulated cold tolerance in the rubber tree (Hevea brasiliensis).

KEY MESSAGE: An ICE-like transcription factor mediates jasmonate-regulated cold tolerance in the rubber tree (Hevea brasiliensis), and confers cold tolerance in transgenic Arabidopsis. The rubber tree (Hevea brasiliensis) is susceptible to low temperatures, and understanding the mechanisms regulating cold stress is of great potential value for enhancing tolerance to this environmental variable. In this study, we find that treatment with exogenous methyl jasmonate (MeJA) could significantly enhance Hevea brasiliensis cold tolerance. In addition, yeast two-hybrid and bimolecular fluorescence complementation (BiFC) experiments show that JASMONATE ZIM-DOMAIN(JAZ) proteins, HbJAZ1 and HbJAZ12, key repressors of JA signaling pathway, interact with HbICE2, a novel ICE (Inducer of CBF Expression)-like protein. HbICE2 was nuclear-localised and bound to the MYC recognition (MYCR) sequence. The transcriptional activation activity of HbICE2 in yeast cells was dependent on the N-terminus, and overexpression of HbICE2 in Arabidopsis resulted in elevated tolerance to chilling stress. Furthermore, dual-luciferase transient assay reveals that HbJAZ1 and HbJAZ12 proteins inhibit the transcriptional function of HbICE2. The expression of C-repeat-binding factor (CBF) signalling pathway genes including HbCBF1, HbCBF2 and HbCOR47 were up-regulated by MeJA. Taken together, our data suggest that the new ICE-like transcription factor HbICE2 is involved in jasmonate-regulated cold tolerance in Hevea brasiliensis.

The relationship between latex metabolism gene expression with rubber yield and related traits in Hevea brasiliensis.

BACKGROUND: Expression patterns of many laticifer-specific gens are closely correlative with rubber yield of Hevea brasiliensis (para rubber tree). To unveil the mechanisms underlying the rubber yield, transcript levels of nine major latex metabolism-related genes, i.e., HMG-CoA synthase (HMGS), HMG-CoA reductase (HMGR), diphosphomevalonate decarboxylase (PMD), farnesyl diphosphate synthase (FPS), cis-prenyltransferase (CPT), rubber elongation factor (REF), small rubber particle protein (SRPP), dihydroxyacid dehydratase (DHAD) and actin depolymerizing factor (ADF), were dertermined, and the relationship between rubber yield with their expression levels was analysed. RESULTS: Except HbHMGR1, HbPMD and HbDHAD, most of these genes were predominantly expressed in latex, and bark tapping markedly elevated the transcript abundance of the analyzed genes, with the 7th tapping producing the greatest expression levels. Both ethephon (ETH) and methyl jasmonate (MeJA) stimulation greatly induced the expression levels of the examined genes, at least at one time point, except HbDHAD, which was unresponsive to MeJA. The genes' expression levels, as well as the rubber yields and two yield characteristics differed significantly among the different genotypes examined. Additionally, the latex and dry rubber yields increased gradually but the dry rubber content did not. Rubber yields and/or yield characteristics were significantly positively correlated with HbCPT, HbFPS, HbHMGS, HbHMGR1 and HbDHAD expression levels, negatively correlated with that of HbREF, but not significantly correlated with HbPMD, HbSRPP and HbADF expression levels. In addition, during rubber production, significantly positive correlations existed between the expression level of HbPMD and the levels of HbREF and HbHMGR1, between HbSRPP and the levels of HbHMGS and HbHMGR1, and between HbADF and HbFPS. CONCLUSIONS: The up-regulation of these genes might be related to the latex production of rubber trees under the stress of bark tapping and latex metabolism. The various correlations among the genes implied that there are differences in their synergic interactions. Thus, these nine genes might be related to rubber yield and yield-related traits in H. brasiliensis, and this work increases our understanding of their complex functions and how they are expressed in both high-and medium-yield rubber tree varieties and low-yield wild rubber tree germplasm.

Calcium-Dependent Protein Kinase Family Genes Involved in Ethylene-Induced Natural Rubber Production in Different Hevea brasiliensis Cultivars.

Natural rubber latex production can be improved by ethylene stimulation in the rubber tree (Hevea brasiliensis). However, the expression levels of most functional proteins for natural rubber biosynthesis are not induced after ethylene application, indicating that post-translational modifications, especially protein phosphorylation, may play important roles in ethylene signaling in Hevea. Here, we performed a comprehensive investigation on evolution, ethylene-induced expression and protein-protein interaction of calcium-dependent protein kinases (CPKs), an important serine/threonine protein kinase family, in Hevea. Nine duplication events were determined in the 30 identified HbCPK genes. Expression profiling of HbCPKs in three rubber tree cultivars with low, medium and high ethylene sensitivity showed that HbCPK6, 17, 20, 22, 24, 28 and 30 are induced by ethylene in at least one cultivar. Evolution rate analysis suggested accelerated evolution rates in two paralogue pairs, HbCPK9/18 and HbCPK19/20. Analysis of proteomic data for rubber latex after ethylene treatment showed that seven HbCPK proteins could be detected, including six ethylene-induced ones. Protein-protein interaction analysis of the 493 different abundant proteins revealed that protein kinases, especially calcium-dependent protein kinases, possess most key nodes of the interaction network, indicating that protein kinase and protein phosphorylation play important roles in ethylene signaling in latex of Hevea. In summary, our data revealed the expression patterns of HbCPK family members and functional divergence of two HbCPK paralogue pairs, as well as the potential important roles of HbCPKs in ethylene-induced rubber production improvement in Hevea.

Salicylic Acid Induces Resistance in Rubber Tree against Phytophthora palmivora.

Induced resistance by elicitors is considered to be an eco-friendly strategy to stimulate plant defense against pathogen attack. In this study, we elucidated the effect of salicylic acid (SA) on induced resistance in rubber tree against Phytophthora palmivora and evaluated the possible defense mechanisms that were involved. For SA pretreatment, rubber tree exhibited a significant reduction in disease severity by 41%. Consistent with the occurrence of induced resistance, the pronounced increase in H2O2 level, catalase (CAT) and peroxidase (POD) activities were observed. For defense reactions, exogenous SA promoted the increases of H2O2, CAT, POD and phenylalanine ammonia lyase (PAL) activities, including lignin, endogenous SA and scopoletin (Scp) contents. However, SA had different effects on the activity of each CAT isoform in the particular rubber tree organs. Besides, three partial cDNAs encoding CAT (HbCAT1, HbCAT2 and HbCAT3) and a partial cDNA encoding PAL (HbPAL) were isolated from rubber tree. Moreover, the expressions of HbCAT1, HbPAL and HbPR1 were induced by SA. Our findings suggested that, upon SA priming, the elevated H2O2, CAT, POD and PAL activities, lignin, endogenous SA and Scp contents, including the up-regulated HbCAT1, HbPAL and HbPR1 expressions could potentiate the resistance in rubber tree against P. palmivora.

Comparative Proteomics of Rubber Latex Revealed Multiple Protein Species of REF/SRPP Family Respond Diversely to Ethylene Stimulation among Different Rubber Tree Clones.

Rubber elongation factor (REF) and small rubber particle protein (SRPP) are two key factors for natural rubber biosynthesis. To further understand the roles of these proteins in rubber formation, six different genes for latex abundant REF or SRPP proteins, including REF138,175,258 and SRPP117,204,243, were characterized from Hevea brasiliensis Reyan (RY) 7-33-97. Sequence analysis showed that REFs have a variable and long N-terminal, whereas SRPPs have a variable and long C-terminal beyond the REF domain, and REF258 has a ß subunit of ATPase in its N-terminal. Through two-dimensional electrophoresis (2-DE), each REF/SRPP protein was separated into multiple protein spots on 2-DE gels, indicating they have multiple protein species. The abundance of REF/SRPP proteins was compared between ethylene and control treatments or among rubber tree clones with different levels of latex productivity by analyzing 2-DE gels. The total abundance of each REF/SRPP protein decreased or changed a little upon ethylene stimulation, whereas the abundance of multiple protein species of the same REF/SRPP changed diversely. Among the three rubber tree clones, the abundance of the protein species also differed significantly. Especially, two protein species of REF175 or REF258 were ethylene-responsive only in the high latex productivity clone RY 8-79 instead of in RY 7-33-97 and PR 107. Some individual protein species were positively related to ethylene stimulation and latex productivity. These results suggested that the specific protein species could be more important than others for rubber production and post-translational modifications might play important roles in rubber biosynthesis.

Genome-wide identification of rubber tree (Hevea brasiliensis Muell. Arg.) aquaporin genes and their response to ethephon stimulation in the laticifer, a rubber-producing tissue.

BACKGROUND: Natural rubber, an important industrial raw material, is specifically synthesized in laticifers located inside the rubber tree (Hevea brasiliensis Muell. Arg.) trunk. Due to the absence of plasmodesmata, the laticifer water balance is mediated by aquaporins (AQPs). However, to date, the characterization of H. brasiliensis AQPs (HbAQPs) is still in its infancy. RESULTS: In this study, 51 full-length AQP genes were identified from the rubber tree genome. The phylogenetic analysis assigned these AQPs to five subfamilies, including 15 plasma membrane intrinsic proteins (PIPs), 17 tonoplast intrinsic proteins (TIPs), 9 NOD26-like intrinsic proteins (NIPs), 4 small basic intrinsic proteins (SIPs) and 6 X intrinsic proteins (XIPs). Functional prediction based on the analysis of the aromatic/arginine (ar/R) selectivity filter, Froger's positions and specificity-determining positions (SDPs) showed a remarkable difference in substrate specificity among subfamilies. Homology analysis supported the expression of 44 HbAQP genes in at least one of the examined tissues. Furthermore, deep sequencing of the laticifer transcriptome in the form of latex revealed a key role of several PIP subfamily members in the laticifer water balance, and qRT-PCR analysis showed diverse expression patterns of laticifer-expressed HbAQP genes upon ethephon treatment, a widely-used practice for the stimulation of latex yield. CONCLUSIONS: This study provides an important genetic resource of HbAQP genes, which will be useful to improve the water use efficiency and latex yield of Hevea.

Molecular characterization of HbCZF1, a Hevea brasiliensis CCCH-type zinc finger protein that regulates hmg1.

KEY MESSAGE: The HbCZF1 protein binds to the hmg1 promoter in yeast and this interaction was confirmed in vitro. The hmg1 promoter was activated in transgenic plants by HbCZF1. Biosynthesis of natural rubber is known to be based on the mevalonate pathway in Hevea brasiliensis. The final step in the mevalonate production is catalyzed by the branch point enzyme, 3-hydroxy-3-methyl-glutaryl coenzyme A reductase (HMGR), which shunts HMG-CoA into the isoprenoid pathway, leading to the synthesis of natural rubber. However, molecular regulation of HMGR expression is not known. To study the transcriptional regulation of HMGR, the yeast one-hybrid experiment was performed to screen the latex cDNA library using the hmg1 (one of the three HMGR in H. brasiliensis) promoter as bait. One cDNA that encodes the CCCH-type zinc finger protein, designated as HbCZF1, was isolated from H. brasiliensis. HbCZF1 interacted with the hmg1 promoter in yeast one-hybrid system and in vitro. HbCZF1 contains a 1110 bp open reading frame that encodes 369 amino acids. The deduced HbCZF1 protein was predicted to possess a typical C-X7-C-X5-C3-H CCCH motif and RNA recognition motif. HbCZF1 was predominant in the latex, but little expression was detected in the leaves, barks, and roots. Furthermore, in transgenic tobacco plants, over-expression of HbCZF1 highly activated the hmg1 promoter. These results suggested that HbCZF1 may participate in the regulation of natural rubber biosynthesis in H. brasiliensis.

Characterization of HbEREBP1, a wound-responsive transcription factor gene in laticifers of Hevea brasiliensis Muell. Arg.

AP2/ERF transcription factors play an important role in regulation of the cross-talk between ethylene and jasmonate signaling pathways mediating defense responses of plants to biotic and abiotic stresses. In this study, an AP2/ERF transcription factor gene was isolated and characterized from laticifers of rubber tree by using RACE and real time PCR. The full length cDNA, referred to as HbEREBP1, was 1,095 bp in length and contained a 732 bp open reading frame encoding a putative protein of 243 amino acid residues. The molecular mass of the putative protein is 26.4 kDa with a pI of 9.46. The deduced amino acid sequence had a specific domain of AP2 superfamily and an ethylene-responsive element binding factor-associated amphiphilic repression motif, sharing 42.4, 39.1, and 38.0% identity with that of AtERF11, AtERF4, and AtERF8 in Arabidopsis, respectively. HbEREBP1 expression was down-regulated by tapping and mechanical wounding in the laticifers of adult trees. It was also down-regulated at early stage while up-regulated at late stage upon treatment with exogenous ethephon or methyl jasmonate, which was reverse to the case of defense genes in laticifers of epicormic shoots of rubber tree. Our results suggest that HbEREBP1 may be a negative regulator of defense genes in laticifers.

Molecular characterization of a novel 14-3-3 protein gene (Hb14-3-3c) from Hevea brasiliensis.

The cDNA encoding a 14-3-3 protein, designated as Hb14-3-3c, was isolated from Hevea brasiliensis. Hb14-3-3c was 1,269 bp long containing a 795 bp open reading frame encoding a putative protein of 264 amino acids, flanked by a 146 bp 5'UTR and a 328 bp 3' UTR. The predicted molecular mass of Hb14-3-3c is 29.67 kDa, with an isoelectric point of 4.52 and the deduced protein showed high similarity to the 14-3-3 protein from other plant species. Expression analysis revealed more significant accumulation of Hb14-3-3c transcripts in latex than in leaves, buds and flowers. The transcription of Hb14-3-3c in latex was induced by jasmonate and ethephon. Overproduction of recombinant Hb14-3-3c protein gave the Escherichia coli cells more tolerance on Co(2+), Cu(2+) and Zn(2+). Through yeast two-hybrid screening, 11 interaction partners of the Hb14-3-3c, which are involved in rubber biosynthesis, stress-related responses, defence etc., were identified in rubber tree latex. Taking these data together, it is proposed that the Hb14-3-3c may participate in regulation of rubber biosynthesis. Thus, the results of this study provide novel insights into the 14-3-3 signaling related to rubber biosynthesis, stress-related responses in rubber tree.

Three MADS-box genes of Hevea brasiliensis expressed during somatic embryogenesis and in the laticifer cells.

Three MADS-box genes, designated HbMADS1, HbMADS2 and HbMADS3, were isolated from Hevea brasiliensis. HbMADS1, HbMADS2 and HbMADS3 encode polypetides consisting of 245, 217 and 239 amino acids, respectively, containing conserved MADS-box motifs at N-terminus. Transcription pattern analysis revealed that three MADS-box genes had highly transcription in the laticifer cells. The transcriptions of HbMADS1and HbMADS3 were induced in the laticifer cells by jamonic acid, while HbMADS2 was not induction by jamonic acid. Ethephone is not effective in inducing their expression. The three genes were differentially expressed during somatic embryogenesis of rubber tree. Characterization of HbMADSs will attribute to understand their possible function in rubber tree.

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.

Ethylene stimulation of latex yield depends on the expression of a sucrose transporter (HbSUT1B) in rubber tree (Hevea brasiliensis).

Hevea brasiliensis is an important industrial crop for natural rubber production. Latex biosynthesis occurs in the cytoplasm of highly specialized latex cells and requires sucrose as the unique precursor. Ethylene stimulation of latex production results in high sugar flow from the surrounding cells of inner bark towards the latex cells. The aim of this work was to understand the role of seven sucrose transporters (HbSUTs) and one hexose transporter (HbHXT1) in this process. Two Hevea clones were used: PB217 and PB260, respectively described as high and low yielding clones. The expression pattern of these sugar transporters (HbSUTs and HbHXT1) was monitored under different physiological conditions and found to be maximal in latex cells. HbSUT1, one of the most abundant isoforms, displayed the greatest response to ethylene treatment. In clone PB217, ethylene treatment led to a higher accumulation of HbSUT1B in latex cells than in the inner bark tissues. Conversely, stronger expression of HbSUT1B was observed in inner bark tissues than in latex cells of PB260. A positive correlation with HbSUT1B transcript accumulation and increased latex production was further supported by its lower expression in latex cells of the virgin clone PB217.

In vitro transient expression system of latex C-serum was used for analysis of hevein promoter in response to abscisic acid in Hevea brasiliensis.

Hevein has been found to be an essential element in coagulation of rubber particles in latex of rubber trees. In a previous study, we cloned a 1,241-bp fragment of a 5' upstream region of the hevein gene by genome walking. This fragment was analyzed by a 5' end nested deletion method in the present study, fused with a uidA (gus) gene to produce a series of tested constructs, which were transferred into C-serum of latex and the Gus activities were detected. Results showed that the fragment from -749 to -292 was sufficient for expression of gus gene in latex, and the fragment from -292 to -168 was crucial in response to abscisic acid inducement. In a transient transgenic test of rubber leaf with particle bombardment, construct Hev749 conferred gus-specific expression in veins, in which the latex tubes mainly distributed. This implies that the fragment from -749 to -292 was laticiferous-specific.

Purification and characterization of cassiicolin, the toxin produced by Corynespora cassiicola, causal agent of the leaf fall disease of rubber tree.

Cassiicolin, a phytotoxin produced by the necrotrophic fungus Corynespora cassiicola, was purified to homogeneity from a rubber tree isolate. The optimized protocol involves reverse phase chromatography followed by size exclusion chromatography, with monitoring of the toxicity on detached rubber tree leaves. Cassiicolin appeared to be a peptide composed of 27 amino acids, glycosylated on the second residue, with a N-terminal pyroglutamic acid and 6 cysteines involved in disulfide bonds. Its molecular mass was estimated to be 2885 Da. No significant sequence homology with other proteins could be found. The availability of pure toxin in sufficient amount is a prerequisite for its structure determination, which is a key step in the understanding of the aggression mechanism.

Ethrel-stimulated prolongation of latex flow in the rubber tree (Hevea brasiliensis Muell. Arg.): an Hev b 7-like protein acts as a universal antagonist of rubber particle aggregating factors from lutoids and C-serum.

Ethrel is the most effective stimuli in prolonging the latex flow that consequently increases yield per tapping. This effect is largely ascribed to the enhanced lutoid stability, which is associated with the decreased release of initiators of rubber particle (RP) aggregation from lutoid bursting. However, the increase in both the bursting index of lutoids and the duration of latex flow after applying ethrel or ethylene gas in high concentrations suggests that a new mechanism needs to be introduced. In this study, a latex allergen Hev b 7-like protein in C-serum was identified by matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI TOF MS). In vitro analysis showed that the protein acted as a universal antagonist of RP aggregating factors from lutoids and C-serum. Ethrel treatment obviously weakened the effect of C-serum on RP aggregation, which was closely associated with the increase in the level of the Hev b 7-like protein and the decrease in the level of the 37 kDa protein, as revealed by sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE), western blotting analysis and antibody neutralization. Thus, the increase of the Hev b 7-like protein level or the ratio of the Hev b 7-like protein to the 37 kDa protein in C-serum should be primarily ascribed to the ethrel-stimulated prolongation of latex flow duration.

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

The secondary laticifer in rubber tree (Hevea brasiliensis Muell. Arg.) is a specific tissue within the secondary phloem. This tissue differentiates from the vascular cambia, and its function is natural rubber biosynthesis and storage. Given that jasmonates play a pivotal role in secondary laticifer differentiation, we established an experimental system with jasmonate (JA) mimic coronatine (COR) for studying the secondary laticifer differentiation: in this system, differentiation occurs within five days of the treatment of epicormic shoots with COR. In the present study, the experimental system was used to perform transcriptome sequencing and gene expression analysis. A total of 67,873 unigenes were assembled, and 50,548 unigenes were mapped at least in one public database. Of these being annotated unigenes, 15,780 unigenes were differentially expressed early after COR treatment, and 19,824 unigenes were differentially expressed late after COR treatment. At the early stage, 8,646 unigenes were up-regulated, while 7,134 unigenes were down-regulated. At the late stage, the numbers of up- and down-regulated unigenes were 7,711 and 12,113, respectively. The annotation data and gene expression analysis of the differentially expressed unigenes suggest that JA-mediated signalling, Ca2+ signal transduction and the CLAVATA-MAPK-WOX signalling pathway may be involved in regulating secondary laticifer differentiation in rubber trees.

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

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

Identification of differentially expressed genes and signalling pathways in bark of Hevea brasiliensis seedlings associated with secondary laticifer differentiation using gene expression microarray.

The natural rubber of Para rubber tree, Hevea brasiliensis, is the main crop involved in industrial rubber production due to its superior quality. The Hevea bark is commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. The laticifer is well defined in the aspect of morphology; however, only some genes associated with its development have been reported. We successfully induced secondary laticifer in the jasmonic acid (JA)-treated and linolenic acid (LA)-treated Hevea bark but secondary laticifer is not observed in the ethephon (ET)-treated and untreated Hevea bark. In this study, we analysed 27,195 gene models using NimbleGen microarrays based on the Hevea draft genome. 491 filtered differentially expressed (FDE) transcripts that are common to both JA- and LA-treated bark samples but not ET-treated bark samples were identified. In the Eukaryotic Orthologous Group (KOG) analysis, 491 FDE transcripts belong to different functional categories that reflect the diverse processes and pathways involved in laticifer differentiation. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and KOG analysis, the profile of the FDE transcripts suggest that JA- and LA-treated bark samples have a sufficient molecular basis for secondary laticifer differentiation, especially regarding secondary metabolites metabolism. FDE genes in this category are from the cytochrome (CYP) P450 family, ATP-binding cassette (ABC) transporter family, short-chain dehydrogenase/reductase (SDR) family, or cinnamyl alcohol dehydrogenase (CAD) family. The data includes many genes involved in cell division, cell wall synthesis, and cell differentiation. The most abundant transcript in FDE list was SDR65C, reflecting its importance in laticifer differentiation. Using the Basic Local Alignment Search Tool (BLAST) as part of annotation and functional prediction, several characterised as well as uncharacterized transcription factors and genes were found in the dataset. Hence, the further characterization of these genes is necessary to unveil their role in laticifer differentiation. This study provides a platform for the further characterization and identification of the key genes involved in secondary laticifer differentiation.