Hormonal treatment of the bark of rubber trees (Hevea brasiliensis) increases latex yield through latex dilution in relation with the differential expression of two aquaporin genes

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.

Transcriptome analysis reveals the molecular mechanisms of rubber biosynthesis and laticifer differentiation during rubber seed germination

The molecular mechanisms underlying the initiation of natural rubber synthesis and laticifer differentiation have not been fully elucidated. In this study, we conducted a time-series transcriptome analysis of five rubber tree tissues at four stages of seed germination. A total of 161,199 DEGs were identified between the two groups, including most 16,673 DEGs (A3 vs B3 and A3 vs C3) and lest 1,210 DEGs (C2 vs D2). We found that the maturation of the seed is accompanied by the formation of laticifer cells in cotyledon. Meanwhile, the analysis of hormones related genes expression may provide effective clues for us to promote the differentiation of laticifer cells in seeds by hormones in the future. In this study, hormone-related gene enrichment analyses revealed that IAA, GA, and CTK were activated in laticifer containing tissues. Similarly, GO and GEGG analysis showed that hormone pathways, especially the auxin pathway, are enriched. Gene expression clustering was analyzed using the short time-series expression miner (STEM), and the analysis revealed four distinct trends in the gene expression profiles. Moreover, we enriched transcription factor (TF) enrichment in cotyledon and embryonic axis tissues, and the MYB type exhibited the most significant difference. Furthermore, our findings revealed that genes related to rubber synthesis exhibited tissue-specific expression patterns during seed germination. Notably, key genes associated with rubber biosynthesis, specifically small rubber particle protein (SRPP) and cis-prenyltransferase (CPT), exhibited significant changes in expression in cotyledon and embryonic axis tissues, suggesting synchronous rubber synthesis with seed germination. Our staining results reveled that laticifer cells were exits in the cotyledon before seed imbibition stage. In conclusion, these results lay the foundation for exploring the molecular mechanisms underlying laticifer differentiation and rubber synthesis during seed germination, deepening our understanding of the initiation stages of rubber biosynthesis and laticifer differentiation.

HbTGA1, a TGA Transcription Factor From Hevea brasiliensis, Regulates the Expression of Multiple Natural Rubber Biosynthesis Genes

The TGA transcription factors are known to modulate the biosynthesis of secondary metabolites in plants. However, their regulatory function in natural rubber (NR) biosynthesis was not revealed in the rubber tree (Hevea brasiliensis). Here, 14 genes encoding TGA transcription factors (name HbTGA1-HbTGA14) were identified in the rubber tree. HbTGAs were differentially expressed in different tissues. HbTGA1 was expressed at its highest level in latex. We found specific in vitro and in vivo binding of the HbTGA1 protein with promoters of multiple NR biosynthesis genes (HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2). The activation of the promoters of HbHMGS2 and HbCPT6 was significantly suppressed by HbTGA1, while the activities of promoters of HbHMGR2, HbCPT8, and HbSRPP2 were increased by HbTGA1. The promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were significantly increased by HbTGA1 under jasmonate stress, while the promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were also significantly increased by HbTGA1 under salicylic acid stress. The present study provides insights into the role of TGA transcription factors in regulating the expression of NR biosynthesis genes from H. brasiliensis.

Analysis of reduced and oxidized antioxidants in Hevea brasiliensis latex reveals new insights into the regulation of antioxidants in response to harvesting stress and tapping panel dryness

Latex diagnosis (LD) is applied to optimize the natural rubber production and prevent tapping panel dryness (TPD), a physiological syndrome affecting latex production in Hevea brasiliensis. The reduced thiol content (RSH) is one of the biochemical parameters associated with the risk of TPD. However, RSH is difficult to interpret because of the influence of the environment. In order to better understand the regulation of antioxidants and to better interpret RSH, a key parameter of LD, this study analysed in latex both oxidised and reduced forms of ascorbic acid (AsA) and glutathione, and their cofactors as well as other latex diagnosis parameters in response to harvesting stress (tapping and ethephon stimulation) and TPD occurrence. The content of antioxidants in latex had a high variability among five rubber clones. The concentration in AsA was about ten times higher than GSH in laticifer, GSH accounting for about 50% of RSH. For short-term harvesting stress, RSH increased with tapping frequency and ethephon stimulation. TPD is associated with high latex viscosity and bursting of lysosomal particles called lutoids, as well as for several rubber clones with lower RSH and GSH contents. These results suggest that a high level of RSH shows the capacity of laticifer metabolism to cope with harvesting stress, while a drop in RSH is the sign of long stress related to lower metabolic activity and TPD occurrence. RSH remains an essential physiological parameter to prevent TPD when associated with reference data under low and high harvesting stress. This study paves the way to understand the role of AsA and GSH, and carry out genetic studies of antioxidants.

