Screening of valid reference genes for real-time RT-PCR data normalization in Hevea brasiliensis and expression validation of a sucrose transporter gene HbSUT3

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.

A novel rubber tree PR-10 protein involved in host-defense response against the white root rot fungus Rigidoporus microporus

BACKGROUND

White root rot disease in rubber trees, caused by the pathogenic fungi Rigidoporus microporus, is currently considered a major problem in rubber tree plantations worldwide. Only a few reports have mentioned the response of rubber trees occurring at the non-infection sites, which is crucial for the disease understanding and protecting the yield losses.

RESULTS

Through a comparative proteomic study using the two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) technique, the present study reveals some distal-responsive proteins in rubber tree leaves during the plant-fungal pathogen interaction. From a total of 12 selected differentially expressed protein spots, several defense-related proteins such as molecular chaperones and ROS-detoxifying enzymes were identified. The expression of 6 candidate proteins was investigated at the transcript level by Reverse Transcription Quantitative PCR (RT-qPCR). In silico, a highly-expressed uncharacterized protein LOC110648447 found in rubber trees was predicted to be a protein in the pathogenesis-related protein 10 (PR-10) class. In silico promoter analysis and structural-related characterization of this novel PR-10 protein suggest that it plays a potential role in defending rubber trees against R. microporus infection. The promoter contains WRKY-, MYB-, and other defense-related cis-acting elements. The structural model of the novel PR-10 protein predicted by I-TASSER showed a topology of the Bet v 1 protein family, including a conserved active site and a ligand-binding hydrophobic cavity.

CONCLUSIONS

A novel protein in the PR-10 group increased sharply in rubber tree leaves during interaction with the white root rot pathogen, potentially contributing to host defense. The results of this study provide information useful for white root rot disease management of rubber trees in the future.

Functional characterization of a bark-specific monoterpene synthase potentially involved in wounding- and methyl jasmonate-induced linalool emission in rubber (Hevea brasiliensis)

Rubber (Hevea brasiliensis) is a latex-producing plant that often encounters mechanical wounding, as well as pathogen and pest attacks through wound sites during and after tapping. Terpenoids play an important role in the ecological interactions of many plant species, and their diversity is mainly generated by enzymes known as terpene synthases (TPS). In this study, one cDNA sequence encoding a putative terpene synthase, HbTPS20, was obtained from the bark tissues of H. brasiliensis. The encoded protein contains 610 amino acids with a putative N-terminal plastid transit peptide of approximately 70 residues. It belongs to the TPS-b subfamily. Further phylogenetic analysis showed that HbTPS20 formed a separate branch that diverged from the progenitor of all other potentially functional terpene synthases of the rubber TPS-b subfamily. The truncated HbTPS20 without the signal peptide coding sequence was successfully expressed in E. coli and in vitro enzymatic assays with geranyl diphosphate (GPP) or neryl diphosphate (NPP) as a substrate defined HbTPS20 as an active linalool synthase (HbLIS) with the ability to produce linalool as the principal product. RT-qPCR analysis showed that the highest transcript levels of HbTPS20 were found in barks, and this gene was expressed at 2.26- and 250-fold greater levels in the bark tissues of wounded and MeJA-treated plants, respectively, than in those of the control plants. This indicates that this gene may be involved in the induced stress responses of rubber.

The rubber tree RALF peptide hormone and its receptor protein kinase FER implicates in rubber production

RAPID ALKALINIZATION FACTORs (RALFs), which are secreted peptides serving as extracellular signals transduced to the inside of the cell, interact with the receptor-like kinase FERONIA (FER) and participates in various biological pathways. Here, we identified 23 RALF and 2 FER genes in Hevea brasiliensis (para rubber tree), and characterized their expression patterns in different tissues, across the process of leaf development, and in response to the rubber yield-stimulating treatments of tapping and ethylene. Four Hevea latex (the cytoplasm of rubber-producing laticifers)-abundant RALF isoforms, HbRALF19, HbRALF3, HbRALF22, and HbRALF16 were listed with descending expression levels. Of the four HbRALFs, expressions of HbRALF3 were markedly regulated in an opposite way by the treatments of tapping (depression) and ethylene (stimulation). All of the four latex-abundant RALFs specifically interacted with the extracellular domain of HbFER1. Transgenic Arabidopsis plants overexpressing these HbRALFs displayed phenotypes similar to those reported for AtRALFs, such as shorter roots, smaller plant architecture, and delayed flowering. The application of HbRALF3 and HbRALF19 recombinant proteins significantly reduced the pH of Hevea latex, an important factor regulating latex metabolism. An in vitro rubber biosynthesis assay in a mixture of latex cytosol (C-serum) revealed a positive role of HbFER1 in rubber biosynthesis. Taken together, these data provide evidence for the participation of the HbRALF-FER module in rubber production.

Trehalose 6-Phosphate/SnRK1 Signaling Participates in Harvesting-Stimulated Rubber Production in the Hevea Tree

Trehalose 6-phosphate (T6P), the intermediate of trehalose biosynthesis and a signaling molecule, affects crop yield via targeting sucrose allocation and utilization. As there have been no reports of T6P signaling affecting secondary metabolism in a crop plant, the rubber tree Hevea brasiliensis serves as an ideal model in this regard. Sucrose metabolism critically influences the productivity of natural rubber, a secondary metabolite of industrial importance. Here, we report on the characterization of the T6P synthase (TPS) gene family and the T6P/SNF1-related protein kinase1 (T6P/SnRK1) signaling components in Hevea laticifers under tapping (rubber harvesting), an agronomic manipulation that itself stimulates rubber production. A total of fourteen TPS genes were identified, among which a class II TPS gene, HbTPS5, seemed to have evolved with a function specialized in laticifers. T6P and trehalose increased when the trees were tapped, this being consistent with the observed enhanced activities of TPS and T6P phosphatase (TPP) and expression of an active TPS-encoding gene, HbTPS1. On the other hand, SnRK1 activities decreased, suggesting the inhibition of elevated T6P on SnRK1. Expression profiles of the SnRK1 marker genes coincided with elevated T6P and depressed SnRK1. Interestingly, HbTPS5 expression decreased significantly with the onset of tapping, suggesting a regulatory function in the T6P pathway associated with latex production in laticifers. In brief, transcriptional, enzymatic, and metabolic evidence supports the participation of T6P/SnRK1 signaling in rubber formation, thus providing a possible avenue to increasing the yield of a valuable secondary metabolite by targeting T6P in specific cells.

Identification and Functional Evaluation of Three Polyubiquitin Promoters from Hevea brasiliensis

Hevea brasiliensis is an economically important tree species that provides the only commercial source of natural rubber. The replacement of the CaMV35S promoter by endogenous polyubiquitin promoters may be a viable way to improve the genetic transformation of this species. However, no endogenous polyubiquitin promoters in Hevea have been reported yet. Here, we identified three Hevea polyubiquitin genes HbUBI10.1, HbUBI10.2 and HbUBI10.3, which encode ubiquitin monomers having nearly identical amino acid sequences to that of AtUBQ10. The genomic fragments upstream of these HbUBI genes, including the signature leading introns, were amplified as putative HbUBI promoters. In silico analysis showed that a number of cis-acting elements which are conserved within strong constitutive polyubiquitin promoters were presented in these HbUBI promoters. Transcriptomic data revealed that HbUBI10.1 and HbUBI10.2 had a constitutive expression in Hevea plants. Semi-quantitative RT-PCR showed that these three HbUBI genes were expressed higher than the GUS gene driven by CaMV35S in transgenic Hevea leaves. All three HbUBI promoters exhibited the capability to direct GFP expression in both transient and stable transformation assays, although they produced lower protoplast transformation efficiencies than the CaMV35S promoter. These HbUBI promoters will expand the availability of promoters for driving the transgene expression in Hevea genetic engineering.

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.

The Calcium Sensor Calcineurin B-Like Proteins -Calcineurin B-Like Interacting Protein Kinases Is Involved in Leaf Development and Stress Responses Related to Latex Flow in Hevea brasiliensis

