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.

BackgroundRubber tree (Hevea brasiliensis) laticifers are the source of natural rubber. Rubber production depends on endogenous and exogenous ethylene (ethephon). AP2/ERF transcription factors, and especially Ethylene-Response Factors, play a crucial role in plant development and response to biotic and abiotic stresses. This study set out to sequence transcript expressed in various tissues using next-generation sequencing and to identify AP2/ERF superfamily in the rubber tree.ResultsThe 454 sequencing technique was used to produce five tissue-type transcript libraries (leaf, bark, latex, embryogenic tissues and root). Reads from all libraries were pooled and reassembled to improve mRNA lengths and produce a global library. One hundred and seventy-three AP2/ERF contigs were identified by in silico analysis based on the amino acid sequence of the conserved AP2 domain from the global library. The 142 contigs with the full AP2 domain were classified into three main families (20 AP2 members, 115 ERF members divided into 11 groups, and 4 RAV members) and 3 soloist members. Fifty-nine AP2/ERF transcripts were found in latex. Alongside the microRNA172 already described in plants, eleven additional microRNAs were predicted to inhibit Hevea AP2/ERF transcripts.ConclusionsHevea has a similar number of AP2/ERF genes to that of other dicot species. We adapted the alignment and classification methods to data from next-generation sequencing techniques to provide reliable information. We observed several specific features for the ERF family. Three HbSoloist members form a group in Hevea. Several AP2/ERF genes highly expressed in latex suggest they have a specific function in Hevea. The analysis of AP2/ERF transcripts in Hevea presented here provides the basis for studying the molecular regulation of latex production in response to abiotic stresses and latex cell differentiation.

BackgroundEthylene production and signalling play an important role in somatic embryogenesis, especially for species that are recalcitrant in in vitro culture. The AP2/ERF superfamily has been identified and classified in Hevea brasiliensis. This superfamily includes the ERFs involved in response to ethylene. The relative transcript abundance of ethylene biosynthesis genes and of AP2/ERF genes was analysed during somatic embryogenesis for callus lines with different regeneration potential, in order to identify genes regulated during that process.ResultsThe analysis of relative transcript abundance was carried out by real-time RT-PCR for 142 genes. The transcripts of ERFs from group I, VII and VIII were abundant at all stages of the somatic embryogenesis process. Forty genetic expression markers for callus regeneration capacity were identified. Fourteen markers were found for proliferating calli and 35 markers for calli at the end of the embryogenesis induction phase. Sixteen markers discriminated between normal and abnormal embryos and, lastly, there were 36 markers of conversion into plantlets. A phylogenetic analysis comparing the sequences of the AP2 domains of Hevea and Arabidopsis genes enabled us to predict the function of 13 expression marker genes.ConclusionsThis first characterization of the AP2/ERF superfamily in Hevea revealed dramatic regulation of the expression of AP2/ERF genes during the somatic embryogenesis process. The gene expression markers of proliferating callus capacity to regenerate plants by somatic embryogenesis should make it possible to predict callus lines suitable to be used for multiplication. Further functional characterization of these markers opens up prospects for discovering specific AP2/ERF functions in the Hevea species for which somatic embryogenesis is difficult.