The Chromosome-Based Rubber Tree Genome Provides New Insights into Spurge Genome Evolution and Rubber Biosynthesis.

The rubber tree, Hevea brasiliensis, produces natural rubber that serves as an essential industrial raw material. Here, we present a high-quality reference genome for a rubber tree cultivar GT1 using single-molecule real-time sequencing (SMRT) and Hi-C technologies to anchor the ∼1.47-Gb genome assembly into 18 pseudochromosomes. The chromosome-based genome analysis enabled us to establish a model of spurge chromosome evolution, since the common paleopolyploid event occurred before the split of Hevea and Manihot. We show recent and rapid bursts of the three Hevea-specific LTR-retrotransposon families during the last 10 million years, leading to the massive expansion by ∼65.88% (∼970 Mbp) of the whole rubber tree genome since the divergence from Manihot. We identify large-scale expansion of genes associated with whole rubber biosynthesis processes, such as basal metabolic processes, ethylene biosynthesis, and the activation of polysaccharide and glycoprotein lectin, which are important properties for latex production. A map of genomic variation between the cultivated and wild rubber trees was obtained, which contains ∼15.7 million high-quality single-nucleotide polymorphisms. We identified hundreds of candidate domestication genes with drastically lowered genomic diversity in the cultivated but not wild rubber trees despite a relatively short domestication history of rubber tree, some of which are involved in rubber biosynthesis. This genome assembly represents key resources for future rubber tree research and breeding, providing novel targets for improving plant biotic and abiotic tolerance and rubber production.

Distinguishing CPT gene family members and vetting the sequence structure of a putative rubber synthesizing variant in Hevea brasiliensis.

cis-Prenyltransferases (cis-PTs) constitute a large family of enzymes conserved during evolution and present in all domains of life. cis-PTs catalyze the cis-1,4-polymerization of isoprene units to generate isoprenoids with carbon skeletons varying from C10 (neryl pyrophosphate) to C > 10,000 (natural rubber). Though the previously reported CPTs in Hevea are designated based on sequence variations, their classification was done mostly by phylogenetic analysis using a mixture of partial as well as full length sequences often excluding the UTRs. In this context an attempt was made to reclassify the CPTs strictly based on their sequence similarity and distinguish the members putatively associated with rubber biosynthesis from the others. Extensive computational analysis was carried out on CPT sequences obtained from public resources and whole genome assemblies of Hevea. Based on the results from BLAST analysis, multiple sequence alignments of protein, nucleotide and untranslated regions, open reading frame analysis, gene prediction analysis and sequence length variations, we conclude that there exists mainly three CPTs namely RubCPT1, RubCPT2 and RubCPT3 putatively associated with rubber biosynthesis in Hevea brasiliensis. The rest were categorised as variants of dehydrodolichyl diphosphate synthase (DHDDS) involved in the synthesis of dolichols having short chain isoprenoids. Analysis of the sequence structure of the most highly expressed RubCPT1 in latex revealed the allele richness and diversity of this important variant prevailing in the popular rubber clones. Haplotypes consisting of SNPs with high degree of heterozygosity were also identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor towards the generation of allelic diversity rather than point mutations. Alternatively, gene expression analysis indicated the possibility of association between specific haplotypes and RubCPT1 expression in Hevea clones which may have downstream impact up to the level of rubber production. The conclusions from this study may pave way for the identification and better understanding of CPTs directly involved with natural rubber biosynthesis in Hevea and the SNP data generated may aid in the development of molecular markers putatively associated with yield in rubber.

The 14-3-3 protein HbGF14a interacts with a RING zinc finger protein to regulate expression of the rubber transferase gene in Hevea brasiliensis.

Hevea brasiliensis is a key commercial source of natural rubber (cis 1,4-polyisoprene). In H. brasiliensis, rubber transferase is responsible for cis-1,4-polymerization of isoprene units from isopentenyl diphosphate and thus affects the yield of rubber. Little is known about the regulatory mechanisms of the rubber transferase gene at a molecular level. In this study we show that the 5'UTR intron of the promoter of the rubber transferase gene (HRT2) suppresses the expression of HRT2. A H. brasiliensis RING zinc finger protein (designated as HbRZFP1) was able to interact specifically with the HRT2 promoter to down-regulate its transcription in vivo. A 14-3-3 protein (named as HbGF14a) was identified as interacting with HbRZFP1, both in yeast and in planta. Transient co-expression of HbGF14a and HbRZFP1-encoding cDNAs resulted in HbRZFP1-mediated HRT2 transcription inhibition being relieved. HbGF14a repressed the protein-DNA binding of HbRZFP1 with the HRT2 promoter in yeast. We propose a regulatory mechanism by which the binding of HbGF14a to HbRZFP1 interferes with the interaction of HbRZFP1 with the HRT2 promoter, thereby repressing the protein-DNA binding between them. This study provides new insights into the role of HbGF14a in mediating expression of the rubber transferase gene in Hevea brasiliensis.

