Comparative proteomics of primary and secondary lutoids reveals that chitinase and glucanase play a crucial combined role in rubber particle aggregation in Hevea brasiliensis

Molecular Identification and Characterization of Hevein Antimicrobial Peptide Genes in Two-Row and Six-Row Cultivars of Barley (Hordeum vulgare L.)

Heveins are one of the most important groups of plant antimicrobial peptides. So far, various roles in plant growth and development and in response to biotic and abiotic stresses have reported for heveins. The present study aimed to identify and characterize the hevein genes in two-row and six-row cultivars of barley. In total, thirteen hevein genes were identified in the genome of two-row and six-row cultivars of barley. The identified heveins were identical in two-row and six-row cultivars of barley and showed a high similarity with heveins from other plant species. The hevein coding sequences produced open reading frames (ORFs) ranged from 342 to 1002 bp. Most of the identified hevein genes were intronless, and the others had only one intron. The hevein ORFs produced proteins ranged from 113 to 333 amino acids. Search for conserved functional domains showed CBD and LYZ domains in barley heveins. All barley heveins comprised extracellular signal peptides ranged from 19 to 35 amino acids. The phylogenetic analysis divided barley heveins into two groups. The promoter analysis showed regulatory elements with different frequencies between two-row and six-row cultivars. These cis-acting elements included elements related to growth and development, hormone response, and environmental stresses. The expression analysis showed high expression level of heveins in root and reproductive organs of both two-row and six-row cultivars. The expression analysis also showed that barley heveins is induced by both biotic and abiotic stresses. The results of antimicrobial activity prediction showed the highest antimicrobial activity in CBD domain of barley heveins. The findings of the current study can improve our knowledge about the role of hevein genes in plant and can be used for future studies.

Construction of the first high-density SNP genetic map and identification of QTLs for the natural rubber content in Taraxacum kok-saghyz Rodin

BACKGROUND

Taraxacum kok-saghyz Rodin (TKS) is a promising commercial alternative natural rubber (NR) yielding plant. Cultivating TKS with a high NR content is an important breeding target, and developing molecular markers related to NR content can effectively accelerate the breeding process of TKS.

RESULTS

To construct a high-density SNP genetic map and uncover genomic regions related to the NR content in TKS, an F₁ mapping population of TKS was constructed by crossing two parents (l66 and X51) with significant differences in NR contents. The NR content of the F₁ plants ranged from 0.30 to 15.14% and was distributed normally with a coefficient of variation of 47.61%, indicating quantitative trait inheritance. Then, employing whole-genome resequencing (WGR), a TKS genetic linkage map of 12,680 bin markers comprising 322,439 SNPs was generated. Based on the genetic map and NR content of the F₁ population, six quantitative trait loci (QTLs) for NR content with LOD > 4.0 were identified on LG01/Chr01 and LG06/Chr06. Of them, the 2.17 Mb genomic region between qHRC-C6-1 and qHRC-C6-2 on ChrA06, with 65.62% PVE in total, was the major QTL region. In addition, the six QTLs have significant additive genetic effects on NR content and could be used to develop markers for marker-assisted selection (MAS) in TKS with a high NR content.

CONCLUSION

This work constructed the first high-density TKS genetic map and identified the QTLs and genomic regions controlling the NR content, which provides useful information for fine mapping, map-based cloning, and MAS in TKS.

Identification of cis conformation natural rubber and proteins in Ficus altissima Blume latex

