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

Abstract

[...]rubber production of the Para rubber tree has reached the limit due to the strict climatic requirements for its planting areas [10], little genetic variability, labor cost, and threats of fatal fungal plant diseases [1]. [...]it is critical to find an alternative source and a model plant for NR production. Morphological Observation of Rubber Particles in the Main Roots of TKS Natural rubber in TKS is produced from the latex system of its enlarged roots, and the biomass of the main roots is important for rubber yield. [...]the mature 6-month TKS plants with flowers were selected to perform morphological analysis, and the enlarged main roots are highlighted (Figure 1A). Establishment of a Visual Two-Dimensional Gel Electrophoresis (2-DE) Proteome Map and Mass Spectrometry (MS) Identification of High Abundance Proteins in Taraxacum Kok-saghyz (TKS) Roots To further obtain the protein accumulation profiling, the total proteins from the 6M TKS roots were extracted and performed 2-DE to separate the total proteins. [...]a high-resolution 2-DE gel was obtained and a visual proteome map was established (Figure 2).

1. Introduction

Natural rubber (NR, cis-1,4-polyisoprene) is a biopolymer with high economic value and it is widely used as a strategic raw material in more than 40,000 products [1,2]. More than 2500 plant species can biosynthesize NR [3,4]. However, high-quality natural rubber in viable quantities is only observed in a few plant species, such as the Para rubber tree Hevea brasiliensis [5,6], Russian dandelion Taraxacum kok-saghyz Rodin (TKS, also called the rubber grass) or its close species Taraxacum brevicorniculatum [7], Eucommia ulmoides Oliver and a guayule shrub Parthenium argentatum [8,9]. Currently, the Para rubber tree is nearly the only commercial plant to cultivate the exclusive source of NR [5,6]. However, rubber production of the Para rubber tree has reached the limit due to the strict climatic requirements for its planting areas [10], little genetic variability, labor cost, and threats of fatal fungal plant diseases [1]. Therefore, it is critical to find an alternative source and a model plant for NR production.

TKS has drawn special attention since the 1940s. It belongs to the Compositae family and originates from the Tekes River Basin near Tianshan Mountain border between Kazakhstan and China [11]. TKS root contains NR ranging from 3% to 28% of total dry weight [1] and its rubber property and molecular characteristics are similar with NR in the Para rubber tree [12,13]. TKS grows widely in temperate and cold areas, and its root also contains about 28% of inulin, which is an important material for bioethanol and the food industry [14,15]. With its advantages of high rubber content and quality, wide planting area, a relatively short life cycle, ease of transformation and harvesting method [14,15], TKS is becoming a promising crop for natural rubber production. In addition, TKS is a perennial herb with a relatively simple genome (1.29 Gb), containing 46,731 protein-coding genes [16]. It is easily manipulated for genetic transformation and can be used as an ideal model plant for studying the rubber biosynthesis mechanism [17].

NR is biosynthesized by two pathways: the mevalonate acid (MVA) pathway [18,19] and the methylerythritol phosphate (MEP) pathway [20]. Recently, to determine the protein-based regulation mechanism of NRB, several proteomics have been conducted on the total latex [21,22,23,24,25] and rubber particles [26,27,28,29,30,31] of the rubber tree H. brasiliensis. But no proteomic analysis has been performed on TKS, and only one study on determining proteins from the rubber phase and pellet phase of the latex from T. brevicorniculatum (TBR) was reported and 278 unique proteins from the one- and two-dimensional gel electrophoresis (1-DE and 2-DE) gels were identified [32]. In order to gain further insight into the protein-based regulation mechanism of rubber biosynthesis, we performed a comprehensive proteomics analysis of the mature roots of TKS at six months (6M) by 2-DE and mass spectrometry (MS), and then identified thousands of proteins by a large-scale shotgun method. These results may deepen our understanding of the root proteome and provide valuable gene candidates for genetic improvement of TKS breeding as a commercial alternative crop.

2. Results 2.1. Morphological Observation of Rubber Particles in the Main Roots of TKS

Natural rubber in TKS is produced from the latex system of its enlarged roots, and the biomass of the main roots is important for rubber yield. Therefore, the mature 6-month TKS plants with flowers were selected to perform morphological analysis, and the enlarged main roots are highlighted (Figure 1A). The light microscope results revealed that the concentric rings of laticifers as specialized tubular vessels in roots can be detected, and the total number of laticifer vessels is about 137 ± 32 (n = 20) in the main roots of 6M TKS. The brown rubber latex granules, known as the specific cytoplasm of laticifer cells, can be clearly detected in the areas of endodermis and cortex in phloem tissues from both the crosscutting and slitting slices of the main roots. All the laticifer cells were considered to originate from the initial cells of vascular cambium to secondary laticifers in paralleling rings (Figure 1B), and most of these laticifer cells seem to connect into tubular vessels in the main roots (Figure 1C). Ultrastructural investigation of rubber particles in TKS root at 6 months was conducted under a transmission electron microscope (TEM). Many spherical or ovoid-ellipsoid rubber particles, ranging from 50 nm to 2000 nm, were examined in the laticifer cells (Figure 1D). We noticed that most rubber particles have a diameter less than 200 nm; they are termed small rubber particles (SRPs). However, rubber particles with diameter larger than 400 nm are fewer in number; they are traditionally called large rubber particles (LRPs). These particles, including SRPs and LRPs, are surrounded by a monolayer membrane (Figure 1D), in which many enzymes related to NRB are anchored or combined with each other [33].

2.2. Establishment of a Visual Two-Dimensional Gel Electrophoresis (2-DE) Proteome Map and Mass Spectrometry (MS) Identification of High Abundance Proteins in Taraxacum Kok-saghyz (TKS) Roots

To further obtain the protein accumulation profiling, the total proteins from the 6M TKS roots were extracted and performed 2-DE to separate the total proteins. Finally, a high-resolution 2-DE gel was obtained and a visual proteome map was established (Figure 2). In this reference 2-DE gel, 428 ± 15 protein spots, ranged from pH 4 to pH 7, were detected from three biological replicates (Figure S1). We manually excised all abundant protein spots from the 2-DE gel, and finally identified 371 protein species by MALDI TOF/TOF MS. These proteins occupied 84.38% volume of all the detected spots in the 2-DE gel and belonged to 231 gene products or named unique proteins (Table S1).

We compared the theoretical and experimental ratios of molecular weight (Mr) and isoelectric point (pI) of the 371 proteins identified from the 2-DE gel, and presented their ratios as radar axis labels (Mr for radial value; pI ratio for annular value) in radial chart (Figure 3A). Most proteins showed similar radial and annular values, and their ratios for theoretical and experimental Mr and pI are near line 1.0. However, some proteins demonstrated differentially experimental Mr and pI values in the 2-DE gel (Figure 2). In the proteomics study based on 2-DE gel, they are called different protein isoforms or protein species. In this study, 90 unique proteins were identified from at least two different spots, and these proteins contained 229 spots in the 2-DE gel. Among them, 56 proteins were identified from 112 spots, and they contained two protein isoforms in the 2-DE gel. Twenty-two proteins had three isoforms, and nine proteins were determined from 36 spots with four isoforms (Table S2). There are three unique proteins containing five isoforms; they are aldolase (spots 305, 313, 314, 403 and 405), enolase (spots 94, 100, 106, 150 and 158) and cytosolic glyceraldehyde-3-phosphate dehydrogenase (GAPC, spots 12, 163, 176, 272 and 275).

It is noteworthy that these protein isoforms show the same theoretical Mr and pI values. Most of their Mr values are ranged from 20 kDa to 60 kDa, and their pI points are mainly distributed from pH 5.0 to 7.0 (Figure 3B). However, their experimental Mr and pI values are different to each other, and most protein isoforms showed a different pI value in the 2-DE gel (Figure 1). These results indicated that these proteins could have post-modifications that resulted in protein forms/species in the 2-DE gel-based proteomics.