Latex flow in Hevea brasiliensis (the Para rubber tree), the sole commercial source of natural rubber (cis-1,4-polyisoprene, NR), renders it uniquely suited for the study of plant stress responses. Calcineurin B-like interacting protein kinases (CIPK) serving as calcium-sensor protein kinases react with calcineurin B-like proteins (CBL) to play crucial roles in hormone signaling transduction and response to abiotic stress in plant developmental processes. However, little is known about their functions in Hevea. In this study, a total of twelve CBL (HbCBL) and thirty CIPK (HbCIPK) genes were identified from the Hevea genome. Structure and phylogenetic analysis assigned these CIPKs to five groups and CBLs to four groups, and mapped onto fourteen of the eighteen Hevea chromosomes. RNA-seq and qPCR analysis showed that the expressions of HbCBL and HbCIPK genes varied in the seven Hevea tissues examined, i.e., latex (cytoplasm of rubber-producing laticifers), bark, leaf, root, seed, female flower, and male flower. The expressions of two HbCBL and sixteen HbCIPK genes showed upward trends during leaf development. Following ethylene yield stimulation and the latex tapping treatment, both practices invoking stress, the expression levels of most latex-expressed genes were significantly altered. Yeast two-hybrid test revealed interactions for multiple combinations of HbCBLs and HbCIPKs with substantial gene expression in latex or other Hevea tissues. However, all the HbCBL-HbCIPK complexes examined did not recruit HbSOS1 or AtSOS1 to form functional salt tolerance SOS pathway in yeast cells. Taken together, the results suggested a role of the Hevea CBL-CIPK network as a point of convergence for several different signaling pathways in growth, development, and stress responses in relation to latex production.

Transcript and Protein Profiling Provides Insights Into the Molecular Mechanisms of Harvesting-Induced Latex Production in Rubber Tree

Natural rubber, an important industrial raw material with wide applications, is harvested in the form of latex (cytoplasm of rubber-producing laticifers) from Hevea brasiliensis (para rubber tree) by the way of tapping. Conspicuous stimulation on latex production is observed for the first few tappings conducted on virgin (untapped before) or resting (tapped before but no tapping for a period) rubber trees. To understand the underlying mechanisms, an integrative analysis of the latex transcriptome and proteome was conducted on virgin or resting Hevea trees for the first five tappings. A total of 505 non-redundant differentially expressed (DE) transcript-derived fragments (TDFs) were identified by silver-staining cDNA-AFLP, with 217 exhibiting patterns of upregulated, 180 downregulated and 108 irregularly-regulated. Meanwhile, 117 two dimensional gel electrophoresis DE-protein spots were isolated and subjected to mass spectrometry analysis, with 89 and 57 being successfully identified by MALDI-TOF and MALDI-TOF/TOF, respectively. About 72.5% DE-TDFs and 76.1% DE-proteins were functionally annotated and categorized. Noteworthily, most of the DE-TDFs implicated in sugar transport and metabolism as well as rubber biosynthesis were upregulated by the tapping treatment. The importance of sugar metabolism in harvesting-induced latex production was reinforced by the identification of abundant relevant DE-protein spots. About 83.8% of the randomly selected DE-TDFs were validated for expression patterns by semi-quantitative RT-PCR, and an 89.7% consistency for the 29 latex regeneration-related DE-TDFs examined by quantitative RT-PCR analysis. In brief, our results reveal extensive physiological and molecular changes in Hevea laticifers incurred by the tapping treatment, and the vast number of DE genes and proteins identified here contribute to unraveling the gene regulatory network of tapping-stimulated latex production.

Evaluation of reference genes and characterization of the MYBs in xylem radial change of Chinese fir stem

The radial change (RC) of tree stem is the process of heartwood formation involved in complex molecular mechanism. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), an evergreen species, is an important fast-growing timber tree in southern China. In this study, the top four stable genes (IDH, UBC2, RCA and H2B) were selected in RC tissues of 15 years old Chinese fir stem (RC15) and the genes (H2B, 18S, TIP41 and GAPDH) were selected in RC tissues of 30 years old Chinese fir stem (RC30). The stability of the reference genes is higher in RC30 than in RC15. Sixty-one MYB transcripts were obtained on the PacBio Sequel platform from woody tissues of one 30 years old Chinese fir stem. Based on the number of MYB DNA-binding domain and phylogenetic relationships, the ClMYB transcripts contained 21 transcripts of MYB-related proteins (1R-MYB), 39 transcripts of R2R3-MYB proteins (2R-MYB), one transcript of R1R2R3-MYB protein (3R-MYB) belonged to 18 function-annotated clades and two function-unknown clades. In RC woody tissues of 30 years old Chinese fir stem, ClMYB22 was the transcript with the greatest fold change detected by both RNA-seq and qRT-PCR. Reference genes selected in this study will be helpful for further verification of transcript abundance patterns during the heartwood formation of Chinese fir.

RNA Sequencing Reveals Phenylpropanoid Biosynthesis Genes and Transcription Factors for Hevea brasiliensis Reaction Wood Formation

Given the importance of wood in many industrial applications, much research has focused on wood formation, especially lignin biosynthesis. However, the mechanisms governing the regulation of lignin biosynthesis in the rubber tree (Hevea brasiliensis) remain to be elucidated. Here, we gained insight into the mechanisms of rubber tree lignin biosynthesis using reaction wood (wood with abnormal tissue structure induced by gravity or artificial mechanical treatment) as an experimental model. We performed transcriptome analysis of rubber tree mature xylem from tension wood (TW), opposite wood (OW), and normal wood (NW) using RNA sequencing (RNA-seq). A total of 214, 1,280, and 32 differentially expressed genes (DEGs) were identified in TW vs. NW, OW vs. NW, and TW vs. OW, respectively. GO and KEGG enrichment analysis of DEGs from different comparison groups showed that zeatin biosynthesis, plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interaction pathways may play important roles in reaction wood formation. Sixteen transcripts involved in phenylpropanoid biosynthesis and 129 transcripts encoding transcription factors (TFs) were used to construct a TF-gene regulatory network for rubber tree lignin biosynthesis. Among them, MYB, C2H2, and NAC TFs could regulate all the DEGs involved in phenylpropanoid biosynthesis. Overall, this study identified candidate genes and TFs likely involved in phenylpropanoid biosynthesis and provides novel insights into the mechanisms regulating rubber tree lignin biosynthesis.

Ethylene response factors regulate expression of HbSUT3, the sucrose influx carrier in laticifers of Hevea brasiliensis

Natural rubber is an important industrial raw material and is commercially produced by rubber trees (Hevea brasiliensis). The sucrose transporter HbSUT3 plays an essential role in rubber production. Its expression in latex (cytoplasm of rubber-producing laticifers) is induced by bark treatment with Ethrel, an ethylene releaser, and the inducing effect correlates well with Ethrel-stimulated rubber yield increase. However, the mechanisms of ethylene induction on HbSUT3 expression are not known. Here, five Ethylene Response Factor (ERF) genes were identified from the cDNA library of Hevea latex by yeast one-hybrid screening with the promoter of HbSUT3 gene as bait. As revealed in a tobacco (Nicotiana tabacum) protoplast transient expression system, these HbERFs were mainly localized in the nucleus and four of them exhibited apparent transactivation activity. Of the five HbERF genes, HbERF-IXc4 was the most frequently screened in yeast one-hybrid, accounting for 65% of the ERF clones obtained. Moreover, among the five HbERFs, HbERF-IXc4 showed the strongest transactivation capacity when expressed in tobacco protoplast, the highest transcript abundance in latex and a close expressional correlation with its target gene, HbSUT3, in response to the Ethrel treatment. Taken together, our results indicate that ERFs, especially HbERF-IXc4, are critically involved in the activation of HbSUT3 expression in latex after Ethrel treatment on Hevea bark, and thus the stimulated latex yield.

Ultrafast and Real-Time Nanoplasmonic On-Chip Polymerase Chain Reaction for Rapid and Quantitative Molecular Diagnostics