Ficus altissima Blume, a horticultural plant in tropical and warm subtropical regions, can produce natural rubber with comparable molecular weight to the Hevea brasiliensis rubber. The F. altissima latex has an acidic pH (about 4.89). The rubber particle size distribution is a unimodal profile, and the peak frequency is at a size of 4.5 μm. The natural rubber of F. altissima was determined to be a cis conformation via ¹³C NMR. The Mp (molecular weight of the peak maxima) of the deproteinized F. altissima rubber was 9.34 × 10⁵ Da. LC-MS was used to identify the proteins of rubber particles and serum. The most abundant protein of the creamy rubber particle layer is an acid phosphatase, while the most abundant proteins of serum were an (R)-mandelonitrilelyase and a polygalacturonase inhibitor. Pharmaceutical proteins (ficins) or enzymes related to the biosynthesis of natural medicines (a cannabidiolic acid synthase and two lupeol synthase) were identified in F. altissima latex. The data of this study may be helpful for research on the functions of latex in latex-borne plants and the biosynthesis mechanism of natural rubber.

Hevea brasiliensis latex proteomics: a review of analytical methods and the way forward

Natural rubber or latex from the Hevea brasiliensis is an important commodity in various economic sectors in today's modern society. Proteins have been detected in latex since the early twentieth century, and they are known to regulate various biological pathways within the H. brasiliensis trees such as the natural rubber biosynthesis, defence against pathogens, wound healing, and stress tolerance. However, the exact mechanisms of the pathways are still not clear. Proteomic analyses on latex have found various proteins and revealed how they fit into the mechanisms of the biological pathways. In the past three decades, there has been rapid latex protein identification due to the improvement of latex protein extraction methods, as well as the emergence of two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). In this manuscript, we reviewed the methods of latex protein extraction that keeps on improving over the past three decades as well as the results of numerous latex protein identification and quantitation.

Comparison of Morphological Characteristics and Determination of Different Patterns for Rubber Particles in Dandelion and Different Rubber Grass Varieties

Russian dandelion Taraxacum kok-saghyz (TKS) is one promising alternative crop for natural rubber production. However, it is easily confused with other dandelions. In this study, we performed a systematical comparison of the morphological characteristics for different TKS varieties and common dandelion Taraxacum officinale (TO). Our results demonstrated that several obvious differences in morphology can be found between TKS and TO. TO leaf is a pinnate shape, its margin is heavily jagged and its base is cuneate, but TKS leaf is more cuneate and its leaf margin is nearly smooth and round. There are obvious differences for the outer bracts of TO and TKS flower buds. TKS bracts are oblanceolate, apex obtuse, margin smooth and sinuate, and its outer layer of flower buds and faceplate involucre sepal is buckled inward to form a certain angle. TKS is self-incompatible, and its seeds are spindle-shaped achene and show upright plumpness. A large amount of laticifer cells and rubber particles can be detected from many TKS tissues, and dry roots of TKS contain high contents of natural rubber. Laticifer cells and rubber particles can only be examined in the vein, stem, and roots of TKS. Our statical results also revealed that the numbers of laticifer cells and rubber particles have a positive relationship with the rubber content in TKS roots. These morphological features can help us to easily distinguish TKS from common dandelion and approximately estimate the rubber content in the roots of different TKS varieties for TKS breeding in future.

Plant Antimicrobial Peptides: State of the Art, In Silico Prediction and Perspectives in the Omics Era

Even before the perception or interaction with pathogens, plants rely on constitutively guardian molecules, often specific to tissue or stage, with further expression after contact with the pathogen. These guardians include small molecules as antimicrobial peptides (AMPs), generally cysteine-rich, functioning to prevent pathogen establishment. Some of these AMPs are shared among eukaryotes (eg, defensins and cyclotides), others are plant specific (eg, snakins), while some are specific to certain plant families (such as heveins). When compared with other organisms, plants tend to present a higher amount of AMP isoforms due to gene duplications or polyploidy, an occurrence possibly also associated with the sessile habit of plants, which prevents them from evading biotic and environmental stresses. Therefore, plants arise as a rich resource for new AMPs. As these molecules are difficult to retrieve from databases using simple sequence alignments, a description of their characteristics and in silico (bioinformatics) approaches used to retrieve them is provided, considering resources and databases available. The possibilities and applications based on tools versus database approaches are considerable and have been so far underestimated.

