Heterozygosity analysis of spontaneous 2n female gametes and centromere mapping of the diploid Hevea brasiliensis based on full-sib triploid populations

KEY MESSAGE

Unreduced megagametophytes via second-division restitution were confirmed through heterozygosity analysis, and four candidate physical centromeres of rubber were located for the first time. The evaluation of maternal heterozygosity restitution (MHR) is vital in identifying the mechanism of 2n gametogenesis and assessing the utilization value of 2n gametes. In this study, three full-sib triploid populations were employed to evaluate the MHR of 2n female gametes of rubber tree clone GT1 and to confirm their genetic derivation. The 2n female gametes of GT1 were derived from second-division restitution (SDR) and transmitted more than half of the parental heterozygosity. In addition, low recombination frequency markers were developed, and four candidate physical centromeres of rubber tree were located for the first time. The confirmation that 2n female gametes of rubber tree clone GT1 are derived from SDR provides insights into the molecular mechanisms of 2n gametogenesis. In addition, the identified centromere location will aid in the development of centromeric markers for the rapid identification of the 2n gametogenesis mechanism.

Bark transcriptome analyses reveals molecular mechanisms involved in tapping panel dryness occurrence and development in rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) has become the mostimportant limiting factor for increasing natural rubber yield, whereas illuminating the molecular mechanisms underlying TPD is the prerequisite for solving the problem of TPD. However, molecular mechanisms underlying TPD are largely unknown. In this study, healthy and different stages of TPD-affected rubber trees were utilized to analyze TPD for the first time. We found that the changing tendencies of key latex physiological parameters were closely related to TPD occurrence and development. To reveal the molecular mechanisms underlying TPD, we sequenced and compared bark transcriptomes among healthy rubber tree, and TPD-affected ones at initial and advanced stages. In total, 8607 genes were identified as TPD-related genes in contrast to healthy rubber tree. According to gene expression profiles, the five samples were divided into three groups including healthy rubber tree, and TPD-affected rubber tree in the initial and advanced stages, which was consistent with the stages of TPD occurrence and development. Interestingly, only asmall proportionof the TPD-related genes were constantly down- or up-regulated with TPD occurrence and development. The TPD-related genes in KEGG pathways significantly enriched were closely associated with protein metabolism, cell division and differentiation, PCD, stress responses, terpene biosynthesis, and various metabolism processes. Moreover, overexpression of HbAPX2 identified as a TPD-related gene enhanced oxidative stress tolerance in S. cerevisiae. The typical symptoms of TPD, partial or complete dry zone (no latex flow) on tapping panel, might attribute to lower IPP available for rubber biosynthesis, and downregulation of the genes in post-IPP steps of rubber biosynthesis and the genes involved in latex flow. Our results not only provide new insights into molecular mechanisms underlying TPD occurrence and development but also contribute to developing effective measures to control TPD in rubber trees.

The impact of non-deproteinization on physicochemical and biological properties of natural rubber latex for biomedical applications

Latex is a colloidal suspension derived from the Hevea brasiliensis tree, derived from natural rubber, poly(isoprene), and assorted constituents including proteins and phospholipids. These constituents are inherent to both natural rubber and latex serum. This investigation was undertaken to examine the impact of the deproteinization process on chemical and biological dynamics of natural rubber latex. Natural Rubber (NR) extracted from the pure latex (LNCP) was obtained through centrifugation, followed by six rounds of solvent purification (LP6). The structure was characterized using Fourier-transform infrared spectroscopy (FTIR), thermogravimetric analysis (TGA), differential scanning calorimetry (DSC), swelling test, surface zeta potential (ζ), scanning electron microscopy (SEM) and in vitro assay. The results revealed that the LP6 group presented decreased swelling kinetics, reduced cell adhesion and proliferation, and a smoother surface with decreased negative surface charge. Conversely, the LNCP group shown accelerated swelling, heightened adhesion and cellular growth, and a more negatively charged and rougher surface. As such, the attributes of latex serum and proteins have potential usage across numerous biomedical applications.

The Interplay of Protein Hydrolysis and Ammonia in the Stability of Hevea Rubber Latex during Storage

Natural rubber (NR) latex derived from Hevea brasiliensis is a complex colloid comprising mainly rubber hydrocarbons (latex particles) and a multitude of minor non-rubber constituents such as non-rubber particles, proteins, lipids, carbohydrates, and soluble organic and inorganic substances. NR latex is susceptible to enzymatic attack after it leaves the trees. It is usually preserved with ammonia and, to a lesser extent, with other preservatives to enhance its colloidal stability during storage. Despite numerous studies in the literature on the influence of rubber proteins on NR latex stability, issues regarding the effect of protein hydrolysis in the presence of ammonia on latex stability during storage are still far from resolved. The present work aims to elucidate the interplay between protein hydrolysis and ammoniation in NR latex stability. Both high- and low-ammonia (with a secondary preservative) NR latexes were used to monitor the changes in their protein compositions during storage. High-ammonia (FNR-A) latex preserved with 0.6% (v/v) ammonia, a low 0.1% ammonia/TMTD/ZnO (FNR-TZ) latex, and a deproteinized NR (PDNR) latex were labeled with fluorescence agents and observed using confocal laser scanning microscopy to determine their protein composition. Protein hydrolysis was confirmed via sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The results revealed that protein hydrolysis increased with the storage duration. The change in protein composition accompanying hydrolysis also allows the spatial distribution of allergenic proteins to be estimated in the latex. Concurrently, the latex stability increased with the storage duration, as measured by the latex's mechanical stability time (MST) and the zeta potential of the latex particles. As monitored by AFM, the surface roughness of the NR latex film increased markedly during extended storage compared with that of the DPNR latex, which remained smooth. These results underscore the pivotal role of ammonia in bolstering NR latex stability brought on by protein hydrolysis, which greatly impacts latex film's formation behavior. NR latex stability underpins the quality of latex-dipped goods during manufacturing, particularly those for medical gloves.

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.

Antioxidant and photoprotective role of latex C-serum from Hevea brasiliensis during 15-week UVB irradiation in male hairless SKH-1 mice

It is known that UVB radiation induces several adverse skin alterations starting from simple photoaging to skin cancer. In addition, it was demonstrated that reactive oxygen species (ROS) were found to be related to cancer development and progression. The aim of study was to examine whether male hairless (SKH-1) mice (Mus musculus) that were subchronically exposed to UVB radiation presented with actinic keratosis (AK) and squamous cell carcinoma lesions, and that treatment with latex C-serum cream significantly prevented abnormal skin development. Data demonstrated for the first time the photoprotective activity of latex C-serum extracted from the rubber tree Hevea brasiliensis var. subconcolor Ducke. Latex C-serum prevented the progression of AK to squamous cell carcinoma in SKH-1 mice, indicating that mice topically treated with latex C-serum presented only AK lesions and treatment with the highest concentration (10%) significantly reduced epidermal thickness, suggesting diminished cell proliferation. Latex C-serum protected the skin of mice against oxidative stress damage, increasing catalase (CAT) activity, regenerating glutathione (GSH) levels, lowering thiobarbituric acid-reactive species (TBARS) production and regenerating the total antioxidant capacity (TAC) of the skin. Evidence that UV radiation in skin induced systemic alterations and erythrocytic analysis indicated that latex C-serum increased CAT activity and GSH levels. Taken together these data indicate that latex C-serum plays an important antioxidant and photoprotective role, preventing serious damage to the skin following exposure to UVB radiation.

HbBIN2 Functions in Plant Cold Stress Resistance through Modulation of HbICE1 Transcriptional Activity and ROS Homeostasis in Hevea brasiliensis