We further determined the protein abundance in the 2-DE gel and the most abundant 20 protein species are highlighted (Figure 3C). The results demonstrated that, based on the spot volume, the top-20 protein spots occupied almost a half of abundance volume (43.84%) for all proteins in the 2–DE gel of 6M TKS roots. The most abundant protein spots were identified as the protein species or isoforms of catalase (spot 355, 6.57%; spot 2, 3.38%), ferritin (spot 357, 3.77%; spot 359, 2.90%), actin (spots 142 and 354), actin (spots 142 and 354), GAPC2 (spots 4 and 118), and the subunit of ATPase (spots 11 and 344). Among them, the most abundant spot was identified as catalase-like isoform X1 (spot 355), followed by allene oxide cyclase (spot 346) and eukaryotic translation initiation factor 5A (eIF-5A, spot 361) (Table S1).

2.3. High-Throughput Shotgun Proteomic Analysis of the 6M Roots of TKS

As proteins with very low abundance or out of the pH range of immobilized pH gradient (IPG) strips are difficult to separate and identify by 2-DE-based proteomic method, a shotgun analysis of the proteins was further performed. The proteins were extracted from the 6M TKS roots, and three biological repeats (R1, R2 and R3) were conducted to obtain a comprehensive proteomic profile. A total of 5205, 5323 and 5654 unique proteins were successfully identified from R1, R2 and R3, respectively (Table S3). Among them, 7481 proteins were identified from at least one independent shotgun experiment, and 5156 proteins were identified from at least two experiments. There were 3545 shared proteins in the three experiments (Figure 4A), and only these shared proteins from three experiments were considered as the 6M TKS root proteins produced by the shotgun proteomic method in the following study. We found 184 shared proteins in the 2-DE and shotgun proteomic methods. There were 3361 and 47 specific proteins in the shotgun and 2-DE methods, respectively. These results indicated that, although the shotgun method can produce large amounts of proteins, the traditional 2-DE gel-based proteomics also generates a few specific proteins from the 6M TKS roots (Figure 4B; Table S3). Our special interest is focusing on the potential biological functions of these shared proteins in TKS roots.

2.4. Pathway Analysis of the Identified Proteins in 6M TKS Roots

To gain insight into the functional categories of the 3592 proteins (including the 231 and 3545 proteins identified by 2-DE and shotgun, respectively), we performed gene ontology (GO) classification and the enriched outputs of their biological process, cellular component and molecular function are presented (Figure 5). The GO enrichment revealed that many proteins were localized in membrane-bounded organelle, and they were mainly involved in metabolic process with catalytic activity or binding ability (Figure S3; Table S4). At the biological process level, 13 main processes were detected. Among them, the largest amount containing 1520 proteins were involved in the metabolic process, followed by 1033 proteins involved in the cellular process and 695 proteins in the single-organism process. A total of 302 proteins were considered to respond to external stimulus. At subcellular level, nine components were observed. Among them, almost half including 866 proteins were localized to the cell part. Many proteins were localized into organelle membrane and macromolecular complex. In molecular function classification, seven pathways were determined. Among them, the most portion including 1688 proteins were taken part into catalytic activity, followed by 1325 proteins with binding ability. There are 107 proteins with transporter activity and 101 proteins with structural molecule activity (Figure 5; Table S4). These different GO term distributions in TKS roots are related to their biological functions and indicate some important biological processes for secondary metabolite in the 6M TKS roots.

To further investigate the biological functions of the identified 3592 proteins, KEGG pathway analysis was performed using the BLAST2GO program. These proteins were clustered into 19 main pathways (Figure 6A), including translation process (616 proteins), carbohydrate metabolism (583 proteins), folding, sorting and degradation (440 proteins), global and overview metabolism (409 proteins), amino acid metabolism (380 proteins), transport and catabolism (342 proteins), lipid metabolism (308 proteins), etc. These results demonstrated that the main proteins in TKS roots were involved in translation and carbohydrate metabolism. Furthermore, 129 sub-pathways were enriched in KEGG annotation (Table S5), and the top 20 sub-pathways are highlighted (Figure 6B). Among them, the most abundant pathway containing 240 proteins is carbon metabolism. The second contains 224 proteins with ribosome activity. The third sub-pathway has 213 proteins for biosynthesis of amino acids, followed by 183 proteins for endocytosis, 182 proteins in the spliceosome, and 172 proteins in the endoplasmic reticulum. Proteins involved in starch and sucrose metabolism (164 proteins), glycolysis (121 proteins), phenylpropanoid biosynthesis (119 proteins), and amino sugar and nucleotide sugar metabolism (116 proteins), are also enriched in these KEGG sub-pathways (Table S5). These proteins are crucial for biosynthesis of isoprenoids in many plants, especially in natural rubber-producing plants [16].

2.5. The Proteins Involved in Natural Rubber Biosynthesis in TKS Roots

We checked the identified 3592 proteins from the 6M TKS roots and paid special attention to the proteins and enzymes involved in NRB. Our proteomic data demonstrated that a total of 58 unique proteins or gene products were identified at least one time by the three shotgun experiments (Table S6), and four of them, named ACAT7 (spots 157, 183 and 424), HMGR1 (spot 146), DXS3 (spot 257) and DXS8 (spots 199 and 204), were also identified from the 2-DE gel of 6M TKS roots (Figure 1; Table S1). Twenty kinds of proteins were determined from the 58 unique gene products (Table S6), and these proteins are involved in both MVA and MEP pathways for NBR in TKS roots (Figure 7). Genomic sequencing data demonstrated that both cytosolic MVA and plastidic MEP pathways are present for NBR, and a total of 102 NBR-related genes have been determined in the genome of TKS [16]. Our proteomic data demonstrated 22 unique proteins were involved in the MVA pathway and 13 proteins in MEP pathway in TKS roots (Figure 7; Table S6). In the MVA pathway, the identified 22 unique proteins belonged to six kinds of essential enzymes for rubber biosynthesis; they are ACAT, HMGS, HMGR, MVK, PMVK and MVD. Among them, seven members for ACAT (from ACAT2 to ACAT8), six members for mevalonate kinase (MVK 1, 2 and 7–11), and five members for phosphomevalonate kinase (PMVK 1 and 2–6), were positively identified from TKS roots by the shotgun method (Table S6). Our proteomics results also showed that, by contrast with rubber tree latex, TKS contains the MEP pathway for NRB in its roots. In the MEP pathway, seven kinds of proteins were identified from 13 unique proteins, including six family members for deoxyxylulose-5-phosphate synthase (DXS 1–3 and 8–10) and two members for D-xylulose 5-phosphate reductoisomerase (DXR 1 and 2).

In the following step for initiation of the synthesis of isopentenyl pyrophosphate (IPP) to form rubber molecule polymer, which is termed initiator synthesis, four crucial enzymes, including isopentenyl-diphosphate delta-isomerase (IPI 2–4), geranylgeranyl pyrophosphate synthase (GGPS 1–3), geranylgeranyl diphosphate synthase (GPS 3, 5 and 7) and farnesyl diphosphate synthase (FPS 1), were identified from 10 unique proteins.

In the final rubber elongation process, cis-isoprene transferase (CPT) can help GGPP to generate natural rubber hydrocarbons contain different length of carbons. Four kinds of proteins, named CPT, SRPP, REF, and HRBP, which is a Nogo-B receptor as a HRT1-REF bridging protein, are considered to play crucial roles for natural rubber elongation. In this proteomic study, we identified the above four members from 12 unique proteins. Among them, eight SRPP members (SRPP 1–7 and 9) and two REF members were determined from the 6M TKS roots.