Identification and Characterization of Glycoproteins and Their Responsive Patterns upon Ethylene Stimulation in the Rubber Latex

Natural rubber is an important industrial material, which is obtained from the only commercially cultivated rubber tree, Hevea brasiliensis. In rubber latex production, ethylene has been extensively used as a stimulant. Recent research showed that post-translational modifications (PTMs) of latex proteins, such as phosphorylation, glycosylation and ubiquitination, are crucial in natural rubber biosynthesis. In this study, comparative proteomics was performed to identify the glycosylated proteins in rubber latex treated with ethylene for different days. Combined with Pro-Q Glycoprotein gel staining and mass spectrometry techniques, we provided the first visual profiling of glycoproteomics of rubber latex and finally identified 144 glycosylated protein species, including 65 differentially accumulated proteins (DAPs) after treating with ethylene for three and/or five days. Gene Ontology (GO) functional annotation showed that these ethylene-responsive glycoproteins are mainly involved in cell parts, membrane components and metabolism. Pathway analysis demonstrated that these glycosylated rubber latex proteins are mainly involved in carbohydrate metabolism, energy metabolism, degradation function and cellular processes in rubber latex metabolism. Protein-protein interaction analysis revealed that these DAPs are mainly centered on acetyl-CoA acetyltransferase and hydroxymethylglutaryl-CoA synthase (HMGS) in the mevalonate pathway for natural rubber biosynthesis. In our glycoproteomics, three protein isoforms of HMGS2 were identified from rubber latex, and only one HMGS2 isoform was sharply increased in rubber latex by ethylene treatment for five days. Furthermore, the HbHMGS2 gene was over-expressed in a model rubber-producing grass Taraxacum Kok-saghyz and rubber content in the roots of transgenic rubber grass was significantly increased over that in the wild type plant, indicating HMGS2 is the key component for natural rubber production.

Natural rubber biosynthesis in plants, the rubber transferase complex, and metabolic engineering progress and prospects

Natural rubber (NR) is a nonfungible and valuable biopolymer, used to manufacture ~50 000 rubber products, including tires and medical gloves. Current production of NR is derived entirely from the para rubber tree (Hevea brasiliensis). The increasing demand for NR, coupled with limitations and vulnerability of H. brasiliensis production systems, has induced increasing interest among scientists and companies in potential alternative NR crops. Genetic/metabolic pathway engineering approaches, to generate NR-enriched genotypes of alternative NR plants, are of great importance. However, although our knowledge of rubber biochemistry has significantly advanced, our current understanding of NR biosynthesis, the biosynthetic machinery and the molecular mechanisms involved remains incomplete. Two spatially separated metabolic pathways provide precursors for NR biosynthesis in plants and their genes and enzymes/complexes are quite well understood. In contrast, understanding of the proteins and genes involved in the final step(s)-the synthesis of the high molecular weight rubber polymer itself-is only now beginning to emerge. In this review, we provide a critical evaluation of recent research developments in NR biosynthesis, in vitro reconstitution, and the genetic and metabolic pathway engineering advances intended to improve NR content in plants, including H. brasiliensis, two other prospective alternative rubber crops, namely the rubber dandelion and guayule, and model species, such as lettuce. We describe a new model of the rubber transferase complex, which integrates these developments. In addition, we highlight the current challenges in NR biosynthesis research and future perspectives on metabolic pathway engineering of NR to speed alternative rubber crop commercial development.

Proteomic Landscape Has Revealed Small Rubber Particles Are Crucial Rubber Biosynthetic Machines for Ethylene-Stimulation in Natural Rubber Production