Identification of Sieve Element Occlusion Gene (SEOs) Family in Rubber Trees (Hevea brasiliensis Muell. Arg.) Provides Insights to the Mechanism of Laticifer Plugging

P proteins encoded by SEOs (sieve element occlusion) have been shown to be associated with the blockage of sieve tubes after injury in many plants, but the presence of SEO genes and their association with rubber tree laticifer plugging and latex yield remain unclear. Through a systematic identification and analysis, seven SEO genes were identified from the rubber tree genome. The physicochemical properties of their proteins, gene structures, conserved domains, and locations on chromosomes were analyzed. According to their phylogenetic distance, HbSEOs were divided into two clusters. The transcript levels of HbSEO genes varied with tissues, in which HbSEO3 and HbSEO4 were most highly expressed in leaf, bark, and latex. HbSEOs could be induced by ethephon, methyl jasmonate, mechanical injury, and tapping; furthermore, they were highly expressed in trees with short flow duration, suggesting their possible association with rubber tree laticifer plugging and latex yield. To our knowledge, this is the first report of HbSEOs in rubber trees. It provides us with a better understanding of the mechanism of laticifer plugging.

Transcript and Protein Profiling Provides Insights Into the Molecular Mechanisms of Harvesting-Induced Latex Production in Rubber Tree

Natural rubber, an important industrial raw material with wide applications, is harvested in the form of latex (cytoplasm of rubber-producing laticifers) from Hevea brasiliensis (para rubber tree) by the way of tapping. Conspicuous stimulation on latex production is observed for the first few tappings conducted on virgin (untapped before) or resting (tapped before but no tapping for a period) rubber trees. To understand the underlying mechanisms, an integrative analysis of the latex transcriptome and proteome was conducted on virgin or resting Hevea trees for the first five tappings. A total of 505 non-redundant differentially expressed (DE) transcript-derived fragments (TDFs) were identified by silver-staining cDNA-AFLP, with 217 exhibiting patterns of upregulated, 180 downregulated and 108 irregularly-regulated. Meanwhile, 117 two dimensional gel electrophoresis DE-protein spots were isolated and subjected to mass spectrometry analysis, with 89 and 57 being successfully identified by MALDI-TOF and MALDI-TOF/TOF, respectively. About 72.5% DE-TDFs and 76.1% DE-proteins were functionally annotated and categorized. Noteworthily, most of the DE-TDFs implicated in sugar transport and metabolism as well as rubber biosynthesis were upregulated by the tapping treatment. The importance of sugar metabolism in harvesting-induced latex production was reinforced by the identification of abundant relevant DE-protein spots. About 83.8% of the randomly selected DE-TDFs were validated for expression patterns by semi-quantitative RT-PCR, and an 89.7% consistency for the 29 latex regeneration-related DE-TDFs examined by quantitative RT-PCR analysis. In brief, our results reveal extensive physiological and molecular changes in Hevea laticifers incurred by the tapping treatment, and the vast number of DE genes and proteins identified here contribute to unraveling the gene regulatory network of tapping-stimulated latex production.

Recent insights on gene expression studies on Hevea Brasiliensis fatal leaf fall diseases

