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.

Dynamic analysis of Tapping Panel Dryness in Hevea brasiliensis reveals new insights on this physiological syndrome affecting latex production.

Tapping Panel Dryness (TPD) is a physiological disorder affecting natural rubber production in Hevea brasiliensis. TPD is associated with clonal susceptibility and overexploitation of rubber trees. Most studies are based on a binary point view of the absence or presence of TPD. This study sets out to characterize the dynamic of the TPD onset through the monthly monitoring of the dry cut length. This reveals the presence of dry spots on the tapped panel of any trees. The frequency of these dry spots increases dramatically in trees developing high level of TPD. Brown bast is an irreversible form of TPD. Brown bast is correlated to a high level of dry cut length. Application of an intensive harvesting system induces early TPD occurrence, which facilitates the study of TPD. Among latex diagnosis parameters, only sucrose content is significantly associated with TPD. Other parameters are more prone to environmental effects and are not reliable as physiological markers. These findings explain the contradictory conclusions of some papers. This study suggests to use intensive harvesting system and monitor the dry cut length for genetic analysis of TPD.

Corrigendum to "Analysis of reduced and oxidized antioxidants in Hevea brasiliensis latex reveals new insights into the regulation of antioxidants in response to harvesting stress and tapping panel dryness" [Heliyon 8 (7) (June 2022) Article e09840].

[This corrects the article DOI: 10.1016/j.heliyon.2022.e09840.].

Analysis of reduced and oxidized antioxidants in Hevea brasiliensis latex reveals new insights into the regulation of antioxidants in response to harvesting stress and tapping panel dryness.

Latex diagnosis (LD) is applied to optimize the natural rubber production and prevent tapping panel dryness (TPD), a physiological syndrome affecting latex production in Hevea brasiliensis. The reduced thiol content (RSH) is one of the biochemical parameters associated with the risk of TPD. However, RSH is difficult to interpret because of the influence of the environment. In order to better understand the regulation of antioxidants and to better interpret RSH, a key parameter of LD, this study analysed in latex both oxidised and reduced forms of ascorbic acid (AsA) and glutathione, and their cofactors as well as other latex diagnosis parameters in response to harvesting stress (tapping and ethephon stimulation) and TPD occurrence. The content of antioxidants in latex had a high variability among five rubber clones. The concentration in AsA was about ten times higher than GSH in laticifer, GSH accounting for about 50% of RSH. For short-term harvesting stress, RSH increased with tapping frequency and ethephon stimulation. TPD is associated with high latex viscosity and bursting of lysosomal particles called lutoids, as well as for several rubber clones with lower RSH and GSH contents. These results suggest that a high level of RSH shows the capacity of laticifer metabolism to cope with harvesting stress, while a drop in RSH is the sign of long stress related to lower metabolic activity and TPD occurrence. RSH remains an essential physiological parameter to prevent TPD when associated with reference data under low and high harvesting stress. This study paves the way to understand the role of AsA and GSH, and carry out genetic studies of antioxidants.

Transcriptome profiling in susceptible and tolerant rubber tree clones in response to cassiicolin Cas1, a necrotrophic effector from Corynespora cassiicola.

Corynespora cassiicola, a fungal plant pathogen with a large host range, causes important damages in rubber tree (Hevea brasiliensis), in Asia and Africa. A small secreted protein named cassiicolin was previously identified as a necrotrophic effector required for the virulence of C. cassiicola in specific rubber tree clones. The objective of this study was to decipher the cassiicolin-mediated molecular mechanisms involved in this compatible interaction. We comparatively analyzed the RNA-Seq transcriptomic profiles of leaves treated or not with the purified cassiicolin Cas1, in two rubber clones: PB260 (susceptible) and RRIM600 (tolerant). The reads were mapped against a synthetic transcriptome composed of all available transcriptomic references from the two clones. Genes differentially expressed in response to cassiicolin Cas1 were identified, in each clone, at two different time-points. After de novo annotation of the synthetic transcriptome, we analyzed GO enrichment of the differentially expressed genes in order to elucidate the main functional pathways impacted by cassiicolin. Cassiicolin induced qualitatively similar transcriptional modifications in both the susceptible and the tolerant clones, with a strong negative impact on photosynthesis, and the activation of defense responses via redox signaling, production of pathogenesis-related protein, or activation of the secondary metabolism. In the tolerant clone, transcriptional reprogramming occurred earlier but remained moderate. By contrast, the susceptible clone displayed a late but huge transcriptional burst, characterized by massive induction of phosphorylation events and all the features of a hypersensitive response. These results confirm that cassiicolin Cas1 is a necrotrophic effector triggering a hypersensitive response in susceptible rubber clones, in agreement with the necrotrophic-effector-triggered susceptibility model.

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.