Rubber particles are a specific organelle for natural rubber biosynthesis (NRB) and storage. Ethylene can significantly improve rubber latex production by increasing the generation of small rubber particles (SRPs), regulating protein accumulation, and activating many enzyme activities. We conducted a quantitative proteomics study of different SRPs upon ethylene stimulation by differential in-gel electrophoresis (DIGE) and using isobaric tags for relative and absolute quantification (iTRAQ) methods. In DIGE, 79 differentially accumulated proteins (DAPs) were determined as ethylene responsive proteins. Our results show that the abundance of many NRB-related proteins has been sharply induced upon ethylene stimulation. Among them, 23 proteins were identified as rubber elongation factor (REF) and small rubber particle protein (SRPP) family members, including 16 REF and 7 SRPP isoforms. Then, 138 unique phosphorylated peptides, containing 129 phosphorylated amino acids from the 64 REF/SRPP family members, were identified, and most serine and threonine were phosphorylated. Furthermore, we identified 226 DAPs from more than 2000 SRP proteins by iTRAQ. Integrative analysis revealed that almost all NRB-related proteins can be detected in SRPs, and many proteins are positively responsive to ethylene stimulation. These results indicate that ethylene may stimulate latex production by regulating the accumulation of some key proteins. The phosphorylation modification of REF and SRPP isoforms might be crucial for NRB, and SRP may act as a complex natural rubber biosynthetic machine.

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.

Proteomic Landscape of the Mature Roots in a Rubber-Producing Grass Taraxacum Kok-saghyz

The rubber grass Taraxacum kok-saghyz (TKS) contains large amounts of natural rubber (cis-1,4-polyisoprene) in its enlarged roots and it is an alternative crop source of natural rubber. Natural rubber biosynthesis (NRB) and storage in the mature roots of TKS is a cascade process involving many genes, proteins and their cofactors. The TKS genome has just been annotated and many NRB-related genes have been determined. However, there is limited knowledge about the protein regulation mechanism for NRB in TKS roots. We identified 371 protein species from the mature roots of TKS by combining two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). Meanwhile, a large-scale shotgun analysis of proteins in TKS roots at the enlargement stage was performed, and 3545 individual proteins were determined. Subsequently, all identified proteins from 2-DE gel and shotgun MS in TKS roots were subject to gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and most proteins were involved in carbon metabolic process with catalytic activity in membrane-bounded organelles, followed by proteins with binding ability, transportation and phenylpropanoid biosynthesis activities. Fifty-eight NRB-related proteins, including eight small rubber particle protein (SRPP) and two rubber elongation factor(REF) members, were identified from the TKS roots, and these proteins were involved in both mevalonate acid (MVA) and methylerythritol phosphate (MEP) pathways. To our best knowledge, it is the first high-resolution draft proteome map of the mature TKS roots. Our proteomics of TKS roots revealed both MVA and MEP pathways are important for NRB, and SRPP might be more important than REF for NRB in TKS roots. These findings would not only deepen our understanding of the TKS root proteome, but also provide new evidence on the roles of these NRB-related proteins in the mature TKS roots.

The formation and accumulation of protein-networks by physical interactions in the rapid occlusion of laticifer cells in rubber tree undergoing successive mechanical wounding

BACKGROUND

Although the wound response of plants has been extensively studied, little is known of the rapid occlusion of wounded cell itself. The laticifer in rubber tree is a specific type of tissue for natural rubber biosynthesis and storage. In natural rubber production, tapping is used to harvest the latex which flows out from the severed laticifer in the bark. Therefore, study of the rapid wound-occlusion of severed laticifer cells is important for understanding the rubber tree being protected from the continuously mechanical wounding.

RESULTS

Using cytological and biochemical techniques, we revealed a biochemical mechanism for the rapid occlusion of severed laticifer cells. A protein-network appeared rapidly after tapping and accumulated gradually along with the latex loss at the severed site of laticifer cells. Triple immunofluorescence histochemical localization showed that the primary components of the protein-network were chitinase, β-1,3-glucanase and hevein together with pro-hevein (ProH) and its carboxyl-terminal part. Molecular sieve chromatography showed that the physical interactions among these proteins occurred under the condition of neutral pH. The interaction of β-1,3-glucanase respectively with hevein, chitinase and ProH was testified by surface plasmon resonance (SPR). The interaction between actin and β-1,3-glucanase out of the protein inclusions of lutoids was revealed by pull-down. This interaction was pharmacologically verified by cytochalasin B-caused significant prolongation of the duration of latex flow in the field.