We checked the annotated genomes of the rubber tree H. brasiliensis and the rubber grass T. kok-saghyz, and found 85 and 106 NRB-related proteins from H. brasiliensis and TKS genomes, respectively. It is noteworthy that 58 out of 106 NRB-related unique proteins were identified from the 6M TKS roots, and 13 out of the identified 58 unique proteins were crucial members in the MEP pathway (Table S6). These results revealed that both MVA and MEP pathways are important in controlling NRB and rubber production in the mature TKS roots.

3. Discussion 3.1. The First Visual Proteome Based on 2-DE Gel Demonstrated Many Protein Isoforms in the Mature TKS Roots

The ease and availability of 2-DE combined MS is a classical standard for proteomics study, and it is still a valuable tool to provide a visual proteome for plants [34]. Large amounts of proteomics studies have been performed in many plant species, however, only several proteomic research works were reported on different latex components of the Para rubber tree [2,35], and only one proteomic analysis was performed to identify 278 unique proteins from the rubber particle phase and pellet phase of the T. brevicorniculatum latex [32]. Our 2-DE gel-based proteomics of TKS roots resulted in 371 abundant protein species (Figure 2; Table S1), but the protein profiles are very different to that in T. brevicorniculatum latex. Compared with the 2-DE gels for a clear aqueous phase and a yellowish pellet phase of the fractionated samples from T. brevicorniculatum latex [32], our 2-DE gel for the mature TKS roots is much clearer and contains more protein spots (Figure 2). Among the identified proteins, 90 unique proteins contained different protein isoforms, which have different experimental Mr and pI values on the 2-DE gel (Table S1), probably due to the existence of modification protein variants.

In the past decade, only 1208 proteins were identified from all the rubber-producing plants, and only several proteomics studies were performed on the rubber particles [36]. In a previous high throughput proteomics study, a total of 186 proteins were positively identified from rubber particles by shotgun tandem MS [28]. After removal of the protein bands corresponding to the RP-abundant REF and SRPP, 137 protein species including 115 unique proteins were identified in a LC-MS/MS based proteome analysis [31]. By investigating the protein profile between LRPs and SRPs, 53 protein spots corresponding to 22 gene products were detected to have significant difference. Among them, most up-regulated proteins in SRPs were identified as SRPP, HMGS, phospholipase D, ethylene response factor, eukaryotic translation initiation factor, etc., but the most abundant proteins in LRPs were REF, glucanase and several hypothetical proteins [27]. Interaction network analysis of rubber particle proteomics revealed the formation of the protein complex consisting of HRT1, REF and HRBP might play crucial roles as a NRB machinery [31], and these proteins are associated with the endoplasmic reticulum [37].

Recently, more than 1600 proteins were identified from total latex of the Para rubber tree by using isobaric tags for relative and absolute quantitation (iTRAQ) method [2] and 1839 unique proteins were determined by LC-MS/MS from the whole translated draft genome of the rubber tree H. brasiliensis [24]. Our high throughput shotgun proteomics produced 7481 proteins from at least one independent shotgun experiment and 3545 shared proteins in the three experiments from the 6M TKS roots (Figure 4B; Table S3), which covered almost all the previously identified proteins from the latex-producing plants. To our best knowledge, it is the first visual proteome based on the combination of 2-DE gel and high throughput proteomics methods, and the results may help us to deepen our understanding of the roles of protein isoforms in the mature TKS roots.

3.2. Large-Scale Shotgun Proteomics Landscape Revealed both Mevalonate Acid (MVA) and Methylerythritol Phosphate (MEP) Pathways are Important for Natural Rubber Biosynthesis (NRB) in the Mature TKS Roots

Natural rubber is an isoprenoid polymer that is synthesized on the rubber particles in many latex-producing plants [16]. Although more than 2500 plants can biosynthesize natural rubber [31], the Para rubber tree H. brasiliensis is the only plant commercially cultivated to produce natural rubber for industry [2]. Rubber latex is a kind of specialized cytoplasm of laticifer cells in the bark phloem of the rubber tree, and natural rubber, as an elastomer with physical and chemical properties, cannot be fully matched by synthetic rubber [31]. Annotation of the rubber tree genome revealed that although 22 MEP genes can be identified, only two DXS genes (DXS7 and DXS10) show substantial and preferential expression in the rubber latex; on the other hand, at least one gene for the identified 18 enzymes in MVA pathway shows latex-biased abundant expression [10]. Proteomics analysis also supported that MVA pathway is more important than MEP pathway for NRB in the latex of the Para rubber tree [2,24,25,29].

Similar results were observed in the recently assembled TKS genome, which contains 102 rubber biosynthesis-related genes. Among them, 40 genes in the MVA pathway and 23 genes in the MEP pathway were identified, respectively. Comparison of gene expression level proved that, for each process in the MVA pathway in TKS, at least one enzyme shows a predominant gene expression level in latex and roots, but most genes in the MEP pathway have a medium or low expression level in the TKS latex [16]. These results indicated that it is the MVA pathway rather than MEP pathway that is involved in the NBR in the rubber latex of both the rubber tree H. brasiliensis and the rubber grass TKS.

However, our comprehensive proteomics data showed that many proteins involved in both MVA (22 unique proteins) and MEP (unique 13 proteins) pathways have been positively identified from the mature TKS roots (Figure 7; Table S6). Among the identified 13 MEP-pathway proteins, two unique proteins named DXS 3 (spot 257) and DXS 8 (spots 199 and 204) are visible on the 2-DE gel (Figure 1), and their relative abundance is very high (Table S1). These proteomic data suggest that, by contrast with that in the inner bark of the Para rubber tree [10], both the MEP pathway and the MVA pathway might be the main source of IPP, and MVA pathway is also crucial for rubber biosynthesis in the mature TKS roots.

3.3. Almost All NRB-Related Proteins Can Be Identified From the Mature TKS Roots

In the Para rubber tree, NRB is known to begin with IPP synthesis in MVA pathway [38,39], which is a cytosolic pathway for rubber biosynthesis [40]. In the early steps of the MVA pathway, ACAT is important for generating acetoacytyl-CoA, then HMGS and HMGR activate the supply of mevalonate substrates [2,41]. ACAT can catalyze a Claisen-type condensation of two acetyl-CoA units to form acetoacetyl-CoA, which is the first step in the MVA pathway and is also important for providing the malonyl-CoA substrate for the biosynthesis of fatty acids [42]. This enzyme was found to mainly accumulate in the mature TKS roots (Table S1; Table S3), but it could not be detected in total latex of the Hevea rubber tree [2]. Through a series of enzymatic reactions, acetyl-CoA can be catalyzed by ACAT to form IPP [43]. Four ACAT genes were determined in the Hevea rubber tree [10] and eight ACAT genes were found in the TKS genome [16]. Overexpression of AtACAT results in the increased accumulation of ACAT in the latex of T. brevicorniculatum, and ultimately increased the pentacyclic triterpene and sterol levels [44]. Our proteomic results demonstrated that seven (ACAT 2-8) out of the eight ACAT members can be identified by the shotgun proteomic method, and one member (ACAT 7, spots 157, 183 and 424) can also be detected on the 2-DE gel (Figure 2), indicating that ACAT 7 is high-abundant in TKS roots.

