Genome-Wide Identification and Expression Profile Analysis of the Phenylalanine Ammonia-Lyase Gene Family in Hevea brasiliensis.

The majority of the world's natural rubber comes from the rubber tree (Hevea brasiliensis). As a key enzyme for synthesizing phenylpropanoid compounds, phenylalanine ammonia-lyase (PAL) has a critical role in plant satisfactory growth and environmental adaptation. To clarify the characteristics of rubber tree PAL family genes, a genome-wide characterization of rubber tree PALs was conducted in this study. Eight PAL genes (HbPAL1-HbPAL8), which spread over chromosomes 3, 7, 8, 10, 12, 13, 14, 16, and 18, were found to be present in the genome of H. brasiliensis. Phylogenetic analysis classified HbPALs into groups I and II, and the group I HbPALs (HbPAL1-HbPAL6) displayed similar conserved motif compositions and gene architectures. Tissue expression patterns of HbPALs quantified by quantitative real-time PCR (qPCR) proved that distinct HbPALs exhibited varying tissue expression patterns. The HbPAL promoters contained a plethora of cis-acting elements that responded to hormones and stress, and the qPCR analysis demonstrated that abiotic stressors like cold, drought, salt, and H2O2-induced oxidative stress, as well as hormones like salicylic acid, abscisic acid, ethylene, and methyl jasmonate, controlled the expression of HbPALs. The majority of HbPALs were also regulated by powdery mildew, anthracnose, and Corynespora leaf fall disease infection. In addition, HbPAL1, HbPAL4, and HbPAL7 were significantly up-regulated in the bark of tapping panel dryness rubber trees relative to that of healthy trees. Our results provide a thorough comprehension of the characteristics of HbPAL genes and set the groundwork for further investigation of the biological functions of HbPALs in rubber trees.

Construction and analysis of the tapping panel dryness-related lncRNA/circRNA-miRNA-mRNA ceRNA network in latex of Hevea brasiliensis.

Tapping panel dryness (TPD) results in a severe reduction in latex yield in Hevea brasiliensis. However, the molecular regulatory mechanisms of TPD occurrence are still largely unclear. In this study, whole-transcriptome sequencing was carried out on latex from TPD and healthy trees. In total, 7078 long noncoding RNAs (lncRNAs), 3077 circular RNAs (circRNAs), 4956 miRNAs, and 25041 mRNAs were identified in latex, among which 435 lncRNAs, 68 circRNAs, 320 miRNAs, and 1574 mRNAs were differentially expressed in the latex of TPD trees. GO and KEGG analyses indicated that plant hormone signal transduction, MAPK signaling pathway, and ubiquitin-mediated proteolysis were the key pathways associated with TPD onset. Phytohormone profiling revealed significant changes in the contents of 28 hormonal compounds, among which ACC, ABA, IAA, GA, and JA contents were increased, while SA content was reduced in TPD latex, suggesting that hormone homeostasis is disrupted in TPD trees. Furthermore, we constructed a TPD-related competitive endogenous RNA (ceRNA) regulatory network of lncRNA/circRNA-miRNA-mRNA with 561 edges and 434 nodes (188 lncRNAs, 5 circRNAs, 191 miRNAs, and 50 mRNAs) and identified two hub lncRNAs (MSTRG.11908.1 and MSTRG.8791.1) and four hub miRNAs (hbr-miR156, miR156-x, miRf10477-y, and novel-m0452-3p). Notably, the lncRNA-miR156/157-SPL module containing three hubs probably plays a crucial role in TPD onset. The expression of network hubs and the lncRNA-miR156/157-SPL module were further validated by qRT-PCR. Our results reveal the TPD-associated ceRNA regulatory network of lncRNA/circRNA-miRNA-mRNA in latex and lay a foundation for further investigation of molecular regulatory mechanisms for TPD onset in H. brasiliensis.

Downregulation of HbFPS1 affects rubber biosynthesis of Hevea brasiliensis suffering from tapping panel dryness.