Cold stress poses significant limitations on the growth, latex yield, and ecological distribution of rubber trees (Hevea brasiliensis). The GSK3-like kinase plays a significant role in helping plants adapt to different biotic and abiotic stresses. However, the functions of GSK3-like kinase BR-INSENSITIVE 2 (BIN2) in Hevea brasiliensis remain elusive. Here, we identified HbBIN2s of Hevea brasiliensis and deciphered their roles in cold stress resistance. The transcript levels of HbBIN2s are upregulated by cold stress. In addition, HbBIN2s are present in both the nucleus and cytoplasm and have the ability to interact with the INDUCER OF CBF EXPRESSION1(HbICE1) transcription factor, a central component in cold signaling. HbBIN2 overexpression in Arabidopsis displays decreased tolerance to chilling stress with a lower survival rate and proline content but a higher level of electrolyte leakage (EL) and malondialdehyde (MDA) than wild type under cold stress. Meanwhile, HbBIN2 transgenic Arabidopsis treated with cold stress exhibits a significant increase in the accumulation of reactive oxygen species (ROS) and a decrease in the activity of antioxidant enzymes. Further investigation reveals that HbBIN2 inhibits the transcriptional activity of HbICE1, thereby attenuating the expression of C-REPEAT BINDING FACTOR (HbCBF1). Consistent with this, overexpression of HbBIN2 represses the expression of CBF pathway cold-regulated genes under cold stress. In conclusion, our findings indicate that HbBIN2 functions as a suppressor of cold stress resistance by modulating HbICE1 transcriptional activity and ROS homeostasis.

Interactions of REF1 and SRPP1 rubber particle proteins from Hevea brasiliensis with synthetic phospholipids: Effect of charge and size of lipid headgroup

According to the fatty acid and headgroup compositions of the phospholipids (PL) from Hevea brasiliensis latex, three synthetic PL were selected (i.e. POPA: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphate POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and POPG: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol) to investigate the effect of PL headgroup on the interactions with two major proteins of Hevea latex, i.e. Rubber Elongation Factor (REF1) and Small Rubber Particle Protein (SRPP1). Protein/lipid interactions were screened using two models (lipid vesicles in solution or lipid monolayers at air/liquid interface). Calcein leakage, surface pressure, ellipsometry, microscopy and spectroscopy revealed that both REF1 and SRPP1 displayed stronger interactions with anionic POPA and POPG, as compared to zwitterionic POPC. A particular behavior of REF1 was observed when interacting with POPA monolayers (i.e. aggregation + modification of secondary structure from α-helices to β-sheets, characteristic of its amyloid aggregated form), which might be involved in the irreversible coagulation mechanism of Hevea rubber particles.

Performance- and Resistance-Related Early Responses of Colombian Elite Rubber Tree Genotypes under Low Pressure of South American Leaf Blight: Implications for Disease Management in the Amazon

The cultivation of Hevea brasiliensis, the primary commercial source of natural rubber, is strongly impacted by South American leaf blight (SALB) disease, caused by the fungus Pseudocercospora ulei. Various management strategies have been implemented, including the selection of resistant genotypes and the identification of escape zones. This study evaluated the growth, early yield, and resistance to SALB of nine Colombian elite genotypes from the ECC-100 series and IAN 873 clone (control) in a large-scale clone trial in an area with low SALB pressure in the Colombian Amazon during 2017-2020. Favorable early performance was evident, although there was a significant increase in the severity and sporulation of P. ulei over time, especially in the ECC 35, ECC 60, and IAN 873 genotypes. However, these scores indicate low susceptibility. Genotypes with higher resistance to SALB demonstrated greater growth and early yield compared to more highly susceptible genotypes. The ECC 64, ECC 73, ECC 90, ECC 25, and ECC 29 genotypes were more desirable in low SALB pressure zones due to their higher resistance and early performance. It is important to highlight that this research contributes to the selection of new SALB-resistant Colombian genotypes of H. brasiliensis. However, it is also necessary to evaluate the productivity of these selections in the mature stage and long-term resistance to SALB before recommending and promoting their commercial adoption in the Colombian Amazon.

Distribution of Aedes mosquito species along the rural-urban gradient in Lambaréné and its surrounding

BACKGROUND

Aedes albopictus and Aedes aegypti are known for their potential as vectors of dengue (DENV) and chikungunya (CHIKV) viruses. However, entomological surveys are mostly carried out during epidemics. In Gabon where outbreaks of both viruses have occurred, there is no vector control program targeting these arboviruses. Therefore, we assessed the presence of Aedes species along a rural-urban gradient in Lambaréné (Gabon) and its surroundings and determined ecological factors associated to their presence.

METHODS

An entomological survey was conducted in Lambaréné and its surrounding rural areas. Mosquitoes were collected with aspirators around human dwellings, and ecological and environmental data were collected from each study area. Morphological identification keys were used to identify Aedes species. RNA was extracted from pools of female mosquitoes and amplified by RT-qPCR to detect the presence of DENV and CHIKV.

RESULTS

Overall, the most common vector collected was Aedes albopictus (97%, 4236/4367 specimens), followed by Aedes aegypti (3%, 131/4367). Albopictus vectors was more abundant in the rural area (Wilcoxon signed-rank test, Z = 627, P = 0.043) than in the urban area. In the urban area, a higher number of mosquitoes (45%) were recorded in the economic zone (zone 3) than in the historical and administrative zones (zone 1 and 2). In the rural area, the proportions of species numbers were significantly higher along the south rural transect (92%) compared to the north rural transect (Wilcoxon signed-rank test, Z = 43, P ˂ 0.016). We also noted a high abundance of vectors in environments characterized by monocultures of Hevea brasiliensis (Hevea) and Manihot esculenta (cassava) (Kruskal-Wallis H-test, H = 25.7, df = 2, P < 0.001). Finally, no mosquito pools were positive for either DENV or CHIKV.

CONCLUSION

Aedes albopictus was the dominant vector across the study sites due to its high invasiveness capacity. This presence re-affirms the potential for local transmission of both DENV and CHIKV, as indicated previously by serological surveys conducted in our study area, even though no transmission was detected during the current study. These findings underscore the need for regular arbovirus surveillance in the study region, with the aim of supporting vector control efforts in the event of outbreaks.

Rubber latex yield is affected by interactions between antecedent temperature, rubber phenology, and powdery mildew disease

Rubber (Hevea brasiliensis) latex production is crucial to the local economy, yet Xishuangbanna's climate is considered sub-optimal for rubber cultivation. The prevalence of the powdery mildew disease (Oidium heveae) in this region has decreased the annual latex yield by 20%. Rubber latex yield is influenced by several factors, including temperature, disease, other biotic conditions, and plantation management. However, the interrelationships and potential influencing networks between rubber latex yield and these factors are rarely quantitatively assessed, and understanding their impacts on latex yield could inform better management practices. To address this gap, we investigated the effects of temperature, phenology, and powdery mildew disease on rubber latex yield in March using observational data on daily rubber latex yield combined with detailed phenology, powdery mildew, and temperature data from 2004 to 2010 in a state farm in the Xishuangbanna, Yunnan, China. We found that the critical influencing periods of daily temperature difference (or diurnal temperature difference) on the rubber latex yield were during Nov 27-Jan 19 and Jan 21-Mar 17. Partial least square regression analysis and variance partitioning analysis were conducted on the 35 phenological variables, eight powdery mildew-related variables, and two climatic variables. The most influential factors were identified as the factors of the daily temperature differences during Jan-Mar, the duration of leaf flushing phenology, and mean and maximum percentage of leaves infected by powdery mildew. Subsequent canonical correlation analysis and linear regression found that temperature difference directly affected the rubber latex yield and indirectly affected the yield through phenology and powdery mildew disease. Raised daily temperature differences from Jan to Mar had the greatest impact, leading to a higher rubber latex yield. Our comprehensive quantitative assessment revealed the relative importance of antecedent daily temperature differences, phenology, and powdery mildew disease as well as their complex interconnections in influencing rubber latex yield. Our findings are essential to future studies on both powdery mildew disease and rubber latex yield, and also develop rubber latex models.