In the following steps, HMGS and HMGR activate the supply of mevalonate substrates [2,41]. HMGS, the second enzyme in the MVA pathway, catalyzes aldol-type condensation of acetoacetyl-CoA and acetyl-CoA to produce HMG-CoA, and it is located as a cytosolic protein [45,46]. In our proteomic data, one HMGS member (HMGS 1, evm.model.utg 270.25; spot 146) was identified from both high-through and 2-DE gel-based proteomics methods (Figure 7; Table S6). Furthermore, HMG-CoA can be converted to MVA by HMGR [47]. HMGR controls the carbon flow and metabolic reaction rate in MVA [42], and it is also a critical player for the regulation of triterpenoid metabolism [38,48]. It was reported that ethylene cannot activate the enzyme activity of HMGR [49]. Five gene members have been characterized in the genome of the Para rubber tree [10], but only HMGR1 gene was found to be induced by ethylene stimulation [49]. In the TKS genome, 12 HMGR members were determined, and only two out of them (TkHMGR1 and TkHMGR2) were predominantly expressed in roots, with the highest expression level in latex [16]. In this proteomic study, two members of HMGR (HMGR1 and HMGR12) were detected from TKS roots by the shotgun method (Table S3). Based on these genomic and proteomic results, we consider that HMGR1 might be the most important player for NRB in TKS roots.

In the following steps, mevalonate kinase (MVK) converts mevalonate into isopentyl pyrophosphate and phosphomevalonate kinase (PMVK) catalyzes mevalonate-5-phosphoate to form mevalonate-5-pyrophosphate [50]. In Hevea, three MVK genes were determined [10]; and in TKS genome, 11 MVK genes were observed [16]. In our proteomic data, we identified eight MVK protein members from the predicted 11 MVK genes, and five PMVK members from the six PMVK genes (Table S6). Finally, MVD catalyzes mevalonate-5-pyrophosphate to form IPP [49]. Two MVD members were determined from both Hevea and TKS genomes [16], and in this proteomic study, they were identified from the mature TKS roots (Table S6). These results revealed that enzymes involved in the conversion of mevalonate into isopentyl pyrophosphate and diphosphate to form IPP are important for NRB in the mature TKS roots.

The MEP pathway is mainly occurred in the plastid and it is an alternative IPP synthesis pathway [51,52]. Recently, the plastidic MEP pathway has also been considered as a possible route for NRB [53,54,55]. The ¹³C-labelling of Hevea seedlings revealed that MEP pathway mainly contributes IPP for carotenoid biosynthesis [38]. Genomic analysis demonstrated that all genes involved in the MEP pathway can be detected in both Hevea and TKS, but most of the MEP-related genes have a medium or low expression level in the latex of the Hevea rubber tree [10] and TKS [16]. However, our proteomics results demonstrated that 13 unique proteins involved in the MEP pathway were determined in the TKS roots; these proteins are DXS, DXR, CMS, CMK, HDS and HDR (Table S6). These results indicated that the MEP pathway may also contribute IPP for NRB in the mature TKS roots.

In the initiator synthesis process, GPS and IPI catalyze IPP and DMAPP to form GPP, and FPS to produce farnesyl diphosphate by adding IPP onto GPP to form polyisoprene [40,56], then GGPS and GPS catalyze farnesyl diphosphate to form geranylgeranyl diphosphate [31,42]. Our proteomics data showed that four crucial enzymes, named IPI, GGPS, GPS and FPS, can be identified from the TKS roots.

Then, in the natural rubber elongation process, cis-isoprene transferase (CPT), also named Hevea rubber transferase, catalyzes multiple isoprene units (C5H8) and polymerizes into IPP long chain molecules and determine the size of the rubber molecule [46,57]. It is closely associated with the rubber particle membrane and is widely known to be a key enzyme in rubber biosynthesis [28,31]. Eleven CPT members were determined in the Para rubber tree genome [10], but only CPT2 shows rubber transferase activity [40], and a rubber transferase activator (RTA) or CPT-like protein (CPTL) can interact with CPTs to activate NRB in the Hevea rubber latex [57]. In the TKS genome, eight CPT members and two CPT-like proteins (CPTL1 and CPTL2) have been characterized, and TkCPT1, TkCPT2 and TkCPTL1 genes play a critical role for the elongation of rubber polymers in the latex and roots of TKS [16]. However, our proteomics data only identified CPT1 in the TKS roots (Table S6), indicated CPT1 may be the crucial player for NRB in TKS roots.

In the final elongation process, three rubber particle membrane binding proteins, named REF, SRPP and a Nogo-B receptor or HRT1-REF bridging protein (HRBP), are widely known to play crucial roles [35], and they were positively identified from the TKS roots by our proteomics study (Figure 7). Among them, REF is anchored inside rubber particle membrane by its auto-assembly ability [42], whereas SRPP largely covers rubber particle surface in an oriented anisotropic manner [58]. REF is also suggested to be a conservative regulator of lipid droplets and it is related to form an intracellular structure [59,60]. However, proteomics of the detergent-washed rubber particles identified none lipid droplet proteins [28,31], indicated that the maturation of rubber particles might be different from that in lipid droplets. Histochemical localization indicated that rubber biosynthesis capability in Hevea laticifer is mostly concentrated in SRPs, and SRPP, a rubber biosynthesis-related protein that accumulated mainly in SRPs, is more important than REF for rubber biosynthesis in the Hevea rubber tree [48]. Five SRPP genes (TbSRPP 1-5) were found in TBR [7]. Among them, except for TbSRPP 2 [61], four members can enhance NRB. But the average molecular weight of rubber hydrocarbons has not been significantly affected by these genes [62]. Overexpression of TbSRPP3 only resulted in a slight change in the content of rubber in TBR root [63]. In this proteomics of TKS roots, eight out of 10 SRPP members were identified (Table S6). Different with that in Hevea genome, which contains eight REF members, only two REF members have been characterized in the TKS genome [16]. Our proteomics data also identified the two REF members in the mature TKS roots (Table S6). These published results, as well as our proteomics of the 6M TKS roots, revealed that SRPP might be more important than REF for rubber biosynthesis in the mature roots of TKS, and more attention should be paid on determination of the detail roles of SRPP in TKS roots in future.

4. Materials and Methods 4.1. Plant Growth Conditions and Root Collection

TKS seeds were collected from the Tekes River basin in Xinjiang, China. Following germination, seedlings were transplanted into the pots in a greenhouse of Shihezi University (Shihezi, China), containing nutritive soil and vermiculite with a ratio of 5/3, and a half strength of Hoagland liquid medium was irrigated with 55% relative humidity at 22 °C in the light and 18 °C in the dark. After growing for 6 months (termed as 6M roots), the main roots from plants were collected, rinsed with tapped water, blotted dry on a filter paper, and the middle parts were dissected into approximately 0.5 cm-thick slices, and then frozen in liquid nitrogen for further study.

4.2. Morphological Observation of Rubber Particles in TKS Roots

For histochemical staining, a 0.5 cm long section from the middle parts of the fresh main roots of TKS were fixed in 4% glutaraldehyde (0.1 M phosphate buffer, pH 7.5) at room temperature, then treated with bromine and iodine in glacial acetic acid as described [64]. Part of the sample was dehydrated through a graded series of ethanol, and embedded in paraffin. Sections (12 μm thickness and area of 5 mm × 5 mm) were cut with microtome (Leica Microsystems, Bannockburn IL, Germany) stained with mercury–bromophenol blue, which is effective in showing the latex protein [65]. Then, the slides were examined under a leica DMLN electron microscope (leica, Wetzlar, Germany). TKS laticifer cells in the root sections could be recognized due to iodine-bromine treatment of the rubber in the laticifers became deep brown or dark color. Another part of the sample was dehydrated in a range of acetone (30%, 50%, 70%, 90%, 1 h for once and pure acetone 2 h for twice). Then it was permeated with acetone and Epon 812 embedding agent and pure Epon 812 embedding agent permeates overnight as described [33]. After that, the samples were polymerized in the oven, flattened and cut into semi-thin slices. Observed under an optical microscope, the laticifers could be recognized by tracing rubber inclusions, which is brown with iodine-bromine staining [65]. For electron microscopy, samples were cut into a smaller size, immediately fixed in the glutaraldehyde solution at 4 °C for 24 h, and then fixed in phosphate buffer (pH 7.2, 2% OsO4 in 0.1 M) for 6 h at room temperature. Ultra-thin sections, double stained with uranium acetate and lead citrate again. A transmission electron microscope (TEM, Hitachi JEM-1230, JEOL, Japan)), country was used to examine these stained slices.

