Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Development of Novel Markers for Yield in Hevea brasiliensis Muell. Arg. 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.

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 relationship between latex metabolism gene expression with rubber yield and related traits in Hevea brasiliensis

Background

Expression patterns of many laticifer-specific gens are closely correlative with rubber yield of Hevea brasiliensis (para rubber tree). To unveil the mechanisms underlying the rubber yield, transcript levels of nine major latex metabolism-related genes, i.e., HMG-CoA synthase (HMGS), HMG-CoA reductase (HMGR), diphosphomevalonate decarboxylase (PMD), farnesyl diphosphate synthase (FPS), cis-prenyltransferase (CPT), rubber elongation factor (REF), small rubber particle protein (SRPP), dihydroxyacid dehydratase (DHAD) and actin depolymerizing factor (ADF), were dertermined, and the relationship between rubber yield with their expression levels was analysed.

Results

Except HbHMGR1, HbPMD and HbDHAD, most of these genes were predominantly expressed in latex, and bark tapping markedly elevated the transcript abundance of the analyzed genes, with the 7th tapping producing the greatest expression levels. Both ethephon (ETH) and methyl jasmonate (MeJA) stimulation greatly induced the expression levels of the examined genes, at least at one time point, except HbDHAD, which was unresponsive to MeJA. The genes’ expression levels, as well as the rubber yields and two yield characteristics differed significantly among the different genotypes examined. Additionally, the latex and dry rubber yields increased gradually but the dry rubber content did not. Rubber yields and/or yield characteristics were significantly positively correlated with HbCPT, HbFPS, HbHMGS, HbHMGR1 and HbDHAD expression levels, negatively correlated with that of HbREF, but not significantly correlated with HbPMD, HbSRPP and HbADF expression levels. In addition, during rubber production, significantly positive correlations existed between the expression level of HbPMD and the levels of HbREF and HbHMGR1, between HbSRPP and the levels of HbHMGS and HbHMGR1, and between HbADF and HbFPS.

Conclusions

The up-regulation of these genes might be related to the latex production of rubber trees under the stress of bark tapping and latex metabolism. The various correlations among the genes implied that there are differences in their synergic interactions. Thus, these nine genes might be related to rubber yield and yield-related traits in H. brasiliensis, and this work increases our understanding of their complex functions and how they are expressed in both high-and medium-yield rubber tree varieties and low-yield wild rubber tree germplasm.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5242-4) contains supplementary material, which is available to authorized users.