Construction of Pará rubber tree genome and multi-transcriptome database accelerates rubber researches.

BACKGROUND: Natural rubber is an economically important material. Currently the Pará rubber tree, Hevea brasiliensis is the main commercial source. Little is known about rubber biosynthesis at the molecular level. Next-generation sequencing (NGS) technologies brought draft genomes of three rubber cultivars and a variety of RNA sequencing (RNA-seq) data. However, no current genome or transcriptome databases (DB) are organized by gene. RESULTS: A gene-oriented database is a valuable support for rubber research. Based on our original draft genome sequence of H. brasiliensis RRIM600, we constructed a rubber tree genome and transcriptome DB. Our DB provides genome information including gene functional annotations and multi-transcriptome data of RNA-seq, full-length cDNAs including PacBio Isoform sequencing (Iso-Seq), ESTs and genome wide transcription start sites (TSSs) derived from CAGE technology. Using our original and publically available RNA-seq data, we calculated co-expressed genes for identifying functionally related gene sets and/or genes regulated by the same transcription factor (TF). Users can access multi-transcriptome data through both a gene-oriented web page and a genome browser. For the gene searching system, we provide keyword search, sequence homology search and gene expression search; users can also select their expression threshold easily. CONCLUSION: The rubber genome and transcriptome DB provides rubber tree genome sequence and multi-transcriptomics data. This DB is useful for comprehensive understanding of the rubber transcriptome. This will assist both industrial and academic researchers for rubber and economically important close relatives such as R. communis, M. esculenta and J. curcas. The Rubber Transcriptome DB release 2017.03 is accessible at http://matsui-lab.riken.jp/rubber/ .

Genome-wide association mapping of latex yield and girth in Amazonian accessions of Hevea brasiliensis grown in a suboptimal climate zone.

Latex yield and growth are the key complex traits in commercial rubber production. The present study is the first to report genome-wide association mapping of latex yield and girth, for 170 Amazonian accessions grown in a suboptimal area characterized by limited rainfall and a lengthy dry season. Targeted sequence enrichment to capture gene transcripts generated 14,155 high quality filtered single nucleotide polymorphisms (SNPs) of which 94.3% resided in coding regions. The rapid decay of linkage disequilibrium over physical and genetic distance found in the accessions was comparable to those previously reported for several outcrossing species. A mixed linear model detected three significant SNPs in three candidate genes involved in plant adaptation to drought stress, individually explaining 12.7-15.7% of the phenotypic variance. The SNPs identified in the study will help to extend understanding, and to support genetic improvement of rubber trees grown in drought-affected regions.

Identification and characterization of the abscisic acid (ABA) receptor gene family and its expression in response to hormones in the rubber tree.

Abscisic acid (ABA) is an essential phytohormone involved in diverse physiological processes. Although genome-wide analyses of the ABA receptor PYR/PYL/RCAR (PYL) protein/gene family have been performed in certain plant species, little is known about the ABA receptor protein/gene family in the rubber tree (Hevea brasiliensis). In this study, we identified 14 ABA receptor PYL proteins/genes (designated HbPYL1 through HbPYL14) in the most recent rubber tree genome. A phylogenetic tree was constructed, which demonstrated that HbPYLs can be divided into three subfamilies that correlate well with the corresponding Arabidopsis subfamilies. Eight HbPYLs are highly expressed in laticifers. Five of the eight genes are simultaneously regulated by ABA, jasmonic acid (JA) and ethylene (ET). The identification and characterization of HbPYLs should enable us to further understand the role of ABA signal in the rubber tree.

Survey of the rubber tree genome reveals a high number of cysteine protease-encoding genes homologous to Arabidopsis SAG12.

Arabidopsis thaliana SAG12, a senescence-specific gene encoding a cysteine protease, is widely used as a molecular marker for the study of leaf senescence. To date, its potential orthologues have been isolated from several plant species such as Brassica napus and Nicotiana tabacum. However, little information is available in rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study presents the identification of SAG12-like genes from the rubber tree genome. Results showed that an unexpected high number of 17 rubber orthologues with a single intron were found, contrasting the single copy with two introns in Arabidopsis. The gene expansion was also observed in another two Euphorbiaceae plants, castor bean (Ricinus communis) and physic nut (Jatropha curcas), both of which contain 8 orthologues. In accordance with no occurrence of recent whole-genome duplication (WGD) events, most duplicates in castor and physic nut were resulted from tandem duplications. In contrast, the duplicated HbSAG12H genes were derived from tandem duplications as well as the recent WGD. Expression analysis showed that most HbSAG12H genes were lowly expressed in examined tissues except for root and male flower. Furthermore, HbSAG12H1 exhibits a strictly senescence-associated expression pattern in rubber tree leaves, and thus can be used as a marker gene for the study of senescence mechanism in Hevea.

Regulation of HbPIP2;3, a Latex-Abundant Water Transporter, Is Associated with Latex Dilution and Yield in the Rubber Tree (Hevea brasiliensis Muell. Arg.).