4.3. Protein Extraction and Two-Dimensional Gel Electrophoresis

Total proteins were extracted from TKS roots using a modified BPP phenol method as described [22]. The washed protein pellets were air-dried for 10 min and dissolved with the Lysis buffer (30 mM Tris-HCl, pH 8.5, 2 M thiourea, 7 M 17 urea, 4% 3-[(3-Cholamidopropyl) dimethylammonio] propanesulfonate (CHAPS)). Protein concentration was determined by the Bradford assaying, and BSA was used as the protein standard. For 2-DE assay, 1300 μg of proteins were diluted to 450 μL with lysis buffer and loaded onto a 24 cm, linear pH gradient 4–7 IPG strips (GE Healthcare, Uppsala, Sweden). Then, they were rehydrated for 24 h at 22 °C. Isoelectric focusing system and gel electrophoresis were performed as described [26]. The gels were stained with Coomassie Brilliant Blue G-250, and analyzed by the ImageMaster 2-D Platinum Software as described [2]. Biological variation analysis module was employed to identify the accumulated protein spots in roots samples with statistically significant (confidence above 95%, p < 0.05). Three biological replicates for each separation were conducted and the three typical 2-DE gels are presented in the supplementary Figure S1.

4.4. Identification of Protein Spots via MALDI-TOF/TOF-MS

