Comparative transcriptome analysis of latex from rubber tree clone CATAS8-79 and PR107 reveals new cues for the regulation of latex regeneration and duration of latex flow

Integrative analysis of transcriptome and metabolome reveals how ethylene increases natural rubber yield in Hevea brasiliensis

Hevea brasiliensis is an important cash crop with the product named natural rubber (NR) for markets. Ethylene (ET) is the most effective yield stimulant in NR production but the molecular mechanism remains incomplete. Here, latex properties analysis, transcriptome analysis, and metabolic profiling were performed to investigate the mechanism of NR yield increase in four consecutive tappings after ET stimulation. The results revealed that sucrose and inorganic phosphate content correlated positively with dry-rubber yield and were induced upon ET stimulation. Stimulation with ET also led to significant changes in gene expression and metabolite content. Genes involved in phytohormone biosynthesis and general signal transduction as well as 51 transcription factors potentially involved in the ET response were also identified. Additionally, KEGG annotation of differentially accumulated metabolites suggested that metabolites involved in secondary metabolites, amino-acid biosynthesis, ABC transporters, and galactose metabolism were accumulated in response to ET. Integrative analysis of the data collected by transcriptomics and metabolomics identified those differentially expressed genes and differentially accumulated metabolites are mainly involved in amino-acid biosynthesis and carbohydrate metabolism. Correlation analysis of genes and metabolites showed a strong correlation between amino-acid biosynthesis during ET stimulation. These findings provide new insights into the molecular mechanism underlying the ET-induced increase in rubber yield and further our understanding of the regulatory mechanism of ethylene signaling in rubber biosynthesis.

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.

Genomic insight into domestication of rubber tree

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

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

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

Comparative analysis of latex transcriptomes reveals the potential mechanisms underlying rubber molecular weight variations between the Hevea brasiliensis clones RRIM600 and Reyan7-33-97

BACKGROUND

The processabilities and mechanical properties of natural rubber depend greatly on its molecular weight (MW) and molecular weight distribution (MWD). However, the mechanisms underlying the regulation of molecular weight during rubber biosynthesis remain unclear.

RESULTS

In the present study, we determined the MW and particle size of latex from 1-year-old virgin trees and 30-year-old regularly tapped trees of the Hevea clones Reyan7-33-97 and RRIM600. The results showed that both the MW and the particle size of latex varied between these two clones and increased with tree age. Latex from RRIM600 trees had a smaller average particle size than that from Reyan7-33-97 trees of the same age. In 1-year-old trees, the Reyan7-33-97 latex displayed a slightly higher MW than that of RRIM600, whereas in 30-year-old trees, the RRIM600 latex had a significantly higher MW than the Reyan7-33-97 latex. Comparative analysis of the transcriptome profiles indicated that the average rubber particle size is negatively correlated with the expression levels of rubber particle associated proteins, and that the high-MW traits of latex are closely correlated with the enhanced expression of isopentenyl pyrophosphate (IPP) monomer-generating pathway genes and downstream allylic diphosphate (APP) initiator-consuming non-rubber pathways. By bioinformatics analysis, we further identified a group of transcription factors that potentially regulate the biosynthesis of IPP.

CONCLUSIONS

Altogether, our results revealed the potential regulatory mechanisms involving gene expression variations in IPP-generating pathways and the non-rubber isoprenoid pathways, which affect the ratios and contents of IPP and APP initiators, resulting in significant rubber MW variations among same-aged trees of the Hevea clones Reyan7-33-97 and RRIM600. Our findings provide a better understanding of rubber biosynthesis and lay the foundation for genetic improvement of rubber quality in H. brasiliensis.

Genome-wide identification and expression analysis of the phosphatase 2A family in rubber tree (Hevea brasiliensis)