Advent and fast spread of pandemic diseases draw worldwide attention to rapid, prompt, and accurate molecular diagnostics with technical development of ultrafast polymerase chain reaction (PCR). Microfluidic on-chip PCR platforms provide highly efficient and small-volume bioassay for point-of-care diagnostic applications. Here we report ultrafast, real-time, and on-chip nanoplasmonic PCR for rapid and quantitative molecular diagnostics at point-of-care level. The plasmofluidic PCR chip comprises glass nanopillar arrays with Au nanoislands and gas-permeable microfluidic channels, which contain reaction microchamber arrays, a precharged vacuum cell, and a vapor barrier. The on-chip configuration allows both spontaneous sample loading and microbubble-free PCR reaction during which the plasmonic nanopillar arrays result in ultrafast photothermal cycling. After rapid sample loading less than 3 min, two-step PCR results for 40 cycles show rapid amplification in 264 s for lambda-DNA, and 306 s for plasmids expressing SARS-CoV-2 envelope protein. In addition, the in situ cyclic real-time quantification of amplicons clearly demonstrates the amplification efficiencies of more than 91%. This PCR platform can provide rapid point-of-care molecular diagnostics in helping slow the fast-spreading pandemic.

A Heterogeneity Radiomic Nomogram for Preoperative Differentiation of Primary Gastric Lymphoma From Borrmann Type IV Gastric Cancer

OBJECTIVE

This study aimed to preoperatively differentiate primary gastric lymphoma from Borrmann type IV gastric cancer by heterogeneity nomogram based on routine contrast-enhanced computed tomographic images.

METHODS

We enrolled 189 patients from 2 hospitals (90 in the training cohort and 99 in the validation cohort). Subjective findings, including high-enhanced mucosal sign, high-enhanced serosa sign, nodular or an irregular outer layer of the gastric wall, and perigastric fat infiltration, were assessed to construct a subjective finding model. A deep learning model was developed to segment tumor areas, from which 1680 three-dimensional heterogeneity radiomic parameters, including first-order entropy, second-order entropy, and texture complexity, were extracted to build a heterogeneity signature by least absolute shrinkage and selection operator logistic regression. A nomogram that integrates heterogeneity signature and subjective findings was developed by multivariate logistic regression. The diagnostic performance of the nomogram was assessed by discrimination and clinical usefulness.

RESULTS

High-enhanced serosa sign and nodular or an irregular outer layer of the gastric wall were identified as independent predictors for building the subjective finding model. High-enhanced serosa sign and heterogeneity signature were significant predictors for differentiating the 2 groups (all, P < 0.05). The area under the curve with heterogeneity nomogram was 0.932 (95% confidence interval, 0.863-0.973) in the validation cohort. Decision curve analysis and stratified analysis confirmed the clinical utility of the heterogeneity nomogram.

CONCLUSIONS

The proposed heterogeneity radiomic nomogram on contrast-enhanced computed tomographic images may help differentiate primary gastric lymphoma from Borrmann type IV gastric cancer preoperatively.

HbWRKY40 plays an important role in the regulation of pathogen resistance in Hevea brasiliensis

KEY MESSAGE

Overexpression of HbWRKY40 induces ROS burst in tobacco and increases disease resistance in Arabidopsis; RNA-seq and ChIP assays revealed the regulatory network of HbWRKY40 in plant defense. WRKY, a family of plant transcription factors, are involved in the regulation of numerous biological processes. In rubber tree Hevea brasiliensis, the roles of WRKYs remain poorly understood. In the present study, a total of 111 genes encoding putative HbWRKY proteins were identified in the H. brasiliensis genome. Among these genes, HbWRKY40 transcripts were significantly induced by Colletotrichum gloeosporioides and salicylic acid. To assess its roles in plant defense, HbWRKY40 was over-expressed in Nicotiana benthamiana and Arabidopsis thaliana. The results showed that HbWRKY40 significantly induced reactive oxygen species burst in N. benthamiana and increased resistance of Arabidopsis against Botrytis cinerea. Transient expression in mesophyll cell protoplasts of H. brasiliensis showed that HbWRKY40 localizes at nuclei. In addition, transcripts of 145 genes were significantly up-regulated and 6 genes were down-regulated in the protoplasts over-expressing HbWRKY40 based on the RNA-seq analysis. Among these potential downstream targets, 12 genes contain potential WRKY-binding sites at the promoter regions. Further analysis through chromatin immunoprecipitation revealed that 10 of these 12 genes were the downstream targets of HbWRKY40. Taken together, our findings indicate that HbWRKY40 plays an important role in the disease resistance by regulating defense-associated genes in H. brasiliensis.

Identification and evaluation of reliable reference genes for quantitative real-time PCR analysis in tea plants under differential biotic stresses

The selection of reliable reference genes (RGs) for normalization under given experimental conditions is necessary to develop an accurate qRT-PCR assay. To the best of our knowledge, only a small number of RGs have been rigorously identified and used in tea plants (Camellia sinensis (L.) O. Kuntze) under abiotic stresses, but no critical RG identification has been performed for tea plants under any biotic stresses till now. In the present study, we measured the mRNA transcriptional levels of ten candidate RGs under five experimental conditions; these genes have been identified as stable RGs in tea plants. By using the ΔCt method, geNorm, NormFinder and BestKeeper, CLATHRIN1 and UBC1, TUA1 and SAND1, or SAND1 and UBC1 were identified as the best combination for normalizing diurnal gene expression in leaves, stems and roots individually; CLATHRIN1 and GAPDH1 were identified as the best combination for jasmonic acid treatment; ACTIN1 and UBC1 were identified as the best combination for Toxoptera aurantii-infested leaves; UBC1 and GAPDH1 were identified as the best combination for Empoasca onukii-infested leaves; and SAND1 and TBP1 were identified as the best combination for Ectropis obliqua regurgitant-treated leaves. Furthermore, our results suggest that if the processing time of the treatment was long, the best RGs for normalization should be recommended according to the stability of the proposed RGs in different time intervals when intragroup differences were compared, which would strongly increase the accuracy and sensitivity of target gene expression in tea plants under biotic stresses. However, when the differences of intergroup were compared, the RGs for normalization should keep consistent across different time points. The results of this study provide a technical guidance for further study of the molecular mechanisms of tea plants under different biotic stresses.

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.

Characterization of a vacuolar sucrose transporter, HbSUT5, from Hevea brasiliensis: involvement in latex production through regulation of intracellular sucrose transport in the bark and laticifers

BACKGROUND

Sucrose (Suc), as the precursor molecule for rubber biosynthesis in Hevea brasiliensis, is transported via phloem-mediated long-distance transport from leaves to laticifers in trunk bark, where latex (cytoplasm of laticifers) is tapped for rubber. In our previous report, six Suc transporter (SUT) genes have been cloned in Hevea tree, among which HbSUT3 is verified to play an active role in Suc loading to the laticifers. In this study, another latex-abundant SUT isoform, HbSUT5, with expressions only inferior to HbSUT3 was characterized especially for its roles in latex production.

RESULTS

Both phylogenetic analysis and subcellular localization identify HbSUT5 as a tonoplast-localized SUT protein under the SUT4-clade (=type III). Suc uptake assay in baker's yeast reveals HbSUT5 to be a typical Suc-H⁺ symporter, but its high affinity for Suc (Km = 2.03 mM at pH 5.5) and the similar efficiency in transporting both Suc and maltose making it a peculiar SUT under the SUT4-clade. At the transcript level, HbSUT5 is abundantly and preferentially expressed in Hevea barks. The transcripts of HbSUT5 are conspicuously decreased both in Hevea latex and bark by two yield-stimulating treatments of tapping and ethephon, the patterns of which are contrary to HbSUT3. Under the ethephon treatment, the Suc level in latex cytosol decreases significantly, but that in latex lutoids (polydispersed vacuoles) changes little, suggesting a role of the decreased HbSUT5 expression in Suc compartmentalization in the lutoids and thus enhancing the Suc sink strength in laticifers.

CONCLUSIONS

Our findings provide insights into the roles of a vacuolar sucrose transporter, HbSUT5, in Suc exchange between lutoids and cytosol in rubber-producing laticifers.

Genome-Wide Identification and Characterization of Heat-Shock Transcription Factors in Rubber Tree