The main protein spots on the 2-DE gel of 6M roots were manually excised from 2-DE gels and subjected to perform an in-gel digestion with the modified bovine trypsin (Cat. 11418025001, Roche, Basel, Switzerland) as described [64]. After digestion, mass spectra of trypsin-digested peptides map fingerprinting were obtained on an AB 5800 MALDI-TOF/TOF mass spectrometry (AB SCIEX, Foster City, MA, USA) with a laser wavelength of 349 nm as described [2]. Ten main spectra were selected for MS/MS analysis and searched for a self-constructed database derived from the original TKS genome. The nucleotide and protein sequences were deposited in the Genome Warehouse (GWH; http://bigd.big.ac.cn/gwh/) under the accession number PRJCA000437 (https://academic.oup.com/nsr), which includes 46,731 gene sequences [16]. Protein matches were considered as a positive identification with a Mascot score higher than 75 and at least two matched peptides, meanwhile individual ion score > 31 (p < 0.05). Detailed information for the identified proteins is provided as the supplementary data (Table S1 and Figure S2).

4.5. Large-Scale Shotgun Proteomics of TKS Roots

Approximately 300 μg proteins from the main roots of TKS were reduced and alkylated by dithiothreitol and iodoacetamide, then digested by trypsin. Enzymic hydrolysates were separated using a Pepmap C18 column (Thermo Fisher, USA) following a 65 min 5%–35% organic gradient by an UltiMate 3000 instrument (Thermo Scientific, Rockford, USA). Fifty washed fractions were collected in 1.5 mL/min tubes. Fractions were dehydrated and merged into 15 fractions, and these samples were subjected to a high-resolution mass spectrometer Triple TOF 6600 system (AB SCIEX). Peptides were automatically selected by ProGroupTM algorithm and ProteinPilot™ software V5.0 (AB SCIEX) to calculate the error factor (EF), reporter peak area and p value as described [30]. A protein database was established by the target protein sequences (https://academic.oup.com/nsr) for the corresponding species. Protein with an unused score > 1.3 (confidence ≥ 95%) was considered as a positive identification. Finally, an in-house BlastP search of the UniProt database was performed for each protein to identify its homologues and potential functions.

4.6. Functional Classification, Hierarchical Clustering and Pathway Analysis of All the Identified Proteins From TKS Roots

The sequences of all the identified proteins were search against the UniProt database (http://www.uniprot.org/) to confirm their functional annotations. The identified proteins were then classified by Gene Ontology (GO, http://www.geneontology.org/) as described [66]. The GO terms were further determined by the WEGO software (http://wego.genomics.org.cn) on biological process, molecular functions and cellular components. Finally, these proteins were performed KEGG pathway analysis (http://www.genome.jp/kegg) to determine their molecular reaction networks.

5. Conclusions

In this study, we obtained the first visualization proteome profiles of mature TKS roots. Combination of 2-DE and shotgun proteomics leads to more than 3,000 unique proteins in TKS roots. Functional category of all TKS proteins revealed that most proteins were involved in carbon metabolic process with catalytic activity in membrane-bounded organelle, followed by proteins with binding ability, transportation activity, starch and sucrose metabolism, and phenylpropanoid biosynthesis. They are important for biosynthesis of isoprenoids in natural rubber-producing plants. Fifty-eight NRB-related proteins were identified, and they are involved in both MVA and MEP pathways. Our results indicated the MEP pathway is also important by contributing IPP for rubber biosynthesis in TKS roots. Almost all NRB-related proteins can be identified, and several members of many key proteins, including ACAT, HMGR, MVK, PMVK, DXS, DXR, GPS, FPS, IPI, REF and SRPP, were determined from the mature TKS roots. Future work should pay more attention to determination of the detailed roles of different members of these NRB-related proteins in the mature TKS roots.

Figure 1. Morphological observation of the laticifer cells and rubber particles in the main roots of Taraxacum kok-saghyz (TKS). A six months old TKS plant and its main roots are highlighted (A). The laticifer cells containing brown rubber can be examined under the crosscutting (B) and slitting (C) slices of the 6M TKS roots. Black arrows indicate the positions of typical laticifer cells. Both small rubber particles (SRP) and large rubber particles (LRP) in the main roots of 6M TKS can be observed under a transmission electron microscope (TEM) (D).

Figure 2. The reference proteome map and MS identification of abundant proteins in 2-DE gels of the 6M TKS roots. The proteins extracted from the 6M mature TKS roots were examined with two-dimensional gel electrophoresis (2-DE) to produce a high-resolution reference proteome map. The 371 major protein spots marked with numbers were positively identified by MALDI TOF/TOF MS and their positions are indicated by white arrows. The detailed identities of these proteins are listed in Table S1 and Figure S2.

Figure 3. Distribution of the identified proteins on the reference 2-DE gel of 6M TKS roots. The theoretical and experimental ratios of Mr and pI for the identified 371 are presented in radial charts with radial and annular radar axes (A). The theoretical distribution of Mr and pI for the protein isoforms is highlighted (B). These protein isoforms were independently identified from two (black circle), three (blue circle), four (green circle), and five (red circle) spots on the reference 2-DE gels of the 6M TKS roots, respectively. The most abundant 20 proteins in the 2-DE gel of 6M roots are highlighted (C). The detailed identities for these protein isoforms are listed in Table S2.

Figure 4. High-throughput proteomics analysis of the 6M TKS roots. Three shotgun proteomics were performed to determine the large-scale protein profiles in the 6M roots, and finally identify 3545 shared proteins in the three experiments (A) The blue, yellow, and pink circles represent for the first, second, and third biological replicates. Among these proteins, 184 proteins could also be detected from the 2-DE gel (B). The pink area represents for the 184 shared proteins that were identified by both the 2-DE and shotgun based proteomic methods.

Figure 5. GO classification of the identified proteins. The 3592 proteins identified from the 6M roots by both shotgun and 2-DE methods were performed GO analysis and these proteins were classified into 3 main categories including biological process, cellular component and molecular function. The number of proteins denotes those with GO annotations.

Figure 6. KEGG pathway analysis of all the identified proteins. KEGG pathways for the identified 3592 proteins are presented (A). The number of proteins enriched in the KEGG database is marked in the right side. Then, the most abundant 20 KEGG pathways are highlighted (B). The size of the circular ring, ranging from 83 to 240, stands for the number of proteins involved in these pathways.

Figure 7. Schematic diagram of natural rubber biosynthesis (NRB)-related proteins in the 6M TKS roots. A schematic diagram of the identified proteins involved in the mevalonate acid (MVA, left) and methylerythritol phosphate (MEP, right) pathways of NRB is highlighted. In this schematic diagram, the identified proteins involved in NRB from the 6M TKS roots were marked with different color. The name of proteins identified from both 2-DE and shotgun methods was marked in red color. The proteins only determined by shotgun analysis are marked in blue color. The proteins identified from the first (6M-1), second (6M-2) and third (6M-3) shotgun experiments are marked with red, blue and black stars, respectively.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/20/10/2596/s1.

Author Contributions

X.W. and H.L. conceived and designed the experiments; Q.X., L.Z., G.D., D.W. and Y.S. performed the experiments; Q.X. wrote the original manuscript. X.W. revised this manuscript. Q.X., L.Z., B.Y., D.W., Y.S., X.G. and L.Y. analyzed the data. All authors have read and approved the final manuscript.

Funding

This research was funded by the Scientific Research Innovation Team Program of Hainan Province, grant number 2018CXTD341, the National Natural Science Foundation of China, grant numbers 31570301 and 31860224, and International scientific and technological cooperation to advance the project by the Shihezi University, grand number GJHZ201708.

Conflicts of Interest

The authors declare that they have no conflict of interests.

Abbreviations

ACL ATP citrate lyase

ACAT acetoacetyl-CoA thiolase

CHAPS 3-[(3-Cholamidopropyl) dimethylammonio] propanesulfonate

CPT cis-isoprene transferase

CPTL cis-isoprene transferase like

DMAPP dimethylallyl diphosphate

DXS 1-deoxy-D-xylulose-5-phosphate synthase

eIF eukaryotic translation initiation factor

FPS farnesyl diphosphate synthase

GAPC glyceraldehyde-3-phosphate dehydrogenase

GGPPS geranylgeranyl diphosphate synthase

GPP geranyl pyrophosphate

GPS geranylgeranyl diphosphate synthase

GO Gene Ontology

KEGG Kyoto Encyclopedia of Genes and Genomes

HMGS 3-hydroxy-glutaryl-CoA reductase

HMGR hydroxy-methyl-glutaryl-CoA reductase

IPI isopentenyl-diphosphate delta-isomerase

IPG immobilized pH gradient

iTRAQ isobaric tags for relative and absolute quantification

LRP large rubber particle protein

MEP methyl-D-erythritol-4-phosphate

MVA mevalonate acid

MVD diphosphomevalonate decarboxylase

MVK mevalonate kinase

MS mass spectrometry

NR natural rubber

NRB natural rubber biosynthesis

PMVK mevalonate disphosphate decarboxylase

REF rubber elongation factor

RTA rubber transferase activating factor

SRPP small rubber particle protein

2-DE two-dimensional gel electrophoresis

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

References

1. Mooibroek, H.; Cornish, K. Alternative sources of natural rubber. Appl. Microbiol. Biotechnol. 2000, 53, 355-365.

2. Wang, X.C.; Wang, D.; Sun, Y.; Yang, Q.; Chang, L.L.; Wang, L.M.; Meng, X.R.; Huang, Q.X.; Jin, X.; Tong, Z. Comprehensive proteomics analysis of laticifer latex reveals new insights into ethylene stimulation of natural rubber production. Sci. Rep. 2015, 5, 13778.

3. Schmidt, T.; Hillebrand, A.; Wurbs, D.; Wahler, D.; Lenders, M.; Gronover, C.; Prüfer, D. Molecular cloning and characterization of rubber biosynthetic genes from Taraxacum kok-saghyz. Plant Mol. Biol. Rep. 2010, 28, 277-284.

4. Cornish, K.; Kopicky, S.L.; Mcnulty, S.K.; Amstutz, N.; Chanon, A.M.; Walker, S.; Kleinhenz, M.D.; Miller, A.R.; Streeter, J.G. Temporal diversity of Taraxacum kok-saghyz plants reveals high rubber yield phenotypes. Biodiversitas 2016, 17, 847-856.

5. Edathil, T.T. South American leaf blight: a potential threat to the natural rubber industry in Asia and Africa. Trop. Pest Manag. 1986, 4, 296-303.

6. Ramirez-Cadavid, D.A.; Cornish, K.; Michel, J. Taraxacum kok-saghyz (TK): compositional analysis of a feedstock for natural rubber and other bioproducts. Ind. Crop. Prod. 2017, 107, 624-640.

7. Collins-Silva, J.; Nural, A.T.; Skaggs, A.; Scott, D.; Hathwaik, U.; Woolsey, R.; Schegg, K.; McMahan, C.; Whalen, M.; Cornish, K.; et al. Altered levels of the Taraxacum kok-saghyz (Russian dandelion) small rubber particle protein, TKSRPP3, result in qualitative and quantitative changes in rubber metabolism. Phytochemistry 2012, 79, 46-56.

8. Kim, I.J.; Ryu, S.B.; Kwak, Y.S.; Kang, H. A novel cDNA from Parthenium argentatum Gray enhances the rubber biosynthetic activity in vitro. J. Exp. Bot. 2004, 396, 377-385.

9. Ikeda, Y.; Junkong, P.; Ohashi, T.; Phakkeeree, T.; Sakaki, Y.; Tohsan, A.; Kohjiya, S.; Cornish, K. Strain-induced crystallization behaviour of natural rubbers from guayule and rubber dandelion revealed by simultaneous time-resolved WAXD/tensile measurements: indispensable function for sustainable resources. RSC Adv. 2016, 6, 95601-95610.

10. Tang, C.R.; Yang, M.; Fang, Y.J.; Luo, Y.F.; Gao, S.H.; Xiao, X.H.; An, Z.W.; Zhou, B.H.; Zhang, B.; Tan, X.Y.; et al. The rubber tree genome reveals new insights into rubber production and species adaptation. Nat. Plants 2016, 6, 16073.

11. Whaley, W.G.; Bowen, J.S. Russian dandelion (kok-saghyz). An emergency source of natural rubber. Misc. Publ. US Dept. Agric. 1947, 618, 1-212.

12. van Beilen, J.B.; Poirier, Y. Establishment of new crops for the production of natural rubber. Trends Biotechnol. 2007, 25, 522-529.

13. Kirschner, J.; Stepanek, J.; Cerny, T.; Heer, P.; Dijk, P.J. Available ex situ germplasm of the potential rubber crop Taraxacum kok-saghyz belongs to a poor rubber producer, T. brevicorniculatum (Compositae-Crepidinae). Genet. Resour. Crop Evol. 2013, 2, 455-471.

14. Kreuzberger, M.; Hahn, T.; Zibek, S.; Schiemann, J.; Thiele, K. Seasonal pattern of biomass and rubber and inulin of wild Russian dandelion (Taraxacum kok-saghyz L. Rodin) under experimental field conditions. Eur. J. Agron. 2016, 80, 66-77.

15. Stolze, A.; Wanke, A.; van Deenen, N.; Geyer, R.; Prüfer, D.; Gronover, C.S. Development of rubber-enriched dandelion varieties by metabolic engineering of the inulin pathway. Plant Biotechnol. J. 2017, 6, 740-753.

16. Lin, T.; Xu, X.; Ruan, J.; Liu, S.Z.; Wu, S.G.; Shao, X.J.; Wang, X.B.; Gan, L.; Qin, B.; Yang, Y.S.; et al. Genome analysis of Taraxacum kok-saghyz Rodin provides new insights into rubber biosynthesis. Natl. Sci. Rev. 2017, 5, 78-87.

17. Ramirez-Cadavid, D.A.; Valles-Ramirez, S.; Cornish, K.; Michel, F.C. Simultaneous quantification of rubber, inulin, and resins in Taraxacum kok-saghyz (TK) roots by sequential solvent extraction. Plant Biotechnol. J. 2018, 122, 647-656.

18. Bach, T.J. Some new aspects of isoprenoid biosynthesis in plants a review. Lipids 1995, 3, 191-202.

19. Walther, T.C.; Farese, R.V. The life of lipid droplets. Biochim. Biophys. Acta 2009, 6, 459-466.

20. Seetang-Nun, Y.; Sharkey, T.D.; Suvachittanont, W. Molecular cloning and characterization of two cDNAs encoding 1-deoxy-D-xylulose 5-phosphate reductoisomerase from Hevea brasiliensis. J. Plant Physiol. 2007, 9, 991-1002.

21. Yagami, T.; Haishima, Y.; Tsuchiya, T.; Tomitaka-Yagami, A.; Kano, H.; Matsunaga, K. Proteomic analysis of putative latex allergens. Int. Arch. Allergy Immunol. 2004, 135, 3-11.

22. Wang, X.C.; Shi, M.J.; Lu, X.L.; Ma, R.F.; Wu, C.G.; Guo, A.P.; Peng, M.; Tian, W.M. A method for protein extraction from different subcellular fractions of laticifer latex in Hevea brasiliensis compatible with 2-DE and MS. Proteome Sci. 2010, 8, 35.

23. D'Amato, A.; Bachi, A.; Fasoli, E.; Boschetti, E.; Peltre, G.; Sénéchal, H.; Sutrad, J.P.; Citteriob, A.; Righettib, P.G. In-depth exploration of Hevea brasiliensis latex proteome and hidden allergens via combinatorial peptide ligand libraries. J. Proteom. 2010, 73, 1368-1380.

24. Habib, M.H.; Yuen, G.C.; Othman, F.; Zainudin, N.N.; Latiff, A.A.; Ismail, M.N. Proteomics analysis of latex from Hevea brasiliensis (clone RRIM 600). Biochem. Cell Biol. 2017, 95, 232-242.

25. Tong, Z.; Wang, D.; Sun, Y.; Yang, Q.; Meng, X.R.; Wang, L.M.; Feng, W.Q.; Li, L.; Wurtele, E.S.; Wang, X.C. Comparative proteomics of rubber latex revealed multiple protein species of REF/SRPP family respond diversely to ethylene stimulation among different rubber tree clones. Int. J. Mol. Sci. 2017, 18, 958.

26. Rojruthai, P.; Sakdapipanich, J.T.; Takahashi, S.; Hyegin, L.; Noike, M.; Koyama, T.; Tanaka, Y. In vitro synthesis of high molecular weight rubber by Hevea small rubber particles. J. Biosci. Bioeng. 2010, 109, 107-114.

27. Xiang, Q.L.; Xia, K.C.; Dai, L.J.; Kang, G.J.; Li, Y.; Nie, Z.Y.; Duan, C.F.; Zeng, R.Z. Proteome analysis of the large and the small rubber particles of Hevea brasiliensis using 2D-DIGE. Plant Phys. Biochem. 2012, 60, 207-213.

28. Dai, L.J.; Kang, G.J.; Li, Y.; Nie, Z.Y.; Duan, C.F.; Zeng, R.Z. In-depth proteome analysis of the rubber particle of Hevea brasiliensis (para rubber tree). Plant Mol. Biol. 2013, 82, 155-168.

29. Dai, L.J.; Kang, G.J.; Nie, Z.Y.; Li, Y.; Zeng, R.Z. Comparative proteomic analysis of latex from Hevea brasiliensis treated with ethrel and methyl jasmonate using iTRAQ-coupled two-dimensional LC-MS/MS. J. Proteom. 2016, 132, 167-175.

30. Wang, D.; Sun, Y.; Tong, Z.; Yang, Q.; Chang, L.L.; Meng, X.R.; Wang, L.M.; Tian, W.M.; Wang, X.C. A protein extraction method for low protein concentration solutions compatible with the proteomic analysis of rubber particles. Electrophoresis 2016, 37, 2930-2939.

31. Yamashita, S.; Yamaguchi, H.; Waki, T.; Aoki, Y.; Mizuno, M.; Yanbe, F.; Ishii, T.; Funaki, A.; Tozawa, Y.; Miyagi-Inoue, Y.; et al. Identification and reconstitution of the rubber biosynthetic machinery on rubber particles from Hevea brasiliensis. eLife 2016, 5, e19022.

32. Wahler, D.; Colby, T.; Kowalski, N.A.; Harzen, A.; Wotzka, S.Y.; Hillebrand, A.; Fischer, R.; Helsper, J.; Schmidt, J.; Gronover, C.S.; et al. Proteomic analysis of latex from the rubber-producing plant Taraxacum brevicorniculatum. Proteomics 2012, 12, 901-905.

33. Singh, A.P.; Wi, S.G.; Chung, G.C.; Kim, Y.S.; Kang, H. The micromorphology and protein characterization of rubber particles in Ficus carica, Ficus benghalensis and Hevea brasiliensis. J. Exp. Bot. 2003, 54, 985-992.

34. Jorrin-Novo, V.J.; Pascual, J.; Sanchez-Lucas, R.; Romero-Rodriguez, M.C.; Rodriguez-Ortega, M.J.; Lenz, C.; Valledor, L. Fourteen years of plant proteomics reflected in proteomics: moving from model species and 2DE-based approaches to orphan species and gel-free platforms. Proteomics 2015, 15, 1089-1112.

35. Wang, D.; Sun, Y.; Chang, L.L.; Tong, Z.; Xie, Q.L.; Jin, X.; Zhu, L.P.; He, P.; Li, H.B.; Wang, X.C. Subcellular proteome profiles of different latex fractions revealed washed solutions from rubber particles contain crucial enzymes for natural rubber biosynthesis. J. Proteom. 2018, 182, 53-64.

36. Cho, W.K.; Jo, Y.; Chu, H.; Park, S.H.; Kim, K.H. Integration of latex protein sequence data provides comprehensive functional overview of latex proteins. Mol. Biol. Rep. 2014, 41, 1469-1481.

37. Brown, D.; Feeney, M.; Ahmadi, M.; Lonoce, C.; Sajari, R.; Di Cola, A.; Frigerio, L. Subcellular localization and interactions among rubber particle proteins from Hevea brasiliensis. J. Exp. Bot. 2017, 68, 5045-5055.

38. Sando, T.; Takeno, S.; Watanabe, N.; Okumoto, H.; Kuzuyama, T.; Yamashita, A.; Hattori, M.; Ogasawara, N.; Fukusaki, E.; Kobayashi, A. Cloning and characterization of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway genes of a natural-rubber producing plant, Hevea brasiliensis. Biosci. Biotechnol. Biochem. 2008, 11, 2903-2917.

39. Wu, C.T.; Li, Y.; Nie, Z.Y.; Dai, L.J.; Kang, G.J.; Zeng, R.Z. Molecular cloning and expression analysis of the mevalonate diphosphate decarboxylase gene from the latex of Hevea brasiliensis. Tree Genet. Genomes 2017, 13, 22.

40. Lau, N.S.; Makita, Y.; Kawashima, M.; Taylor, T.D.; Kondo, S.; Othman, A.S.; Shu-Chien, A.C.; Matsui, M. The rubber tree genome shows expansion of gene family associated with rubber biosynthesis. Sci. Rep. 2016, 6, 28594.

41. Puskas, J.E.; Gautriaud, E.; Deffieux, A.; Kennedy, J.P. Natural rubber biosynthesis-a living carbocationic polymerization. Prog. Polym. Sci. 2006, 31, 533-548.

42. Cornish, K. Similarities and differences in rubber biochemistry among plant species. Phytochemistry 2001, 57, 1123-1134.

43. Kirby, J.; Keasling, J.D. Biosynthesis of plant isoprenoids perspectives for microbial engineering. Annu. Rev. Plant Biol. 2009, 60, 335-355.

44. Putter, K.M.; van Deenen, N.; Unland, K.; Prufer, D.; Gronover, C.S. Isoprenoid biosynthesis in dandelion latex is enhanced by the overexpression of three key enzymes involved in the mevalonate pathway. BMC Plant Biol. 2017, 17, 88.

45. Sirinupong, N.; Suwanmanee, P.; Doolittle, R.F.; Suvachitanont, W. Molecular cloning of a new cDNA and expression of 3-hydroxy-3-methylglutaryl-CoA synthase gene from Hevea brasiliensis. Planta 2005, 221, 502-512.

46. Grabinska, K.A.; Park, E.J.; Sessa, W.C. cis-Prenyltransferase: new insights into protein glycosylation, rubber synthesis, and human diseases. J. Biol. Chem. 2016, 291, 18582-18590.

47. Schaller, H.; Crausem, B.; Benveniste, P.; Chye, M.L.; Tan, C.T.; Song, Y.H.; Chua, N.H. Expression of the Hevea brasiliensis Müll. Arg. 3-Hydroxy-3-Methylglutaryl-Coenzyme A Reductase 1 in tobacco results in sterol overproduction. Plant Physiol. 1995, 109, 761-770.

48. Sando, T.; Hayashi, T.; Takeda, T.; Akiyama, Y.; Nakazawa, Y.; Fukusaki, E.; Kobayashi, A. Histochemical study of detailed laticifer structure and rubber biosynthesis-related protein localization in Hevea brasiliensis using spectral confocal laser scanning microscopy. Planta 2009, 230, 215-225.

49. Chye, M.L.; Tan, C.T.; Chua, N.H. Three genes encode 3-hydroxy-3-methylglutaryl-coenzyme A reductase in Hevea brasiliensis: hmg1 and hmg3 are differentially expressed. Plant Mol. Biol. 1992, 19, 473-484.

50. Williamson, B.Y.; Kekwick, R.G. The enzymes forming isopentenyl pyrophosphate from 5-phosphomevalonate (mevalonate 5-phosphate) in the latex of Hevea brasiliensis. Biochem. J. 1971, 124, 407-417.

51. Ko, J.H.; Chow, K.S.; Han, K.H. Transcriptome analysis reveals novel features of the molecular events occurring in the laticifers of Hevea brasiliensis (para rubber tree). Plant Mol. Biol. 2003, 53, 479-492.

52. Swiezewska, E.; Danikiewicz, W. Polyisoprenoids: Structure, biosynthesis and function. Prog. Lipid Res. 2005, 44, 235-258.

53. Chow, K.S.; Wan, K.L.; Isa, M.N.; Bahari, A.; Tan, S.H.; Harikrishna, K.; Yeang, H.Y. Insights into rubber biosynthesis from transcriptome analysis of Hevea brasiliensis latex. J. Exp. Bot. 2007, 58, 2429-2440.

54. Chow, K.S.; Mat-Isa, M.N.; Bahari, A.; Ghazali, A.K.; Alias, H.; Mohd-Zainuddin, Z.; Hoh, C.C.; Wan, K.L. Metabolic routes affecting rubber biosynthesis in Hevea brasiliensis latex. J. Exp. Bot. 2012, 63, 1863-1871.

55. Makita, Y.; Ng, K.K.; Veera Singham, G.; Kawashima, M.; Hirakawa, H.; Sato, S.; Othman, A.S.; Matsui, M. Large-scale collection of full-length cDNA and transcriptome analysis in Hevea brasiliensis. DNA Res. 2017, 24, 159-167.

56. Uthup, T.K.; Saha, T.; Ravindran, M.; Bini, K. Impact of an intragenic retrotransposon on the structural integrity and evolution of a major isoprenoid biosynthesis pathway gene in Hevea brasiliensis. Plant Physiol. Biochem. 2013, 73, 176-188.

57. Epping, J.; Deenen, N.V.; Niephaus, E.; Stolze, A.; Fricke, J.; Huber, C.; Eisenreich, W.; Twyman, R.M.; Prufer, D.; Gronover, C.S. A rubber transferase activator is necessary for natural rubber biosynthesis in dandelion. Nat. Plants 2015, 1, 15048.

58. Berthelot, K.; Lecomte, S.; Estevez, Y.; Peruch, F. Hevea brasiliensis REF (Hev b 1) and SRPP (Hev b 3): An overview on rubber particle proteins. Biochimie 2014, 106, 1-9.

59. Gidda, S.K.; Park, S.; Pyc, M.; Yurchenko, O.; Cai, Y.; Wu, P.; Andrews, D.W.; Chapman, K.D.; Dyer, J.M.; Mullen, R.T. Lipid drop let associated proteins (LDAPs) are required for the dynamic regulation of neutral lipid compartmentation in plant cells. Plant Physiol. 2016, 170, 2052-2071.

60. Kim, E.Y.; Park, K.Y.; Seo, Y.S.; Kim, W.T. Arabidopsis Small rubber particle protein homolog SRPs play dual roles as positive factors for tissue growth and development and in drought stress responses. Plant Physiol. 2016, 4, 2494-2510.

61. Post, J.; van Deenen, N.; Fricke, J.; Kowalski, N.; Wurbs, D.; Schaller, H.; Eisenreich, W.; Huber, C.; Twyman, R.M.; Prüfer, D.; et al. Laticifer-specific cis-prenyltransferase silencing affects the rubber, triterpene, and inulin content of Taraxacum brevicorniculatum. Plant Physiol. 2012, 158, 1406-1417.

62. Hillebrand, A.; Post, J.J.; Wurbs, D.; Wahler, D.; Lenders, M.; Krzyzanek, V.; Prufer, D.; Gronover, C.S. Down-regulation of small rubber particle protein expression affects integrity of rubber particles and rubber content in Taraxacum brevicorniculatum. PLoS ONE 2012, 7, e41874.

63. Tata, S.K.; Choi, J.Y.; Jung, J.Y.; Lim, K.Y.; Shin, J.S.; Ryu, S.B. Laticifer tissue-specific activation of the Hevea SRPP promoter in Taraxacum brevicorniculatum and its regulation by light, tapping and cold stress. Ind. Crop. Prod. 2012, 40, 219-224.

64. Wang, X.C.; Shi, M.; Wang, D.; Chen, Y.; Cai, F.; Zhang, S.; Wang, L.; Tong, Z.; Tian, W.M. Comparative proteomics of primary and secondary lutoids reveals that chitinase and glucanase play a crucial combined role in rubber particle aggregation in Hevea brasiliensis. J. Proteome Res. 2013, 12, 5146-5159.

65. Hao, B.; Wu, J. Laticifer differentiation in Hevea brasiliensis induction by exogenous jasmonic acid and linolenic acid. Ann. Bot. 2000, 85, 37-43.

66. Conesa, A.; Gotz, S.; García-Gomez, J.M.; Terol, J.; Talon, M.; Robles, M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674-3676.