Hevea brasiliensis is one of the most important agricultural commodities globally, heavily cultivated in Southeast Asia. Fatal leaf fall diseases cause aggressive leaf defoliation, linked to lower latex yield and death of crops before maturity. Due to the significant consequences of the disease to H. brasiliensis, the recent gene expression studies from four fall leaf diseases of H. brasiliensis were gathered; South American leaf blight, powdery mildew, Corynespora cassiicola and Phytophthora leaf fall disease. The differential analysis observed the pattern of commonly expressed genes upon fungi triggers using RT-PCR, DDRT-PCR, Real-time qRT-PCR and RNA-Seq. We have observed that RNA-Seq is the best tool to seek novel genes. Among the identified genes with defence-against fungi were pathogenesis-related genes such as β-1,3-glucanase and chitinase, the reactive oxygen species, and the phytoalexin biosynthesis. This manuscript also provided functional elaboration on the responsive genes and predicted possible biosynthetic pathways to identify and characterise novel genes in the future. At the end of the manuscript, the PCR methods and proteomic approaches were presented for future molecular and biochemical studies in the related diseases to H. brasiliensis.

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.

Transcriptome analysis of Pará rubber tree (H. brasiliensis) seedlings under ethylene stimulation

BACKGROUND

Natural rubber (cis-1,4-polyioprene, NR) is an indispensable industrial raw material obtained from the Pará rubber tree (H. brasiliensis). Natural rubber cannot be replaced by synthetic rubber compounds because of the superior resilience, elasticity, abrasion resistance, efficient heat dispersion, and impact resistance of NR. In NR production, latex is harvested by periodical tapping of the trunk bark. Ethylene enhances and prolongs latex flow and latex regeneration. Ethephon, which is an ethylene-releasing compound, applied to the trunk before tapping usually results in a 1.5- to 2-fold increase in latex yield. However, intense mechanical damage to bark tissues by excessive tapping and/or over-stimulation with ethephon induces severe oxidative stress in laticifer cells, which often causes tapping panel dryness (TPD) syndrome. To enhance NR production without causing TPD, an improved understanding of the molecular mechanism of the ethylene response in the Pará rubber tree is required. Therefore, we investigated gene expression in response to ethephon treatment using Pará rubber tree seedlings as a model system.

RESULTS

After ethephon treatment, 3270 genes showed significant differences in expression compared with the mock treatment. Genes associated with carotenoids, flavonoids, and abscisic acid biosynthesis were significantly upregulated by ethephon treatment, which might contribute to an increase in latex flow. Genes associated with secondary cell wall formation were downregulated, which might be because of the reduced sugar supply. Given that sucrose is an important molecule for NR production, a trade-off may arise between NR production and cell wall formation for plant growth and for wound healing at the tapping panel.

CONCLUSIONS

Dynamic changes in gene expression occur specifically in response to ethephon treatment. Certain genes identified may potentially contribute to latex production or TPD suppression. These data provide valuable information to understand the mechanism of ethylene stimulation, and will contribute to improved management practices and/or molecular breeding to attain higher yields of latex from Pará rubber trees.

The jasmonic acid pathway, rather than abscisic acid, may partly explain contrasting stomatal responses in two strawberry cultivars under osmotic stress

Drought is a major threat in agriculture and horticulture, including commercial strawberry production. Here, we compare hormonal regulation of a first-line drought stress response, namely stomatal closure, in two Fragaria x ananassa cultivars, known to differ in their drought stress phenotype. We show that the observed difference in xylem abscisic acid accumulation cannot explain the different stomatal responses under osmotic stress. Foliar abscisic acid accumulation cannot fully account for the stomatal behavior in one of both cultivars either. An indirect effect of abscisic acid on stomatal conductance via an impact on leaf hydraulic conductance, possibly mediated via aquaporins, as is recently proposed in literature, was not observed here. Next, we show that these two cultivars respond differently to jasmonic acid and one of its precursors. This difference in sensitivity of the jasmonates pathway between both cultivars may partly explain the different stomatal response. This study contributes to the understanding of the regulation of an important drought stress response in an economically important crop prone to water deficit stress.

Plant Aquaporins: Diversity, Evolution and Biotechnological Applications