Structural and Functional Annotation of Transposable Elements Revealed a Potential Regulation of Genes Involved in Rubber Biosynthesis by TE-Derived siRNA Interference in Hevea brasiliensis.

The natural rubber biosynthetic pathway is well described in Hevea, although the final stages of rubber elongation are still poorly understood. Small Rubber Particle Proteins and Rubber Elongation Factors (SRPPs and REFs) are proteins with major function in rubber particle formation and stabilization. Their corresponding genes are clustered on a scaffold1222 of the reference genomic sequence of the Hevea brasiliensis genome. Apart from gene expression by transcriptomic analyses, to date, no deep analyses have been carried out for the genomic environment of SRPPs and REFs loci. By integrative analyses on transposable element annotation, small RNAs production and gene expression, we analysed their role in the control of the transcription of rubber biosynthetic genes. The first in-depth annotation of TEs (Transposable Elements) and their capacity to produce TE-derived siRNAs (small interfering RNAs) is presented, only possible in the Hevea brasiliensis clone PB 260 for which all data are available. We observed that 11% of genes are located near TEs and their presence may interfere in their transcription at both genetic and epigenetic level. We hypothesized that the genomic environment of rubber biosynthesis genes has been shaped by TE and TE-derived siRNAs with possible transcriptional interference on their gene expression. We discussed possible functionalization of TEs as enhancers and as donors of alternative transcription start sites in promoter sequences, possibly through the modelling of genetic and epigenetic landscapes.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis.

BACKGROUND: Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea. METHODS: Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20-22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators. RESULTS: Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Proteomic Landscape Has Revealed Small Rubber Particles Are Crucial Rubber Biosynthetic Machines for Ethylene-Stimulation in Natural Rubber Production.

Rubber particles are a specific organelle for natural rubber biosynthesis (NRB) and storage. Ethylene can significantly improve rubber latex production by increasing the generation of small rubber particles (SRPs), regulating protein accumulation, and activating many enzyme activities. We conducted a quantitative proteomics study of different SRPs upon ethylene stimulation by differential in-gel electrophoresis (DIGE) and using isobaric tags for relative and absolute quantification (iTRAQ) methods. In DIGE, 79 differentially accumulated proteins (DAPs) were determined as ethylene responsive proteins. Our results show that the abundance of many NRB-related proteins has been sharply induced upon ethylene stimulation. Among them, 23 proteins were identified as rubber elongation factor (REF) and small rubber particle protein (SRPP) family members, including 16 REF and 7 SRPP isoforms. Then, 138 unique phosphorylated peptides, containing 129 phosphorylated amino acids from the 64 REF/SRPP family members, were identified, and most serine and threonine were phosphorylated. Furthermore, we identified 226 DAPs from more than 2000 SRP proteins by iTRAQ. Integrative analysis revealed that almost all NRB-related proteins can be detected in SRPs, and many proteins are positively responsive to ethylene stimulation. These results indicate that ethylene may stimulate latex production by regulating the accumulation of some key proteins. The phosphorylation modification of REF and SRPP isoforms might be crucial for NRB, and SRP may act as a complex natural rubber biosynthetic machine.

Genome-wide analysis in Hevea brasiliensis laticifers revealed species-specific post-transcriptional regulations of several redox-related genes.