The protein phosphatase 2As (PP2As) play a key role in manipulating protein phosphorylation. Although a number of proteins in the latex of laticifers are phosphorylated during latex regeneration in rubber tree, information about the PP2A family is limited. In the present study, 36 members of the HbPP2A family were genome-wide identified. They were clustered into five subgroups: the subgroup HbPP2AA (4), HbPP2AB' (14), HbPP2AB'' (6), HbPP2AB55 (4), and HbPP2AC (8). The members within the same subgroup shared highly conserved gene structures and protein motifs. Most of HbPP2As possessed ethylene- and wounding-responsive cis-acting elements. The transcripts of 29 genes could be detected in latex by using published high-throughput sequencing data. Of the 29 genes, seventeen genes were significantly down-regulated while HbPP2AA1-1 and HbPP2AB55α/Bα-1were up-regulated by tapping. Of the 17 genes, 14 genes were further significantly down-regulated by ethrel application. The down-regulated expression of a large number of HbPP2As may attribute to the enhanced phosphorylation of the proteins in latex from the tapped trees and the trees treated with ethrel application.

Response of bacterial communities in rubber plantations to different fertilizer treatments

In the present study, the effects of chemical fertilizer (CF) and organic fertilizer plus chemical fertilizer application (OF-CF) on natural rubber yield, soil properties, and soil bacterial community were systematically investigated in rubber plantations. The rubber dry yield was 26.3% more in the OF treatment group than in the CF treatment group. The contents of total nitrogen (TN), available nitrogen (AN), available phosphorus (AP), and available potassium (AK) as well as soil organic matter (SOM) and pH value were higher following OF-CF treatment. Using Illumina sequencing, a total of 927 operational taxonomic units (OTUs) were obtained following CF treatment, while 955 OTUs were obtained after OF-CF treatment. Relative abundance analysis showed the relative abundances of four phyla (Acidobacteria, Proteobacteria, Actinobacteria, Gemmatimonadetes) were different between the two treatment groups. Correlation analysis revealed Acidobacteria, Bacteroidetes, Thaumarchaeota, Elusimicrobia, Verrucomicrobia were the key taxa that determined the soil properties. Additionally, five OTUs (OTU_506, OTU_391, OTU_189, OTU_278, OTU_1057) were thought to be related to the biodegradation of natural rubber. Taken together, these results improve our understanding of the OF-mediated improvement in soil fertility and contribute to the identification of rubber-degrading bacteria in rubber plantations.

Comparative Proteomic Analysis of the Hevea brasiliensis Latex under Ethylene and Calcium Stimulation

BACKGROUND

Calcium ions usually act as a second messenger in the signal transmission process and a major element required by plants. In Hevea, calcium ion could alleviate the negative effects of long-term ethylene application to a certain extent. However, the molecular mechanisms remain unclear.

METHODS

Two-dimensional electrophoresis was used to determine the pattern of protein changes in latex after treatments with calcium and/or ethylene. Quantitative real-time polymerase chain reaction and Western blotting were used to determine the expression levels of some proteins and genes. STRING software was used to determine the protein-protein interaction network of the identified proteins.

RESULTS

Comparative proteomics identified 145 differentially expressed proteins, which represented 103 unique proteins. The abundance change patterns of some proteins involved in signal transduction, rubber particle aggregation, and natural rubber biosynthesis were altered upon calcium stimulation. Quantitative real-time polymerase chain reaction analysis of 29 proteins showed that gene expression did not always maintain the same trend as protein expression. The increased enzyme activities of superoxide dismutase, ascorbate peroxidase, and glutathione reductase suggested that calcium can enhance the antistress ability of plants by increasing the activity of their antioxidant enzyme systems.

CONCLUSION

These results supplement the rubber latex proteome, and provide evidence for investigating the molecular mechanisms by which calcium alleviates the negative effects of ethylene stimulation.

Genomic technologies for Hevea breeding