The plasma membrane forms a permeable barrier that separates the cytoplasm from the external environment, defining the physical and chemical limits in each cell in all organisms. The movement of molecules and ions into and out of cells is controlled by the plasma membrane as a critical process for cell stability and survival, maintaining essential differences between the composition of the extracellular fluid and the cytosol. In this process aquaporins (AQPs) figure as important actors, comprising highly conserved membrane proteins that carry water, glycerol and other hydrophilic molecules through biomembranes, including the cell wall and membranes of cytoplasmic organelles. While mammals have 15 types of AQPs described so far (displaying 18 paralogs), a single plant species can present more than 120 isoforms, providing transport of different types of solutes. Such aquaporins may be present in the whole plant or can be associated with different tissues or situations, including biotic and especially abiotic stresses, such as drought, salinity or tolerance to soils rich in heavy metals, for instance. The present review addresses several aspects of plant aquaporins, from their structure, classification, and function, to in silico methodologies for their analysis and identification in transcriptomes and genomes. Aspects of evolution and diversification of AQPs (with a focus on plants) are approached for the first time with the aid of the LCA (Last Common Ancestor) analysis. Finally, the main practical applications involving the use of AQPs are discussed, including patents and future perspectives involving this important protein family.

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

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

Identification of the targets of HbEIN3/EILs in genomic wide in Hevea brasiliensis

EIN3/EILs are key regulators in ET signaling pathway. In this work, 4 members of EIN3/EILs of Hevea brasiliensis (HbEIN3/EILs) showed interaction with two F box proteins, HbEBF1 and HbEBF2. HbEIN3 located in nucleus and exhibited strong transcriptional activity. HbEIN3 was induced by ET treatment in C-serum, but not in B-serum of latex. HbEIN3/EILs bound to G-box cis-element. To globally search the potential targets of HbEIN3/EILs, genomic sequences of H. brasiliensis was re-annotated and an HCES (Hevea Cis-Elements Scanning) program was developed ( www.h-brasiliensis.com ). HCES scanning results showed that ET- and JA- responsive cis-elements distribute overlapping in gene promoters. 3146 genes containing G-box in promoters are potential targets of HbEIN3, including 41 genes involved in biosynthesis and drainage of latex, of which 7 rate-limiting genes of latex production were regulated by both ET and JA, suggesting that ET and JA signaling pathways coordinated the latex biosynthesis and drainage in H. brasiliensis. Abbreviations: ABRE: ABA responsive elements; bHLH: basic helix-loop-helix; COG: Orthologous Groups; DRE: dehydration response element; ERE: ethylene responsive element; ET: Ethylene; GO: Gene Ontology; HCES: Hevea Cis-Elements Scanning; JA: jasmonates; JRE: Jasmonate-responsive element; KEGG: Kyoto Encyclopedia of Genes and Genomes; NR: non-redundant database; PLACE: Plant Cis-acting Regulatory DNA Elements; qRT-PCR: quantitative real-time RT-PCR.

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.

Identification of JAZ-interacting MYC transcription factors involved in latex drainage in Hevea brasiliensis

Hevea brasiliensis Müll. Arg. is one of the most frequently wounded plants worldwide. Expelling latex upon mechanical injury is a wound response of rubber trees. However, JA-mediated wound responses in rubber trees are not well documented. In this work, three JAZ-interacting MYC transcription factors of H. brasiliensis (termed HbMYC2/3/4) were identified by yeast two-hybrid screening. HbMYC2/3/4 each showed specific interaction profiles with HbJAZs. HbMYC2/3/4 each localized in the nucleus and exhibited strong transcriptional activity. To identify the target genes potentially regulated by HbMYC2/3/4, cis-elements interacting with HbMYC2/3/4 were first screened by yeast one-hybrid assays; the results indicated that HbMYC2/3/4 each could bind G-box elements. Additional analysis confirmed that HbMYC2/3/4 bound the HbPIP2;1 promoter, which contains five G-box cis-elements, and regulated the expression of reporter genes in yeast cells and in planta. HbMYC2/3/4 were induced by exogenous JA treatment but suppressed by ethylene (ET) treatment; in contrast, HbPIP2;1 was positively regulated by ET but negatively regulated by JA treatment. Given that HbPIP2;1 is involved in latex drainage, it could be proposed that HbMYC2/3/4 are involved in the regulation of HbPIP2;1 expression as well as latex drainage, both of which are coordinated by the JA and ET signalling pathways.