Tapping panel dryness (TPD) is a century-old problem that has plagued the natural rubber production of Hevea brasiliensis. TPD may result from self-protective mechanisms of H. brasiliensis in response to stresses such as excessive hormone stimulation and mechanical wounding (bark tapping). It has been hypothesized that TPD impairs rubber biosynthesis; however, the underlying mechanisms remain poorly understood. In the present study, we firstly verified that TPD-affected rubber trees exhibited lower rubber biosynthesis activity and greater rubber molecular weight compared to healthy rubber trees. We then demonstrated that HbFPS1, a key gene of rubber biosynthesis, and its expression products were downregulated in the latex of TPD-affected rubber trees, as revealed by transcriptome sequencing and iTRAQ-based proteome analysis. We further discovered that the farnesyl diphosphate synthase HbFPS1 could be recruited to small rubber particles by HbSRPP1 through protein-protein interactions to catalyze farnesyl diphosphate (FPP) synthesis and facilitate rubber biosynthesis initiation. FPP content in the latex of TPD-affected rubber trees was significantly decreased with the downregulation of HbFPS1, ultimately resulting in abnormal development of rubber particles, decreased rubber biosynthesis activity, and increased rubber molecular weight. Upstream regulator assays indicated that a novel regulator, MYB2-like, may be an important regulator of downregulation of HbFPS1 in the latex of TPD-affected rubber trees. Our findings not only provide new directions for studying the molecular events involved in rubber biosynthesis and TPD syndrome and contribute to rubber management strategies, but also broaden our knowledge of plant isoprenoid metabolism and its regulatory networks.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis).

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Involvement of HbMC1-mediated cell death in tapping panel dryness of rubber tree (Hevea brasiliensis).

Tapping panel dryness (TPD) causes a significant reduction in the latex yield of rubber tree (Hevea brasiliensis Muell. Arg.). It is reported that TPD is a typical programmed cell death (PCD) process. Although PCD plays a vital role in TPD occurrence, there is a lack of detailed and systematic study. Metacaspases are key regulators of diverse PCD in plants. Based on our previous result that HbMC1 was associated with TPD, we further elucidate the roles of HbMC1 on rubber tree TPD in this study. HbMC1 was up-regulated by TPD-inducing factors including wounding, ethephon and H2O2. Moreover, the expression level of HbMC1 was increased along with TPD severity in rubber tree, suggesting a positive correlation between HbMC1 expression and TPD severity. To investigate its biological function, HbMC1 was overexpressed in yeast (Saccharomyces cerevisiae) and tobacco (Nicotiana benthamiana). Transgenic yeast and tobacco overexpressing HbMC1 showed growth retardation compared with controls under H2O2-induced oxidative stress. In addition, overexpression of HbMC1 in yeast and tobacco reduced cell survival after high-concentration H2O2 treatment and resulted in enhanced H2O2-induced leaf cell death, respectively. A total of 11 proteins, rbcL, TM9SF2-like, COX3, ATP9, DRP, HbREF/Hevb1, MSSP2-like, SRC2, GATL8, CIPK14-like and STK, were identified and confirmed to interact with HbMC1 by yeast two-hybrid screening and co-transformation in yeast. The 11 proteins mentioned above are associated with many biological processes, including rubber biosynthesis, stress response, autophagy, carbohydrate metabolism, signal transduction, etc. Taken together, our results suggest that HbMC1-mediated PCD plays an important role in rubber tree TPD, and the identified HbMC1-interacting proteins provide valuable information for further understanding the molecular mechanism of HbMC1-mediated TPD in rubber tree.

Over-expression of 3-hydroxy-3- methylglutaryl-coenzyme A reductase 1 (hmgr1) gene under super-promoter for enhanced latex biosynthesis in rubber tree (Hevea brasiliensis Muell. Arg.).

Natural rubber (cis-1, 4-polyisoprene) is being produced from bark laticifer cells of Hevea brasiliensis and the popular high latex yielding Indian rubber clones are easily prone to onset of tapping panel dryness syndrome (TPD) which is considered as a physiological syndrome affecting latex production either partially or completely. This report describes an efficient protocol for development of transgenic rubber plants by over-expression of 3-hydroxy 3-methylglutaryl Co-enzyme A reductase 1 (hmgr1) gene which is considered as rate limiting factor for latex biosynthesis via Agrobacterium-mediated transformation. The pBIB plasmid vector containing hmgr1 gene cloned under the control of a super-promoter was used for genetic transformation using embryogenic callus. Putatively transgenic cell lines were obtained on selection medium and produced plantlets with 44% regeneration efficiency. Transgene integration was confirmed by PCR amplification of 1.8 kb hmgr1 and 0.6 kb hpt genes from all putatively transformed callus lines as well as transgenic plants. Southern blot analysis showed the stable integration and presence of transgene in the transgenic plants. Over expression of hmgr1 transgene was determined by Northern blot hybridization, semi-quantitative PCR and real-time PCR (qRT-PCR) analysis. Accumulation of hmgr1 mRNA transcripts was more abundant in transgenic plants than control. Increased level of photosynthetic pigments, protein contents and HMGR enzyme activity was also noticed in transgenic plants over control. Interestingly, the latex yield was significantly enhanced in all transgenic plants compared to the control. The qRT-PCR results exhibit that the hmgr1 mRNA transcript levels was 160-fold more abundance in transgenic plants over untransformed control. These results altogether suggest that there is a positive correlation between latex yield and accumulation of mRNA transcripts level as well as HMGR enzyme activity in transgenic rubber plants. It is presumed that there is a possibility for enhanced level of latex biosynthesis in transgenic plants as the level of mRNA transcripts and HMGR enzyme activity is directly correlated with latex yield in rubber tree. Further, the present results clearly suggest that the quantification of HMGR enzyme activity in young seedlings will be highly beneficial for early selection of high latex yielding plants in rubber breeding programs.