Latex C-serum from Hevea brasiliensis induces apoptotic cell death in a leukemic cell line

BACKGROUND

Hevea brasiliensis latex is generally cultivated for the use of rubber particles. Previous studies have shown that the antiproliferative activity of C-serum in hepatocellular carcinoma is not induced through the classical apoptotic signaling pathway. However, in a leukemic cell line, the anti-proliferation effect of latex C serum remained unclear.

METHODS

Leukemic cell lines (K562 and U937) and human peripheral blood mononuclear cells (PBMCs) were examined for cell viability using the MTT assay. Flow cytometry was used for apoptotic cell detection by annexin V/PI staining. The expression levels of proapoptotic and antiapoptotic marker genes were measured by qRT‒PCR. Moreover, the caspase activities of the extrinsic and intrinsic apoptotic pathways were detected by enzymatic activities.

RESULTS

Latex C-serum inhibited cell proliferation in the K562 and U937 leukemic cell lines but did not affect human PBMCs. Latex C-serum significantly induced the percentage of early and late apoptotic cells in the leukemic cell line. The expression levels of the pro-apoptotic marker genes BAD, BAX, and CASPASE3 significantly increased in the leukemic cell line after post-latex C-serum leukemic cell treatment. The extrinsic, intrinsic and common apoptotic pathways were also studied through caspase-8, -9, and -3 activities. Latex C-serum treatment significantly induced caspase-8, -9, and -3 activation in the K562 cell line and U937 cell line compared to the untreated cells.

CONCLUSIONS

These results indicate that latex C-serum enhanced anti-proliferation in leukemic cell lines by inducing apoptosis and caspase activation.

[Hevea brasiliensis como fuente alergénica: revisión bibliográfica]

Abstract

La planta Hevea brasiliensis se utiliza ampliamente en la industria como fuente de extracción de caucho, un elemento empleado en diversas áreas comerciales y médicas. Los estudios inmunológicos de esta especie indican que es una fuente alergénica importante, que puede provocar sensibilización y alergia. Se han identificado diferentes componentes alergénicos de esta planta, con diversas propiedades inmunitarias y bioquímicas, y estudiado más de diez tipos diferentes de alérgenos, cada uno con distinta capacidad de inducir síntomas alérgicos. En esta revisión informamos los avances actuales en el estudio de Hevea brasiliensis.

Multiomics studies with co-transformation reveal microRNAs via miRNA-TF-mRNA network participating in wood formation in Hevea brasiliensis

Introduction

MicroRNAs (miRNAs) are small endogenous non-coding RNAs that play an important role in wood formation in plants. However, the significance of the link between miRNAs and their target transcripts in wood formation remains unclear in rubber tree (Hevea brasiliensis).

Methods

In this study, we induced the formation of reaction wood by artificially bending rubber trees for 300 days and performed small RNA sequencing and transcriptome deep sequencing (RNA-seq) to describe the complement of miRNAs and their targets contributing to this process.

Results and discussion

We identified 5, 11, and 2 differentially abundant miRNAs in normal wood (NW) compared to tension wood (TW), in NW relative to opposite wood (OW), and between TW and OW, respectively. We also identified 12 novel miRNAs and 39 potential miRNA-mRNA pairs with different accumulation patterns in NW, TW, and OW. We noticed that many miRNAs targeted transcription factor genes, which were enriched in KEGG pathways associated with phenylpropanoid biosynthesis, phenylalanine metabolism, and pyruvate metabolism. Thus, miRNA-TF-mRNA network involved in wood formation via tension wood model were constructed. We validated the differential accumulation of miRNAs and their targets by RT-qPCR analysis and overexpressed miRNA in Nicotiana benthamiana with its potential target gene. These results will provide a reference for a deep exploration of growth and development in rubber tree.

Single-cell transcriptomic analyses reveal cellular and molecular patterns of rubber tree response to early powdery mildew infection

As a perennial woody plant, the rubber tree (Hevea brasiliensis) must adapt to various environmental challenges through gene expression in multiple cell types. It is still unclear how genes in this species are expressed at the cellular level and the precise mechanisms by which cells respond transcriptionally to environmental stimuli, especially in the case of pathogen infection. Here, we characterized the transcriptomes in Hevea leaves during early powdery mildew infection using single-cell RNA sequencing. We identified 10 cell types and constructed the first single-cell atlas of Hevea leaves. Distinct gene expression patterns of the cell clusters were observed under powdery mildew infection, which was especially significant in the epidermal cells. Most of the genes involved in host-pathogen interactions in epidermal cells exhibited a pattern of dramatically increased expression with increasing pseudotime. Interestingly, we found that the HbCNL2 gene, encoding a nucleotide-binding leucine-rich repeat protein, positively modulated the defence of rubber leaves against powdery mildew. Overexpression of the HbCNL2 gene triggered a typical cell death phenotype in tobacco leaves and a higher level of reactive oxygen species in the protoplasts of Hevea leaves. The HbCNL2 protein was located in the cytomembrane and nucleus, and its leucine-rich repeat domain interacted with the histidine kinase-like ATPase domain of the molecular chaperone HbHSP90 in the nucleus. Collectively, our results provide the first observation of the cellular and molecular responses of Hevea leaves to biotrophic pathogen infection and can guide the identification of disease-resistance genes in this important tree species.

Synthesis of epoxidized natural rubber grafted with hyaluronic acid for the development of biomaterials

Natural Rubber (NR), extracted from Hevea brasiliensis rubber trees, is a biocompatible biopolymer with properties that support in the tissue repair process. However, its biomedical applications are limited due to the presence of allergenic proteins, hydrophobicity, and unsaturated bonds. To overcome these limitations and contribute to the development of new biomaterials, this study aims to deproteinize, epoxidize, and subject NR to copolymerization by grafting with hyaluronic acid (HA), which is widely recognized for its bioactive properties in the medical field. The deproteinization, epoxidation, and graft copolymerization through the esterification reaction were confirmed by Fourier Transform Infrared Spectroscopy and Hydrogen Nuclear Magnetic Resonance Spectroscopy analysis. Thermogravimetry and Differential Scanning Calorimetry demonstrated that the grafted sample exhibited a lower degradation rate and a higher glass transition temperature, indicating strong intermolecular interactions. Moreover, contact angle measurement revealed that the grafted NR exhibited a high hydrophilic character. The results obtained suggest the formation of a novel material with great potential for application in biomaterials involved in tissue repair processes.

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker-trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.

Induction of Axillary Bud Swelling of Hevea brasiliensis to Regenerate Plants through Somatic Embryogenesis and Analysis of Genetic Stability

To overcome rubber tree (RT) tissue culture explant source limitations, the current study aimed to establish a new Hevea brasiliensis somatic embryogenesis (SE) system, laying the technical foundation for the establishment of an axillary-bud-based seedling regeneration system. In this study, in vitro plantlets of Hevea brasiliensis Chinese Academy of Tropical Agricultural Sciences 917 (CATAS 917) were used as the experimental materials. Firstly, the optimum conditions for axillary bud swelling were studied; then, the effects of phenology, the swelling time of axillary buds (ABs), and medium of embryogenic callus induction were studied. Plantlets were obtained through somatic embryogenesis. Flow cytometry, inter-simple sequence repeat (ISSR molecular marker) and chromosome karyotype analysis were used to study the genetic stability of regenerated plants along with budding seedlings (BSs) and secondary somatic embryo seedlings (SSESs) as the control. The results show that the rubber tree's phenology period was mature, and the axillary bud induction rate was the highest in the 2 mg/L 6-benzyladenine (6-BA) medium (up to 85.83%). Later, 3-day-old swelling axillary buds were used as explants for callogenesis and somatic embryogenesis. The callus induction rate was optimum in MH (Medium in Hevea) + 1.5 mg/L 2,4-dichlorophenoxyacetic acid (2,4-D) + 1.5 mg/L 1-naphthalene acetic acid (NAA) + 1.5 mg/L Kinetin (KT) + 70 g/L sucrose (56.55%). The regenerated plants were obtained after the 175-day culture of explants through callus induction, embryogenic callus induction, somatic embryo development, and plant regeneration. Compared with the secondary somatic embryo seedling control, axillary bud regeneration plants (ABRPs) were normal diploid plants at the cellular and molecular level, with a variation rate of 7.74%.

