Genomic technologies for Hevea breeding

Genome-Wide Identification of the Geranylgeranyl Pyrophosphate Synthase (GGPS) Gene Family Associated with Natural Rubber Synthesis in Taraxacum kok-saghyz L. Rodin

Taraxacum kok-saghyz Rodin (TKS) is a recognized alternative source of natural rubber comparable to the rubber tree. The geranylgeranyl pyrophosphate synthase (GGPS) catalyzed the synthesis of geranylgeranyl pyrophosphate (GGPP), which is an important enzyme in the secondary metabolism pathway. In this study, we present the first analysis of the GGPS gene family in TKS, where a total of seven TkGGPS family members were identified. Their core motifs, conserved structural domains, gene structures, and cis-acting elements were described. In addition, two phylogenetic trees were constructed based on the Neighbor-Joining and Maximum-Likelihood methods, and the TkGGPSs were highly conserved and exhibited good collinearity with the other species. Transcriptome data showed that seven TkGGPS gene members were expressed in all the 12 tissues measured, and TkGGPS1, TkGGPS3, and TkGGPS6 were highly expressed in latex, suggesting that they may be associated with natural rubber synthesis. Meanwhile, quantitative real-time PCR (qRT-PCR) showed that the expression levels of the TkGGPS genes were regulated by the ethylene and methyl jasmonate (MeJA) pathways. Subcellular localization results indicated that all the TkGGPS proteins were also located in chloroplasts involved in photosynthesis in plants. This study will provide valuable insights into the selection of candidate genes for molecular breeding and natural rubber biosynthesis in TKS.

Overexpression of HbGRF4 or HbGRF4-HbGIF1 Chimera Improves the Efficiency of Somatic Embryogenesis in Hevea brasiliensis

Transgenic technology is a crucial tool for gene functional analysis and targeted genetic modification in the para rubber tree (Hevea brasiliensis). However, low efficiency of plant regeneration via somatic embryogenesis remains a bottleneck of successful genetic transformation in H. brasiliensis. Enhancing expression of GROWTH-REGULATING FACTOR 4 (GRF4)-GRF-INTERACTING FACTOR 1 (GIF1) has been reported to significantly improve shoot and embryo regeneration in multiple crops. Here, we identified endogenous HbGRF4 and HbGIF1 from the rubber clone Reyan7-33-97, the expressions of which dramatically increased along with somatic embryo (SE) production. Intriguingly, overexpression of HbGRF4 or HbGRF4-HbGIF1 markedly enhanced the efficiency of embryogenesis in two H. brasiliensis callus lines with contrasting rates of SE production. Transcriptional profiling revealed that the genes involved in jasmonic acid response were up-regulated, whereas those in ethylene biosynthesis and response as well as the S-adenosylmethionine-dependent methyltransferase activity were down-regulated in HbGRF4- and HbGRF4-HbGIF1-overexpressing H. brasiliensis embryos. These findings open up a new avenue for improving SE production in rubber tree, and help to unravel the underlying mechanisms of HbGRF4-enhanced somatic embryogenesis.

Comprehensive Analysis of the Catalase (CAT) Gene Family and Expression Patterns in Rubber Tree (Hevea brasiliensis) under Various Abiotic Stresses and Multiple Hormone Treatments

Catalase (CAT) is one of the key enzymes involved in antioxidant defense systems and mainly scavenges H2O2 and plays a vital role in plant growth, development, and various adverse stresses. To date, a systematic study of the CAT gene family in rubber tree has not been reported. In this study, five HbCAT gene family members were identified from the rubber tree genome, and these were mainly clustered into two subfamilies. Gene structure and motif analysis showed that exon-intron and motif patterns were conserved across different plant species. Sequence analysis revealed that HbCAT proteins contain one active catalytic site, one heme-ligand signature sequence, three conserved amino acid residues (His, Tyr, and Asn), and one peroxisome-targeting signal 1 (PTS1) sequence. Fragment duplication is a selection pressure for the evolution of the HbCAT family based on Ka/Ks values. Analysis of cis-acting elements in the promoters indicated that HbCAT gene expression might be regulated by abscisic acid (ABA), salicylic acid (SA), and MYB transcription factors; furthermore, these genes might be involved in plant growth, development, and abiotic stress responses. A tissue-specific expression analysis showed that HbCATs gradually increased with leaf development and were highly expressed in mature leaves. Gene expression profiling exhibited the differential expression of the HbCATs under cold, heat, drought, and NaCl stresses. Our results provide comprehensive information about the HbCAT gene family, laying the foundation for further research on its function in rubber tree.

Genomic insight into domestication of rubber tree