Heat-shock transcription factors (Hsfs) play a pivotal role in the response of plants to various stresses. The present study aimed to characterize the Hsf genes in the rubber tree, a primary global source of natural rubber. In this study, 30 Hsf genes were identified in the rubber tree using genome-wide analysis. They possessed a structurally conserved DNA-binding domain and an oligomerization domain. On the basis of the length of the insert region between HR-A and HR-B in the oligomerization domain, the 30 members were clustered into three classes, Classes A (18), B (10), and C (2). Members within the same class shared highly conserved gene structures and protein motifs. The background expression levels of 11 genes in cold-tolerant rubber-tree clone 93-14 were significantly higher than those in cold-sensitive rubber-tree clone Reken501, while four genes exhibited inverse expression patterns. Upon cold stress, 20 genes were significantly upregulated in 93-114. Of the upregulated genes, HbHsfA2b, HbHsfA3a, and HbHsfA7a were also significantly upregulated in three other cold-tolerant rubber-tree clones at one or more time intervals upon cold stress. Their nuclear localization was verified, and the protein–protein interaction network was predicted. This study provides a basis for dissecting Hsf function in the enhanced cold tolerance of the rubber tree.

Identification and evaluation of reference genes for quantitative real-time PCR analysis in Polygonum cuspidatum based on transcriptome data

BACKGROUND

Polygonum cuspidatum of the Polygonaceae family is a traditional medicinal plant with many bioactive compounds that play important roles in human health and stress responses. Research has attempted to identify biosynthesis genes and metabolic pathways in this species, and quantitative real-time PCR (RT-qPCR) has commonly been used to detect gene expression because of its speed, sensitivity, and specificity. However, no P. cuspidatum reference genes have been identified, which hinders gene expression studies. Here, we aimed to identify suitable reference genes for accurate and reliable normalization of P. cuspidatum RT-qPCR data.

RESULTS

Twelve candidate reference genes, including nine common (ACT, TUA, TUB, GAPDH, EF-1γ, UBQ, UBC, 60SrRNA, and eIF6A) and three novel (SKD1, YLS8, and NDUFA13), were analyzed in different tissues (root, stem, and leaf) without treatment and in leaves under abiotic stresses (salt, ultraviolet [UV], cold, heat, and drought) and hormone stimuli (abscisic acid [ABA], ethylene [ETH], gibberellin [GA₃], methyl jasmonate [MeJA], and salicylic acid [SA]). Expression stability in 65 samples was calculated using the △CT method, geNorm, NormFinder, BestKeeper, and RefFinder. Two reference genes (NDUFA13 and EF-1γ) were sufficient to normalize gene expression across all sample sets. They were also the two most stable genes for abiotic stresses and different tissues, whereas NDUFA13 and SKD1 were the top two choices for hormone stimuli. Considering individual experimental sets, GAPDH was the top-ranked gene under ABA, ETH, and GA₃ treatments, while 60SrRNA showed good stability under MeJA and cold treatments. ACT, UBC, and TUB were suitable genes for drought, UV, and ABA treatments, respectively. TUA was not suitable because of its considerable variation in expression under different conditions. The expression patterns of PcPAL, PcSTS, and PcMYB4 under UV and SA treatments and in different tissues normalized by stable and unstable reference genes demonstrated the suitability of the optimal reference genes.

CONCLUSIONS

We propose NDUFA13 and EF-1γ as reference genes to normalize P. cuspidatum expression data. To our knowledge, this is the first systematic study of reference genes in P. cuspidatum which could help advance molecular biology research in P. cuspidatum and allied species.

Cloning and functional analysis of five TERMINAL FLOWER 1/CENTRORADIALIS-like genes from Hevea brasiliensis

Five TERMINAL FLOWER 1 (TFL1)/CENTRORADIALIS (CEN)-like genes were isolated and characterized from rubber tree (Hevea brasiliensis). All genes, except HbCEN1, were found to have conserved genomic organization, characteristic of the phosphatidyl ethanolamine-binding protein (PEBP) family. Overexpression of all of them delayed flowering and altered flower architecture compared with the wild-type (wt) counterpart. In addition, as premature-flowering of the terminal bud was successfully overcome in the tfl1-1 mutant of Arabidopsis, all these genes have a potential function similar to TFL1. Quantitative reverse transcriptase-polymerase chain reaction analysis showed higher expressions of HbCEN1 and HbCEN2 in the shoot apices and stems of both immature and mature rubber trees than in reproductive organs. HbTFL1-1 and HbTFL1-2 expression was confined to roots of 3-month-old seedlings and HbTFL1-3 was significantly higher in the shoot apices of these seedlings. These results suggested that HbCEN1 and HbCEN2 could be associated with the development of vegetative growth, whereas HbTFL1-1, HbTFL1-2 and HbTFL1-3 seem to be mainly related with maintenance of juvenility. In addition, four of the five genes displayed variable diurnal expression, HbTFL1-1 and HbTFL1-3 being mainly expressed during the night whereas HbCEN1 and HbCEN2 showed irregular diurnal rhythms.

Asymmetric birth and death of type I and type II MADS-box gene subfamilies in the rubber tree facilitating laticifer development

The rubber tree (Hevea brasiliensis Muell. Arg.) is a rubber producing crop and contains specialized laticifers. MADS-box genes are a family of transcription factor genes that regulate plant development, especially floral organ and gametophyte development. 97 MADS-box genes were identified in the rubber tree through transcriptomes and genome mining. 93.8% of the genes were mapped onto the genome scaffolds in correspondence to the coverage (93.8%) of current version of sequenced genome. Phylogenetic analysis indicates that type II MADS-box genes have been more actively duplicated than their orthologous genes in Arabidopsis and rice, so that most (70, 72.2%) of the MADS-box genes in the rubber tree belong to type II subfamily. This is a high percentage compared to those in Arabidopsis (43.7%) and rice (56.8%). Moreover, 69 out of 70 type II genes in the rubber tree are transcribed, and they are mostly predominantly expressed in flowers, but some genes are predominantly expressed in laticifers, suggesting their roles in both flower and laticifer development. The number of type I genes in the rubber tree is only 27 (27.8%), a much smaller number compared to their orthologous genes in Arabidopsis (56.3%) and rice (43.2%). At the same time, most of the type I genes (55.6%, 15) in the rubber tree are silent and are probably pseudogenes. The high birth rate and low death rate of type II genes and low birth rate and high death rate of type I genes may corresponds to special developmental requirements in the rubber tree, e.g. the development of laticifer system for biosynthesis of cis-polyisoprene, the rubber. Moreover, atypical MIKC* factors (e.g. HbMADS1 in S-clade, and HbMADS20 in P-clade) are identified. These genes are diverged to typical MIKC* genes in sequences and facilitate functions required in laticifer development and rubber biosynthesis, which is not necessary in Arabidopsis and rice.

Sulphated Polysaccharide from Acanthophora spicifera Induced Hevea brasiliensis Defense Responses Against Phytophthora palmivora Infection

Elicitors from seaweeds are considered an alternative stimulant of plant defenses against pathogenic infection. Finding new sources of elicitors and exploring their effects on plant defenses is a significant undertaking. In this study, we extracted crude polysaccharide (CPS) from Acanthophora spicifera (a red alga) and tested the effects of the compound on rubber tree (Hevea brasiliensis) defense responses. Accumulations of salicylic acid (SA) and scopoletin (Scp) were measured by HPLC. The expression of SA- and Jasmonic acid (JA)-responsive genes was analyzed by semi-qRT-PCR. Strong anion exchange chromatography and Fourier-transform infrared (FTIR) spectroscopy were used for purification and functional characterization of CPS, respectively. The extracted CPS enhanced rubber tree defenses against Phytophthora palmivora infection. It induced SA and Scp accumulations and SA-responsive gene expression, but suppressed JA-responsive gene expression. We successfully separated the non-sulphated polysaccharide (F1) from the sulphated polysaccharides (SPS). Both peaks of SPS (F2 and F3) were identified as lambda (λ)-carrageenan. The F3 fraction showed greater elicitor activity on tobacco leaves. It induced SA and Scp accumulations and peroxidase activity but suppressed catalase activity. Furthermore, the purified λ-carrageenan did not cause cell death in tobacco or rubber tree leaves. Therefore, the elicitor from A. spicifera could be an alternative plant stimulant.

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.

Evaluation of reference genes for RT-qPCR analysis in wild and cultivated Cannabis