CONCLUSIONS

The formation of protein-network by interactions of the proteins with anti-pathogen activity released from lutoids and accumulation of protein-network by binding to the cytoskeleton are crucial for the rapid occlusion of laticifer cells in rubber tree. The protein-network at the wounded site of laticifer cells provides not only a physical barrier but also a biochemical barrier to protect the wounded laticifer cells from pathogen invasion.

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.

Subcellular proteome profiles of different latex fractions revealed washed solutions from rubber particles contain crucial enzymes for natural rubber biosynthesis

Rubber particle (RP) is a specific organelle for natural rubber biosynthesis (NRB) and storage in rubber tree Hevea brasiliensis. NRB is processed by RP membrane-localized proteins, which were traditionally purified by repeated washing. However, we noticed many proteins in the discarded washing solutions (WS) from RP. Here, we compared the proteome profiles of WS, C-serum (CS) and RP by 2-DE, and identified 233 abundant proteins from WS by mass spectrometry. Many spots on 2-DE gels were identified as different protein species. We further performed shotgun analysis of CS, WS and RP and identified 1837, 1799 and 1020 unique proteins, respectively. Together with 2-DE, we finally identified 1825 proteins from WS, 246 were WS-specific. These WS-specific proteins were annotated in Gene Ontology, indicating most abundant pathways are organic substance metabolic process, protein degradation, primary metabolic process, and energy metabolism. Protein-protein interaction analysis revealed these WS-specific proteins are mainly involved in ribosomal metabolism, proteasome system, vacuolar protein sorting and endocytosis. Label free and Western blotting revealed many WS-specific proteins and protein complexes are crucial for NRB initiation. These findings not only deepen our understanding of WS proteome, but also provide new evidences on the roles of RP membrane proteins in NRB.

SIGNIFICANCE

Natural rubber is stored in rubber particle from the rubber tree. Rubber particles were traditionally purified by repeated washing, but many proteins were identified from the washing solutions (WS). We obtained the first visualization proteome profiles with 1825 proteins from WS, including 246 WS-specific ones. These WS proteins contain almost all enzymes for polyisoprene initiation and may play important roles in rubber biosynthesis.

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.

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.

Expression Profiling of Mitogen-Activated Protein Kinase Genes Reveals Their Evolutionary and Functional Diversity in Different Rubber Tree (Hevea brasiliensis) Cultivars

Rubber tree (Heveabrasiliensis) is the only commercially cultivated plant for producing natural rubber, one of the most essential industrial raw materials. Knowledge of the evolutionary and functional characteristics of kinases in H. brasiliensis is limited because of the long growth period and lack of well annotated genome information. Here, we reported mitogen-activated protein kinases in H. brasiliensis (HbMPKs) by manually checking and correcting the rubber tree genome. Of the 20 identified HbMPKs, four members were validated by proteomic data. Protein motif and phylogenetic analyses classified these members into four known groups comprising Thr-Glu-Tyr (TEY) and Thr-Asp-Tyr (TDY) domains, respectively. Evolutionary and syntenic analyses suggested four duplication events: HbMPK3/HbMPK6, HbMPK8/HbMPK9/HbMPK15, HbMPK10/HbMPK12 and HbMPK11/HbMPK16/HbMPK19. Expression profiling of the identified HbMPKs in roots, stems, leaves and latex obtained from three cultivars with different latex yield ability revealed tissue- and variety-expression specificity of HbMPK paralogues. Gene expression patterns under osmotic, oxidative, salt and cold stresses, combined with cis-element distribution analyses, indicated different regulation patterns of HbMPK paralogues. Further, Ka/Ks and Tajima analyses suggested an accelerated evolutionary rate in paralogues HbMPK10/12. These results revealed HbMPKs have diverse functions in natural rubber biosynthesis, and highlighted the potential possibility of using MPKs to improve stress tolerance in future rubber tree breeding.

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.