Identification and characterization of the abscisic acid (ABA) receptor gene family and its expression in response to hormones in the rubber tree

Abscisic acid (ABA) is an essential phytohormone involved in diverse physiological processes. Although genome-wide analyses of the ABA receptor PYR/PYL/RCAR (PYL) protein/gene family have been performed in certain plant species, little is known about the ABA receptor protein/gene family in the rubber tree (Hevea brasiliensis). In this study, we identified 14 ABA receptor PYL proteins/genes (designated HbPYL1 through HbPYL14) in the most recent rubber tree genome. A phylogenetic tree was constructed, which demonstrated that HbPYLs can be divided into three subfamilies that correlate well with the corresponding Arabidopsis subfamilies. Eight HbPYLs are highly expressed in laticifers. Five of the eight genes are simultaneously regulated by ABA, jasmonic acid (JA) and ethylene (ET). The identification and characterization of HbPYLs should enable us to further understand the role of ABA signal in the rubber tree.

The Hevea brasiliensis XIP aquaporin subfamily: genomic, structural and functional characterizations with relevance to intensive latex harvesting

X-Intrinsic Proteins (XIP) were recently identified in a narrow range of plants as a full clade within the aquaporins. These channels reportedly facilitate the transport of a wide range of hydrophobic solutes. The functional roles of XIP in planta remain poorly identified. In this study, we found three XIP genes (HbXIP1;1, HbXIP2;1 and HbXIP3;1) in the Hevea brasiliensis genome. Comprehensive bioinformatics, biochemical and structural analyses were used to acquire a better understanding of this AQP subfamily. Phylogenetic analysis revealed that HbXIPs clustered into two major groups, each distributed in a specific lineage of the order Malpighiales. Tissue-specific expression profiles showed that only HbXIP2;1 was expressed in all the vegetative tissues tested (leaves, stem, bark, xylem and latex), suggesting that HbXIP2;1 could take part in a wide range of cellular processes. This is particularly relevant to the rubber-producing laticiferous system, where this isoform was found to be up-regulated during tapping and ethylene treatments. Furthermore, the XIP transcriptional pattern is significantly correlated to latex production level. Structural comparison with SoPIP2;1 from Spinacia oleracea species provides new insights into the possible role of structural checkpoints by which HbXIP2;1 ensures glycerol transfer across the membrane. From these results, we discuss the physiological involvement of glycerol and HbXIP2;1 in water homeostasis and carbon stream of challenged laticifers. The characterization of HbXIP2;1 during rubber tree tapping lends new insights into molecular and physiological response processes of laticifer metabolism in the context of latex exploitation.

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.

Quantitative and Functional Phosphoproteomic Analysis Reveals that Ethylene Regulates Water Transport via the C-Terminal Phosphorylation of Aquaporin PIP2;1 in Arabidopsis

Ethylene participates in the regulation of numerous cellular events and biological processes, including water loss, during leaf and flower petal wilting. The diverse ethylene responses may be regulated via dynamic interplays between protein phosphorylation/dephosphorylation and ubiquitin/26S proteasome-mediated protein degradation and protease cleavage. To address how ethylene alters protein phosphorylation through multi-furcated signaling pathways, we performed a (15)N stable isotope labelling-based, differential, and quantitative phosphoproteomics study on air- and ethylene-treated ethylene-insensitive Arabidopsis double loss-of-function mutant ein3-1/eil1-1. Among 535 non-redundant phosphopeptides identified, two and four phosphopeptides were up- and downregulated by ethylene, respectively. Ethylene-regulated phosphorylation of aquaporin PIP2;1 is positively correlated with the water flux rate and water loss in leaf. Genetic studies in combination with quantitative proteomics, immunoblot analysis, protoplast swelling/shrinking experiments, and leaf water loss assays on the transgenic plants expressing both the wild-type and S280A/S283A-mutated PIP2;1 in the both Col-0 and ein3eil1 genetic backgrounds suggest that ethylene increases water transport rate in Arabidopsis cells by enhancing S280/S283 phosphorylation at the C terminus of PIP2;1. Unknown kinase and/or phosphatase activities may participate in the initial up-regulation independent of the cellular functions of EIN3/EIL1. This finding contributes to our understanding of ethylene-regulated leaf wilting that is commonly observed during post-harvest storage of plant organs.

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.