The commercial production of high quality natural rubber (NR) solely depends on Hevea brasiliensis Muell. Arg, (Para rubber tree) and accounts for >98% of total production worldwide. NR with its unique properties is an essential commodity for the automobile industry and its synthetic counterparts are in no way substitute to it. The rubber tree genome is very complex and plays an important role in delivering the unique properties of Hevea. But a lack of knowledge on the molecular mechanisms of rubber biosynthesis, disease resistance, etc., in elite clones of rubber still persists. Marker-assisted selection and transgenic techniques were proved to be advantageous in improving the breeding efficiency for latex yield, disease resistance, etc. The suppression subtractive hybridization (SSH), in the form of subtracted cDNA libraries and microarrays, can assist in searching the functions of expressed genes (candidate gene approach). Expressed sequence tags (ESTs) related to various metabolic aspects are well utilized to create EST banks that broadly represent the genes expressed in one tissue, such as latex cells, that assists in the study of gene function and regulation. Transcriptome analysis and gene mapping have been accomplished in Hevea at various stages. However, a selection criterion to delineate high yielding genotypes at the juvenile stage has not been accomplished so far. This is the main pit fall for rubber breeding apart from stock-scion interactions leading to yield differences among a clonally multiplied population. At least four draft genome sequences have been published on Hevea rubber, and all give different genome size and contig lengths-a comprehensive and acceptable genomic map remains unfulfilled. The progress made in molecular markers, latex biosynthesis genes, transcriptome analysis, chloroplast and mitochondrial DNA diversity, paternity identification through Breeding without Breeding (BwB), stimulated latex production and its molecular intricacies, molecular biology of tapping panel dryness, genomics for changed climates and genome mapping are discussed in this review. These information can be utilized to improvise the molecular breeding programs of Hevea in future.

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

Biosynthesis of Natural Rubber: Current State and Perspectives

Natural rubber is a kind of indispensable biopolymers with great use and strategic importance in human society. However, its production relies almost exclusively on rubber-producing plants Hevea brasiliensis, which have high requirements for growth conditions, and the mechanism of natural rubber biosynthesis remains largely unknown. In the past two decades, details of the rubber chain polymerization and proteins involved in natural rubber biosynthesis have been investigated intensively. Meanwhile, omics and other advanced biotechnologies bring new insight into rubber production and development of new rubber-producing plants. This review summarizes the achievements of the past two decades in understanding the biosynthesis of natural rubber, especially the massive information obtained from the omics analyses. Possibilities of natural rubber biosynthesis in vitro or in genetically engineered microorganisms are also discussed.

Comparative transcriptome analysis reveals an early gene expression profile that contributes to cold resistance in Hevea brasiliensis (the Para rubber tree)

The rubber tree (Hevea brasiliensis Muell. Arg) is a tropical, perennial, woody plant that is susceptible to cold stress. In China, cold stress has been found to severely damage rubber plants in plantations in past decades. Although several Hevea clones that are resistant to cold have been developed, their cold hardiness mechanism has yet to be elucidated. For the study reported herein, we subjected the cold-resistant clone CATAS93-114 and the cold-sensitive clone Reken501 to chilling stress, and characterized their transcriptomes at 0, 2, 8 and 24 h after the start of chilling. We found that 7870 genes were differentially expressed in the transcriptomes of the two clones. In CATAS93-114, a greater number of genes were found to be up- or downregulated between 2 h and 8 h than in Reken501, which indicated a more rapid and intensive response by CATAS93-114 than by Reken501. The differentially expressed genes were grouped into seven major clusters, according to their Gene Ontology terms. The expression profiles for genes involved in abscisic acid metabolism and signaling, in an abscisic acid-independent pathway, and in early signal perception were found to have distinct expression patterns for the transcriptomes of the two clones. The differential expression of 22 genes that appeared to have central roles in response to chilling was confirmed by quantitative real-time PCR.

Physiological and Proteomic Analyses of Molybdenum- and Ethylene-Responsive Mechanisms in Rubber Latex