Identification and functional characterization of HbOsmotin from Hevea brasiliensis

Latex in the laticiferous cell network of Hevea brasiliensis tree is composed of cytoplasm that synthesizes natural rubber. Ethylene stimulation of the tree bark enhances latex production partly by prolonging the duration of latex flow during the tapping process. Here, we identified an osmotin-like cDNA sequence (HbOsmotin) from H. brasiliensis that belongs to the pathogenesis-related 5 (PR-5) gene family. The HbOsmotin protein is present in the lutoids of latex in H. brasiliensis, whereas in onion epidermal cells, this protein is predominantly distributed around the cell wall, suggesting that it may be secreted from the cytoplasm. We investigated the effects of exogenous ethylene on HbOsmotin transcription and protein accumulation in rubber latex, and further determined the protein function after osmotic stress in Arabidopsis. In regularly tapped trees, HbOsmotin expression was drastically inhibited in rubber latex after tapping, although the expression was subsequently recovered by ethylene stimulation. However, in virgin plants that had never been tapped, exogenous ethylene application slightly decreased HbOsmotin expression. HbOsmotin overexpression in Arabidopsis showed that HbOsmotin reduced the osmotic stress tolerance of the plant, which likely occurred by raising the water potential. These data indicated that HbOsmotin may contribute to osmotic regulation in laticiferous cells.

Transcriptome and proteome analyses provide insight into laticifer's defense of Euphorbia tirucalli against pests

The cytoplasm of laticifers, which are plant cells specialized for rubber production and defense against microbes and herbivores, is a latex. Although laticifers share common functions, the protein constituents of latexes are highly variable among plant species and even among organs. In this study, transcriptomic and proteomic analyses of Euphorbia tirucalli's (Euphorbiaceae) latex were conducted to determine the molecular basis of the laticifer's functions in this plant. The hybrid de novo assembly of Illumina mRNA-seq and expressed sequence tags obtained by Sanger's sequencing revealed 26,447 unigenes. A unigene similar to Arabidopsis embryo-specific protein 3 (AT5G62200), which is a PLAT domain-containing protein, and rubber elongation factor showed the highest expression levels. The proteome analysis, studied by liquid chromatography-mass spectrometry with the de novo assembled unigenes as the database, revealed 161 proteins in the latex, 107 of which were not detected in the stem. A gene ontology analysis indicated that the laticifer's proteome was enriched with proteins related to proteolysis, phosphatase, defense against various environmental stresses and lipid metabolisms. D-mannose-binding lectin, ricin (which lacked the N-terminal conserved ribosome-inactivating protein domain), chitinase and peroxidase were highly accumulated, as confirmed by two-dimensional polyacrylamide gel electrophoresis. Thus, the lectins and chitinase may be the major defensive proteins against pests, and the other defense-related proteins and transcripts detected in latex may work in coordination with them. Highly expressing unigenes with unknown functions are candidate novel defense- or rubber production-related genes.

A protein extraction method for low protein concentration solutions compatible with the proteomic analysis of rubber particles

The extraction of high-purity proteins from the washing solution (WS) of rubber particles (also termed latex-producing organelles) from laticifer cells in rubber tree for proteomic analysis is challenging due to the low concentration of proteins in the WS. Recent studies have revealed that proteins in the WS might play crucial roles in natural rubber biosynthesis. To further examine the involvement of these proteins in natural rubber biosynthesis, we designed an efficiency method to extract high-purity WS proteins. We improved our current borax and phenol-based method by adding reextraction steps with phenol (REP) to improve the yield from low protein concentration samples. With this new method, we extracted WS proteins that were suitable for proteomics. Indeed, compared to the original borax and phenol-based method, the REP method improved both the quality and quantity of isolated proteins. By repeatedly extracting from low protein concentration solutions using the same small amount of phenol, the REP method yielded enough protein of sufficiently high-quality from starting samples containing less than 0.02 mg of proteins per milliliter. This method was successfully applied to extract the rubber particle proteins from the WS of natural rubber latex samples. The REP-extracted WS proteins were resolved by 2DE, and 28 proteins were positively identified by MS. This method has the potential to become widely used for the extraction of proteins from low protein concentration solutions for proteomic analysis.

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.