Comparative proteomic analysis of latex from Hevea brasiliensis treated with Ethrel and methyl jasmonate using iTRAQ-coupled two-dimensional LC-MS/MS

Ethrel (ET) is an effective and widely used latex yield stimulant of Hevea brasiliensis (Pará rubber tree), and jasmonate (JA) is a key inducer of laticifer differentiation in this plant. To examine variations in the latex proteome caused by these phytohormones, ET and methyl jasmonate (MeJA) were applied to Reyan 7-33-97 rubber tree clones, and comparative proteomic analyses were conducted. On the basis of a transcriptome shotgun assembly (TSA) sequence database and an iTRAQ-coupled two-dimensional LC-MS/MS approach, 1499 latex proteins belonging to 1078 clusters were identified. With a 1.5-fold cut-off value to determine up- and down-regulated proteins, a total of 101 latex proteins were determined to be regulated by ET and/or MeJA via pairwise comparisons among the three exposure durations (0 h, 6 h, and 48 h). Proteins associated with latex regeneration, including phosphoenolpyruvate carboxylase and acetyl-CoA C-acetyltransferase, and those associated with latex flow, such as chitinase and a sieve element occlusion protein, were affected by the application of ET. Chitinase and polyphenol oxidase were also found to be regulated by MeJA. The findings of this study may provide new insight into the roles of phytohormones in latex yield and the causative mechanisms of laticifer differentiation in rubber trees.

SIGNIFICANCE

On the basis of a transcriptome shotgun assembly (TSA) sequence database and an iTRAQ-coupled two-dimensional LC-MS/MS approach, the most comprehensive proteome of the latex was profiled, and the ethylene-/jasmonate-responsive proteins were identified in the latex of H. brasiliensis. The findings of this study may provide new insight into the role of phytohormones in latex yield and the causative mechanisms of laticifer differentiation in rubber trees.

Functional characterization of FaNIP1;1 gene, a ripening-related and receptacle-specific aquaporin in strawberry fruit

Strawberry fruit (Fragaria × ananassa) is a soft fruit with high water content at ripe stage (more than 90% of its fresh weight). Aquaporins play an important role in plant water homeostasis, through the facilitation of water transport and solutes. We report the role played by FaNIP1;1 in the receptacle ripening process. The analysis by qRT-PCR of FaNIP1;1 showed that this gene is mainly expressed in fruit receptacle and has a ripening-related expression pattern that was accompanied by an increase in both the abscisic acid and water content of the receptacle throughout fruit ripening. Moreover, FaNIP1;1 was induced in situations of water deficit. Additionally, we show that FaNIP1;1 expression was positively regulated by abscisic acid and negatively regulated by auxins. The water transport capacity of FaNIP1;1 was determined by a stopped-flow spectroscopy in yeast over-expressing FaNIP1;1. Glycerol, H2O2 and boron transport were also demonstrated in yeast. On the other hand, GFP-FaNIP1;1 fusion protein was located in plasma membrane. In conclusion, FaNIP1;1 seems to play an important role increasing the plasma membrane permeability, that allows the water accumulation in the strawberry fruit receptacle throughout the ripening process.

Regulation of HbPIP2;3, a Latex-Abundant Water Transporter, Is Associated with Latex Dilution and Yield in the Rubber Tree (Hevea brasiliensis Muell. Arg.)