Molybdenum (Mo) is an essential micronutrient in many plants. In the rubber tree Hevea brasiliensis, Mo application can reduce the shrinkage of the tapping line, decrease tapping panel dryness, and finally increase rubber latex yield. After combined Mo with ethylene (Eth), these effects become more obvious. However, the molecular mechanism remains unclear. Here, we compared the changed patterns of physiological parameters and protein accumulation in rubber latex after treated with Mo and/or Eth. Our results demonstrated that both Eth and Mo can improve the contents of thiol, sucrose, and dry yield in rubber latex. However, lutoid bursting is significantly inhibited by Mo. Comparative proteomics identified 169 differentially expressed proteins, including 114 unique proteins, which are mainly involved in posttranslational modification, carbohydrate metabolism, and energy production. The abundances of several proteins involved in rubber particle aggregation are decreased upon Mo stimulation, while many enzymes related to natural rubber biosynthesis are increased. Comparison of the accumulation patterns of 25 proteins revealed that a large portion of proteins have different changed patterns with their gene expression levels. Activity assays of six enzymes revealed that Mo stimulation can increase latex yield by improving the activity of some Mo-responsive enzymes. These results not only deepen our understanding of the rubber latex proteome but also provide new insights into the molecular mechanism of Mo-stimulated rubber latex yield.

Transcriptomic analyses reveal molecular mechanisms underlying growth heterosis and weakness of rubber tree seedlings

BACKGROUND

Breeding rubber tree seedling with growth heterosis is vital for natural rubber production. It is the prerequisites for effectively utilizing growth heterosis to elucidate its molecular mechanisms, but the molecular mechanisms remain poorly understood in rubber tree. To elucidate seedling growth heterosis, we conducted comparative transcriptomic analyses between the two hybrids and their parents.

RESULTS

By identifying and comparing differently expressed genes (DEGs), we found that the hybrids (BT 3410 and WC 11) show significantly differential expression profiles from their parents (PR 107 and RRIM 600). In BT 3410-parent triad, 1092 (49.95%) and 1094 (50.05%) DEGs indicated clear underdominance or overdominance, respectively. Whereas in WC 11-parent triad, most DEGs (78.2%, 721) showed low- or high-parent dominance; 160 (17.35%) exhibited expression patterns that are not statistically distinguishable from additivity, and 8 (0.87%) and 33 (3.58%) DEGs exhibited underdominance and overdominance, respectively. Furthermore, some biological processes are differentially regulated between two hybrids. Interestingly, the pathway in response to stimulus is significantly downregulated and metabolic pathways are more highly regulated in BT 3410.

CONCLUSIONS

Taken together, the genotypes, transcriptomes and biological pathways (especially, carbohydrate metabolism) are highly divergent between two hybrids, which may be associated with growth heterosis and weakness. Analyzing gene action models in hybrid-parent triads, we propose that overdominance may play important roles on growth heterosis, whereas dominance on hybrid weakness in rubber tree seedlings. These findings bring new insights into our understanding of growth heterosis of rubber tree seedling.

Comparison of Gene Expression Profiles in Nonmodel Eukaryotic Organisms with RNA-Seq

With recent advances of next-generation sequencing technology, RNA-Sequencing (RNA-Seq) has emerged as a powerful approach for the transcriptomic profiling. RNA-Seq has been used in almost every field of biological studies, and has greatly extended our view of transcriptomic complexity in different species. In particular, for nonmodel organisms which are usually without high-quality reference genomes, the de novo transcriptome assembly from RNA-Seq data provides a solution for their comparative transcriptomic study. In this chapter, we focus on the comparative transcriptomic analysis of nonmodel organisms. Two analysis strategies (without or with reference genome) are described step-by-step, with the differentially expressed genes explored.

Transcript Profiling of Hevea brasiliensis during Latex Flow

Latex exploitation enhances latex regeneration in rubber trees. The latex exploitation-caused latex flow lasts from 10 min to a few hours, which is convenient for exploring the transcript profiling of latex metabolism-related genes at the different stages of latex flow. In the present study, the expression pattern of 62 latex metabolism-related genes involved in water transportation, carbohydrate metabolism, natural rubber biosynthesis, hormone signaling, ROS generation and scavenging, and latex coagulum across three stages of latex flow between rubber tree clones CATAS7-33-97 and CATAS8-79 were comparatively analyzed by quantitative real-time PCR. The two clones show differences in latex regeneration and have a different duration of latex flow. The results showed that the expression levels of 38 genes were significantly higher in CATAS8-79 latex than in CATAS7-33-97 during latex regeneration, while 45 genes had a notably higher expression level in CATAS8-79 latex during latex flow. Together with the activation of the MEP pathway and jasmonate pathway in CATAS8-79 latex, HbPIP1;3, HbPIP1;4, HbSUT3, HbSus3, HbHMGS1-2, HbMK should contribute to the high latex regeneration ability. The up-regulation of ethylene signaling and Hb44KD and the down-regulation of latex coagulation-related genes in CATAS8-79 latex might contribute to its longer latex flow duration. This study provides some cues for revealing the regulation of latex metabolism in rubber trees.

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

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