RT-qPCR has been widely used for gene expression analysis in recent years. The accuracy of this technique largely depends on the selection of suitable reference genes. In order to facilitate gene expression analysis in wild and cultivated Cannabis, the expression stability of seven candidate reference genes (ACT2, 18S rRNA, GAPDH, UBQ, TUB, PP2A and EF1α) were assessed in leaves samples of different development stages and different organs of both wild and cultivated Cannabis in the present study. Their expression stabilities were evaluated through three software packages (GeNorm, Normfinder and Bestkeeper). Results showed that UBQ and EF1α were the highly ranked genes in different leaves samples, and PP2A was the most stable reference gene in different organs, while GAPDH was the least stable one. And the validation of the reference genes selected was further confirmed by the expression patterns of MDS and OLS.

Salicylic Acid Induces Resistance in Rubber Tree against Phytophthora palmivora

Induced resistance by elicitors is considered to be an eco-friendly strategy to stimulate plant defense against pathogen attack. In this study, we elucidated the effect of salicylic acid (SA) on induced resistance in rubber tree against Phytophthora palmivora and evaluated the possible defense mechanisms that were involved. For SA pretreatment, rubber tree exhibited a significant reduction in disease severity by 41%. Consistent with the occurrence of induced resistance, the pronounced increase in H₂O₂ level, catalase (CAT) and peroxidase (POD) activities were observed. For defense reactions, exogenous SA promoted the increases of H₂O₂, CAT, POD and phenylalanine ammonia lyase (PAL) activities, including lignin, endogenous SA and scopoletin (Scp) contents. However, SA had different effects on the activity of each CAT isoform in the particular rubber tree organs. Besides, three partial cDNAs encoding CAT (HbCAT1, HbCAT2 and HbCAT3) and a partial cDNA encoding PAL (HbPAL) were isolated from rubber tree. Moreover, the expressions of HbCAT1, HbPAL and HbPR1 were induced by SA. Our findings suggested that, upon SA priming, the elevated H₂O₂, CAT, POD and PAL activities, lignin, endogenous SA and Scp contents, including the up-regulated HbCAT1, HbPAL and HbPR1 expressions could potentiate the resistance in rubber tree against P. palmivora.

Characterization of the rubber tree metallothionein family reveals a role in mitigating the effects of reactive oxygen species associated with physiological stress

Metallothioneins (MTs) as reactive oxygen species (ROS) scavengers play important roles in stress response and heavy metal homeostasis. In Hevea brasiliensis (the para rubber tree that is the source of commercial natural rubber) and in other trees, the functions of MTs are not well understood. Latex exudes when the rubber tree is tapped. The flow of latex and its regeneration can be enhanced by tapping, wounding and ethylene treatment, all of which produce ROS as a by-product. Here, we show the presence of four MT genes in H. brasiliensis, comprising three Type 2 (HbMT2, -2a and -2b) and one Type 3 (HbMT3L) isoforms, representing one of the smallest MT gene families among angiosperms. The four HbMTs exhibited distinct tissue expression patterns: HbMT2 and HbMT3L mainly in leaves, HbMT2a specifically in flowers and HbMT2b in diverse tissues. The expression of HbMT2b, an isoform present in latex, decreased significantly in the latex following the stress-inducing treatments of tapping, wounding and ethephon (an ethylene generator). The expressions of the leaf-abundant isoforms, HbMT2 and -3L were up-regulated following pathogenic fungus infection and high-temperature stress, but down-regulated by low-temperature stress. These reactions were consistent with multiple defense- and hormone-responsive cis-acting elements in the HbMT promoters. Nine transcription factors were shown to implicate in the high-temperature responsiveness of HbMT2 and -3L in leaves. Overexpression of HbMT2 in Escherichia coli enhanced the bacterium's tolerance to heavy metals and ROS, consistent with its predicted role as an ROS scavenger. Taken together, our results, along with other relevant studies, suggest an important role of HbMTs in latex regeneration as well as species adaptation via the regulation of ROS homeostasis.

Chilling-induced DNA Demethylation is associated with the cold tolerance of Hevea brasiliensis

BACKGROUND

Low temperature influences the development and latex production of rubber trees (Hevea brasiliensis) when extension to suboptimal high-latitude areas. The successful extension of Hevea brasiliensis cultivation to high-latitude areas has long believed to benefit from the breeding of cold-tolerant cultivars. A puzzling incongruity is the variation in cold tolerance among the cultivated clones despite their similar genetic make-up.

RESULTS

To investigate this, we first transferred cultivar Reyan 7-33-97 to short-term cold treatment, and showed that cold-related genes (such as HbICE1 and HbCBF2), cold-responsive (COR) genes, and DNA-methylation related genes (such as HbMET1) were induced by cold treatment. Furthermore, long-term cold treatment not only elevated the transcriptional activities of the HbICE1, HbCBF2, and HbMET, but also induced DNA demethylation of their promoters. Cold treatment increased the transcriptional activities of demethylation-related genes such as the HbDME, HbROS, and HbDML genes, but did not alter the promoter methylation status. Furthermore, the HbICE1 and HbMET promoters showed hypomethylation status in samples collected at the end of winter from 12 different cultivars grown in four geographical locations, but switched to hypermethylation status at the end of summer. Expression of COR was correlated with the low temperature. Given that little genetic diversity exists in the HbICE1 and HbMET promoters among different cultivars, the DNA demethylation induced by cold was highly correlated with low temperature, but not with the genetic backgrounds of cultivars.

CONCLUSION

Cold-induced epigenetic modification might play an important role in cold tolerance of H. brasiliensis.

Comparative transcriptome analysis reveals phytohormone signalings, heat shock module and ROS scavenger mediate the cold-tolerance of rubber tree

Two contrasting cold response rubber tree clones, the cold-resistant ‘93-114’ and cold-sensitive ‘Reken501’, were subject to a global transcriptome response assessing via high-throughput RNA-seq technique and comprehensive bioinformatics analysis using the referenced rubber tree genome with the purpose of exploring the potential molecular cues underlying the tolerance of rubber trees to cold stress. As a result, a total of 1919 genes had significantly higher expression, while 2929 genes had significantly lower expression in ‘93–114’ than in ‘Reken501’ without cold stress. Upon cold stress, the numbers of genes with significantly higher expression decreased to 1501 at 1 h treatment and to 1285 at 24 h treatment in ‘93–114’ than that of ‘Reken501’, conversely, the numbers of genes with significantly lower expression increased to 7567 at 1 h treatment and to 5482 at 24 h treatment. Functional annotation of the differentially expressed genes between ‘93–114’ and ‘Reken501’ suggests that down-regulation of auxin and ethylene signaling and activation of heat shock module and ROS scavengers is a primary strategy for H. brasiliensis to cope with cold stress. Our identified vital differentially expressed genes may be beneficial for elucidation of the molecular mechanisms underlying cold tolerance and for genetic improvement of H. brasiliensis clones.

Characterization of Sugar Contents and Sucrose Metabolizing Enzymes in Developing Leaves of Hevea brasiliensis

Sucrose-metabolizing enzymes in plant leaves have hitherto been investigated mainly in temperate plants, and rarely conducted in tandem with gene expression and sugar analysis. Here, we investigated the sugar content, gene expression, and the activity of sucrose-metabolizing enzymes in the leaves of Hevea brasiliensis, a tropical tree widely cultivated for natural rubber. Sucrose, fructose and glucose were the major sugars detected in Hevea leaves at four developmental stages (I to IV), with starch and quebrachitol as minor saccharides. Fructose and glucose contents increased until stage III, but decreased strongly at stage IV (mature leaves). On the other hand, sucrose increased continuously throughout leaf development. Activities of all sucrose-cleaving enzymes decreased markedly at maturation, consistent with transcript decline for most of their encoding genes. Activity of sucrose phosphate synthase (SPS) was low in spite of its high transcript levels at maturation. Hence, the high sucrose content in mature leaves was not due to increased sucrose-synthesizing activity, but more to the decline in sucrose cleavage. Gene expression and activities of sucrose-metabolizing enzymes in Hevea leaves showed striking differences compared with other plants. Unlike in most other species where vacuolar invertase predominates in sucrose cleavage in developing leaves, cytoplasmic invertase and sucrose synthase (cleavage direction) also featured prominently in Hevea. Whereas SPS is normally responsible for sucrose synthesis in plant leaves, sucrose synthase (synthesis direction) was comparable or higher than that of SPS in Hevea leaves. Mature Hevea leaves had an unusually high sucrose:starch ratio of about 11, the highest reported to date in plants.

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

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