Rubber tree (Hevea brasiliensis) latex, the source of natural rubber, is synthesised in the cytoplasm of laticifers. Efficient water inflow into laticifers is crucial for latex flow and production since it is the determinant of the total solid content of latex and its fluidity after tapping. As the mature laticifer vessel rings are devoid of plasmodesmata, water exchange between laticifers and surrounding cells is believed to be governed by plasma membrane intrinsic proteins (PIPs). To identify the most important PIP aquaporin in the water balance of laticifers, the transcriptional profiles of ten-latex-expressed PIPs were analysed. One of the most abundant transcripts, designated HbPIP2;3, was characterised in this study. When tested in Xenopus laevis oocytes HbPIP2;3 showed a high efficiency in increasing plasmalemma water conductance. Expression analysis indicated that the HbPIP2;3 gene was preferentially expressed in latex, and the transcripts were up-regulated by both wounding and exogenously applied Ethrel (a commonly-used ethylene releaser). Although regular tapping up-regulated the expression of HbPIP2;3 during the first few tappings of the virginal rubber trees, the transcriptional kinetics of HbPIP2;3 to Ethrel stimulation in the regularly tapped tree exhibited a similar pattern to that of the previously reported HbPIP2;1 in the virginal rubber trees. Furthermore, the mRNA level of HbPIP2;3 was associated with clonal yield potential and the Ethrel stimulation response. Together, these results have revealed the central regulatory role of HbPIP2;3 in laticifer water balance and ethylene stimulation of latex production in Hevea.

Identification of laticifer-specific genes and their promoter regions from a natural rubber producing plant Hevea brasiliensis

Latex, the milky cytoplasm of highly differentiated cells called laticifers, from Hevea brasiliensis is a key source of commercial natural rubber production. One way to enhance natural rubber production would be to express genes involved in natural rubber biosynthesis by a laticifer-specific overexpression system. As a first step to identify promoters which could regulate the laticifer-specific expression, we identified random clones from a cDNA library of H. brasiliensis latex, resulting in 4325 expressed sequence tags (ESTs) assembled into 1308 unigenes (692 contigs and 617 singletons). Quantitative analyses of the transcription levels of high redundancy clones in the ESTs revealed genes highly and predominantly expressed in laticifers, such as Rubber Elongation Factor (REF), Small Rubber Particle Protein and putative protease inhibitor proteins. HRT1 and HRT2, cis-prenyltransferases involved in rubber biosynthesis, was also expressed predominantly in laticifers, although these transcript levels were 80-fold lower than that of REF. The 5'-upstream regions of these laticifer-specific genes were cloned and analyzed in silico, revealing seven common motifs consisting of eight bases. Furthermore, transcription factors specifically expressed in laticifers were also identified. The common motifs in the laticifer-specific genes and the laticifer-specific transcription factors are potentially involved in the regulation of gene expression in laticifers.

The symbiosis with the arbuscular mycorrhizal fungus Rhizophagus irregularis drives root water transport in flooded tomato plants

It is known that the presence of arbuscular mycorrhizal fungi within the plant roots enhances the tolerance of the host plant to different environmental stresses, although the positive effect of the fungi in plants under waterlogged conditions has not been well studied. Tolerance of plants to flooding can be achieved through different molecular, physiological and anatomical adaptations, which will affect their water uptake capacity and therefore their root hydraulic properties. Here, we investigated the root hydraulic properties under non-flooded and flooded conditions in non-mycorrhizal tomato plants and plants inoculated with the arbuscular mycorrhizal fungus Rhizophagus irregularis. Only flooded mycorrhizal plants increased their root hydraulic conductivity, and this effect was correlated with a higher expression of the plant aquaporin SlPIP1;7 and the fungal aquaporin GintAQP1. There was also a higher abundance of the PIP2 protein phoshorylated at Ser280 in mycorrhizal flooded plants. The role of plant hormones (ethylene, ABA and IAA) in root hydraulic properties was also taken into consideration, and it was concluded that, in mycorrhizal flooded plants, ethylene has a secondary role regulating root hydraulic conductivity whereas IAA may be the key hormone that allows the enhancement of root hydraulic conductivity in mycorrhizal plants under low oxygen conditions.

Comparative analysis of latex transcriptome reveals putative molecular mechanisms underlying super productivity of Hevea brasiliensis