The SWEET gene family in Hevea brasiliensis - its evolution and expression compared with four other plant species

SWEET proteins play an indispensable role as a sugar efflux transporter in plant development and stress responses. The SWEET genes have previously been characterized in several plants. Here, we present a comprehensive analysis of this gene family in the rubber tree, Hevea brasiliensis. There are 36 members of the SWEET gene family in this species, making it one of the largest families in plant genomes sequenced so far. Structure and phylogeny analyses of these genes in Hevea and in other species demonstrated broad evolutionary conservation. RNA‐seq analyses revealed that SWEET2, 16, and 17 might represent the main evolutionary direction of SWEET genes in plants. Our results in Hevea suggested the involvement of HbSWEET1a, 2e, 2f, and 3b in phloem loading, HbSWEET10a and 16b in laticifer sugar transport, and HbSWEET9a in nectary‐specific sugar transport. Parallel studies of RNA‐seq analyses extended to three other plant species (Manihot esculenta, Populus trichocarpa, and Arabidopsis thaliana) produced findings which implicated MeSWEET10a, 3a, and 15b in M. esculenta storage root development, and the involvement of PtSWEET16b and PtSWEET16d in P. trichocarpa xylem development. RT‐qPCR results further revealed that HbSWEET10a, 16b, and 1a play important roles in phloem sugar transport. The results from this study provide a foundation not only for further investigation into the functionality of the SWEET gene family in Hevea, especially in its sugar transport for latex production, but also for related studies of this gene family in the plant kingdom.

Comparative morphology and transcriptome analysis reveals distinct functions of the primary and secondary laticifer cells in the rubber tree

Laticifers are highly specialized cells that synthesize and store natural rubber. Rubber trees (Hevea brasiliensis Muell. Arg.) contain both primary and secondary laticifers. Morphological and functional differences between the two types of laticifers are largely unknown, but such information is important for breeding and cultivation practices. Morphological comparison using paraffin sections revealed only distribution differences: the primary laticifers were distributed randomly, while the secondary laticifers were distributed in concentric rings. Using isolated laticifer networks, the primary laticifers were shown to develop via intrusive "budding" and formed necklace-like morphology, while the secondary laticifers developed straight and smooth cell walls. Comparative transcriptome analysis indicated that genes involved in cell wall modification, such as pectin esterase, lignin metabolic enzymes, and expansins, were highly up-regulated in the primary laticifers and correspond to its necklace-like morphology. Genes involved in defense against biotic stresses and rubber biosynthesis were highly up-regulated in the primary laticifers, whereas genes involved in abiotic stresses and dormancy were up-regulated in the secondary laticifers, suggesting that the primary laticifers are more adequately prepared to defend against biotic stresses, while the secondary laticifers are more adequately prepared to defend against abiotic stresses. Therefore, the two types of laticifers are morphologically and functionally distinct.

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 REF_138,175,258 and SRPP_117,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 REF₂₅₈ 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 REF₁₇₅ or REF₂₅₈ 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.

Reactive oxygen species in Hevea brasiliensis latex and relevance to Tapping Panel Dryness

Environmental stress can lead to oxidative stress resulting from an increase in reactive oxygen species (ROS) and involves redox adjustments. Natural rubber is synthesized in laticifers, which is a non-photosynthetic tissue particularly prone to oxidative stress. This paper reviews the current state of knowledge on the ROS production and ROS-scavenging systems in laticifers. These regulations have been the subject of intense research into a physiological syndrome, called Tapping Panel Dryness (TPD), affecting latex production in Hevea brasiliensis. In order to prevent TPD occurrence, monitoring thiol content appeared to be a crucial factor of latex diagnosis. Thiols, ascorbate and γ-tocotrienol are the major antioxidants in latex. They are involved in membrane protection from ROS and likely have an effect on the quality of raw rubber. Some transcription factors might play a role in the redox regulatory network in Hevea, in particular ethylene response factors, which have been the most intensively studied given the role of ethylene on rubber production. Current challenges for rubber research and development with regard to redox systems will involve improving antioxidant capacity using natural genetic variability.