Understanding the genetic basis of rubber tree (Hevea brasiliensis) domestication is crucial for further improving natural rubber production to meet its increasing demand worldwide. Here we provide a high-quality H. brasiliensis genome assembly (1.58 Gb, contig N50 of 11.21 megabases), present a map of genome variations by resequencing 335 accessions and reveal domestication-related molecular signals and a major domestication trait, the higher number of laticifer rings. We further show that HbPSK5, encoding the small-peptide hormone phytosulfokine (PSK), is a key domestication gene and closely correlated with the major domestication trait. The transcriptional activation of HbPSK5 by myelocytomatosis (MYC) members links PSK signaling to jasmonates in regulating the laticifer differentiation in rubber tree. Heterologous overexpression of HbPSK5 in Russian dandelion (Taraxacum kok-saghyz) can increase rubber content by promoting laticifer formation. Our results provide an insight into target genes for improving rubber tree and accelerating the domestication of other rubber-producing plants.

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.

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.

Genome-Wide Identification and Comparative Analysis of WOX Genes in Four Euphorbiaceae Species and Their Expression Patterns in Jatropha curcas

The WUSCHEL-related homeobox (WOX) proteins are widely distributed in plants and play important regulatory roles in growth and development processes such as embryonic development and organ development. Here, series of bioinformatics methods were utilized to unravel the structural basis and genetic hierarchy of WOX genes, followed by regulation of the WOX genes in four Euphorbiaceae species. A genome-wide survey identified 59 WOX genes in Hevea brasiliensis (H. brasiliensis: 20 genes), Jatropha curcas (J. curcas: 10 genes), Manihot esculenta (M. esculenta: 18 genes), and Ricinus communis (R. communis: 11 genes). The phylogenetic analysis revealed that these WOX members could be clustered into three close proximal clades, such as namely ancient, intermediate and modern/WUS clades. In addition, gene structures and conserved motif analyses further validated that the WOX genes were conserved within each phylogenetic clade. These results suggested the relationships among WOX members in the four Euphorbiaceae species. We found that WOX genes in H. brasiliensis and M. esculenta exhibit close genetic relationship with J. curcas and R. communis. Additionally, the presence of various cis-acting regulatory elements in the promoter of J. curcas WOX genes (JcWOXs) reflected distinct functions. These speculations were further validated with the differential expression profiles of various JcWOXs in seeds, reflecting the importance of two JcWOX genes (JcWOX6 and JcWOX13) during plant growth and development. Our quantitative real-time PCR (qRT-PCR) analysis demonstrated that the JcWOX11 gene plays an indispensable role in regulating plant callus. Taken together, the present study reports the comprehensive characteristics and relationships of WOX genes in four Euphorbiaceae species, providing new insights into their characterization.

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.

Complete chloroplast genome of the medicinal plant Evolvulus alsinoides: comparative analysis, identification of mutational hotspots and evolutionary dynamics with species of Solanales

Evolvulus alsinoides, belonging to the family Convolvulaceae, is an important medicinal plant widely used as a nootropic in the Indian traditional medicine system. In the genus Evolvulus, no research on the chloroplast genome has been published. Hence, the present study focuses on annotation, characterization, identification of mutational hotspots, and phylogenetic analysis in the complete chloroplast genome (cp) of E. alsinoides. Genome comparison and evolutionary dynamics were performed with the species of Solanales. The cp genome has 114 genes (80 protein-coding genes, 30 transfer RNA, and 4 ribosomal RNA genes) that were unique with total genome size of 157,015 bp. The cp genome possesses 69 RNA editing sites and 44 simple sequence repeats (SSRs). Predicted SSRs were randomly selected and validated experimentally. Six divergent hotspots such as trnQ-UUG, trnF-GAA, psaI, clpP, ndhF, and ycf1 were discovered from the cp genome. These microsatellites and divergent hot spot sequences of the Taxa 'Evolvulus' could be employed as molecular markers for species identification and genetic divergence investigations. The LSC area was found to be more conserved than the SSC and IR region in genome comparison. The IR contraction and expansion studies show that nine genes rpl2, rpl23, ycf1, ycf2, ycf1, ndhF, ndhA, matK, and psbK were present in the IR-LSC and IR-SSC boundaries of the cp genome. Fifty-four protein-coding genes in the cp genome were under negative selection pressure, indicating that they were well conserved and were undergoing purifying selection. The phylogenetic analysis reveals that E. alsinoides is closely related to the genus Cressa with some divergence from the genus Ipomoea. This is the first time the chloroplast genome of the genus Evolvulus has been published. The findings of the present study and chloroplast genome data could be a valuable resource for future studies in population genetics, genetic diversity, and evolutionary relationship of the family Convolvulaceae.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01051-w.