Increasing demand for natural rubber prompts studies into the mechanisms governing the productivity of rubber tree (Heveabrasiliensis). It is very interesting to notice that a rubber tree of clone PR107 in Yunnan, China is reported to yield more than 20 times higher than the average rubber tree. This super-high-yielding (SHY) rubber tree (designated as SY107), produced 4.12 kg of latex (cytoplasm of rubber producing laticifers, containing about 30% of rubber) per tapping, more than 7-fold higher than that of the control. This rubber tree is therefore a good material to study how the rubber production is regulated at a molecular aspect. A comprehensive cDNA-AFLP transcript profiling was performed on the latex of SY107 and its average counterparts by using the 384 selective primer pairs for two restriction enzyme combinations (ApoI/MseI and TaqI/MseI). A total of 746 differentially expressed (DE) transcript-derived fragments (TDFs) were identified, of which the expression patterns of 453 TDFs were further confirmed by RT-PCR. These RT-PCR confirmed TDFs represented 352 non-redundant genes, of which 215 had known or partially known functions and were grouped into 10 functional categories. The top three largest categories were transcription and protein synthesis (representing 24.7% of the total genes), defense and stress (15.3%), and primary and secondary metabolism (14.0%). Detailed analysis of the DE-genes suggests notable characteristics of SHY phenotype in improved sucrose loading capability, rubber biosynthesis-preferred sugar utilization, enhanced general metabolism and timely stress alleviation. However, the SHY phenotype has little correlation with rubber-biosynthesis pathway genes.

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.

Transcriptome and gene expression analysis in cold-acclimated guayule (Parthenium argentatum) rubber-producing tissue

Natural rubber biosynthesis in guayule (Parthenium argentatum Gray) is associated with moderately cold night temperatures. To begin to dissect the molecular events triggered by cold temperatures that govern rubber synthesis induction in guayule, the transcriptome of bark tissue, where rubber is produced, was investigated. A total of 11,748 quality expressed sequence tags (ESTs) were obtained. The vast majority of ESTs encoded proteins that are similar to stress-related proteins, whereas those encoding rubber biosynthesis-related proteins comprised just over one percent of the ESTs. Sequence information derived from the ESTs was used to design primers for quantitative analysis of the expression of genes that encode selected enzymes and proteins with potential impact on rubber biosynthesis in field-grown guayule plants, including 3-hydroxy-3-methylglutaryl-CoA synthase, 3-hydroxy-3-methylglutaryl-CoA reductase, farnesyl pyrophosphate synthase, squalene synthase, small rubber particle protein, allene oxide synthase, and cis-prenyl transferase. Gene expression was studied for field-grown plants during the normal course of seasonal variation in temperature (monthly average maximum 41.7 °C to minimum 0 °C, from November 2005 through March 2007) and rubber transferase enzymatic activity was also evaluated. Levels of gene expression did not correlate with air temperatures nor with rubber transferase activity. Interestingly, a sudden increase in night temperature 10 days before harvest took place in advance of the highest CPT gene expression level.

Molecular characterization and expression analysis of cDNAs encoding four Rab and two Arf GTPases in the latex of Hevea brasiliensis

In plants, Rab and Arf GTPases are key regulators of vesicle trafficking. To investigate whether these small GTPases (SG) play a role in the regulation of the regeneration of latex (the cytoplasm of the rubber-producing laticifer cell) in Hevea brasiliensis (Hevea hereafter), full-length cDNAs that encode four HbRab and two HbArf GTPases were cloned. The four HbRab proteins showed specific GTP-binding activity when expressed in Escherichia coli. Transcript expression of the six SG genes was investigated by real-time RT-PCR. All genes revealed to be expressed in each of the six Hevea tissues examined, but the expression patterns were different. Four genes, HbArf1, HbRab2, HbRab3 and HbRab4, displayed a preferential expression in latex. The expression of all genes was upregulated by the act of latex exploitation (tapping), and HbRab1 had the highest level of upregulation. Wounding markedly upregulated the expression of two SG genes (HbRab1 and HbArf2), and exogenous methyl jasmonate upregulated all six SG genes. Wounding might upregulate the expression of HbRab1 and HbArf2 through a jasmonic acid-mediated signaling pathway. None of the genes were markedly upregulated by Ethrel (an ethylene releaser and latex stimulator); instead, HbArf2 and HbRab4 were downregulated significantly after a 24 h treatment with Ethrel. This paper gives the first description of Rab and Arf GTPases in Hevea and provides clues for their involvement, HbRab1 in particular, in latex regeneration.