Comparative transcriptome investigation of global gene expression changes caused by miR156 overexpression in Medicago sativa

Background

Medicago sativa (alfalfa) is a low-input forage and potential bioenergy crop, and improving its yield and quality has always been a focus of the alfalfa breeding industry. Transgenic alfalfa plants overexpressing a precursor of alfalfa microRNA156 (MsmiR156) were recently generated by our group. These plants (miR156OE) showed enhanced biomass yield, reduced internodal length, increased shoot branching and trichome density, and a delay in flowering time. Transcripts of three SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes (MsSPL6, MsSPL12, and MsSPL13) were found to be targeted for cleavage by MsmiR156 in alfalfa.

Results

To further illustrate the molecular mechanisms underlying the effects of miR156 in alfalfa, two miR156OE genotypes (A11a and A17) were subjected to Next Generation RNA Sequencing with Illumina HiSeq. More than 1.11 billion clean reads were obtained from our available sequenced samples. A total of 160,472 transcripts were generated using Trinity de novo assembly and 4,985 significantly differentially expressed genes were detected in miR156OE plants A11a and A17 using the Medicago truncatula genome as reference. A total of 17 genes (including upregulated, downregulated, and unchanged) were selected for quantitative real-time PCR (qRT-PCR) validation, which showed that gene expression levels were largely consistent between qRT-PCR and RNA-Seq data. In addition to the established SPL genes MsSPL6, MsSPL12 and MsSPL13, four new SPLs; MsSPL2, MsSPL3, MsSPL4 and MsSPL9 were also down-regulated significantly in both miR156OE plants. These seven SPL genes belong to genes phylogeny clades VI, IV, VIII, V and VII, which have been reported to be targeted by miR156 in Arabidopsis thaliana. The gene ontology terms characterized electron transporter, starch synthase activity, sucrose transport, sucrose-phosphate synthase activity, chitin binding, sexual reproduction, flavonoid biosynthesis and lignin catabolism correlate well to the phenotypes of miR156OE alfalfa plants.

Conclusions

This is the first report of changes in global gene expression in response to miR156 overexpression in alfalfa. The discovered miR156-targeted SPL genes belonging to different clades indicate miR156 plays fundamental and multifunctional roles in regulating alfalfa plant development.

Electronic supplementary material

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

Evaluation of Reference Genes for Quantitative Real-Time PCR Analysis of the Gene Expression in Laticifers on the Basis of Latex Flow in Rubber Tree (Hevea brasiliensis Muell. Arg.)

Latex exploitation-caused latex flow is effective in enhancing latex regeneration in laticifer cells of rubber tree. It should be suitable for screening appropriate reference gene for analysis of the expression of latex regeneration-related genes by quantitative real-time PCR (qRT-PCR). In the present study, the expression stability of 23 candidate reference genes was evaluated on the basis of latex flow by using geNorm and NormFinder algorithms. Ubiquitin-protein ligase 2a (UBC2a) and ubiquitin-protein ligase 2b (UBC2b) were the two most stable genes among the selected candidate references in rubber tree clones with differential duration of latex flow. The two genes were also high-ranked in previous reference gene screening across different tissues and experimental conditions. By contrast, the transcripts of latex regeneration-related genes fluctuated significantly during latex flow. The results suggest that screening reference gene during latex flow should be an efficient and effective clue for selection of reference genes in qRT-PCR.

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.

The rubber tree genome reveals new insights into rubber production and species adaptation