Cloning and Expression Analysis of HbPR-1b and HbPR-3 in Hevea brasiliensis During Inoculation with Rigidoporus microporus

BACKGROUND AND OBJECTIVE

Pathogenesis-related (PR) proteins are dramatically accumulated after pathogen infection. Innate defense response through increasing PR-proteins is important for rubber rootstock selection that is tolerant to the white root disease caused by Rigidoporus microporus. This study was aimed to investigate the expression levels of PR-1 and PR-3 genes in tolerant (PB5/51) and susceptible (BPM24 and RRIM600) rubber clones after R. microporus infection.

MATERIALS AND METHODS

The mRNA of HbPR-1b and HbPR-3 was isolated and characterized from rubber leaves. Gene expression levels of HbPR-1b and HbPR-3 were compared among three rubber clones (PB5/51, BPM24 and RRIM600) after R. microporus infection at 0, 12, 24, 48, 72 and 96 h using quantitative real-time PCR. The relative transcript abundances between inoculated and control plants were compared using the means of gene expression between time points and by Tukey's HSD test. A probability value (p<0.05) was used to determine the significance of difference between time points.

RESULTS

The open reading frame of HbPR-1b is 492 bp with deduced 163 amino acid residues and the phylogenetic analysis showed it shared significant evolutionary history and clustering into group I of PR-protein. Moreover, the partial HbPR-3 was isolated with 390 bp. Gene expression levels of HbPR-1b and HbPR-3 showed marked differences in both transcripts depending on the rubber clones. Two genes demonstrated up-regulation of both tolerance and susceptibility in response to attack by R. microporus. The highest expression levels were found in seedlings of PB5/51 after inoculation. In RRIM600, low expression levels of HbPR-1b and HbPR-3 were initially observed but gradually increased at 24 h post inoculation. The transcription profile of HbPR-1b was stable expression in BPM24.

CONCLUSION

The results demonstrated that the level ofHbPR-1b and HbPR-3 transcription can distinguish between tolerant and susceptible clones. The candidate defense genes to the white root disease were observed in PB5/51 seedlings, particularly HbPR-1b.

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

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

The plant defense and pathogen counterdefense mediated by Hevea brasiliensis serine protease HbSPA and Phytophthora palmivora extracellular protease inhibitor PpEPI10

Rubber tree (Hevea brasiliensis Muell. Arg) is an important economic crop in Thailand. Leaf fall and black stripe diseases caused by the aggressive oomycete pathogen Phytophthora palmivora, cause deleterious damage on rubber tree growth leading to decrease of latex production. To gain insights into the molecular function of H. brasiliensis subtilisin-like serine proteases, the HbSPA, HbSPB, and HbSPC genes were transiently expressed in Nicotiana benthamiana via agroinfiltration. A functional protease encoded by HbSPA was successfully expressed in the apoplast of N. benthamiana leaves. Transient expression of HbSPA in N. benthamiana leaves enhanced resistance to P. palmivora, suggesting that HbSPA plays an important role in plant defense. P. palmivora Kazal-like extracellular protease inhibitor 10 (PpEPI10), an apoplastic effector, has been implicated in pathogenicity through the suppression of H. brasiliensis protease. Semi-quantitative RT-PCR revealed that the PpEPI10 gene was significantly up-regulated during colonization of rubber tree by P. palmivora. Concurrently, the HbSPA gene was highly expressed during infection. To investigate a possible interaction between HbSPA and PpEPI10, the recombinant PpEPI10 protein (rPpEPI10) was expressed in Escherichia coli and purified using affinity chromatography. In-gel zymogram and co-immunoprecipitation (co-IP) assays demonstrated that rPpEPI10 specifically inhibited and interacted with HbSPA. The targeting of HbSPA by PpEPI10 revealed a defense-counterdefense mechanism, which is mediated by plant protease and pathogen protease inhibitor, in H. brasiliensis-P. palmivora interactions.

Reliable reference genes for normalization of gene expression data in tea plants (Camellia sinensis) exposed to metal stresses

Tea plants [Camellia sinensis (L.) O. Kuntze] are an important leaf-type crop that are widely used for the production of non-alcoholic beverages in the world. Exposure to excessive amounts of heavy metals adversely affects the quality and yield of tea leaves. To analyze the molecular responses of tea plants to heavy metals, a reliable quantification of gene expression is important and of major importance herein is the normalization of the measured expression levels for the target genes. Ideally, stably expressed reference genes should be evaluated in all experimental systems. In this study, 12 candidate reference genes (i.e., 18S rRNA, Actin, CYP, EF-1α, eIF-4α, GAPDH, MON1, PP2AA3, TBP, TIP41, TUA, and UBC) were cloned from tea plants, and the stability of their expression was examined systematically in 60 samples exposed to diverse heavy metals (i.e., manganese, aluminum, copper, iron, and zinc). Three Excel-based algorithms (geNorm, NormFinder, and BestKeeper) were used to evaluate the expression stability of these genes. PP2AA3 and 18S rRNA were the most stably expressed genes, even though their expression profiles exhibited some variability. Moreover, commonly used reference genes (i.e., GAPDH and TBP) were the least appropriate reference genes for most samples. To further validate the suitability of the analyzed reference genes, the expression level of a phytochelatin synthase gene (i.e., CsPCS1) was determined using the putative reference genes for data normalizations. Our results may be beneficial for future studies involving the quantification of relative gene expression levels in tea plants.

Effect of HbDHN1 and HbDHN2 Genes on Abiotic Stress Responses in Arabidopsis

Dehydrin is a type of late embryogenesis abundant (LEA) protein. The dehydrin genes, HbDHN1 and HbDHN2, in Hevea brasiliensis were previously found to be induced at the wounding site of epicormic shoots, with local tissue dehydration identified as the key signal for laticifer differentiation. However, the exact role of the HbDHNs remains unknown. In this study, HbDHN1 and HbDHN2 expression was examined under multiple abiotic stresses; namely, cold, salt, drought, wounding, abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA) treatment. Although, both HbDHNs were defined as SK2-type dehydrin, they showed different cellular localizations. Overexpression of the HbDHNs in Arabidopsis thaliana further revealed a significant increase in tolerance to salt, drought and osmotic stresses. Increased accumulation of proline and a reduction in electrolyte leakage were also observed under salt and drought stress, and a higher water content was indicated under osmotic stress. The transgenic plants also showed higher activity levels of ascorbate peroxidase (APX), superoxide dismutase (SOD) and catalase, and accumulated less hydrogen peroxide (H₂O₂) and superoxide ([Formula: see text]). Given that reactive oxygen species (ROS) are thought to be a key signal for laticifer differentiation, these findings suggest that HbDHNs act as ROS scavengers, directly or indirectly affecting laticifer differentiation. Both HbDHNs therefore influence physiological processes, improving plant tolerance to multiple abiotic stresses.

The calcium-dependent protein kinase (CDPK) and CDPK-related kinase gene families in Hevea brasiliensis-comparison with five other plant species in structure, evolution, and expression

Calcium‐dependent protein kinases (CDPKs or CPKs) play important roles in various physiological processes of plants, including growth and development, stress responses and hormone signaling. Although the CDPK gene family has been characterized in several model plants, little is known about this gene family in Hevea brasiliensis (the Para rubber tree). Here, we characterize the entire H. brasiliensis CDPK and CDPK‐related kinase (CRK) gene families comprising 30 CDPK genes (HbCPK1 to 30) and nine CRK genes (HbCRK1 to 9). Structure and phylogeny analyses of these CDPK and CRK genes demonstrate evolutionary conservation in these gene families across H. brasiliensis and other plant species. The expression of HbCPK and HbCRK genes was investigated via Solexa sequencing in a range of experimental conditions (different tissues, phases of leaf development, ethylene treatment, and various abiotic stresses). The results suggest that HbCPK and HbCRK genes are important components in growth, development, and stress responses of H. brasiliensis. Parallel studies on the CDPK and CRK gene families were also extended to five other plant species (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Manihot esculenta, and Ricinus communis). The CDPK and CRK genes from different plant species that exhibit similar expression patterns tend to cluster together, suggesting a coevolution of gene structure and expression behavior in higher plants. The results serve as a foundation to further functional studies of these gene families in H. brasiliensis as well as in the whole plant kingdom.