Gene expression analysis and SNP/InDel discovery to investigate yield heterosis of two rubber tree F1 hybrids

As an important industrial material, natural rubber is mainly harvested from the rubber tree. Rubber tree breeding is inefficient, expensive and time-consuming, whereas marker-assisted selection is a feasible method for early selection of high-yield hybrids. We thus sequenced and analyzed the transcriptomes of two parent rubber trees (RRIM 600 and PR 107) and their most productive hybrids (RY 7-33-97 and RY 7-20-59) to understand their gene expression patterns and genetic variations including single nucleotide polymorphisms (SNPs) and small insertions/deletions (InDels). We discovered >31,000 genetic variations in 112,702 assembled unigenes. Our results showed that the higher yield in F₁ hybrids was positively associated with their higher genome heterozygosity, which was further confirmed by genotyping 10 SNPs in 20 other varieties. We also showed that RY 7-33-97 and RY 7-20-59 were genetically closer to RRIM 600 and PR 107, respectively, in agreement with both their phenotypic similarities and gene expression profiles. After identifying ethylene- and jasmonic acid–responsive genes at the transcription level, we compared and analyzed the genetic variations underlying rubber biosynthesis and the jasmonic acid and ethylene pathways in detail. Our results suggest that genome-wide genetic variations play a substantive role in maintaining rubber tree heterosis.

Genome-Wide Identification of Jatropha curcas Aquaporin Genes and the Comparative Analysis Provides Insights into the Gene Family Expansion and Evolution in Hevea brasiliensis

Aquaporins (AQPs) are channel-forming integral membrane proteins that transport water and other small solutes across biological membranes. Despite the vital role of AQPs, to date, little is known in physic nut (Jatropha curcas L., Euphorbiaceae), an important non-edible oilseed crop with great potential for the production of biodiesel. In this study, 32 AQP genes were identified from the physic nut genome and the family number is relatively small in comparison to 51 in another Euphorbiaceae plant, rubber tree (Hevea brasiliensis Muell. Arg.). Based on the phylogenetic analysis, the JcAQPs were assigned to five subfamilies, i.e., nine plasma membrane intrinsic proteins (PIPs), nine tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and four small basic intrinsic proteins (SIPs). Like rubber tree and other plant species, functional prediction based on the aromatic/arginine selectivity filter, Froger's positions, and specificity-determining positions showed a remarkable difference in substrate specificity among subfamilies of JcAQPs. Genome-wide comparative analysis revealed the specific expansion of PIP and TIP subfamilies in rubber tree and the specific gene loss of the XIP subfamily in physic nut. Furthermore, by analyzing deep transcriptome sequencing data, the expression evolution especially the expression divergence of duplicated HbAQP genes was also investigated and discussed. Results obtained from this study not only provide valuable information for future functional analysis and utilization of Jc/HbAQP genes, but also provide a useful reference to survey the gene family expansion and evolution in Euphorbiaceae plants and other plant species.

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.

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.

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

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

Comprehensive Proteomics Analysis of Laticifer Latex Reveals New Insights into Ethylene Stimulation of Natural Rubber Production

