<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2024.1444693</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Integrative analysis of transcriptome and metabolome reveals how ethylene increases natural rubber yield in <italic>Hevea brasiliensis</italic>
</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Hong</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/333209"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Dai</surname>
<given-names>Longjun</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Mingyang</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2830239"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fan</surname>
<given-names>Xiaokang</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2830219"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Lu</surname>
<given-names>Liangruinan</given-names>
</name>
<uri xlink:href="https://loop.frontiersin.org/people/2830250"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Bingbing</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Zhenhui</given-names>
</name>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Wang</surname>
<given-names>Lifeng</given-names>
</name>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1713163"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<institution>Rubber Research Institute, Key Laboratory of Biology and Genetic Resources of Rubber Tree, Ministry of Agriculture and Rural Affairs, Chinese Academy of Tropical Agricultural Sciences (CATAS)</institution>, <addr-line>Haikou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Damian Gruszka, University of Silesia in Katowice, Poland</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Quaid Hussain, Shenzhen University, China</p>
<p>Xiyang Zhao, Jilin Agricultural University, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Zhenhui Wang, <email xlink:href="mailto:wzh-36@163.com">wzh-36@163.com</email>; Lifeng Wang, <email xlink:href="mailto:lfwang@catas.cn">lfwang@catas.cn</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>09</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1444693</elocation-id>
<history>
<date date-type="received">
<day>06</day>
<month>06</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>08</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Yang, Dai, Liu, Fan, Lu, Guo, Wang and Wang</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Yang, Dai, Liu, Fan, Lu, Guo, Wang and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>
<italic>Hevea brasiliensis</italic> is an important cash crop with the product named natural rubber (NR) for markets. Ethylene (ET) is the most effective yield stimulant in NR production but the molecular mechanism remains incomplete. Here, latex properties analysis, transcriptome analysis, and metabolic profiling were performed to investigate the mechanism of NR yield increase in four consecutive tappings after ET stimulation. The results revealed that sucrose and inorganic phosphate content correlated positively with dry-rubber yield and were induced upon ET stimulation. Stimulation with ET also led to significant changes in gene expression and metabolite content. Genes involved in phytohormone biosynthesis and general signal transduction as well as 51 transcription factors potentially involved in the ET response were also identified. Additionally, KEGG annotation of differentially accumulated metabolites suggested that metabolites involved in secondary metabolites, amino-acid biosynthesis, ABC transporters, and galactose metabolism were accumulated in response to ET. Integrative analysis of the data collected by transcriptomics and metabolomics identified those differentially expressed genes and differentially accumulated metabolites are mainly involved in amino-acid biosynthesis and carbohydrate metabolism. Correlation analysis of genes and metabolites showed a strong correlation between amino-acid biosynthesis during ET stimulation. These findings provide new insights into the molecular mechanism underlying the ET-induced increase in rubber yield and further our understanding of the regulatory mechanism of ethylene signaling in rubber biosynthesis.</p>
</abstract>
<kwd-group>
<kwd>natural rubber</kwd>
<kwd>ethylene</kwd>
<kwd>latex yield</kwd>
<kwd>transcriptome</kwd>
<kwd>metabolome</kwd>
</kwd-group>
<counts>
<fig-count count="9"/>
<table-count count="1"/>
<equation-count count="0"/>
<ref-count count="53"/>
<page-count count="15"/>
<word-count count="4932"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Crop and Product Physiology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The rubber tree (<italic>Hevea brasiliensis</italic> Muell. Arg.) is cultivated in almost 30 tropical and subtropical countries to produce natural rubber (NR) (<xref ref-type="bibr" rid="B34">Qi et&#xa0;al., 2023</xref>). NR, the main constituent of which is cis-1,4-polyisoprene, is widely used owing to its unique characteristics, e.g. superior resilience, elasticity, efficient heat dispersion, abrasion resistance, and impact resistance (<xref ref-type="bibr" rid="B52">Zhang et&#xa0;al., 2023</xref>). Most importantly, it cannot be replaced by synthetic rubbers (<xref ref-type="bibr" rid="B9">Cornish, 2001</xref>). NR is obtained from the milky latex synthesized in highly specialized rubber-tree cells in the phloem, known as laticifers. The cytoplasmic content of these specialized laticifers is expelled in the form of a milky liquid when the bark is tapped, which is known as natural raw rubber latex (<xref ref-type="bibr" rid="B1">Adiwilaga and Kush, 1996</xref>; <xref ref-type="bibr" rid="B35">Qin et&#xa0;al., 2022</xref>). It is composed of water, rubber hydrocarbon, proteins, enzymes, nucleic acids, lipids, esters, and inorganic salts (<xref ref-type="bibr" rid="B40">Suksup et&#xa0;al., 2017</xref>). Rubber hydrocarbon, which accounts for 25% to 41% of the latex, is a polymer of isoprene units (<xref ref-type="bibr" rid="B24">Lehman et&#xa0;al., 2022</xref>). The other non-rubber substances also greatly affect the unique chemical and physical properties of NR (<xref ref-type="bibr" rid="B24">Lehman et&#xa0;al., 2022</xref>). Meanwhile, the capacity of rubber trees to regenerate latex components after tapping is a crucial factor that determines the efficiency of rubber yield (<xref ref-type="bibr" rid="B53">Zhu and Zhang, 2009</xref>). During the process of latex regeneration, the cells of the rubber tree undergo rapid metabolism to replenish the lost cytoplasmic contents. This process involves synthesizing and transporting large amounts of biomolecules, such as proteins, enzymes, nucleic acids, carbohydrates, lipids, and inorganic ions including potassium, inorganic phosphate, magnesium, etc.</p>
<p>The biosynthesis of NR, like many other secondary metabolites, is influenced by various phytohormones (<xref ref-type="bibr" rid="B53">Zhu and Zhang, 2009</xref>). As one of the major phytohormones, ET is widely used in agricultural production practices including natural raw latex production (<xref ref-type="bibr" rid="B30">Montoro et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B3">An et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B14">Fern&#xe1;ndez-Cancelo et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B38">Shu et&#xa0;al., 2023</xref>). ET releaser Ethephon or Ethrel are commonly used stimulators of rubber production worldwide and usually increase latex yield by 1.5- to 2-fold (<xref ref-type="bibr" rid="B10">Coup&#xe9; and Chrestin, 1989</xref>; <xref ref-type="bibr" rid="B53">Zhu and Zhang, 2009</xref>). A substantial amount of research in recent decades has focused on deciphering the mechanism of ethylene action. Concerning the physiology and biochemistry aspects of rubber trees, ET affects the permeability of cellular membranes, leading to increased and prolonged latex flow and general regenerative metabolism (<xref ref-type="bibr" rid="B10">Coup&#xe9; and Chrestin, 1989</xref>; <xref ref-type="bibr" rid="B29">Lu et&#xa0;al., 2023</xref>). As a side effect, ET treatment causes senescence of rubber tree seedlings, such as yellowing and abscised leaves (<xref ref-type="bibr" rid="B31">Nakano et&#xa0;al., 2021</xref>). Simultaneously, the enzyme activities of glutamine synthetase (GS) and chitinase are upregulated, and many metabolites accumulate in laticifers (<xref ref-type="bibr" rid="B10">Coup&#xe9; and Chrestin, 1989</xref>; <xref ref-type="bibr" rid="B2">Amalou et&#xa0;al., 1991</xref>; <xref ref-type="bibr" rid="B33">Pujade-Renaud et&#xa0;al., 1994</xref>). Paradoxically, stimulation with ET does not affect the transcription level of certain genes involved in rubber biosynthesis or even downregulated (<xref ref-type="bibr" rid="B8">Chye et&#xa0;al., 1992</xref>; <xref ref-type="bibr" rid="B39">Sirinupong et&#xa0;al., 2005</xref>). Indeed, it is generally accepted that ET does not directly promote rubber biosynthesis, instead many antioxidant genes are upregulated upon exposure of rubber trees to ET. To reveal the molecular mechanisms underlying ET-induced stimulation of NR yield, much work has been performed with the latex or bark of rubber trees based on signal transcriptomics or proteomics analysis (<xref ref-type="bibr" rid="B48">Wang et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B11">Dai et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B25">Li et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B31">Nakano et&#xa0;al., 2021</xref>). Although these studies have identified numerous genes and proteins involved in response to ET, much information remains to be delineated to achieve a complete understanding of the mechanism by which ET improves NR production.</p>
<p>Owing to the complexity of the ET stimulation of rubber latex yield, a multi-omics analysis might help to uncover the mystery of ET in rubber trees. Therefore, we analyzed the physiological indicators of latex yield and generated metabolomes and transcriptomes for latex or bark from four consecutive tapping with or without ET stimulation. The results provide new insights into the molecular mechanism of the ET-induced increase in rubber yield.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Plant materials and ET treatments</title>
<p>CATAS73397, an elite rubber tree cultivar was planted at the experimental farm of the Chinese Academy of Tropical Agricultural Sciences in Danzhou City, Hainan Province, (19&#xb0;51&#x2019; 51N; 109&#xb0;55&#x2019; 63E). All the trees were planted in 2007 and had their first tapping in 2014 with the latex harvesting system S/2 d3 (half-spiral tapping with a tapping frequency of 3 days), and trees with the same physical growth status were selected for treatment with ET. Before treatment, all trees were tapped to collect latex and bark samples for the negative control (ET0), and then the stimulation treatments were performed with 1.5% aqueous ethephon applied on the tapping panel by brush 3 days later. After 24 h, these trees were tapped to collect latex and bark samples (ET24, 24 h post-exposure to ET). After an additional 3 days, the trees were again tapped to collect latex and bark samples (ET1T, 96 h post-exposure to ET). After an additional 3 days, the trees were tapped again (ET2T, 168 h post-exposure to ET). All the treatments comprised three independent biological replicates, and every fifth tree served as a biological replicate. The latex samples were collected by a tumble placed on ice after tapping for the physiological character test and metabolomics analysis immediately, and the bark samples were collected and kept in liquid nitrogen for the following assays.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Latex yield and physiological indices test</title>
<p>Fresh latex (1 mL) was mixed with 9 mL of 2.5% trichloroacetic acid (TCA), and the mixture was centrifuged for 10 min at 4000 &#xd7; <italic>g</italic>. The resultant rubber clump and TCA-treated serum were separated, and then the treated serum was used for further analysis after filtration through qualitative filter paper. Latex physiological index, i.e., dry-rubber content, sucrose content (mM), and inorganic phosphate content (mM) were measured following the method of <xref ref-type="bibr" rid="B23">Koryati and Elfiana (2019)</xref>. Latex yield was converted into dry-rubber yield (gram tree<sup>&#x2212;1</sup> tap<sup>&#x2212;1</sup>) to remove the effect of water contents.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Metabolite extraction and metabolomics analysis</title>
<p>Samples of fresh latex were centrifuged at 15,000 &#xd7; g for 45 min at 4&#xb0;C, and serum was collected for widely targeted metabolomics analysis. The extracted serum samples were sent to Shanghai Biotree Biotech for metabolomic sequencing and analysis. Briefly, SIMCA software (V16.0.2, Sartorius Stedim Data Analytics AB, Umea, Sweden) was used for the principal component analysis (PCA) and orthogonal projections to latent structures-discriminate analysis&#x2019; (OPLS-DA). Differentially accumulated metabolites (DAMs) were identified based on a Student&#x2019;s t-test p-value of less than 0.05 and &#x2018;variable importance in the projection&#x2019; of the OPLS-DA model greater than 1. DAMs were used to generate hierarchical clustering plots based on the Euclidean distance method, and complete linkage and heat maps were created using the Pheatmap package in R Studio. Additionally, commercial databases including KEGG (<ext-link ext-link-type="uri" xlink:href="http://www.genome.jp/kegg/">http://www.genome.jp/kegg/</ext-link>) and MetaboAnalyst (<ext-link ext-link-type="uri" xlink:href="http://www.metaboanalyst.ca/">http://www.metaboanalyst.ca/</ext-link>) were used for pathway enrichment analysis (<xref ref-type="bibr" rid="B19">Kanehisa et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B49">Xia et&#xa0;al., 2015</xref>).</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>RNA preparation and transcriptome sequencing</title>
<p>The bark samples taken at ET0, ET24, ET1T, and ET2T were ground in liquid nitrogen, and total RNA was extracted according to <xref ref-type="bibr" rid="B50">Yang et&#xa0;al. (2020)</xref>. Sangon Biotech Corporation (Shanghai, China) constructed the sequencing library and carried out Illumina sequencing. After filtering out reads containing the adaptor sequence and deleting low-quality reads, the clean reads were mapped to the <italic>H. brasiliensis</italic> genome using HISAT2.2.4 tools, and statistical analysis of the alignment results was conducted by RseQC (<xref ref-type="bibr" rid="B47">Wang et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B20">Kim et&#xa0;al., 2015</xref>). Differentially expressed genes (DEGs) were analyzed using the DESeq program of the R software package with parameters q-value &lt; 0.05 and |log2 (fold change)| &#x2265; 1 (<xref ref-type="bibr" rid="B28">Love et&#xa0;al., 2014</xref>). Functional analysis of the obtained DEGs was performed per our previous study (<xref ref-type="bibr" rid="B50">Yang et&#xa0;al., 2020</xref>).</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>Quantitative PCR analysis</title>
<p>A total of 6 unigenes were selected to test the results of RNA-Seq by real-time quantitative PCR (RT-qPCR). <italic>HbActin</italic> was used as an internal control. All RT-qPCR primers were designed using Primer Premier 5 software (PREMIER Biosoft, Palo Alto, CA, USA) and synthesized by Sangon. Primer pairs are listed in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;1</bold>
</xref>. RNA was reverse transcribed using the PrimeScript&#x2122; RT reagent kit with gDNA Eraser (Perfect Real Time; Takara, Da Lian, China). The amplified PCR products were monitored with the CFX96 Real-Time PCR detection System (Bio-Rad). All RT-qPCR experiments were carried out three times in independent runs for all references and selected genes. Ct values represent the mean &#xb1; SD of three biological replicates, and relative gene expression was evaluated using the 2<sup>&#x2013;&#x394;&#x394;Ct</sup> method (<xref ref-type="bibr" rid="B27">Livak and Schmittgen, 2001</xref>).</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Statistical analysis</title>
<p>Statistical analysis was performed with SPSS statistics 27 (IBM Corporation, Armonk, NY, USA). One-way analysis of variance and the Student-Newman-Keuls test were performed for comparisons between multiple groups.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Latex physiological changes of the trees under ET stimulation</title>
<p>The physiological parameters of latex usually reflect the latex metabolic status of rubber trees. To study how stimulation with ethephon affects these physiological parameters, we designed ET stimulation experiments and the dry-rubber yield and characteristics of physiological parameters such as sucrose and inorganic phosphate at multiple time points were evaluated (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Compared with ET0 (before ET exposure, control), treatment with ET significantly increased the dry-rubber yield, especially at ET24 and ET1T (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>). The sucrose content increased continuously up to ET2T (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>). Slightly different from sucrose, a sustained increase of inorganic phosphate was observed starting at ET24 to ET2T (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1D</bold>
</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Experimental design and physiological indicators of rubber tree latex. <bold>(A)</bold> Experimental design for ET stimulation and sample collection. <bold>(B)</bold> The dry-rubber yield of rubber trees in response to ET. <bold>(C)</bold> The characteristics of sucrose (Suc) in rubber-tree latex during ET stimulation. <bold>(D)</bold> The characteristics of inorganic phosphate (Pi) in rubber-tree latex during ET stimulation. * indicated the significant difference at p&lt;0.05 and ** indicated the significant difference at p&lt;0.01.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g001.tif"/>
</fig>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Identification of DEGs involved in the ET response</title>
<p>To identify DEGs involved in the ET response of <italic>H. brasiliensis</italic>, rubber-tree bark was subjected to RNA sequencing (RNA-seq), and the raw data were deposited into the Genome Sequence Archive in the National Genomics Data Center, China National Center for Bioinformation, Chinese Academy of Sciences (accession number CRA012198). To verify the reliability of the RNA-seq data, six DEGs from different pathways were randomly selected for real-time quantitative PCR (RT-qPCR). The results showed that all the selected DEGs showed similar expression profiles in RNA-seq and RT-qPCR analyses, indicating that the RNA-seq data were accurate, credible, and highly reproducible (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Further, the number of DEGs was compared between each treatment group. Compared with ET0, 3400 DEGs were identified for ET24 (709 upregulated, 2691 downregulated), 2595 for ET1T (445 upregulated, 2150 downregulated), and 683 for ET2T (170 upregulated, 513 downregulated) (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;2</bold>
</xref>). These results indicated that the transcription of most DEGs was suppressed by ET, especially at ET24 and ET1T (<xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A, B</bold>
</xref>). Based on K-means cluster analysis of their expression patterns, the DEGs were divided into six clusters (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). Compared with ET0, 866 DEGs in cluster 1 displayed a negative response to ET (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3C</bold>
</xref>). Moreover, 261 and 30 DEGs were downregulated or upregulated at all three time points relative to ET0 (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Figure&#xa0;1</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;3</bold>
</xref>). Among the common ET-upregulated DEGs, four (LOC110661531, LOC110638059, LOC110652110, LOC110668984) encode chitinase and one encodes pathogenesis-related protein 1 (LOC110662812) (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;3</bold>
</xref>). KEGG analysis of DEGs from comparisons ET24/ET0, ET1T/ET0 and ET2T/ET0 revealed that most DEGs were enriched in the pathways for starch and sucrose metabolism (ko00500), phytohormone signal transduction (ko04075), and amino-acid metabolism (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Relative expression level of six DEGs obtained by RT-qPCR and RNA-seq. The columns and lines represent log 2 (Fold change) determined by RT-qPCR and RNA-seq.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>DEGs involved in the ET response. <bold>(A)</bold> Number of DEGs identified between different treatments of rubber trees (horizontal axis). Up, upregulated; Down, downregulated. <bold>(B)</bold> Heat map of DEGs identified upon stimulation with ET (ET0, control). <bold>(C)</bold> K-means cluster analysis of DEGs.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Analysis of KEGG terms enriched among the DEGs. The bubble diagram shows the significant enrichment KEGG pathways of the DEGs from ET24/ET0 <bold>(A)</bold>, ET1T/ET0 <bold>(B)</bold>, and ET2T/ET0 <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g004.tif"/>
</fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Identification of phytohormone-signaling genes and transcription-factor genes among the candidate DEGs</title>
<p>Exogenous ET application may alter and rebalance the cellular abundance of other major phytohormones such as auxin (indole acetic acid, IAA), cytokinin, gibberellin, abscisic acid, ET, brassinosteroid, jasmonate, and salicylic acid. To investigate the crosstalk between the ET signaling pathway and pathways governed by other phytohormones, the expression patterns involved in their biosynthesis and signal transduction pathways were analyzed (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A&#x2013;H</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>). A total of 134 target genes were identified, including 36 involved in IAA signaling, 11 in salicylic acid and cytokinin signaling, 10 in gibberellin signaling, 22 in abscisic acid signaling, 17 in brassinosteroid signaling, 14 in ET signaling, and 13 in jasmonate signaling (<xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A&#x2013;H</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). Based on expression relative to ET0, these genes were divided into four clusters (<xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). The genes in cluster 1 were upregulated in response to ET at ET24, whereas these genes were downregulated at ET1T and ET2T, including one TIR1 and four SAUR for IAA, one GA2ox and one PIF4 for gibberellin, four PYR/PYL and one PP2C for abscisic acid, two ET response factors (ERF1/2), three S-adenosylmethionine synthases (SAMS) and one ACO for ET, one fadA and one JAR1 for jasmonate, and finally two 4CL and three PR-1 for salicylic acid (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). The genes in clusters 2 and 4 were downregulated in response to ET at ET24, whereas these genes were upregulated at ET1T and ET2T, including one TIR1, two SAUR, four GH3, seventeen AUX/IAA, four ARF, and two YUC for IAA, two CRE1, two B-ARR, six A-ARR, and one AHP for cytokinin, two GA2ox and one PIF4 for gibberellin, one AAO3, three NCED, one CYP70A, one SnRK2, two PP2C, and two ABF for abscisic acid, two ERF1/2, four ACS, and one ACO for ET, one CPS, four TCH4, five CYCD3, two BRI1, one BKI1, and one BAK1 for brassinosteroid, three LOX, one AOS, one OPR, and four JAZ for jasmonate, and finally one TGA, two PR-1, and three NPR1 for salicylic acid (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). The remaining genes were in cluster 3, in which genes were upregulated at ET24 and ET2T but downregulated at ET1T. Interestingly, the expression of two key ET biosynthesis genes (ACS and ACO) was suppressed by exogenous ET, whereas the expression of two downstream regulatory genes (EBF1/2 and ERF1/2) of the ET signaling pathway was increased in response to ET (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF2">
<bold>Supplementary Figure&#xa0;2</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>).</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Phytohormone signaling and transcription factor genes involved in ET stimulation. Expression of genes involved in the biosynthesis of IAA <bold>(A)</bold>, SA <bold>(B)</bold>, CK <bold>(C)</bold>, GA <bold>(D)</bold>, ABA <bold>(E)</bold>, BR <bold>(F)</bold>, ET <bold>(G)</bold>, and JA <bold>(H)</bold> and heat map of transcription factors <bold>(I)</bold> involved in the ET response. The heat map was drawn with Tbtools software (<xref ref-type="bibr" rid="B7">Chen et&#xa0;al., 2020</xref>). Transcript level data from A to H were normalized using Z-score calculation. IAA, indole acetic acid; SA, salicylic acid; CK, cytokinin; GA, gibberellin; ABA, abscisic acid; BR, brassinosteroid; ET, ethylene; JA, jasmonic acid.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g005.tif"/>
</fig>
<p>To screen for transcription factors that are potentially involved in the ET response, we identified putative transcription factors among the DEGs, and a total of 101 transcription factors from 18 transcription-factor families were detected. The MYB superfamily was the largest member with 32 genes, followed by HD_ZIP with 16, GATA with 10, HSF and MEIS1 with 8, bHLH with 7, bZIP_ATF with 3, CREB/ATF with 2, UAF with 2, MOT2 with 2, and HLTF, TFIIIC, UBA_NAC, E2F/DP, bZIP_CPRF2, and trihelix with 1 (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5I</bold>
</xref>). With the parameters q-value &lt; 0.05 and |log2 (fold change)| &#x2265; 2, 51 transcription factors were identified as major ET-responsive transcription factors, including 3 MADS, 14 MYB, 10 HD-ZIP, 5 HSF, 2 MEIS1, 1 HLTF, 3 bHLH, 7 GATA, 1 UAF, 1 E2F/DP, 1 bZIP_HY5, and 3 bZIP_ATF6 (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;5</bold>
</xref>). These results implied that these genes likely participate in metabolic processes that increase rubber yield under ET stimulation.</p>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Metabolomic response of rubber trees under ET stimulation</title>
<p>Widely targeted metabolomics was performed to investigate changes in metabolites during ET treatment. Then, we facilitated statistical analysis of the number of common and unique DAMs in each contrast combination (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>; <xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>). In the comparison combination of ET24/ET0, there were a total of 175 DAMs (<xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>). Similarly, there were 104 DAMs in the ET1T/ET0 comparison combination and 189 in the ET2T/ET0 comparison combination (<xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>). Additionally, there were 16 DAMs shared between ET24/ET0, ET1T/ET0 and ET2T/ET0. Furthermore, there were 39, 18, and 10 compounds shared between ET24/ET0 and ET1T/ET0, ET24/ET0 and ET24/ET0, ET1T/ET0 and ET2T/ET0, respectively (<xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>). The unique DAMs in the comparison combination of ET24 <italic>vs</italic> ET0, ET1T <italic>vs </italic>ET0 and ET2T <italic>vs</italic> ET0 were 102, 39, and 145, respectively (<xref ref-type="supplementary-material" rid="SF3">
<bold>Supplementary Figure&#xa0;3</bold>
</xref>). To further analyze the composition and enrichment of DAMs in each contrast combination, we facilitated statistical analysis of the number of major classes of the detected metabolites and visualized by heat maps. For contrast combination ET24/ET0, the 175 DAMs were mainly enriched in nucleosides, nucleotides, lipids, organic acids, and organic oxygen compounds (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6A</bold>
</xref>). For contrast combination ET1T/ET0, the 104 DAMs were mainly enriched in lipids, organic oxygen compounds, and organic acids (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6B</bold>
</xref>). For contrast combination ET2T/ET0, the 189 DAMs were mainly enriched in organic acids (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6C</bold>
</xref>). These results indicated that the metabolic levels of organic acids and organic oxygen compounds might be highly relevant to the ET-induced stimulation of latex production in rubber trees.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Number and heat map of DAMs involved in the ET response. The number of major classes of the detected metabolites and heat maps for ET24/ET0 <bold>(A)</bold>, ET1T/ET0 <bold>(B)</bold>, and ET2T/ET0 <bold>(C)</bold>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g006.tif"/>
</fig>
<p>We next carried out KEGG annotation of the DAMs identified at ET24/ET0, ET1T/ET0, and ET2T/ET0. For ET24, 50 pathways were identified, with notable enrichment in metabolic pathways (hbr01100), biosynthesis of secondary metabolites (hbr01110), biosynthesis of amino acids (hbr01230), ABC transporters (hbr02010), and galactose metabolism (hbr00052) (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). ET1T had 21 pathways, with notable enrichment in metabolic pathways (hbr01100), biosynthesis of secondary metabolites (hbr01110), ABC transporters (hbr02010), biosynthesis of amino acids (hbr01230), and galactose metabolism (hbr00052) (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>). ET2T had 29 pathways, with notable enrichment in metabolic pathways (hbr01100), biosynthesis of secondary metabolites (hbr01110), ABC transporters (hbr02010), biosynthesis of amino acids (hbr01230), purine metabolism (hbr00230), and galactose metabolism (hbr00052) (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>).</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Key KEGG pathways in the metabolome.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="center">Pathway ID</th>
<th valign="middle" rowspan="2" align="center">Description</th>
<th valign="middle" colspan="3" align="center">Number of DAMs</th>
</tr>
<tr>
<th valign="middle" align="center">ET24/ET0</th>
<th valign="middle" align="center">ET1T/ET0</th>
<th valign="middle" align="center">ET2T/ET0</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">hbr01100</td>
<td valign="middle" align="center">Metabolic pathways</td>
<td valign="middle" align="center">59</td>
<td valign="middle" align="center">27</td>
<td valign="middle" align="center">38</td>
</tr>
<tr>
<td valign="middle" align="center">hbr01110</td>
<td valign="middle" align="center">Biosynthesis of secondary metabolites</td>
<td valign="middle" align="center">24</td>
<td valign="middle" align="center">11</td>
<td valign="middle" align="center">14</td>
</tr>
<tr>
<td valign="middle" align="center">hbr01230</td>
<td valign="middle" align="center">Biosynthesis of amino acids</td>
<td valign="middle" align="center">14</td>
<td valign="middle" align="center">6</td>
<td valign="middle" align="center">7</td>
</tr>
<tr>
<td valign="middle" align="center">hbr02010</td>
<td valign="middle" align="center">ABC transporters</td>
<td valign="middle" align="center">14</td>
<td valign="middle" align="center">8</td>
<td valign="middle" align="center">9</td>
</tr>
<tr>
<td valign="middle" align="center">hbr00052</td>
<td valign="middle" align="center">Galactose metabolism</td>
<td valign="middle" align="center">9</td>
<td valign="middle" align="center">5</td>
<td valign="middle" align="center">4</td>
</tr>
<tr>
<td valign="middle" align="center">hbr00230</td>
<td valign="middle" align="center">Purine metabolism</td>
<td valign="middle" align="center">5</td>
<td valign="middle" align="center">&#x2013;</td>
<td valign="middle" align="center">6</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>-represents undetected.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_5">
<label>3.5</label>
<title>Integrated analysis of the transcriptome and metabolome involved in the ET response of rubber trees</title>
<p>To investigate how treatment with ET increases latex yield, we performed an integrated analysis of the metabolome and transcriptome of ET-treated rubber trees. As shown in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>, amino-acid biosynthesis and carbohydrate metabolism were identified as the key pathways related to the ET-induced stimulation of latex production.</p>
<p>The biosynthesis of amino acids is essential for plant growth and development, including responses to biotic and abiotic stresses (<xref ref-type="bibr" rid="B51">Zeier, 2013</xref>; <xref ref-type="bibr" rid="B15">Galili et&#xa0;al., 2016</xref>). At the metabolism level, 21 types of metabolites related to amino-acid biosynthesis were identified (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;6</bold>
</xref>). For ET24/ET0, the following compounds were found to have accumulated: <sc>l</sc>-lysine, <sc>l</sc>-aspartic acid, phosphoenolpyruvic acid, <sc>l</sc>-ornithine, <sc>l</sc>-cysteine, dihydroxyacetone phosphate, <sc>l</sc>-leucine, <sc>l</sc>-histidine, <sc>l</sc>-valine, <sc>l</sc>-threonine, <sc>l</sc>-homoserine, <sc>d</sc>-erythrose 4-phosphate, and o-phosphoserine 2, whereas <sc>l</sc>-tyrosine and <sc>l</sc>-cysteine levels decreased (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). However, the abundance of most of these metabolites was decreased, except serine, cysteine, o-phosphoserine 2, alanine, and proline (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>). Expression was altered for 54 genes related to amino-acid biosynthesis (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;7</bold>
</xref>). The expression of 27 genes was significantly upregulated in response to ET signaling, especially early during stimulation (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;7</bold>
</xref>). The expression of 16 genes was significantly downregulated in response to ET signaling, especially early during stimulation (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;7</bold>
</xref>). Therefore, these genes and metabolites may be involved in the early ET response by regulating amino-acid biosynthesis. The results of the Pearson correlation analysis between the expression levels of the candidate genes and the accumulation of amino acids revealed that 23 genes were significantly positively correlated to at least 9 metabolites related to amino-acid biosynthesis (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref>). Thus, these genes may play key roles in regulating amino-acid biosynthesis under ET stimulation in rubber trees.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Expression profiles and Pearson correlation analysis of DEGs and DAMs relevant to amino-acid biosynthesis. <bold>(A)</bold> Schematic diagram of major genes and metabolites involved in amino-acid biosynthesis. Transcript level data were normalized using Z-score calculation. <bold>(B)</bold> Pearson correlation analysis between genes and metabolites involved in amino-acid biosynthesis. Genes that were significantly upregulated or downregulated in response to ET are indicated with *(p &lt; 0.05) or **(p &lt; 0.01).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g007.tif"/>
</fig>
<p>Rubber biosynthesis requires a substantial amount of metabolic energy due to its high energy consumption, and carbohydrate metabolism provides energy and carbon skeletons (for use in biosynthesis) for this process. (<xref ref-type="bibr" rid="B41">Tang et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B6">Chao et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B13">Fang et&#xa0;al., 2021</xref>) At the metabolomic level, the content of each of raffinose, UDP-galactose, melibiose, and &#x3b1;-galactose 1-phosphate was increased, whereas the content of sucrose was decreased (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;8</bold>
</xref>). However, the expression of seven DEGs related to carbohydrate metabolism was downregulated at ET24 and even more so at ET1T and ET2T (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;9</bold>
</xref>).</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Expression profiles of DEGs and DAMs involved in carbohydrate metabolism. Transcript level data were normalized using Z-score calculation. GOLS, inositol 3-alpha-galactosyltransferase; UGE, UDP-glucose 4-epimerase; RS, raffinose synthase; GLA, alpha-galactosidase; INV, beta-fructofuranosidase.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g008.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>Identifying the key metabolic compounds that contribute to the ET-mediated rubber yield increase and revealing their underlying transcriptomic and transcriptional regulatory networks are important steps toward improving both NR yield and quality. Although much research has been based on the analysis of the transcriptome and/or proteome underlying the cellular signaling pathways of the latex or bark of rubber trees, the molecular mechanism of the ET-mediated enhancement of latex production and the effects of ET on NR quality remain a mystery. Here, through integrative analyses of metabolic and transcriptomic data across four consecutive tappings of rubber trees with or without ET treatment, we revealed the changes in global metabolites and genes that are regulated by ET.</p>
<sec id="s4_1">
<label>4.1</label>
<title>ET-mediated rubber yield increase is a stress or defense response process</title>
<p>ET is a gaseous phytohormone that promotes fruit ripping, leaf senescence, abscission, and stress response in plants (<xref ref-type="bibr" rid="B37">Schaller, 2012</xref>; <xref ref-type="bibr" rid="B18">Iqbal et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B16">Gambhir et&#xa0;al., 2023</xref>). In many plants, the synthesis of ET is associated with biotic and abiotic stresses (<xref ref-type="bibr" rid="B43">Tao&#xa0;et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B44">Thao et&#xa0;al., 2015</xref>). In this respect, ET is often called &#x201c;stress ET&#x201d; (<xref ref-type="bibr" rid="B4">Bakshi et&#xa0;al., 2015</xref>). Over many decades of NR production, ET has been the most effective and commonly used enhancer of latex yield, but the adverse effects of ET should not be ignored, including tapping panel dryness (<xref ref-type="bibr" rid="B32">Osborne and Sargent, 1974</xref>). Like other latex-producing plants, the NR of the rubber tree is synthesized and stored in a specialized cellular organelle called the laticifer, which contributes to the defense strategy of plants (<xref ref-type="bibr" rid="B22">Konno, 2011</xref>; <xref ref-type="bibr" rid="B17">Hua et&#xa0;al., 2017</xref>). In laticifers, many latex proteins, such as peptidases, chitinases, lectins, pathogenesis-related proteins, and protease inhibitors, function as defense molecules to protect plants from stress-inducing agents such as fungi, insects, viruses, and wounding (<xref ref-type="bibr" rid="B17">Hua et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B21">Kitajima et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B36">Ramos et&#xa0;al., 2019</xref>). As such, in rubber trees, laticifers exude latex for defensive purposes, and humans have exploited the resultant increased turgor pressure by tapping trees. Our present results revealed that the expression of numerous genes was repressed by more than 75% in response to ET, e.g., 79% for ET24/ET0, 83% for ET1T/ET0, and 75% for ET2T/ET0 (<xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3A, B</bold>
</xref>). In contrast, the genes encoding chitinase and PR1 were oftentimes upregulated by ET (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;3</bold>
</xref>), suggesting that ET constitutes a potent defensive signal. This may explain why overstimulation with ET or overtapping causes tapping panel dryness, a complex physiological response of the rubber tree that results in the cessation of latex flow during tapping.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Rubber yield increase under ET stimulation is a result of the cross-talk between different hormones</title>
<p>Endogenous phytohormone balance plays a vital role in latex regeneration and the duration of latex flow (<xref ref-type="bibr" rid="B46">Tungngoen et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B5">Chao et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B12">Deng et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B30">Montoro et&#xa0;al., 2018</xref>). However, how different phytohormones work together to regulate these two processes remains unclear. In our study, the biosynthesis and signaling pathways of ET and all the other phytohormones were particularly affected by ET (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table&#xa0;4</bold>
</xref>). For ET biosynthesis, treatment with ET triggered the expression of SAMS, whereas the expression of ACS and ACO were repressed. Although the rubber-tree genome contains 181 genes encoding ERFs, only one <italic>ERF</italic> gene was dramatically upregulated. These results suggest that ET does not promote the biosynthesis of endogenous ET but rather ET signal transduction&#x2014;especially at the early stage of ET stimulation. Similar results have been observed in the latex of rubber-tree clones PB206 and Reyan7-33-97 (<xref ref-type="bibr" rid="B42">Tang&#xa0;et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B30">Montoro et&#xa0;al., 2018</xref>). In addition, ET inhibits the expression of most genes related to phytohormone biosynthesis and signaling pathways. This suggests that the increase in yield stimulated by ET is the result of a combination of many phytohormones. Evidence reported by <xref ref-type="bibr" rid="B46">Tungngoen et&#xa0;al. (2011)</xref> revealed that auxin, salicylic acid, abscisic acid, and ET had identical effects on latex yield and latex total solid content in clone PB217. Together, research to date has demonstrated that treatment with ET increases latex yield by facilitating cross-talk between signaling pathways modulated by various phytohormones.</p>
</sec>
<sec id="s4_3">
<label>4.3</label>
<title>Amino-acid biosynthesis and carbohydrate metabolism involved in latex flow under ET stimulation</title>
<p>The effect of ET stimulation on amino-acid biosynthesis in rubber-tree latex remains unknown. Our results reveal for the first time the regulatory network of amino-acid biosynthesis in latex in response to ET stimulation using integrated metabolomic and transcriptomic analyses. The results show that ET altered the cellular content of each of the 21 amino acids or their derivatives, and 54 genes relevant to amino-acid biosynthesis were either upregulated or downregulated (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>; <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Tables&#xa0;6</bold>
</xref>, <xref ref-type="supplementary-material" rid="SM1">
<bold>7</bold>
</xref>). Pearson correlation analysis suggested that GS, PK, SAMS, SDH, DAHPS, CYSS, IPMS, METH, NAGK, PFK, PHGDH, PSAT, RPIA, and SHMT are the key genes controlling amino-acid biosynthesis under ET stimulation (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref>). In rubber trees, NR biosynthesis requires a substantial amount of metabolic energy and induces a water deficit in laticifer cells owing to latex flow (<xref ref-type="bibr" rid="B45">Tungngoen et&#xa0;al., 2009</xref>; <xref ref-type="bibr" rid="B41">Tang et&#xa0;al., 2010</xref>). In all plants, i.e., including the rubber tree, raffinose-family oligosaccharides are critical for resistance to osmotic stress to assist in resistance to osmotic stress (<xref ref-type="bibr" rid="B26">Lintunen et&#xa0;al., 2016</xref>). In our present study, the cellular content of metabolites and altered expression of genes involved in the biosynthesis of raffinose-family oligosaccharides were altered by ET stimulation (<xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>). Thus, the raffinose-family oligosaccharides would be involved in response to the water deficit regulation induced by ET stimulation during latex flow.</p>
</sec>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusion</title>
<p>This study clarifies the molecular mechanisms of the ethylene-induced increase in latex yield in the rubber tree utilizing a multi-omics approach. We developed a potential underlying model based on our physiological, transcriptomic, and metabolomic analysis to elucidate the ethylene-mediated natural rubber yield increase (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref>). Briefly, exogenous ethylene acts as a stress signal and suppresses the expression of many genes related to primary metabolism by inducing endogenous hormone rebalance. The phytohormone then activates the expression of certain transcription factors, which regulate amino acid biosynthesis and carbohydrate metabolism. Amino acids provide sufficient proteins and enzymes for the NR biosynthesis process, which is essential for increasing rubber production. Ethylene promoted the biosynthesis of sucrose, galactose, raffinose and melibiose, indicating that these carbohydrates provide the rubber tree with high-energy substrates for NR biosynthesis under ethylene stimulation. These findings provide new insights into the molecular mechanism of the ethylene-induced increase in rubber yield.</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Putative model for ethylene-mediated natural rubber yield increase.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1444693-g009.tif"/>
</fig>
</sec>
</body>
<back>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material</bold>
</xref>.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>HY: Conceptualization, Data curation, Funding acquisition, Methodology, Visualization, Writing &#x2013; original draft, Writing &#x2013; review &amp; editing. LD: Data curation, Methodology, Visualization, Writing &#x2013; original draft. ML: Data curation, Methodology, Writing &#x2013; review &amp; editing. XF: Data curation, Validation, Writing &#x2013; original draft. LL: Data curation, Validation, Writing &#x2013; original draft. BG: Validation, Writing &#x2013; original draft. ZW: Funding acquisition, Writing &#x2013; original draft. LW: Supervision, Writing &#x2013; review &amp; editing.</p>
</sec>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Major Science and Technology Plan of Hainan Province (ZDKJ2021004) and the Natural Science Foundation of Hainan Province (320QN350).</p>
</sec>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2024.1444693/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2024.1444693/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Image1.jpeg" id="SF1" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;1</label>
<caption>
<p>Venn diagram for comparison of the numbers of genes commonly downregulated (left) or upregulated (right) in response to ethylene.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image2.jpeg" id="SF2" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;2</label>
<caption>
<p>Cluster analysis of DEGs involved in phytohormone biosynthesis.</p>
</caption>
</supplementary-material>
<supplementary-material xlink:href="Image3.jpeg" id="SF3" mimetype="image/jpeg">
<label>Supplementary Figure&#xa0;3</label>
<caption>
<p>Venn diagram of the number of DAMs.</p>
</caption>
</supplementary-material>
  <supplementary-material xlink:href="Table1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Adiwilaga</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kush</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Cloning and characterization of cDNA encoding farnesyl diphosphate synthase from rubber tree (<italic>Hevea brasiliensis</italic>)</article-title>. <source>Plant Mol. Biol.</source> <volume>30</volume>, <fpage>935</fpage>&#x2013;<lpage>946</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/BF00020805</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amalou</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Bangratz</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chrestin</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>1991</year>). <article-title>Ethrel (ethylene releaser)-induced increases in the adenylate pool and transtonoplast &#x394;pH within <italic>Hevea</italic> latex cells</article-title>. <source>Plant Physiol.</source> <volume>98</volume>, <fpage>1270</fpage>&#x2013;<lpage>1276</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.98.4.1270</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>An</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Althiab</surname> <given-names>A. R.</given-names>
</name>
<name>
<surname>Bouzayen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Zouine</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Chervin</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Auxin and ethylene regulation of fruit set</article-title>. <source>Plant Sci.</source> <volume>292</volume>, <elocation-id>110381</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.plantsci.2019.110381</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bakshi</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Shemansky</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Binder</surname> <given-names>B. M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>History of research on the plant hormone ethylene</article-title>. <source>J. Plant Growth Regu.</source> <volume>34</volume>, <fpage>809</fpage>&#x2013;<lpage>827</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00344-015-9522-9</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Tian</surname> <given-names>W. M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>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</article-title>. <source>BMC Plant Biol.</source> <volume>15</volume>, <fpage>104</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-015-0488-3</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Tian</surname> <given-names>W. M.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Transcript profiling of <italic>Hevea brasiliensis</italic> during latex flow</article-title>. <source>Front. Plant Sci.</source> <volume>8</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2017.01904</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Frank</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>He</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>TBtools: An integrative toolkit developed for interactive analyses of big biological data</article-title>. <source>Mol. Plant</source> <volume>13</volume>, <fpage>1194</fpage>&#x2013;<lpage>1202</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.molp.2020.06.009</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chye</surname> <given-names>M. L.</given-names>
</name>
<name>
<surname>Tan</surname> <given-names>C. T.</given-names>
</name>
<name>
<surname>Chua</surname> <given-names>N. H.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Three genes encode 3-hydroxy-3-methylglutaryl-coenzyme a reductase in <italic>Hevea brasiliensis</italic>: hmg1 and hmg3 are differentially expressed</article-title>. <source>Plant Mol. Biol.</source> <volume>19</volume>, <fpage>473</fpage>&#x2013;<lpage>484</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/bf00023395</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cornish</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Similarities and differences in rubber biochemistry among plant species</article-title>. <source>Phytochemistry</source> <volume>57</volume>, <fpage>1123</fpage>&#x2013;<lpage>1134</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/s0031-9422(01)00097-8</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Coup&#xe9;</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Chrestin</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>1989</year>). &#x201c;<article-title>Physico-chemical and biochemical mechanisms of hormonal (ethylene) stimulation</article-title>,&#x201d; in <source>Physiology of the Rubber Tree Latex</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>d'Auzac</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Jacob</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Chrestin</surname> <given-names>H.</given-names>
</name>
</person-group> (<publisher-name>CRC Press</publisher-name>, <publisher-loc>Boca Raton, FL</publisher-loc>), <fpage>295</fpage>&#x2013;<lpage>319</lpage>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname> <given-names>L. J.</given-names>
</name>
<name>
<surname>Kang</surname> <given-names>G. J.</given-names>
</name>
<name>
<surname>Nie</surname> <given-names>Z. Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>R. Z.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Comparative proteomic analysis of latex from <italic>Hevea brasiliensis</italic> treated with ethrel and methyl jasmonate using iTRAQ-coupled two-dimensional LC&#x2013;MS/MS</article-title>. <source>J. Proteomics</source> <volume>132</volume>, <fpage>167</fpage>&#x2013;<lpage>175</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jprot.2015.11.012</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deng</surname> <given-names>X. M.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>S. H.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>S. G.</given-names>
</name>
<name>
<surname>Chao</surname> <given-names>J. Q.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Comparative transcriptome analysis reveals phytohormone signalings, heat shock module and ROS scavenger mediate the cold-tolerance of rubber tree</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>4931</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-018-23094-y</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Long</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Fang</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A predominant isoform of fructokinase, <italic>HbFRK2</italic>, is involved in <italic>Hevea brasiliensis</italic> (para rubber tree) latex yield and regeneration</article-title>. <source>Plant Physiol. Biochem.</source> <volume>162</volume>, <fpage>211</fpage>&#x2013;<lpage>220</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.plaphy.2021.02.039</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fern&#xe1;ndez-Cancelo</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Mu&#xf1;oz</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Echeverr&#xed;a</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Larrigaudi&#xe8;re</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Teixid&#xf3;</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Munn&#xe9;-Bosch</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Ethylene and abscisic acid play a key role in modulating apple ripening after harvest and after cold-storage</article-title>. <source>Postharvest. Biol. Technol.</source> <volume>188</volume>, <fpage>111902</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.postharvbio.2022.11190</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Galili</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Amir</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Fernie</surname> <given-names>A. R.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The regulation of essential amino acid synthesis and accumulation in plants</article-title>. <source>Annu. Rev. Plant Biol.</source> <volume>67</volume>, <fpage>153</fpage>&#x2013;<lpage>178</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev-arplant-043015-112213</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gambhir</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Raghuvanshi</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Parida</surname> <given-names>A. P.</given-names>
</name>
<name>
<surname>Kujur</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Sharma</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Sopory</surname> <given-names>S. K.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Elevated methylglyoxal levels inhibit tomato fruit ripening by preventing ethylene biosynthesis</article-title>. <source>Plant Physiol.</source> <volume>192</volume>, <fpage>2161</fpage>&#x2013;<lpage>2184</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/plphys/kiad142</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hua</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>C. J.</given-names>
</name>
<name>
<surname>Jing</surname> <given-names>S. X.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>S. H.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S. H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Chemical profile and defensive function of the latex of <italic>Euphorbia peplus</italic>
</article-title>. <source>Phytochemistry</source> <volume>136</volume>, <fpage>56</fpage>&#x2013;<lpage>64</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.phytochem.2016.12.021</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Iqbal</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>N. A.</given-names>
</name>
<name>
<surname>Ferrante</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Trivellini</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Francini</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>M. I. R.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Ethylene role in plant growth, development and senescence: interaction with other phytohormones</article-title>. <source>Front. Plant Sci.</source> <volume>8</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2017.00475</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kanehisa</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Sato</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Kawashima</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Furumichi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Tanabe</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>KEGG as a reference resource for gene and protein annotation</article-title>. <source>Nucleic Acids Res.</source> <volume>44</volume>, <fpage>457</fpage>&#x2013;<lpage>462</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkv1070</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Langmead</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Salzberg</surname> <given-names>S. L.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>HISAT: a fast spliced aligner with low memory requirements</article-title>. <source>Nat. Methods</source> <volume>12</volume>, <fpage>357</fpage>&#x2013;<lpage>360</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nmeth.3317</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kitajima</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Aoki</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Shibata</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Nakajima</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Sakurai</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Yazaki</surname> <given-names>K.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Comparative multi-omics analysis reveals diverse latex-based defense strategies against pests among latex-producing organs of the fig tree (<italic>Ficus carica</italic>)</article-title>. <source>Planta</source> <volume>247</volume>, <fpage>1423</fpage>&#x2013;<lpage>1438</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00425-018-2880-3</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Konno</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Plant latex and other exudates as plant defense systems: roles of various defense chemicals and proteins contained therein</article-title>. <source>Phytochemistry</source> <volume>72</volume>, <fpage>1510</fpage>&#x2013;<lpage>1530</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.phytochem.2011.02.016</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Koryati</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Elfiana</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Production and physiology characteristics of some rubber clones in early tapping with growth regulator</article-title>. <source>IOP. Conf. Ser. Mater. Sci. Engi.</source> <volume>536</volume>, <elocation-id>12141</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1088/1757-899X/536/1/012141</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lehman</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Tuljittraporn</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Songtipya</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Uthaipan</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Sengloyluan</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Johns</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Influence of non-rubber components on the properties of unvulcanized natural rubber from different clones</article-title>. <source>Polymers</source> <volume>14</volume>, <elocation-id>1759</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/polym14091759</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. C.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>He</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Transcriptome analyses reveal molecular mechanism underlying tapping panel dryness of rubber tree (<italic>Hevea brasiliensis</italic>)</article-title>. <source>Sci. Rep.</source> <volume>6</volume>, <elocation-id>23540</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/srep23540</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lintunen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Paljakka</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Jyske</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Peltoniemi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Sterck</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Von Arx</surname> <given-names>G.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Osmolality and non-structural carbohydrate composition in the secondary phloem of trees across a latitudinal gradient in Europe</article-title>. <source>Front. Plant Sci.</source> <volume>7</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2016.00726</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Livak</surname> <given-names>K. J.</given-names>
</name>
<name>
<surname>Schmittgen</surname> <given-names>T. D.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Analysis of relative gene expression data using real-time quantitative PCR and the 2<sup>&#x2013;&#x394;&#x394;Ct</sup> method</article-title>. <source>Methods</source> <volume>25</volume>, <fpage>402</fpage>&#x2013;<lpage>408</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1006/meth.2001.1262</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Love</surname> <given-names>M. I.</given-names>
</name>
<name>
<surname>Huber</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Anders</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2</article-title>. <source>Genome Biol.</source> <volume>15</volume>, <elocation-id>550</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s13059-014-0550-8</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname> <given-names>Y. L.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>K. Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y. Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X. D.</given-names>
</name>
<name>
<surname>Cao</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Expressional correlation of latex-regeneration genes with latex physiological parameters varies with population type and growth stage in <italic>Hevea</italic> tree</article-title>. <source>Trees</source> <volume>37</volume>, <fpage>1107</fpage>&#x2013;<lpage>1119</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00468-023-02410-8</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Montoro</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Favreau</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Herlinawati</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Labrune</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Martin-Magniette</surname> <given-names>M. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Transcriptome analysis in <italic>Hevea brasiliensis</italic> latex revealed changes in hormone signaling pathways during ethephon stimulation and consequent tapping panel dryness</article-title>. <source>Sci. Rep.</source> <volume>8</volume>, <fpage>8483</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-018-26854-y</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakano</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Mitsuda</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Ide</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Mori</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Mira</surname> <given-names>F. R.</given-names>
</name>
<name>
<surname>Rosmalawati</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Transcriptome analysis of Par&#xe1; rubber tree (<italic>H. brasiliensis</italic>) seedlings under ethylene stimulation</article-title>. <source>BMC Plant Biol.</source> <volume>21</volume>, <fpage>420</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-021-03196-y</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Osborne</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Sargent</surname> <given-names>J. A.</given-names>
</name>
</person-group> (<year>1974</year>). <article-title>A model for the mechanism of stimulation of latex flow in <italic>Hevea brasiliensis</italic> by ethylene</article-title>. <source>Ann. Appl. Biol.</source> <volume>78</volume>, <fpage>83</fpage>&#x2013;<lpage>88</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1744-7348.1974.tb01488.x</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pujade-Renaud</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Clement</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Perrot-Rechenmann</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Prevot</surname> <given-names>J. C.</given-names>
</name>
<name>
<surname>Hrestin</surname> <given-names>H. C.</given-names>
</name>
<name>
<surname>Jacob</surname> <given-names>J. L.</given-names>
</name>
<etal/>
</person-group>. (<year>1994</year>). <article-title>Ethylene-induced increase in glutamine synthetase activity and mRNA levels in <italic>Hevea brasiliensis</italic> latex cells</article-title>. <source>Plant Physiol.</source> <volume>105</volume>, <fpage>127</fpage>&#x2013;<lpage>132</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.105.1.127</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qi</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Yun</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Z. X.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>The main service functions and driving forces of rubber (<italic>Hevea brasiliensis</italic>) plantation ecosystem in China</article-title>. <source>J. Rubber. Res.</source> <volume>26</volume>, <fpage>155</fpage>&#x2013;<lpage>164</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s42464-023-00202-w</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qin</surname> <given-names>Y. X.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fang</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>X. Y.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>J. H.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Anaerobic metabolism in <italic>Hevea brasiliensis</italic> laticifers is relevant to rubber synthesis when tapping is initiated</article-title>. <source>Ind. Crop Prod.</source> <volume>178</volume>, <elocation-id>114663</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.indcrop.2022.114663</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramos</surname> <given-names>M. V.</given-names>
</name>
<name>
<surname>Demarco</surname> <given-names>D.</given-names>
</name>
<name>
<surname>da Costa Souza</surname> <given-names>I. C.</given-names>
</name>
<name>
<surname>de Freitas</surname> <given-names>C. D. T.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Laticifers, latex, and their role in plant defense</article-title>. <source>Trends Plant Sci.</source> <volume>24</volume>, <fpage>553</fpage>&#x2013;<lpage>567</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tplants.2019.03.006</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schaller</surname> <given-names>G. E.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Ethylene and the regulation of plant development</article-title>. <source>BMC Biol.</source> <volume>10</volume>, <elocation-id>9</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/1741-7007-10-9</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shu</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Sheng</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>
<italic>SlMAPK3</italic> positively regulates the ethylene production of postharvest tomato fruits and is involved in ethylene-mediated cold tolerance</article-title>. <source>J. Agric. Food Chem.</source> <volume>71</volume>, <fpage>6003</fpage>&#x2013;<lpage>6013</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1021/acs.jafc.2c08723</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sirinupong</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Suwanmanee</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Doolittle</surname> <given-names>R. F.</given-names>
</name>
<name>
<surname>Suvachitanont</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Molecular cloning of a new cDNA and expression of 3-hydroxy-3-methylglutaryl-CoA synthase gene from <italic>Hevea brasiliensis</italic>
</article-title>. <source>Planta</source> <volume>221</volume>, <fpage>502</fpage>&#x2013;<lpage>512</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00425-004-1463-7</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Suksup</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Imkaew</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Smitthipong</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Cream concentrated latex for foam rubber products</article-title>. <source>IOP. Conf. Ser. Mater. Sci. Eng.</source> <volume>272</volume>, <elocation-id>12025</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1088/1757-899X/272/1/012025</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Sakr</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>H. P.</given-names>
</name>
<etal/>
</person-group>. (<year>2010</year>). <article-title>The sucrose transporter <italic>HbSUT3</italic> plays an active role in sucrose loading to laticifer and rubber productivity in exploited trees of <italic>Hevea brasiliensis</italic> (para rubber tree)</article-title>. <source>Plant Cell Environ.</source> <volume>33</volume>, <fpage>1708</fpage>&#x2013;<lpage>1720</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-3040.2010.02175.x</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Fang</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>X.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>The rubber tree genome reveals new insights into rubber production and species adaptation</article-title>. <source>Nat. Plants</source> <volume>2</volume>, <fpage>16073</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nplants.2016.73</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tao</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>H. W.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>W. K.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>S. Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J. S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>The role of ethylene in plants under salinity stress</article-title>. <source>Front. Plant Sci.</source> <volume>6</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2015.01059</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thao</surname> <given-names>N. P.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>M. I. R.</given-names>
</name>
<name>
<surname>Thu</surname> <given-names>N. B. A.</given-names>
</name>
<name>
<surname>Hoang</surname> <given-names>X. L. T.</given-names>
</name>
<name>
<surname>Asgher</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>N. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Role of ethylene and its cross talk with other signaling molecules in plant responses to heavy metal stress</article-title>. <source>Plant Physiol.</source> <volume>169</volume>, <fpage>73</fpage>&#x2013;<lpage>84</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.15.00663</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tungngoen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kongsawadworakul</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Viboonjun</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Katsuhara</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Brunel</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Sakr</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2009</year>). <article-title>Involvement of HbPIP2;1 and HbTIP1;1 aquaporins in ethylene stimulation of latex yield through regulation of water exchanges between inner liber and latex cells in <italic>Hevea brasiliensis</italic>
</article-title>. <source>Plant Physiol.</source> <volume>151</volume>, <fpage>843</fpage>&#x2013;<lpage>856</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1104/pp.109.140228</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tungngoen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Viboonjun</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Kongsawadworakul</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Katsuhara</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Julien</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Sakr</surname> <given-names>S.</given-names>
</name>
<etal/>
</person-group>. (<year>2011</year>). <article-title>Hormonal treatment of the bark of rubber trees (<italic>Hevea brasiliensis</italic>) increases latex yield through latex dilution in relation with the differential expression of two aquaporin genes</article-title>. <source>J. Plant Physiol.</source> <volume>168</volume>, <fpage>253</fpage>&#x2013;<lpage>262</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jplph.2010.06.009</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>RSeQC: quality control of RNA-seq experiments</article-title>. <source>Bioinformatics</source> <volume>28</volume>, <fpage>2184</fpage>&#x2013;<lpage>2185</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/bts356</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Comprehensive proteomics analysis of laticifer latex reveals new insights into ethylene stimulation of natural rubber production</article-title>. <source>Sci. Rep.</source> <volume>5</volume>, <elocation-id>13778</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/srep13778</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Sinelnikov</surname> <given-names>I. V.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Wishart</surname> <given-names>D. S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>MetaboAnalyst 3.0&#x2014;making metabolomics more meaningful</article-title>. <source>Nucleic Acids Res.</source> <volume>43</volume>, <fpage>W251</fpage>&#x2013;<lpage>W257</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/nar/gkv380</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Dai</surname> <given-names>L. J.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>Y. X.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Melatonin enhances salt stress tolerance in rubber tree (<italic>Hevea brasiliensis</italic>) seedlings</article-title>. <source>Ind. Crops Prod.</source> <volume>145</volume>, <elocation-id>111990</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.indcrop.2019.111990</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeier</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>New insights into the regulation of plant immunity by amino acid metabolic pathways</article-title>. <source>Plant Cell Environ.</source> <volume>36</volume>, <fpage>2085</fpage>&#x2013;<lpage>2103</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/pce.12122</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Ou</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Z.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Intermediate-volatility aromatic hydrocarbons from the rubber products industry in China</article-title>. <source>Sci. Total. Environ.</source> <volume>898</volume>, <elocation-id>165583</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.scitotenv.2023.165583</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname> <given-names>J. H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Ethylene stimulation of latex production in <italic>Hevea brasiliensis</italic>
</article-title>. <source>Plant Signaling Behav.</source> <volume>4</volume>, <fpage>1072</fpage>&#x2013;<lpage>1074</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4161/psb.4.11.9738</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>