Rubber tree (Hevea brasiliensis) latex, the source of natural rubber, is synthesised in the cytoplasm of laticifers. Efficient water inflow into laticifers is crucial for latex flow and production since it is the determinant of the total solid content of latex and its fluidity after tapping. As the mature laticifer vessel rings are devoid of plasmodesmata, water exchange between laticifers and surrounding cells is believed to be governed by plasma membrane intrinsic proteins (PIPs). To identify the most important PIP aquaporin in the water balance of laticifers, the transcriptional profiles of ten-latex-expressed PIPs were analysed. One of the most abundant transcripts, designated HbPIP2;3, was characterised in this study. When tested in Xenopus laevis oocytes HbPIP2;3 showed a high efficiency in increasing plasmalemma water conductance. Expression analysis indicated that the HbPIP2;3 gene was preferentially expressed in latex, and the transcripts were up-regulated by both wounding and exogenously applied Ethrel (a commonly-used ethylene releaser). Although regular tapping up-regulated the expression of HbPIP2;3 during the first few tappings of the virginal rubber trees, the transcriptional kinetics of HbPIP2;3 to Ethrel stimulation in the regularly tapped tree exhibited a similar pattern to that of the previously reported HbPIP2;1 in the virginal rubber trees. Furthermore, the mRNA level of HbPIP2;3 was associated with clonal yield potential and the Ethrel stimulation response. Together, these results have revealed the central regulatory role of HbPIP2;3 in laticifer water balance and ethylene stimulation of latex production in Hevea.

Cloning and characterization of HbMT2a, a metallothionein gene from Hevea brasiliensis Muell. Arg differently responds to abiotic stress and heavy metals.

Metallothioneins (MTs) are of low molecular mass, cysteine-rich proteins. They play an important role in the detoxification of heavy metals and homeostasis of intracellular metal ions, and protecting against intracellular oxidative damages. In this study a full-length cDNA of type 2 plant metallothioneins, HbMT2a, was isolated from 25 mM Polyethyleneglycol (PEG) stressed leaves of Hevea brasiliensis by RACE. The HbMT2a was 372bp in length and had a 237bp open reading frame (ORF) encoding for a protein of 78 amino acid residues with molecular mass of 7.772 kDa. The expression of HbMT2a in the detached leaves of rubber tree clone RY7-33-97 was up-regulated by Me-JA, ABA, PEG, H2O2, Cu(2+) and Zn(2+), but down-regulated by water. The role of HbMT2a protein in protecting against metal toxicity was demonstrated in vitro. PET-28a-HbMT2-beared Escherichia coli. Differential expression of HbMT2a upon treatment with 10 °C was observed in the detached leaves of rubber tree clone 93-114 which is cold-resistant and Reken501 which is cold-sensitive. The expression patterns of HbMT2a in the two rubber tree clones may be ascribed to a change in the level of endogenous H2O2.

[Ecological issues and loss risk of cold-resistant rubber germplasm resource in Xishuangbanna].

Based on the field investigation data and related literatures, this paper analyzed the ecological issues in the breeding and utilization of highly cold-resistant and high-yielding rubber germplasm in Xishuangbanna, and the loss risk of this germplasm. The rubber farmers in Xishuangbanna had built a cycle pattern of establishing rubber nursery and managing rubber plantations to develop the local rubber industry, but the unlimited utilization of this germplasm made the rubber plantations expand to the highland with an elevation of 1300 m, resulting in the ecological issues such as the destruction of natural forests and species, decline of soil fertility, regional drought, rubber diseases and insect pests, and human-elephant conflicts. This germplsm and its intellectual property were under transnational loss without any control, making even more potential ecological and economic risks. Great attentions should be paid by the related organizations to take effective measures to conserve this germplsm resource to reduce its potential risks.

[Development and application of EST-SSR markers in Hevea brasiliensis Muell. Arg].

Three thousand and ninety Unigenes were obtained from 10 778 Hevea brasiliensis ESTs. Four hundred and thirty SSRs were distributed in 353 Unigenes, which accounts for 11.42% of the total number of Unigenes. The frequency of SSRs was 1/3.93 kb. Dinucleotide and trinucleotide repeats were the dominant types among the obtained unigenes, accounting for 63.49% and 32.09%, respectively. TC/AG, CT/GA and CTT/GAA, AAG/TTC, and AGA/TCT were the most abundant motifs for dinucleotide and trinucleotide motifs. One hundred and forty-eight primer pairs were designed by PRIMER5.0 and 21 primer pairs were synthesized. Among them, 15 primer pairs can produce clear and stable bands, and the PCR products were screened in denaturing polyacrylamide gel following silver staining. Genetic diversity of 44 rubber clones were investigated with these primer pairs, and a dendrogram of 44 rubber clones was built. The results indicated that it is an effective and feasible way to develop EST-SSR markers from H. brasiliensis EST sequences, and the primers designed in this study can be used in genetic study of rubber tree.