The Para rubber tree (Hevea brasiliensis) is an economically important tropical tree species that produces natural rubber, an essential industrial raw material. Here we present a high-quality genome assembly of this species (1.37 Gb, scaffold N50 = 1.28 Mb) that covers 93.8% of the genome (1.47 Gb) and harbours 43,792 predicted protein-coding genes. A striking expansion of the REF/SRPP (rubber elongation factor/small rubber particle protein) gene family and its divergence into several laticifer-specific isoforms seem crucial for rubber biosynthesis. The REF/SRPP family has isoforms with sizes similar to or larger than SRPP1 (204 amino acids) in 17 other plants examined, but no isoforms with similar sizes to REF1 (138 amino acids), the predominant molecular variant. A pivotal point in Hevea evolution was the emergence of REF1, which is located on the surface of large rubber particles that account for 93% of rubber in the latex (despite constituting only 6% of total rubber particles, large and small). The stringent control of ethylene synthesis under active ethylene signalling and response in laticifers resolves a longstanding mystery of ethylene stimulation in rubber production. Our study, which includes the re-sequencing of five other Hevea cultivars and extensive RNA-seq data, provides a valuable resource for functional genomics and tools for breeding elite Hevea cultivars.

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.

Comparative Transcriptome Analysis of Latex Reveals Molecular Mechanisms Underlying Increased Rubber Yield in Hevea brasiliensis Self-Rooting Juvenile Clones

Rubber tree (Hevea brasiliensis) self-rooting juvenile clones (JCs) are promising planting materials for rubber production. In a comparative trial between self-rooting JCs and donor clones (DCs), self-rooting JCs exhibited better performance in rubber yield. To study the molecular mechanism associated with higher rubber yield in self-rooting JCs, we sequenced and comparatively analyzed the latex of rubber tree self-rooting JCs and DCs at the transcriptome level. Total raw reads of 34,632,012 and 35,913,020 bp were obtained from the library of self-rooting JCs and DCs, respectively, by using Illumina HiSeq 2000 sequencing technology. De novo assemblies yielded 54689 unigenes from the library of self-rooting JCs and DCs. Among 54689 genes, 1716 genes were identified as differentially expressed between self-rooting JCs and DCs via comparative transcript profiling. Functional analysis showed that the genes related to the mass of categories were differentially enriched between the two clones. Several genes involved in carbohydrate metabolism, hormone metabolism and reactive oxygen species scavenging were up-regulated in self-rooting JCs, suggesting that the self-rooting JCs provide sufficient molecular basis for the increased rubber yielding, especially in the aspects of improved latex metabolisms and latex flow. Some genes encoding epigenetic modification enzymes were also differentially expressed between self-rooting JCs and DCs. Epigenetic modifications may lead to gene differential expression between self-rooting JCs and DCs. These data will provide new cues to understand the molecular mechanism underlying the improved rubber yield of H. brasiliensis self-rooting clones.

Cloning, heterologous expression and characterization of ascorbate peroxidase (APX) gene in laticifer cells of rubber tree (Hevea brasiliensis Muell. Arg.)

Ascorbate peroxidases (APXs) are a kind of crucial enzymes for removing reactive oxygen species (ROS) in plant cell. In the present study, a full-length cDNA encoding an APX, designated HbAPX, was isolated from Hevea brasiliensis by the rapid amplification of cDNA ends (RACE) method. HbAPX was 1174-bp in length and contained a 912-bp open reading frame (ORF) encoding a putative protein of 304 amino acids. The predicted molecular mass of HbAPX was 27.6 kDa (kDa) with an isoelectric point (pI) of 6.73. The phylogenetic analysis showed that HbAPX belonged to the cytosolic subgroup and was more relative to PtAPX and MdAPX2. By using PlantCare online analysis, such cis-acting elements as W-box and MRE were detected in the promoter region of HbAPX. Overproduction of recombinant HbAPX protein either in Escherichia coli or yeast enhanced their tolerance to such abiotic stresses as Cu(2+), Zn(2+), Na(2+) and hydrogen peroxide (H2O2). Ethrel application significantly down-regulated the expression of HbAPX and inhibited the activity of HbAPX in vivo. The ethrel-caused down-regulation of HbAPX may disturb the redox homeostasis in laticifer cells of rubber tree.