Genome-wide identification and expression analysis of the metacaspase gene family in Hevea brasiliensis.

Metacaspases, a family of cysteine proteases, have been suggested to play important roles in programmed cell death (PCD) during plant development and stress responses. To date, no systematic characterization of this gene family has been reported in rubber tree (Hevea brasiliensis). In the present study, nine metacaspase genes, designated as HbMC1 to HbMC9, were identified from whole-genome sequence of rubber tree. Multiple sequence alignment and phylogenetic analyses suggested that these genes were divided into two types: type I (HbMC1-HBMC7) and type II (HbMC8 and HbMC9). Gene structure analysis demonstrated that type I and type II HbMCs separately contained four and two introns, indicating the conserved exon-intron organization of HbMCs. Quantitative real-time PCR analysis revealed that HbMCs showed distinct expression patterns in different tissues, suggesting the functional diversity of HbMCs in various tissues during development. Most of the HbMCs were regulated by drought, cold, and salt stress, implying their possible functions in regulating abiotic stress-induced cell death. Of the nine HbMCs, HbMC1, HbMC2, HbMC5, and HbMC8 displayed a significantly higher relative transcript accumulation in barks of tapping panel dryness (TPD) trees compared with healthy trees. In addition, the four genes were up-regulated by ethephon (ET) and methyl jasmonate (MeJA), indicating their potential involvement in TPD resulting from ET- or JA-induced PCD. In summary, this work provides valuable information for further functional characterization of HbMC genes in rubber tree.

Molecular cloning, expression profiles, and characterization of a novel polyphenol oxidase (PPO) gene in Hevea brasiliensis.

The polyphenol oxidase (PPO) is involved in undesirable browning in many plant foods. Although the PPOs have been studied by several researchers, the isolation and expression profiles of PPO gene were not reported in rubber tree. In this study, a new PPO gene, HbPPO, was isolated from Hevea brasiliensis. The sequence alignment showed that HbPPO indicated high identities to plant PPOs and belonged to dicot branch. The cis-acting regulatory elements related to stress/hormone responses were predicted in the promoter region of HbPPO. Real-time RT-PCR analyses showed that HbPPO expression varied widely depending on different tissues and developmental stages of leaves. Besides being associated with tapping panel dryness, the HbPPO transcripts were regulated by ethrel, wounding, H2O2, and methyl jasmonate treatments. Moreover, the correlation between latex coagulation rate and PPO activity was further confirmed in this study. Our results lay the foundation for further analyzing the function of HbPPO in rubber tree.

Molecular characterization and expression analysis of the small GTPase ROP members expressed in laticifers of the rubber tree (Hevea brasiliensis).

ROP (Rho of plants) proteins are plant-specific Rho-type small GTPases which play important roles in cellular processes and stress responses. This study explores the characteristics and possible functions of ROPs that are expressed primarily in laticifers of the rubber tree (Hevea brasiliensis). The work serves as a preliminary step to determining their involvement in latex flow and regeneration, laticifers formation and tapping panel dryness (TPD, a physiological disorder in rubber trees that result in the stoppage of latex flow). In this connection, we (i) identified five HbROPs (HbROP1-HbROP5) by searching latex transcripts database and the genome databases, (ii) characterized molecular and phylogenic aspects of the HbROPs and examined the cis-regulatory elements in their promoter regions; (iii) analyzed by Real-time Quantitative PCR (QPCR) the tissue specificity of the HbROPs and their expression patterns in response to tapping, bark wounding and growth regulator treatments. All five HbROP genes were strongly expressed in the latex, with HbROP1, 3, 4 and 5 showing the highest expression among the six Hevea tissues examined, viz. latex, bud, mature leaf, bark, male flower and seed. When tapping was initiated on previously untapped trees, HbROP3 transcription was substantially down-regulated whereas HbROP5 expression was markedly up-regulated. Transcripts of HbROP3 rose gradually with the development of TPD. Except for the cytokinin 6-benzyl aminopurine that induced a rise in HbROP5 transcripts by more than 2-fold, the other growth regulators tested had little effect on HbROPs expression. The roles of HbROPs in rubber tree are discussed in relation to the diverse functions of ROP homologs reported in other plant species.