Ethylene is a stimulant to increase natural rubber latex. After ethylene application, both fresh yield and dry matter of latex are substantially improved. Moreover, we found that ethylene improves the generation of small rubber particles. However, most genes involved in rubber biosynthesis are inhibited by exogenous ethylene. Therefore, we conducted a proteomics analysis of ethylene-stimulated rubber latex, and identified 287 abundant proteins as well as 143 ethylene responsive latex proteins (ERLPs) with mass spectrometry from the 2-DE and DIGE gels, respectively. In addition, more than 1,600 proteins, including 404 ERLPs, were identified by iTRAQ. Functional classification of ERLPs revealed that enzymes involved in post-translational modification, carbohydrate metabolism, hydrolase activity, and kinase activity were overrepresented. Some enzymes for rubber particle aggregation were inhibited to prolong latex flow, and thus finally improved latex production. Phosphoproteomics analysis identified 59 differential phosphoproteins; notably, specific isoforms of rubber elongation factor and small rubber particle protein that were phosphorylated mainly at serine residues. This post-translational modification and isoform-specific phosphorylation might be important for ethylene-stimulated latex production. These results not only deepen our understanding of the rubber latex proteome but also provide new insights into the use of ethylene to stimulate rubber latex production.

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

BACKGROUND

Rubber tree (Hevea brasiliensis Muell. Arg.) is the primarily commercial source of natural rubber in the world. Latex regeneration and duration of latex flow after tapping are the two factors that determine rubber yield of rubber tree, and exhibit a huge variation between rubber tree clones CATAS8-79 and PR107.

RESULTS

To dissect the molecular mechanism for the regulation of latex regeneration and duration of latex flow, we sequenced and comparatively analyzed latex of rubber tree clone CATAS8-79 and PR107 at transriptome level. More than 26 million clean reads were generated in each pool and 51,829 all-unigenes were totally assembled. A total of 6,726 unigenes with differential expression patterns were detected between CATAS8-79 and PR107. Functional analysis showed that genes related to mass of categories were differentially enriched between the two clones. Expression pattern of genes which were involved in latex regeneration and duration of latex flow upon successive tapping was analyzed by quantitative PCR. Several genes related to rubber biosynthesis, cellulose and lignin biosynthesis and rubber particle aggregation were differentially expressed between CATAS8-79 and PR107.

CONCLUSIONS

This is the first report about probing latex regeneration and duration of latex flow by comparative transcriptome analysis. Among all the suggested factors, it is more important that the level of endogenous jasmonates, carbohydrate metabolism, hydroxymethylglutaryl-CoA reductase (HMGR) and Hevea rubber transferase (HRT) in mevalonate (MVA) parthway for latex regeneration while the level of endogenous ethylene (ETH), lignin content of laticifer cell wall, antioxidants and glucanases for the duration of latex flow. These data will provide new cues for understanding the molecular mechanism for the regulation of latex regeneration and duration of latex flow in rubber tree.

SNP identification from RNA sequencing and linkage map construction of rubber tree for anchoring the draft genome

Hevea brasiliensis, or rubber tree, is an important crop species that accounts for the majority of natural latex production. The rubber tree nuclear genome consists of 18 chromosomes and is roughly 2.15 Gb. The current rubber tree reference genome assembly consists of 1,150,326 scaffolds ranging from 200 to 531,465 bp and totalling 1.1 Gb. Only 143 scaffolds, totalling 7.6 Mb, have been placed into linkage groups. We have performed RNA-seq on 6 varieties of rubber tree to identify SNPs and InDels and used this information to perform target sequence enrichment and high throughput sequencing to genotype a set of SNPs in 149 rubber tree offspring from a cross between RRIM 600 and RRII 105 rubber tree varieties. We used this information to generate a linkage map allowing for the anchoring of 24,424 contigs from 3,009 scaffolds, totalling 115 Mb or 10.4% of the published sequence, into 18 linkage groups. Each linkage group contains between 319 and 1367 SNPs, or 60 to 194 non-redundant marker positions, and ranges from 156 to 336 cM in length. This linkage map includes 20,143 of the 69,300 predicted genes from rubber tree and will be useful for mapping studies and improving the reference genome assembly.