Peptidomic Analysis and Antimicrobial Activity of Serum Peptide from Hevea brasiliensis clone BPM24

BACKGROUND

Hevea brasiliensis is severely affected by the fungal disease caused by Phytophthora spp. Significant loss of rubber yield is widespread and extensive use of chemical fungicides has resulted in health and environmental problems.

OBJECTIVE

This work aims to extract and identify the latex serum peptides from a disease tolerant clone of Hevea brasiliensis, and study the inhibitory efficacy against pathogenic bacteria and fungi.

METHODS

Serum peptides were extracted from H. brasiliensis BPM 24 using mixed lysis solution. Low molecular weight peptides were screened and fractionated by solid-phase extraction and then identified by tandem mass spectrometry. Total and fractionated serum peptides were assayed for bacterial and fungal inhibition using broth microdilution and poisoned food methods. An inhibitory control study in the greenhouse was also performed using susceptible clones for pre- and post-infection with Phytophthora spp.

RESULTS

Forty-three serum peptide sequences were successfully identified. Thirty-four peptides matched with the proteins associated with plant defense response signaling, host resistance, and adverse environmental factors. The inhibitory study of total serum peptides demonstrated anti-bacterial and anti-fungal properties. The greenhouse study exhibited disease inhibitory efficacy of 60% for the treatment of Phytophthora spp. in post-infected plants and 80% for pre-treated samples.

CONCLUSION

Latex serum peptides from disease tolerant H. brasiliensis revealed several proteins and peptides associated with plant defense and disease resistance. The peptides play a vital role for defense against bacteria and fungi pathogens, including Phytophthora spp. Enhanced disease protection can be obtained when the extracted peptides were applied to the susceptible plants before exposure to the fungi. These findings provided an insight and may pave the way for the development of biocontrol peptides from natural resources.

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 kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Study on in vitro NR biosynthesis by rapid quantitative determination of substrate depletion

A convenient and nonradioactive method for quantifying in vitro NR biosynthesis is presented that is based upon the quantitation of substrate depletion by ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS). NR oligomers could be in vitro biosynthesized with the enzyme source from Hevea brasiliensis (Hevea) or Taraxacum kok-saghyz (TKS) by exogenous monomers (IPP) and initiators (FPP). The IPP incorporation rate and FPP consumption rate were 62.24% and 51.14% respectively when the washed rubber particles (WRP) of Hevea was the enzyme source. The IPP incorporation rate and the FPP consumption rate were 74.49% and 95.90% respectively when the sediment bottom fraction (BF) of Hevea was the enzyme source. The in vitro NR biosynthesis can be divided into two stages:(1) the initiation reaction of FPP, which occurs more in BF, and (2) the growth reaction of IPP, which occurs more in WRP. In addition, the IPP incorporation and FPP consumption rates were 59.39% and 34.15% respectively when the BF of TKS was selected as an enzyme source.

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.

Rubber Genotypes with Contrasting Drought Factor Index Revealed Different Mechanisms for Drought Resistance in Hevea brasiliensis

It is predicted that drought will be more frequent and sustained in the future, which may affect the decline of rubber tree production. Therefore, it is critical to research some of the variables related to the drought-resistance mechanism of the rubber tree. As a result, it can be used to guide the selection of new rubber drought-resistance clones. The goal of this study was to identify drought-resistance mechanisms in rubber clones from the high drought factor index (DFI) group using ecophysiological and biochemical variables. The treatments consist of two factors, namely water deficit and contrasting clones based on the DFI variable. The first factor consisted of three levels, namely normal (fraction of transpirable soil water (FTSW) > 0.75), severe water deficit (0.1 < FTSW < 0.20), and recovery condition (FTSW > 0.75 after rewatering). The second factor consisted of seven clones, namely clones G239, GT1 (low DFI), G127, SP 217, PB 260 (moderate DFI), as well as G206 and RRIM 600 (high DFI). RRIM 600 had the highest DFI among the other clones as a drought-tolerance mechanism characteristic. Furthermore, clones RRIM 600, GT1, and G127 had lower stomatal conductance and transpiration rate than drought-sensitive clone PB 260. As a result, as drought avoidance mechanisms, clones RRIM 600, GT1, and G127 consume less water than clone PB 260. These findings indicated that clone RRIM 600 was a drought-resistant clone with drought tolerance and avoidance mechanisms.

Analysis of Genetic Diversity and Resistance to Foliar Pathogens Based on Genotyping-by-Sequencing of a Para Rubber Diversity Panel and Progeny of an Interspecific Cross

Para rubber trees (Hevea brasiliensis) are the largest major source of natural rubber in the world. Its major pathogens are Phytophthora spp., Corynespora cassiicola, and Colletotrichum spp. A rubber diversity panel of 116 clones using over 12,000 single nucleotide polymorphisms (SNPs) from DArTSeq genotyping revealed clear phylogenetic differences in clones that originated from different geographical regions of the world. An integrated linkage map constructed with an F₁ progeny of 86 from an interspecific cross between H. brasiliensis and H. benthamiana using 23,978 markers [10,323 SNPs and 13,655 SilicoDArTs] spanned 3947.83 cM with 0.83 cM average marker-interval. The genome scaffolds that were anchored to the linkage map, covering 1.44 Gb of H. brasiliensis reference genome, revealed a high level of collinearity between the genetic map and reference genome. Association analysis identified 12 SNPs significantly associated with the resistance against Phytophthora, Corynespora, and Colletotrichum in six linkage groups: 2, 6, 12, 14, 17, and 18. Kompetitive Allele-Specific PCR marker assays were developed for those 12 SNPs, screened with 178 individuals, and detected clear separation between two genotypes. Within the proximity to those SNPs, 41 potentially key genes that have previously been reported to associate with plant disease resistance were predicted with high confidence.

Cloning, Expression Analysis, and Functional Characterization of Candidate Oxalate Transporter Genes of HbOT1 and HbOT2 from Rubber Tree (Hevea brasiliensis)

Secretion of oxalic acid from roots is an important aluminum detoxification mechanism for many plants such as Hevea brasiliensis (rubber tree). However, the underlying molecular mechanism and oxalate transporter genes in plants have not yet been reported. In this study, the oxalate transporter candidate genes HbOT1 and HbOT2 from the rubber tree were cloned and preliminarily identified. It was found that HbOT1 had a full length of 1163 bp with CDS size of 792 bp, encoding 263 amino acids, and HbOT2 had a full length of 1647 bp with a CDS region length of 840 bp, encoding 279 amino acid residues. HbOT1 and HbOT2 were both stable hydrophobic proteins with transmembrane structure and SNARE_assoc domains, possibly belonging to the SNARE_assoc subfamily proteins of the SNARE superfamily. qRT-PCR assays revealed that HbOT1 and HbOT2 were constitutively expressed in different tissues, with HbOT1 highly expressed in roots, stems, barks, and latex, while HbOT2 was highly expressed in latex. In addition, the expressions of HbOT1 and HbOT2 were up-regulated in response to aluminum stress, and they were inducible by metals, such as copper and manganese. Heterologous expression of HbOT1 and HbOT2 in the yeast mutant AD12345678 enhanced the tolerance to oxalic acid and high concentration aluminum stress, which was closely correlated with the secretion of oxalic acid. This study is the first report on oxalate transporter genes in plants, which provides a theoretical reference for the study on the molecular mechanism of oxalic acid secretion to relieve aluminum toxicity and on aluminum-tolerance genetic engineering breeding.

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.