Transcriptome Analysis of the Signalling Networks in Coronatine-Induced Secondary Laticifer Differentiation from Vascular Cambia in Rubber Trees

The secondary laticifer in rubber tree (Hevea brasiliensis Muell. Arg.) is a specific tissue within the secondary phloem. This tissue differentiates from the vascular cambia, and its function is natural rubber biosynthesis and storage. Given that jasmonates play a pivotal role in secondary laticifer differentiation, we established an experimental system with jasmonate (JA) mimic coronatine (COR) for studying the secondary laticifer differentiation: in this system, differentiation occurs within five days of the treatment of epicormic shoots with COR. In the present study, the experimental system was used to perform transcriptome sequencing and gene expression analysis. A total of 67,873 unigenes were assembled, and 50,548 unigenes were mapped at least in one public database. Of these being annotated unigenes, 15,780 unigenes were differentially expressed early after COR treatment, and 19,824 unigenes were differentially expressed late after COR treatment. At the early stage, 8,646 unigenes were up-regulated, while 7,134 unigenes were down-regulated. At the late stage, the numbers of up- and down-regulated unigenes were 7,711 and 12,113, respectively. The annotation data and gene expression analysis of the differentially expressed unigenes suggest that JA-mediated signalling, Ca²⁺ signal transduction and the CLAVATA-MAPK-WOX signalling pathway may be involved in regulating secondary laticifer differentiation in rubber trees.

De novo Transcriptome Analysis Reveals Distinct Defense Mechanisms by Young and Mature Leaves of Hevea brasiliensis (Para Rubber Tree)

Along with changes in morphology in the course of maturation, leaves of Hevea brasiliensis become more resistant to leaf diseases, including the South American Leaf Blight (SALB), a devastating fungal disease of this economically important tree species. To understand the underlying mechanisms of this defense, and to identify the candidate genes involved, we sequenced the Hevea leaf transcriptome at four developmental stages (I to IV) by Illumina sequencing. A total of 62.6 million high-quality reads were generated, and assembled into 98,796 unique transcripts. We identified 3,905 differentially expressed genes implicated in leaf development, 67.8% (2,651) of which were during the transition to leaf maturation. The genes involved in cyanogenic metabolism, lignin and anthocyanin biosynthesis were noteworthy for their distinct patterns of expression between developing leaves (stages I to III) and mature leaves (stage IV), and the correlation with the change in resistance to SALB and the Oidium/Colletotrichum leaf fall. The results provide a first profile of the molecular events that relate to the dynamics of leaf morphology and defense strategies during Hevea leaf development. This dataset is beneficial to devising strategies to engineer resistance to leaf diseases as well as other in-depth studies in Hevea tree.

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

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

Expression analysis of ROS producing and scavenging enzyme-encoding genes in rubber tree infected by Pseudocercospora ulei

South American Leaf Blight (SALB), caused by the ascomycete Pseudocercospora ulei, is responsible for the low productivity of rubber trees in Latin America and is a serious threat to rubber plantations in Asia and Africa, where the rubber trees are derived from highly susceptible clones. Three contrasted genotypes were chosen for their levels of resistance to the pathogen: FX2784 (totally resistant), MDF180 (partially resistant) and PB314 (susceptible). Array analyses were previously performed to identify genes differentially expressed in resistant and susceptible genotypes. Twenty-one genes were selected for further gene expression analysis in non-inoculated and inoculated genotypes from 24 to 216 h post infection (hpi). These genes are involved in ROS production (HbRBOHA, HbRBOHB, HbRBOHC, HbRBOHD), ROS-scavenging systems (cytoplasmic and chloroplastic HbCuZnSOD, HbMnSOD, HbCAT, HbAPX1, HbAPX2, HbMDHAR, HbGCL1, HbGCL2, HbOASTL, HbGPX, HbDHAR), and leaf senescence (HbCASP, HbPCYST, HbWRKY2, HbPLY, HbKAT2). First, a genotype-dependent level of expression was observed. The genes HbRBOHA, HbCuZnSOD cyto, HbCAT, HbGCL and HbWRKY2 were constitutively expressed at lower levels in the MDF180 genotype than in the FX2784 and PB314 genotypes. Conversely, the levels of expression of HbDHAR, HbGPX and HbPCYST were higher in the older, non-inoculated leaves of MDF180. Lower production of ROS and efficient regeneration of reduced ascorbate ensure a balanced redox intracellular state in this genotype. Second, inoculation of the leaves induced few modifications in the expression level of the studied genes. In the MDF180 partially resistant genotype, an increase in the expression level of HbRBOHB, HbRBOHD 48 hpi and a decrease in the expression level of HbDHAR 216 hpi were observed. In the FX2784 totally resistant genotype, an increase in the expression level of HbRBOHD and HbCuZnSOD cyto and a decrease in HbCAT were observed 48 hpi. This transitory variation could be associated with the oxidative burst classically observed in hypersensitive response (HR). The increase in the synthesis of reduced glutathione in this genotype could ensure redox balance and consequently cell homeostasis. In the PB314 susceptible genotype, HbROHC, HbCuZnSOD chloro was up-regulated 216 hpi concomitantly with a decrease in the expression level of HbCAT, consequently causing an accumulation of H2O2 and programmed cell death. The level of expression of a transcription factor, HbWRKY2, was also modulated by the P. ulei infection with early transient up-regulation in the FX2784 totally resistant genotype and permanent up-regulation in the MDF180 partially resistant genotype. These results complement studies on genetic determinism of SALB resistance and a recent publication on Hevea glutathione reductase gene.

The Hevea brasiliensis XIP aquaporin subfamily: genomic, structural and functional characterizations with relevance to intensive latex harvesting

X-Intrinsic Proteins (XIP) were recently identified in a narrow range of plants as a full clade within the aquaporins. These channels reportedly facilitate the transport of a wide range of hydrophobic solutes. The functional roles of XIP in planta remain poorly identified. In this study, we found three XIP genes (HbXIP1;1, HbXIP2;1 and HbXIP3;1) in the Hevea brasiliensis genome. Comprehensive bioinformatics, biochemical and structural analyses were used to acquire a better understanding of this AQP subfamily. Phylogenetic analysis revealed that HbXIPs clustered into two major groups, each distributed in a specific lineage of the order Malpighiales. Tissue-specific expression profiles showed that only HbXIP2;1 was expressed in all the vegetative tissues tested (leaves, stem, bark, xylem and latex), suggesting that HbXIP2;1 could take part in a wide range of cellular processes. This is particularly relevant to the rubber-producing laticiferous system, where this isoform was found to be up-regulated during tapping and ethylene treatments. Furthermore, the XIP transcriptional pattern is significantly correlated to latex production level. Structural comparison with SoPIP2;1 from Spinacia oleracea species provides new insights into the possible role of structural checkpoints by which HbXIP2;1 ensures glycerol transfer across the membrane. From these results, we discuss the physiological involvement of glycerol and HbXIP2;1 in water homeostasis and carbon stream of challenged laticifers. The characterization of HbXIP2;1 during rubber tree tapping lends new insights into molecular and physiological response processes of laticifer metabolism in the context of latex exploitation.

Identification and Validation of Reference Genes for RT-qPCR Analysis in Non-Heading Chinese Cabbage Flowers

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis Makino) is an important vegetable member of Brassica rapa crops. It exhibits a typical sporophytic self-incompatibility (SI) system and is an ideal model plant to explore the mechanism of SI. Gene expression research are frequently used to unravel the complex genetic mechanism and in such studies appropriate reference selection is vital. Validation of reference genes have neither been conducted in Brassica rapa flowers nor in SI trait. In this study, 13 candidate reference genes were selected and examined systematically in 96 non-heading Chinese cabbage flower samples that represent four strategic groups in compatible and self-incompatible lines of non-heading Chinese cabbage. Two RT-qPCR analysis software, geNorm and NormFinder, were used to evaluate the expression stability of these genes systematically. Results revealed that best-ranked references genes should be selected according to specific sample subsets. DNAJ, UKN1, and PP2A were identified as the most stable reference genes among all samples. Moreover, our research further revealed that the widely used reference genes, CYP and ACP, were the least suitable reference genes in most non-heading Chinese cabbage flower sample sets. To further validate the suitability of the reference genes identified in this study, the expression level of SRK and Exo70A1 genes which play important roles in regulating interaction between pollen and stigma were studied. Our study presented the first systematic study of reference gene(s) selection for SI study and provided guidelines to obtain more accurate RT-qPCR results in non-heading Chinese cabbage.