MicroRNA-mediated post-transcriptional regulation has been reported on ROS production and scavenging systems. Although microRNAs first appeared highly conserved among plant species, several aspects of biogenesis, function and evolution of microRNAs were shown to differ. High throughput transcriptome and degradome analyses enable to identify small RNAs and their mRNA targets. A non-photosynthetic tissue particularly prone to redox reactions, laticifers from Hevea brasiliensis, revealed species-specific post-transcriptional regulations. This paper sets out to identify the 407 genes of the thirty main redox-related gene families harboured by the Hevea genome. There are 161 redox-related genes expressed in latex. Thirteen of these redox-related genes were targeted by 11 microRNAs. To our knowledge, this is the first report on a mutation in the miR398 binding site of the cytosolic CuZnSOD. A working model was proposed for transcriptional and post-transcriptional regulation with respect to the predicted subcellular localization of deduced proteins.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness.

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

Overexpression of Hevea brasiliensis ethylene response factor HbERF-IXc5 enhances growth and tolerance to abiotic stress and affects laticifer differentiation.

Ethylene response factor 1 (ERF1) is an essential integrator of the jasmonate and ethylene signalling pathways coordinating a large number of genes involved in plant defences. Its orthologue in Hevea brasiliensis, HbERF-IXc5, has been assumed to play a major role in laticifer metabolism and tolerance to harvesting stress for better latex production. This study sets out to establish and characterize rubber transgenic lines overexpressing HbERF-IXc5. Overexpression of HbERF-IXc5 dramatically enhanced plant growth and enabled plants to maintain some ecophysiological parameters in response to abiotic stress such as water deficit, cold and salt treatments. This study revealed that HbERF-IXc5 has rubber-specific functions compared to Arabidopsis ERF1 as transgenic plants overexpressing HbERF-IXc5 accumulated more starch and differentiated more latex cells at the histological level. The role of HbERF-IXc5 in driving the expression of some target genes involved in laticifer differentiation is discussed.

Reactive oxygen species in Hevea brasiliensis latex and relevance to Tapping Panel Dryness.

Environmental stress can lead to oxidative stress resulting from an increase in reactive oxygen species (ROS) and involves redox adjustments. Natural rubber is synthesized in laticifers, which is a non-photosynthetic tissue particularly prone to oxidative stress. This paper reviews the current state of knowledge on the ROS production and ROS-scavenging systems in laticifers. These regulations have been the subject of intense research into a physiological syndrome, called Tapping Panel Dryness (TPD), affecting latex production in Hevea brasiliensis. In order to prevent TPD occurrence, monitoring thiol content appeared to be a crucial factor of latex diagnosis. Thiols, ascorbate and γ-tocotrienol are the major antioxidants in latex. They are involved in membrane protection from ROS and likely have an effect on the quality of raw rubber. Some transcription factors might play a role in the redox regulatory network in Hevea, in particular ethylene response factors, which have been the most intensively studied given the role of ethylene on rubber production. Current challenges for rubber research and development with regard to redox systems will involve improving antioxidant capacity using natural genetic variability.

Expression Profiles, Characterization and Function of HbTCTP in Rubber Tree (Hevea brasiliensis).

As a highly conserved protein, the translationally controlled tumor protein (TCTP) carries out vital roles in various life processes. In rubber tree, two TCTP genes, HbTCTP and HbTCTP1, were cloned, but only HbTCTP1 was studied in details. In this study, cis-acting regulatory elements, expression patterns, subcellular localization, interacting proteins, and antioxidant activity of HbTCTP were systematically analyzed. Besides the common cis-acting regulatory elements, HbTCTP promoter also harbored various known cis-elements that respond to hormone/stresses. Being consistent with the aforementioned results, HbTCTP was regulated by drought, low temperature, high salt, ethylene (ET), wounding, H2O2, and methyl jasmonate (MeJA) treatments. HbTCTP was expressed throughout different tissues and developmental stages of leaves. In addition, HbTCTP was associated with tapping panel dryness (TPD). HbTCTP was localized in the membrane, cytoplasm and the nucleus, and interacted with four proteins rubber elongation factor (REF), 17.5 kDa heat shock family protein, annexin, and REF-like stress related protein 1. Being similar to HbTCTP1, HbTCTP also indicated antioxidant activity in metal-catalyzed oxidation (MCO) system. Our results are useful for further understanding the molecular characterization and expression profiles of HbTCTP, but also lay a solid foundation for elucidating the function of HbTCTP in rubber tree.