Genome-wide identification of the rubber tree superoxide dismutase (SOD) gene family and analysis of its expression under abiotic stress

Background

The rubber tree (Hevea brasiliensis) is the only species capable of producing high-quality natural rubber for commercial use, and is often subjected to various abiotic stresses in non-traditional rubber plantation areas. Superoxide dismutase (SOD) is a vital metalloenzyme translated by a SOD gene family member and acts as a first-line of protection in plant cells by catalysing the disproportionation of reactive oxygen species (ROS) to produce H₂O₂ and O₂. However, the SOD gene family is not reported in rubber trees.

Methods

Here, we used hidden markov model (HMM) and BLASTP methods to identify SOD genes in the H. brasiliensis genome. Phylogenetic tree, conserved motifs, gene structures, cis elements, and gene ontology annotation (GO) analyses were performed using MEGA 6.0, MEME, TBtools, PlantCARE, and eggNOG database, respectively. HbSOD gene expression profiles were analysed using quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results

We identified nine HbSOD genes in the rubber tree genome, including five HbCSDs, two HbFSDs, and two HbMSDs. Phylogenetic relationship analysis classified the SOD proteins from the rubber tree and other related species into three subfamilies. The results of gene structure and conserved motif analysis illustrated that most HbSOD genes have similar exon-intron numbers and conserved motifs in the same evolutionary branch. Five hormone-related, four stress-related, and light-responsive elements were detected in the HbSODs' promoters. HbSODs were expressed in different tissues, gradually increased with leaf development, and were abundantly expressed in mature leaves. HbCSD2 and HbCSD4 was significantly upregulated under low and high temperatures, and salt stress, except for HbCSD2, by heat. Furthermore, most HbSOD genes were significantly upregulated by drought, except HbMSD2. These findings imply that these genes may play vital roles in rubber tree stress resistance. Our results provide a basis for further studies on the functions of HbSOD genes in rubber trees and stress response mechanisms.

Insights to the superoxide dismutase genes and its roles in Hevea brasiliensis under abiotic stress

The superoxide dismutase (SOD) protein significantly influences the development and growth of plants and their reaction to abiotic stresses. However, little is known about the characteristics of rubber tree SOD genes and their expression changes under abiotic stresses. The present study recognized 11 SOD genes in the rubber tree genome, including 7 Cu/ZnSODs, 2 MnSODs, and 2 FeSODs. Except for HbFSD1, SODs were scattered on five chromosomes. The phylogenetic analysis of SOD proteins in rubber trees and a few other plants demonstrated that the SOD proteins contained three major subgroups. Moreover, the genes belonging to the same clade contained similar gene structures, which confirmed their classification further. The extension of the SOD gene family in the rubber tree was mainly induced by the segmental duplication events. The cis-acting components analysis showed that HbSODs were utilized in many biological procedures. The transcriptomics data indicated that the phosphorylation of the C-terminal domain of RNA polymerase II might control the cold response genes through the CBF pathway and activate the SOD system to respond to cold stress. The qRT-PCR results showed that the expression of HbCSD1 was significantly downregulated under drought and salt stresses, which might dominate the adaption capability to different stresses. Additionally, salt promoted the expression levels of HbMSD1 and HbMSD2, exhibiting their indispensable role in the salinity reaction. The study results will provide a theoretical basis for deep research on HbSODs in rubber trees.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03328-7.

Simultaneous quantification of volatile fatty acids and non-volatile organic acids in Hevea brasiliensis latex

The current method used in latex industries to determine the volatile fatty acids (VFAs) contents of Hevea brasiliensis latex is steam distillation. However, the accuracy of the method has been debated for some time. We assessed the accuracy of the method and developed a new, more reliable high performance liquid chromatographic method of determining acids in latex. The VFAs (formic, acetic, propionic, butyric and valeric acids) and non-volatile organic acids (Oxalic, malic, lactic, citric and succinic acids) in latex are directly determined simultaneously for the first time with high sensitivity and without losses during sample preparation. To avoid errors from derivatization, an acid resistant Prevail^TM HPLC column and a gradient mobile phase of 25 mM potassium dihydrogen phosphate (pH 2.5) and acetonitrile were employed. Under optimum conditions, the calibrations of both types of acids demonstrated satisfactory correlation coefficients of  ≥0.990, with limits of detection ranging from 0.02-395 mM. The developed method demonstrated the profiles of acids in field and concentrated latex of the same batch. Moreover, the evolution of the profiles of all studied acids in both types of latex during a three-month period was also revealed. This article is protected by copyright. All rights reserved.

HbTGA1, a TGA Transcription Factor From Hevea brasiliensis, Regulates the Expression of Multiple Natural Rubber Biosynthesis Genes

The TGA transcription factors are known to modulate the biosynthesis of secondary metabolites in plants. However, their regulatory function in natural rubber (NR) biosynthesis was not revealed in the rubber tree (Hevea brasiliensis). Here, 14 genes encoding TGA transcription factors (name HbTGA1-HbTGA14) were identified in the rubber tree. HbTGAs were differentially expressed in different tissues. HbTGA1 was expressed at its highest level in latex. We found specific in vitro and in vivo binding of the HbTGA1 protein with promoters of multiple NR biosynthesis genes (HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2). The activation of the promoters of HbHMGS2 and HbCPT6 was significantly suppressed by HbTGA1, while the activities of promoters of HbHMGR2, HbCPT8, and HbSRPP2 were increased by HbTGA1. The promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were significantly increased by HbTGA1 under jasmonate stress, while the promoter activities of HbHMGS2, HbHMGR2, HbCPT6, HbCPT8, and HbSRPP2 were also significantly increased by HbTGA1 under salicylic acid stress. The present study provides insights into the role of TGA transcription factors in regulating the expression of NR biosynthesis genes from H. brasiliensis.

Natural latex serum: characterization and biocompatibility assessment using Galleria mellonella as an alternative in vivo model

Natural latex serum (NLS) is one of the natural rubber latex fractions from Hevea brasiliensis tree, which is formed by centrifuged serum and is composed of proteins, acids, nucleotides, salts and carbohydrates. The proteins present in NLS have demonstrated several interesting biological properties, including angiogenic, healing, osteogenic, anti-inflammatory, antimicrobial, in addition to inducing neovascularization, bone formation and osseointegration. Thus, we proposed to characterize NLS by physicochemical techniques and to investigate the biocompatibility by toxicological assays and safety test in Galleria mellonella. Infrared spectrum showed vibrational bands characteristic of amide I, II and III that are linked to the protein content, which was confirmed by the High Performance Liquid Chromatography profile and by the Electrophoresis analysis. This material did not exhibit hemolytic (rate <0.5%) and cytotoxic effects (viability >70%) and was able to enhance the proliferation of fibroblasts (>600%) after 3 days. The pronounced proliferative effect observed in fibroblast cells can be explained by the presence of the fibroblast growth factor (FGF) like protein revealed by the Western blot test. Moreover, NLS did not provoke toxic effects (survival ∼ 80%) on the G. mellonella model, indicating that it is a biocompatible and safe material.

Reprogramming of the Hevea brasiliensis Epigenome and Transcriptome in Response to Cold Stress

Low temperature is a key factor limiting the rubber plantation extending to high latitude area. Previous work has shown that cold-induced DNA demethylation was coordinated with the expression of cold-responsive (COR) genes in Hevea brasiliensis. In this work, reduced representation bisulphite sequencing analysis of H. brasiliensis showed that cold treatment induced global genomic DNA demethylation and altered the sequence contexts of methylated cytosines, but the levels of mCG methylation in transposable elements were slightly enhanced by cold treatment. Integrated analysis of the DNA methylome and transcriptome revealed 400 genes whose expression correlated with altered DNA methylation. DNA demethylation in the upstream region of gene seems to correlate with higher gene expression, whereas demethylation in the gene body has less association. Our results suggest that cold treatment globally change the genomic DNA methylation status of the rubber tree, which might coordinate reprogramming of the transcriptome.