Expression Stabilities of Candidate Reference Genes for RT-qPCR in Chinese Jujube (Ziziphus jujuba Mill.) under a Variety of Conditions

Reverse transcription-quantitative real-time polymerase chain reaction (RT-qPCR) is a powerful method for evaluating patterns of gene expression. Jujube whole-genome sequencing has been completed, and analysis of gene function, an important part of any follow-up study, requires the appropriate selection of reference genes. Indeed, suitable reference gene selection for RT-qPCR is critical for accurate normalization of target gene expression. In this study, the software packages geNorm and NormFinder were employed to examine the expression stabilities of nine candidate reference genes under a variety of conditions. Actin-depolymerizing factor 1 (ACT1), Histone-H3 (His3), and Polyadenylate-binding protein-interacting protein (PAIP) were determined to be the most stably expressed genes during five stages of fruit development and ACT1, SiR-Fd, BTF3, and Tubulin alpha chain (TUA) across different tissues/organs. Whereas ACT1, Basic Transcription factor 3 (BTF3), Glyceraldehyde-3-phosphate dehydrogenase (GADPH), and PAIP were the most stable under dark conditions. ACT1, PAIP, BTF3, and Elongation factor 1- gamma (EF1γ) were the most stably expressed genes under phytoplasma infection. Among these genes, SiR-Fd and PAIP are here first reported as stable reference genes. When normalized using these most stable reference genes, the expression patterns of four target genes were found to be in accordance with physiological data, indicating that the reference genes selected in our study are suitable for use in such analyses. This study provides appropriate reference genes and corresponding primers for further RT-qPCR studies in Chinese jujube and emphasizes the importance of validating reference genes for gene expression analysis under variable experimental conditions.

Validation of putative reference genes for gene expression studies in heat stressed and α-MSH treated melanocyte cells of Bos indicus using real-time quantitative PCR

Normalization of cellular mRNA data using internal reference gene (IRG) is an essential step in expression analysis studies. MIQE guidelines ensure that the choice and appropriateness of IRG should be validated for particular tissues or cell types and specific experimental designs. The objective of the present study was to assess 15 IRGs from different functional classes that could serve as best IRGs for Bos indicus (Tharparkar cattle) melanocyte cells under heat stress and hormonal treatment. We implemented the use of geNorm, NormFinder and BestKeeper algorithm to measure the stability of the gene transcript. A total of 15 IRGs (ACTB, BZM, EEF1, GAPDH, GTP, HMBS, HPRT, RPL22, RPL4, RPS15, RPS18, RPS23, RPS9, UBC and UXT) from different functional classes were evaluated. Pair wise comparisons using geNorm revealed that HPRT and RPS23 were the most stable combination of IRGs with M-value of 0.29 followed by UXT (0.30) and RPL4 (0.31). The NormFinder analysis also identified the same set of stably expressed genes (UXT, RPL4, RPS23 and HPRT); however, the rank order was little different. The UXT gene showed lowest crossing point SD and CV values of 0.30 and 1.17, respectively indicating its maximum expression stability through BestKeeper analysis. The present study indicated that, ACTB and HMB were not reliable IRGs for melanocytes cells on account of their lower expression stability. Current study further revealed that UXT, HPRT and RPS23 are the best IRGs for normalization of qPCR data in Bos indicus melanocyte cells under heat stress and hormonal treatment.

Identification, Functional Study, and Promoter Analysis of HbMFT1, a Homolog of MFT from Rubber Tree (Hevea brasiliensis)

A homolog of MOTHER OF FT AND TFL1 (MFT) was isolated from Heveabrasiliensis and its biological function was investigated. Protein multiple sequence alignment and phylogenetic analysis revealed that HbMFT1 conserved critical amino acid residues to distinguish MFT, FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1)-like proteins and showed a closer genetic relationship to the MFT-like group. The accumulation of HbMFT1 was generally detected in various tissues except pericarps, with the highest expression in embryos and relatively higher expression in roots and stems of seedlings, flowering inflorescences, and male and female flowers. HbMFT1 putative promoter analysis showed that tissue-specific, environmental change responsive and hormone-signaling responsive elements were generally present. HbMFT1 was strongly induced under a short-day condition at 28 °C, with the highest expression after the onset of a day. Overexpression of HbMFT1 inhibited seed germination, seedling growth, and flowering in transgenic Arabidopsis. The qRT-PCR further confirmed that APETALA1 (AP1) and FRUITFULL (FUL) were drastically down-regulated in 35S::HbMFT1 plants. A histochemical β-glucuronidase (GUS) assay showed that HbMFT1::GUS activity was mainly detected in stamens and mature seeds coinciding with its original expression and notably induced in rosette leaves and seedlings of transgenic Arabidopsis by exogenous abscisic acid (ABA) due to the presence of ABA cis-elements in HbMFT1 promoter. These results suggested that HbMFT1 was mainly involved in maintenance of seed maturation and stamen development, but negatively controlled germination, growth and development of seedlings and flowering. In addition, the HbMFT1 promoter can be utilized in controlling transgene expression in stamens and seeds of rubber tree or other plant species.

Identification and Characterization of the Glucose-6-Phosphate Dehydrogenase Gene Family in the Para Rubber Tree, Hevea brasiliensis

As a key enzyme in the pentose phosphate pathway (PPP), glucose-6-phosphate dehydrogenase (G6PDH) provides nicotinamide adenine dinucleotide phosphate (NADPH) and intermediary metabolites for rubber biosynthesis, and plays an important role in plant development and stress responses. In this study, four Hevea brasiliensis (Para rubber tree) G6PDH genes (HbG6PDH1 to 4) were identified and cloned using a genome-wide scanning approach. All four HbG6PDH genes encode functional G6PDH enzymes as shown by heterologous expression in E. coli. Phylogeny analysis and subcellular localization prediction show that HbG6PDH3 is a cytosolic isoform, while the other three genes (HbG6PDH1, 2 and 4) are plastidic isoforms. The subcellular locations of HbG6PDH3 and 4, two latex-abundant isoforms were further verified by transient expression in rice protoplasts. Enzyme activity assay and expression analysis showed HbG6PDH3 and 4 were implicated in PPP during latex regeneration, and to influence rubber production positively in rubber tree. The cytosolic HbG6PDH3 is a predominant isoform in latex, implying a principal role for this isoform in controlling carbon flow and NADPH production in the PPP during latex regeneration. The expression pattern of plastidic HbG6PDH4 correlates well with the degree of tapping panel dryness, a physiological disorder that stops the flow of latex from affected rubber trees. In addition, the four HbG6PDHs responded to temperature and drought stresses in root, bark, and leaves, implicating their roles in maintaining redox balance and defending against oxidative stress.

Identification of Powdery Mildew Responsive Genes in Hevea brasiliensis through mRNA Differential Display

Powdery mildew is an important disease of rubber trees caused by Oidium heveae B. A. Steinmann. As far as we know, none of the resistance genes related to powdery mildew have been isolated from the rubber tree. There is little information available at the molecular level regarding how a rubber tree develops defense mechanisms against this pathogen. We have studied rubber tree mRNA transcripts from the resistant RRIC52 cultivar by differential display analysis. Leaves inoculated with the spores of O. heveae were collected from 0 to 120 hpi in order to identify pathogen-regulated genes at different infection stages. We identified 78 rubber tree genes that were differentially expressed during the plant–pathogen interaction. BLAST analysis for these 78 ESTs classified them into seven functional groups: cell wall and membrane pathways, transcription factor and regulatory proteins, transporters, signal transduction, phytoalexin biosynthesis, other metabolism functions, and unknown functions. The gene expression for eight of these genes was validated by qRT-PCR in both RRIC52 and the partially susceptible Reyan 7-33-97 cultivars, revealing the similar or differential changes of gene expressions between these two cultivars. This study has improved our overall understanding of the molecular mechanisms of rubber tree resistance to powdery mildew.