The rubber tree genome reveals new insights into rubber production and species adaptation.

The Para rubber tree (Hevea brasiliensis) is an economically important tropical tree species that produces natural rubber, an essential industrial raw material. Here we present a high-quality genome assembly of this species (1.37 Gb, scaffold N50 = 1.28 Mb) that covers 93.8% of the genome (1.47 Gb) and harbours 43,792 predicted protein-coding genes. A striking expansion of the REF/SRPP (rubber elongation factor/small rubber particle protein) gene family and its divergence into several laticifer-specific isoforms seem crucial for rubber biosynthesis. The REF/SRPP family has isoforms with sizes similar to or larger than SRPP1 (204 amino acids) in 17 other plants examined, but no isoforms with similar sizes to REF1 (138 amino acids), the predominant molecular variant. A pivotal point in Hevea evolution was the emergence of REF1, which is located on the surface of large rubber particles that account for 93% of rubber in the latex (despite constituting only 6% of total rubber particles, large and small). The stringent control of ethylene synthesis under active ethylene signalling and response in laticifers resolves a longstanding mystery of ethylene stimulation in rubber production. Our study, which includes the re-sequencing of five other Hevea cultivars and extensive RNA-seq data, provides a valuable resource for functional genomics and tools for breeding elite Hevea cultivars.

Involvement of Ethylene in the Latex Metabolism and Tapping Panel Dryness of Hevea brasiliensis.

Ethephon, an ethylene releaser, is used to stimulate latex production in Hevea brasiliensis. Ethylene induces many functions in latex cells including the production of reactive oxygen species (ROS). The accumulation of ROS is responsible for the coagulation of rubber particles in latex cells, resulting in the partial or complete stoppage of latex flow. This study set out to assess biochemical and histological changes as well as changes in gene expression in latex and phloem tissues from trees grown under various harvesting systems. The Tapping Panel Dryness (TPD) susceptibility of Hevea clones was found to be related to some biochemical parameters, such as low sucrose and high inorganic phosphorus contents. A high tapping frequency and ethephon stimulation induced early TPD occurrence in a high latex metabolism clone and late occurrence in a low latex metabolism clone. TPD-affected trees had smaller number of laticifer vessels compared to healthy trees, suggesting a modification of cambial activity. The differential transcript abundance was observed for twenty-seven candidate genes related to TPD occurrence in latex and phloem tissues for ROS-scavenging, ethylene biosynthesis and signalling genes. The predicted function for some Ethylene Response Factor genes suggested that these candidate genes should play an important role in regulating susceptibility to TPD.

Development of a new pCAMBIA binary vector using Gateway® technology.

pCAMBIA vectors have become popular for their easy handling, stability and the existence of a range of selection and reporter genes. However, these vectors have yet to integrate the Gateway® cloning system, which has enabled site-specific recombination without the need for restriction enzymes and ligases. This paper sets out to convert the pCambia2300 binary vector into a destination vector with the Gateway® cassette driven by the CaMV35S promoter. The destination vector, pCamway35S, was then evaluated using the uidA reporter gene. Transient and stable transformation experiments were successfully assayed, either by particle bombardment or by Agrobacterium tumefaciens in Allium cepa and Hevea embryogenic calli. After counting the transformation units, the statistical analysis performed on the data showed that the pCamway 35S::uidA vector was as efficient as pCambia2301, a pCAMBIA2300 containing the uidA reporter gene under the CaMV 35S promoter.

Ethylene Response Factors Are Controlled by Multiple Harvesting Stresses in Hevea brasiliensis.