Based on Candidate Genes from Biosynthetic Pathways Associated with Latex Production

Scarcity of functional genetic markers associated with candidate genes (CGs) is a serious constraint for marker-assisted selection in the natural rubber producing tree, Hevea brasiliensis. In order to develop markers associated with rubber yield, five CGs involved in latex biosynthesis were characterized from 16 popular Hevea varieties. Novel SNPs and indels were identified and developed into markers using simple genotyping techniques like allele-specific PCR, CAPS, etc. A progeny population was genotyped using these markers to validate them, to understand their segregation pattern and to map them to a genetic linkage map. Parent-specific maps were constructed using pseudo-test cross strategy with the help of additional markers. The sequence structure information generated will be useful for future studies on gene mapping, functional relevance of coding SNPs and evolution of rubber biosynthesis genes in Hevea. Concurrently, the markers developed may serve as powerful tools for yield-based selection and for genetic diversity and pedigree studies in Hevea. Above all, the marker assays designed for genotyping could be economically carried out in any laboratory having basic molecular biology infrastructure and expertise.

Influence of Particle Size and Zeta Potential in Treating Highly Coloured Old Landfill Leachate by Tin Tetrachloride and Rubber Seed

Old leachate normally has a low organic compound content, poor biodegradability and is hard to biologically treat. The efficacy of tetravalent metal salts as a coagulant and the application of a natural coagulant as a flocculant in landfill leachate treatment is still inconclusive. Hence, this study aimed to evaluate the potential application of tin tetrachloride (SnCl₄) as the main coagulant and the rubber seed (Hevea brasiliensis) (RS) as the natural coagulant aid as the sole treatment in eradicating highly coloured and turbid stabilised landfill leachate present at one of the old local landfills in Malaysia. The standard jar test conducted revealed that SnCl₄ was able to eliminate 99% and 97.3% of suspended solids (SS) and colour, respectively, at pH8, with 10,000 mg/L dosages, an average particle size of 2419 d·nm, and a zeta potential (ZP) of −0.4 mV. However, RS was found to be ineffective as the main coagulant and could only remove 46.7% of SS and 76.5% of colour at pH3 with 6000 mg/L dosages, and also exhibited smaller particles (933 d·nm) with ZP values of −6.3 mV. When used as a coagulant aid, the polymer bridging mechanism in RS helped in reducing the SnCl₄ concentration from 10,000 mg/L to 8000 mg/L by maintaining the same performances. The presence of 1000 mg/L RS as a coagulant aid was able to remove 100% of SS and 97.6% of colour. The study concluded that RS has the potential to be used together with SnCl₄ in treating concentrated leachate with SS and colour.

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.

Differential Analysis of Mycelial Proteins and Metabolites From Rigidoporus Microporus During In Vitro Interaction With Hevea Brasiliensis

Rigidoporus microporus is the fungus accountable for the white root rot disease that is detrimental to the rubber tree, Hevea brasiliensis. The pathogenicity mechanism of R. microporus and the identity of the fungal proteins and metabolites involved during the infection process remain unclear. In this study, the protein and metabolite profiles of two R. microporus isolates, Segamat (SEG) and Ayer Molek (AM), were investigated during an in vitro interaction with H. brasiliensis. The isolates were used to inoculate H. brasiliensis clone RRIM 2025, and mycelia adhering to the roots of the plant were collected for analysis. Transmission electron microscope (TEM) images acquired confirms the hyphae attachment and colonization of the mycelia on the root of the H. brasiliensis clones after 4 days of inoculation. The protein samples were subjected to 2-DE analysis and analyzed using MALDI-ToF MS/MS, while the metabolites were extracted using methanol and analyzed using LC/MS-QTOF. Based on the differential analyses, upregulation of proteins that are essential for fungal evolution such as malate dehydrogenase, fructose 1,6-biphosphate aldolase, and glyceraldehyde-3-phosphate dehydrogenase hints an indirect role in fungal pathogenicity, while metabolomic analysis suggests an increase in acidic compounds which may lead to increased cell wall degrading enzyme activity. Bioinformatics analyses revealed that the carbohydrate and amino acid metabolisms were prominently affected in response to the fungal pathogenicity. In addition to that, other pathways that were significantly affected include "Protein Ubiquitination Pathway," Unfolded Protein Response," "HIFα Signaling," and "Sirtuin Signaling Pathway." The identification of responsive proteins and metabolites from this study promotes a better understanding of mechanisms underlying R. microporus pathogenesis and provides a list of potential biological markers for early recognition of the white root rot disease.

Recent insights on gene expression studies on Hevea Brasiliensis fatal leaf fall diseases

Hevea brasiliensis is one of the most important agricultural commodities globally, heavily cultivated in Southeast Asia. Fatal leaf fall diseases cause aggressive leaf defoliation, linked to lower latex yield and death of crops before maturity. Due to the significant consequences of the disease to H. brasiliensis, the recent gene expression studies from four fall leaf diseases of H. brasiliensis were gathered; South American leaf blight, powdery mildew, Corynespora cassiicola and Phytophthora leaf fall disease. The differential analysis observed the pattern of commonly expressed genes upon fungi triggers using RT-PCR, DDRT-PCR, Real-time qRT-PCR and RNA-Seq. We have observed that RNA-Seq is the best tool to seek novel genes. Among the identified genes with defence-against fungi were pathogenesis-related genes such as β-1,3-glucanase and chitinase, the reactive oxygen species, and the phytoalexin biosynthesis. This manuscript also provided functional elaboration on the responsive genes and predicted possible biosynthetic pathways to identify and characterise novel genes in the future. At the end of the manuscript, the PCR methods and proteomic approaches were presented for future molecular and biochemical studies in the related diseases to H. brasiliensis.

Prediction of South American Leaf Blight and Disease-Induced Photosynthetic Changes in Rubber Tree, Using Machine Learning Techniques on Leaf Hyperspectral Reflectance

The efficiency of visible and near-infrared (VIS/NIR) sensors and predictive modeling for detecting and classifying South American Leaf Blight (SALB) (Pseudocercospora ulei) in rubber trees (Hevea brasiliensis) has been poorly explored. Furthermore, the performance of VIS/NIR analysis combined with machine learning (ML) algorithms for predicting photosynthetic alterations caused by SALB is unknown. Therefore, this study aimed to detect and classify the SALB levels, as well as to predict, for the first time, disease-induced photosynthetic changes in rubber trees. Leaf hyperspectral reflectance combined with five ML techniques (random forest (RF), boosted regression tree (BRT), bagged classification and regression trees (BCART), artificial neural network (ANN), and support vector machine (SVM)) were used. The RF, ANN, and BCART models achieved the best performance for classifying the SALB levels on the training dataset (accuracies of 98.0 to 99.8%), with 10-fold cross-validation repeated five times, and test dataset (accuracies of 97.1 to 100%). The ANN and RF models were better at predicting leaf gas exchange-related traits such as net CO₂ assimilation rate (A) and extrinsic water use efficiency (WUEe) in the training (R² ranged from 0.97 to 0.99) and testing (R² ranged from 0.96 to 0.99) phases. In comparison, lower performances (R² ranged from 0.24 to 0.52) were evidenced for the photochemical traits. This research provides a basis for future designs of a remote monitoring system based on early detection and accurate diagnosis of biotic stress caused by SALB, which is fundamental for more effective rubber crop protection.

Differential Expression of lncRNAs and miRNAs Between Self-Rooting Juvenile and Donor Clones Unveils Novel Insight Into the Molecular Regulation of Rubber Biosynthesis in Hevea brasiliensis

The rubber tree (Hevea brasiliensis Muell. Arg.) is a tropical tree species that produce natural rubber. Self-rooted juvenile clones (SRJCs) are novel rubber tree planting materials developed through primary somatic embryogenesis. SRJCs have a higher rubber yield compared with donor clones (DCs). The molecular basis underlying increased rubber yield in SRJCs remains largely unknown. Here, the latex from SRJCs and DCs were collected for strand-specific and small RNA-seq methods. A total of 196 differentially expressed long noncoding RNAs (DELs), and 11 differentially expressed microRNAs were identified in latex between SRJCs and DCs. Targeted genes of DELs were markedly enriched for various biological pathways related to plant hormone signal transduction, photosynthesis, glutathione metabolism, and amino acids biosynthesis. DELs probably acted as cis-acting regulation was calculated, and these DELs relevant to potentially regulate rubber biosynthesis, reactive oxygen species metabolism, and epigenetic modification. Furthermore, the DELs acting as microRNA targets were studied. The interaction of microRNA and DELs might involve in the regulation of natural rubber biosynthesis.

Estimation of Heterografting Associated DNA Methylation Changes in Tree Crops by MSAP Analysis

The genetic incompatibility of the seedlings which are used as rootstocks (stock-scion interactions) and the mechanical stress induced by grafting are two major factors responsible for the high intraclonal variations observed in tree crops which are propagated through bud grafting. Since stress-induced DNA methylation changes associated with heterografting is a major contributor of such variations in grafted tree crops, a proper assessment of this epigenetic phenomenon is inevitable to devise strategies for the development of more uniform planting materials with minimal intraclonal variations in the future. In order to evaluate and establish the effects of heterografting on the epigenome of plants, availability of ideal plant materials and a standard procedure for testing is very essential. Development of genetically uniform own-rooted seedlings through induction of cleavage polyembryony by a novel technique of half ovulo embryo culture is the first step. Grafting of buds from these genetically and epigenetically uniform plants to genetically divergent rootstock and identification of DNA methylation polymorphism among them forms the second part of the methodology for detecting epigenetic changes associated with grafting in tree crops. Methylation-sensitive amplification polymorphism technique (MSAP), a modified version of AFLP using a pair of methylation-sensitive and insensitive isoschizomers (such as HpaII and MspI), is an ideal methodology to assess DNA methylation polymorphisms on a genomic scale in such plants. Comparative analysis of two sets of restriction digestion products (EcoRI/HpaII and EcoRI/MspI) allows the identification of DNA methylation polymorphisms induced by grafting and will aid in the detection of differentially methylated regions (DMRs) among grafted plants. This chapter describes a detailed protocol for inducing multiple embryos of single zygotic origin and regeneration of seedlings in rubber tree (Hevea brasiliensis), grafting of buds from these genetically uniform own-rooted seedlings to divergent rootstocks, identification of epigenetic changes induced by grafting or stock-scion interactions through MSAP analysis, and locating the differentially methylated genomic region. The methodology described here could be applied to any tree species commercially propagated through grafting for detecting epigenetic changes putatively associated with intraclonal variability.

Unravelling Rubber Tree Growth by Integrating GWAS and Biological Network-Based Approaches

Hevea brasiliensis (rubber tree) is a large tree species of the Euphorbiaceae family with inestimable economic importance. Rubber tree breeding programs currently aim to improve growth and production, and the use of early genotype selection technologies can accelerate such processes, mainly with the incorporation of genomic tools, such as marker-assisted selection (MAS). However, few quantitative trait loci (QTLs) have been used successfully in MAS for complex characteristics. Recent research shows the efficiency of genome-wide association studies (GWAS) for locating QTL regions in different populations. In this way, the integration of GWAS, RNA-sequencing (RNA-Seq) methodologies, coexpression networks and enzyme networks can provide a better understanding of the molecular relationships involved in the definition of the phenotypes of interest, supplying research support for the development of appropriate genomic based strategies for breeding. In this context, this work presents the potential of using combined multiomics to decipher the mechanisms of genotype and phenotype associations involved in the growth of rubber trees. Using GWAS from a genotyping-by-sequencing (GBS) Hevea population, we were able to identify molecular markers in QTL regions with a main effect on rubber tree plant growth under constant water stress. The underlying genes were evaluated and incorporated into a gene coexpression network modelled with an assembled RNA-Seq-based transcriptome of the species, where novel gene relationships were estimated and evaluated through in silico methodologies, including an estimated enzymatic network. From all these analyses, we were able to estimate not only the main genes involved in defining the phenotype but also the interactions between a core of genes related to rubber tree growth at the transcriptional and translational levels. This work was the first to integrate multiomics analysis into the in-depth investigation of rubber tree plant growth, producing useful data for future genetic studies in the species and enhancing the efficiency of the species improvement programs.

HbSnRK2.6 Functions in ABA-Regulated Cold Stress Response by Promoting HbICE2 Transcriptional Activity in Hevea brasiliensis

Low temperature remarkably limits rubber tree (Hevea brasiliensis Muell. Arg.) growth, latex production, and geographical distribution, but the underlying mechanisms of Hevea brasiliensis cold stress response remain elusive. Here, we identified HbSnRK2.6 as a key component in ABA signaling functions in phytohormone abscisic acid (ABA)-regulated cold stress response in Hevea brasiliensis. Exogenous application of ABA enhances Hevea brasiliensis cold tolerance. Cold-regulated (COR) genes in the CBF pathway are upregulated by ABA. Transcript levels of all five HbSnRK2.6 members are significantly induced by cold, while HbSnRK2.6A, HbSnRK2.6B, and HbSnRK2.6C can be further activated by ABA under cold conditions. Additionally, HbSnRK2.6s are localized in the cytoplasm and nucleus, and can physically interact with HbICE2, a crucial positive regulator in the cold signaling pathway. Overexpression of HbSnRK2.6A or HbSnRK2.6B in Arabidopsis extensively enhances plant responses to ABA and expression of COR genes, leading to increased cold stress tolerance. Furthermore, HbSnRK2.6A and HbSnRK2.6B can promote transcriptional activity of HbICE2, thus, increasing the expression of HbCBF1. Taken together, we demonstrate that HbSnRK2.6s are involved in ABA-regulated cold stress response in Hevea brasiliensis by regulating transcriptional activity of HbICE2.

Genome-Wide Identification and Analysis of Small Nucleolar RNAs and Their Roles in Regulating Latex Regeneration in the Rubber Tree (Hevea brasiliensis)

Small nucleolar RNAs (snoRNAs) are a class of conserved nuclear RNAs that play important roles in the modification of ribosomal RNAs (rRNAs) in plants. In rubber trees, rRNAs are run off with latex flow during tapping and need to be regenerated for maintaining the functions of the laticifer cells. SnoRNAs are expected to play essential roles in the regeneration of rRNAs. However, snoRNAs in the rubber tree have not been sufficiently characterized thus far. In this study, we performed nuclear RNA sequencing (RNA-seq) to identify snoRNAs globally and investigate their roles in latex regeneration. We identified a total of 3,626 snoRNAs by computational prediction with nuclear RNA-seq data. Among these snoRNAs, 50 were highly expressed in latex; furthermore, the results of reverse transcription polymerase chain reaction (RT-PCR) showed the abundant expression of 31 of these snoRNAs in latex. The correlation between snoRNA expression and adjusted total solid content (TSC/C) identified 13 positively yield-correlated snoRNAs. To improve the understanding of latex regeneration in rubber trees, we developed a novel insulated tapping system (ITS), which only measures the latex regenerated in specific laticifers. Using this system, a laticifer-abundant snoRNA, HbsnoR28, was found to be highly correlated with latex regeneration. To the best of our knowledge, this is the first report to globally identify snoRNAs that might be involved in latex regeneration regulation and provide new clues for unraveling the mechanisms underlying the regulation of latex regeneration.