Tolerance of recurrent mechanical wounding and exogenous ethylene is a feature of the rubber tree. Latex harvesting involves tapping of the tree bark and ethephon is applied to increase latex flow. Ethylene is an essential element in controlling latex production. The ethylene signalling pathway leads to the activation of Ethylene Response Factor (ERF) transcription factors. This family has been identified in Hevea brasiliensis. This study set out to understand the regulation of ERF genes during latex harvesting in relation to abiotic stress and hormonal treatments. Analyses of the relative transcript abundance were carried out for 35 HbERF genes in latex, in bark from mature trees and in leaves from juvenile plants under multiple abiotic stresses. Twenty-one HbERF genes were regulated by harvesting stress in laticifers, revealing an overrepresentation of genes in group IX. Transcripts of three HbERF-IX genes from HbERF-IXc4, HbERF-IXc5 and HbERF-IXc6 were dramatically accumulated by combining wounding, methyl jasmonate and ethylene treatments. When an ethylene inhibitor was used, the transcript accumulation for these three genes was halted, showing ethylene-dependent induction. Subcellular localization and transactivation experiments confirmed that several members of HbERF-IX are activator-type transcription factors. This study suggested that latex harvesting induces mechanisms developed for the response to abiotic stress. These mechanisms probably depend on various hormonal signalling pathways. Several members of HbERF-IX could be essential integrators of complex hormonal signalling pathways in Hevea.

Transcriptional and post-transcriptional regulation of the jasmonate signalling pathway in response to abiotic and harvesting stress in Hevea brasiliensis.

BACKGROUND: Latex harvesting in Hevea brasiliensis amounts to strong abiotic stress that can cause a halt in production in the most susceptible clones. Although the role of jasmonic acid has been suggested in laticifer differentiation, its role in latex production and in the response to harvesting stress has received very little attention. Only a few key genes acting in the COI-JAZ-MYC module have been isolated and studied at transcriptional level. RESULTS: Use of a reference transcriptome obtained on rubber clone PB 260 covering a large number of tissues under different environmental conditions enabled us to identify 24 contigs implicated in the jasmonate signalling pathway in the rubber tree. An analysis of their expression profile by qPCR, combined with hierarchical clustering, suggested that the jasmonate signalling pathway is highly activated in laticifer cells and, more particularly, in the response to harvesting stress. By comparison with their genomic sequences, the existence of regulation by alternative splicing was discovered for JAZ transcripts in response to harvesting stress. Lastly, positive transcriptional regulation of the HbJAZ_1405 gene by MYC was demonstrated. CONCLUSION: This study led to the identification of all actors of jasmonate signalling pathway and revealed a specific gene expression pattern in latex cells. In-depth analysis of this regulation showed alternative splicing that has been previously shown in Arabidopsis. Interestingly, genotypic variation was observed in Hevea clones with contrasting latex metabolism. This result suggests an involvement of jasmonate signalling pathway in latex production. The data suggest that specific variability of the JA pathway may have some major consequences for resistance to stress. The data support the hypothesis that a better understanding of transcriptional regulations of jasmonate pathway during harvesting stress, along with the use of genotypic diversity in response to such stress, can be used to improve resistance to stress and rubber production in Hevea.

Sequence and expression analyses of ethylene response factors highly expressed in latex cells from Hevea brasiliensis.

The AP2/ERF superfamily encodes transcription factors that play a key role in plant development and responses to abiotic and biotic stress. In Hevea brasiliensis, ERF genes have been identified by RNA sequencing. This study set out to validate the number of HbERF genes, and identify ERF genes involved in the regulation of latex cell metabolism. A comprehensive Hevea transcriptome was improved using additional RNA reads from reproductive tissues. Newly assembled contigs were annotated in the Gene Ontology database and were assigned to 3 main categories. The AP2/ERF superfamily is the third most represented compared with other transcription factor families. A comparison with genomic scaffolds led to an estimation of 114 AP2/ERF genes and 1 soloist in Hevea brasiliensis. Based on a phylogenetic analysis, functions were predicted for 26 HbERF genes. A relative transcript abundance analysis was performed by real-time RT-PCR in various tissues. Transcripts of ERFs from group I and VIII were very abundant in all tissues while those of group VII were highly accumulated in latex cells. Seven of the thirty-five ERF expression marker genes were highly expressed in latex. Subcellular localization and transactivation analyses suggested that HbERF-VII candidate genes encoded functional transcription factors.