Genotypic variation in 9-Cis-Epoxycarotenoid Dioxygenase3 gene expression and abscisic acid accumulation in relation to drought tolerance of Hevea brasiliensis

Abscisic acid (ABA) is a stress-related plant hormone, which is reported to confer drought tolerance. A key enzyme in ABA biosynthesis is 9-cis-epoxycarotenoid dioxygenase. In this study, changes in morphological, physiological response, HbNCED3, and ABA accumulation of RRIM 623 and PB 5/51 rubber clones were observed at different time points of water deficit conditions (0, 3, 5, 7, and 9 days of withholding water). During water deficit, the relative water content (RWC), photosynthetic rate (Pn), and stomatal conductance (Gs) decreased, whereas the electro leakage (EL) increased. The magnitudes of the changes in these parameters were greater for PB 5/51 than for RRIM 623. Therefore, RRIM 623 was designated as representative of drought-tolerant clone and PB 5/51 as a drought-sensitive clone. The HbNCED3 transcription level of RRIM 623 showed lower expression compared with that of PB 5/51, which corresponded to the accumulation of ABA. RRIM 623 accumulated less ABA than PB 5/51. The ABA in RRIM 623 gradually increased, especially on the 7th day of withholding water, whereas that in PB 5/51 rapidly increased during the early periods of drought conditions. Additionally, the sensitivity of stomatal response to ABA showed that RRIM 623 had a higher sensitivity than PB 5/51. These results demonstrate that the drought-tolerant rubber clone, RRIM 623, was characterized by lower ABA accumulation during drought stress than the drought-sensitive clone, PB 5/51. The drought tolerance mechanism of the RRIM 623 might be associated with stomatal sensitivity to ABA accumulation under drought stress.

Changes in anti-nutrient, phytochemical, and micronutrient contents of different processed rubber (Hevea brasiliensis) seed meals

Rubber (Hevea brasiliensis) is a perennial plant crop grown in many parts of Africa, South East Asia, and South America, especially within the hot and humid climatic regions. Rubber seed, either as feed or food, is a useful raw material to produce edible oil and protein. Despite the huge quantity of rubber seeds produced in Nigeria and its potential as a protein source, rubber seeds still appear neglected and under-utilised as feed/food given its perception as inedible and toxic due to the high concentration of cyanogenic glycoside. Therefore, the quest for effective processing technique(s) that would enhance its food use application is very fitting. This current study was directed to determine the changes in anti-nutrient, phytochemical, and micronutrient contents of different processed rubber seed meals. Specifically, the rubber seeds underwent processing, which employed boiling and the combined action of boiling and fermentation methods that brought about three seed meal flour groups, i.e., raw (RRSM), boiled (BRSM), and fermented (FRSM) seed meals. These were subsequently analysed for anti-nutrient/phytochemical (oxalate, phytate, tannin, phenols, saponin, hydrogen cyanide (HCN), alkaloids, flavonoids, and trypsin inhibitors), and micronutrient (which involved minerals (magnesium, phosphorus, calcium, iron, zinc, potassium, sodium, manganese, lead, and selenium) and vitamin (vitamin B1, B2, B3, C, E, and beta carotene)) contents. The results showed that the processing methods used to achieve the RRSM, BRSM, and FRSM, reduced the anti-nutrients (phytate, tannin, and oxalate) below the acceptable limits, and the HCN below the toxic levels. Importantly, the processing methods herein have not yet succeeded in removing HCN in the (processed) rubber seed meals, but can be seen to be heading toward the right direction. The FRSM obtained significantly lower (p < 0.05) anti-nutrient/phytochemical, but significantly higher (p < 0.05) mineral contents, compared with the other groups (RRSM and BRSM), except for flavonoids that obtained a 30% increase over the BRSM. Some mineral and vitamin contents could be lost in the BRSM compared to the others (RRSM and FRSM) in this study. Additionally, the FRSM obtained higher vitamin contents, after those of RRSM. Overall, the combined action of boiling and fermentation should be recommended for the proper utilisation of rubber seed as food/feed.

Structural and Biochemical Aspects Related to Resistance and Susceptibility of Rubber Tree Clones to Anthracnose

The aim of the present study was to evaluate structural and biochemical aspects related to the interaction of resistant (RRIM 937, IAC 502 and 507) and susceptible (RRIM 600) rubber tree clones with C. tamarillo. For such analysis, ultrathin sections of the leaf limb were embedded in historesin and differently stained to verify structural alterations and presence of starch grains, arginine, lipids, tannins and lignins. The total proteins and activity of the enzymes peroxidase and (PAL) were quantified. Stomatal density was also analyzed under a scanning electron microscope. Data indicated alterations in the cell content of resistant clones inoculated with the pathogen, as well as greater lignin and lipid accumulation in these samples. For tannins, there was no difference between inoculated and non-inoculated clones. Arginine was found at greater quantities in IAC 502 and 507. Starch grains were not detected in any of the analyzed samples. Protein level and stomatal density were lower in resistant clones. Peroxidase activity was more expressive in resistant clones. PAL activity, there was no significant difference between clones. The lignin and lipids, total protein, peroxidase activity and stomatal density may be related to the resistance of rubber tree clones to anthracnose.

Transcriptomes analysis reveals novel insight into the molecular mechanisms of somatic embryogenesis in Hevea brasiliensis

BACKGROUND

Somatic embryogenesis (SE) is a promising technology for plant vegetative propagation, which has an important role in tree breeding. Though rubber tree (Hevea brasiliensis Muell. Arg.) SE has been founded, few late SE-related genes have been identified and the molecular regulation mechanisms of late SE are still not well understood.

RESULTS

In this study, the transcriptomes of embryogenic callus (EC), primary embryo (PE), cotyledonary embryo (CE), abnormal embryo (AE), mature cotyledonary embryo (MCE) and withered abnormal embryo (WAE) were analyzed. A total of 887,852,416 clean reads were generated, 85.92% of them were mapped to the rubber tree genome. The de novo assembly generated 36,937 unigenes. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of CE vs. AE and MCE vs. WAE, respectively. The specific common DEGs were mainly involved in the phytohormones signaling pathway, biosynthesis of phenylpropanoid and starch and sucrose metabolism. Among them, hormone signal transduction related genes were significantly enriched, especially the auxin signaling factors (AUX-like1, GH3.1, SAUR32-like, IAA9-like, IAA14-like, IAA27-like, IAA28-like and ARF5-like). The transcription factors including WRKY40, WRKY70, MYBS3-like, MYB1R1-like, AIL6 and bHLH93-like were characterized as molecular markers for rubber tree late SE. CML13, CML36, CAM-7, SERK1 and LEAD-29-like were also related to rubber tree late SE. In addition, histone modification had crucial roles during rubber tree late SE.

CONCLUSIONS

This study provides important information to elucidate the molecular regulation during rubber tree late SE.

Tobacco rattle virus-induced gene silencing in Hevea brasiliensis

Virus-induced gene silencing (VIGS) is a powerful gene-silencing tool that has been intensively applied in plants. To data, the application of VIGS in rubber tree has not yet been reported. In this study, we described the efficient gene silencing in rubber tree by VIGS. The gene encoding Hevea brasiliensis phytoene desaturase (HbPDS) was identified in rubber tree genome. Small interfering RNAs from HbPDS and the silencing gene fragment were predicted and a length of 399 bp was selected to be tested. We showed that the tobacco rattle virus (TRV)-VIGS could induce effective HbPDS silencing in rubber tree. This study was the first to report VIGS in rubber tree. The present TRV-VIGS method could be used to perform reverse genetic approaches to identify unknown gene functions and might be further applied to produce gene silenced rubber tree plants, to advance functional gene of rubber tree.