De novo assembly and Transcriptome characterization of an endemic species of Vietnam, Panax vietnamensis Ha et Grushv., including the development of EST-SSR markers for population genetics

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Phylogeographic analysis reveals extensive genetic variation of native grass Elymus nutans (Poaceae) on the Qinghai-Tibetan plateau

Introduction

Elymus nutans holds ecological and pastoral significance due to its adaptability and nutritional value, the Qinghai-Tibet Plateau (QTP) is a key hub for its genetic diversity. To conserve and harness its genetic resources in highland ecosystems, a thorough assessment is vital. However, a comprehensive phylogeographic exploration of E. nutans is lacking. The objective of this study was to unravel the genetic diversity, adaptation, and phylogenetics of E. nutans populations.

Methods

Encompassing 361 individuals across 35 populations, the species' genetic landscape and dynamic responses to diverse environments were decoded by using four chloroplast DNA (cpDNA) sequences and nine microsatellite markers derived from the transcriptome.

Results and discussion

This study unveiled a notable degree of genetic diversity in E. nutans populations at nuclear (I = 0.46, He = 0.32) and plastid DNA levels (Hd = 0.805, π = 0.67). Analysis via AMOVA highlighted genetic variation predominantly within populations. Despite limited isolation by distance (IBD), the Mekong-Salween Divide (MSD) emerged as a significant factor influencing genetic differentiation and conserving diversity. Furthermore, correlations were established between external environmental factors and effective alleles of three EST-SSRs (EN5, EN57 and EN80), potentially linked to glutathione S-transferases T1 or hypothetical proteins, affecting adaptation. This study deepens the understanding of the intricate relationship between genetic diversity, adaptation, and environmental factors within E. nutans populations on the QTP. The findings shed light on the species' evolutionary responses to diverse ecological conditions and contribute to a broader comprehension of plant adaptation mechanisms.

Suspension cell cultures of Panax vietnamensis as a biotechnological source of ginsenosides: growth, cytology, and ginsenoside profile assessment

Introduction

Panax vietnamensis is a valuable medicinal plant and a source of a broad spectrum of biologically active ginsenosides of different structural groups. Overexploitation and low adaptability to planation cultivation have made this species vulnerable to human pressure and prompted the development of cell cultivation in vitro as a sustainable alternative to harvesting wild plants for their bioactive components. Despite high interest in biotechnological production, little is known about the main factors affecting cell growth and ginsenoside biosynthesis of this species under in vitro conditions. In this study, the potential of cell cultures of P. vietnamensis as a biotechnological source of ginsenosides was was assessed.

Methods

Six suspension cell lines that were developed from different sections of a single rhizome through a multi-step culture optimization process and maintained for over 3 years on media with different mineral salt base and varying contents of auxins and cytokinins. These cell lines were evaluated for productivity parameters and cytological characteristics. Ginsenoside profiles were assessed using a combination of the reversed-phase ultra-high-performance liquid chromatography-Orbitrap-tandem mass spectrometry (UHPLC-Orbitrap-MS/MS) and ultra-performance liquid chromatography-time of flight-mass spectrometry (UPLC-TOF-MS).

Results

All lines demonstrated good growth with a specific growth rate of 0.1-0.2 day-1, economic coefficient of 0.31-0.70, productivity on dry weight (DW) of 0.30-0.83 gDW (L·day)-1, and maximum biomass accumulation varying from 10 to 22 gDW L-1. Ginsenosides of the protopanaxadiol (Rb1, Rb2/Rb3, malonyl-Rb1, and malonyl-Rb2/Rb3), oleanolic acid (R0 and chikusetsusaponin IV), and ocotillol (vinaginsenoside R1) groups and their isomers were identified in cell biomass extracts. Chikusetsusaponin IV was identified in P. vietnamensis cell culture for the first time.

Discussion

These results suggest that suspension cell cultures of Vietnamese ginseng have a high potential for the biotechnological production of biomass containing ginsenosides, particularly of the oleanolic acid and ocotillol groups.

Genetic population structure of the Vietnamese ginseng (Panax vietnamensis Ha et Grushv.) detected by microsatellite analysis

Abstract Panax vietnamensis Ha et Grushv. is a precious medicinal species native to the tropical forests of Vietnam. Due to habitat loss and over-harvesting, this species is endangered in Vietnam. To conserve the species, we investigated genetic variability and population structure using nine microsatellites for 148 individuals from seven populations across the current distribution range of P. vietnamensis in Vietnam. We determined a moderate genetic diversity within populations (HO = 0.367, HE = 0.437) and relatively low population differentiation (the Weir and Cockerham index of 0.172 and the Hedrick index of 0.254) and showed significant differentiation (P < 0.05), which suggested fragmented habitats, over-utilization and over-harvesting of P. vietnamensis. Different clustering methods revealed that individuals were grouped into two major clusters, which were associated with gene flow across the geographical range of P. vietnamensis. This study also detected that ginseng populations can have undergone a recent bottleneck. We recommend measures in future P. vietnamensis conservation and breeding programs.

Development of unigene-derived SSR markers from RNA-seq data of Uraria lagopodioides (Fabaceae) and their application in the genus Uraria Desv. (Fabaceae)

BACKGROUND

Uraria Desv. belongs to the tribe Desmodieae (Fabaceae), a group of legume plants, some of which have medicinal properties. However, due to a lack of genomic information, the interspecific relationships, genetic diversity, population genetics, and identification of functional genes within Uraria species are still unclear.

RESULTS

Using RNA-Seq, a total of 66,026 Uraria lagopodioides unigenes with a total sequence content of 52,171,904 bp were obtained via de novo assembly and annotated using GO, KEGG, and KOG databases. 17,740 SSRs were identified from a set of 66,026 unigenes. Cross-species amplification showed that 54 out of 150 potential unigene-derived SSRs were transferable in Uraria, of which 19 polymorphic SSRs were developed. Cluster analysis based on polymorphisms successfully distinguished seven Uraria species and revealed their interspecific relationships. Seventeen samples of seven Uraria species were clustered into two monophyletic clades, and phylogenetic relationships of Uraria species based on unigene-derived SSRs were consistent with classifications based on morphological characteristics.

CONCLUSIONS

Unigenes annotated in the present study will provide new insights into the functional genomics of Uraria species. Meanwhile, the unigene-derived SSR markers developed here will be invaluable for assessing the genetic diversity and evolutionary history of Uraria and relatives.

Molecular Cloning and Functional Characterization of Oxidosqualene Cyclases from Panax vietnamensis

Panax vietnamensis is a valuable medicinal resource with promising preclinical applications. Ginsenosides, which are triterpenoids, are the primary active components in P. vietnamensis. Oxidosqualene cyclases (OSCs) catalyze the formation of the basic skeleton of triterpenes from 2,3-oxidosqualene, which is a crucial step in the biosynthesis of triterpenoids. The OSCs involved in triterpenoid biosynthesis in P. vietnamensis have not yet been characterized. Four OSC genes (PvOSC1-4) were cloned from P. vietnamensis and functionally characterized via heterologous expression in yeast. Transgenic yeast expressing PvOSC1, PvOSC3, and PvOSC4 produced the corresponding products β-amyrin, cycloartenol, and dammarenediol-II, respectively. PvOSC1, PvOSC3, and PvOSC4 are monofunctional OSCs. In this study, we characterized three PvOSC genes, providing a better understanding of the biosynthesis of triterpenoids in P. vietnamensis and the multiple choices of plant OSCs for metabolic engineering in yeast and other hosts.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

Upcoming progress of transcriptomics studies on plants: An overview

Transcriptome sequencing or RNA-Sequencing is a high-resolution, sensitive and high-throughput next-generation sequencing (NGS) approach used to study non-model plants and other organisms. In other words, it is an assembly of RNA transcripts from individual or whole samples of functional and developmental stages. RNA-Seq is a significant technique for identifying gene predictions and mining functional analysis that improves gene ontology understanding mechanisms of biological processes, molecular functions, and cellular components, but there is limited information available on this topic. Transcriptomics research on different types of plants can assist researchers to understand functional genes in better ways and regulatory processes to improve breeding selection and cultivation practices. In recent years, several advancements in RNA-Seq technology have been made for the characterization of the transcriptomes of distinct cell types in biological tissues in an efficient manner. RNA-Seq technologies are briefly introduced and examined in terms of their scientific applications. In a nutshell, it introduces all transcriptome sequencing and analysis techniques, as well as their applications in plant biology research. This review will focus on numerous existing and forthcoming strategies for improving transcriptome sequencing technologies for functional gene mining in various plants using RNA- Seq technology, based on the principles, development, and applications.

Molecular Characterization of Tinospora cordifolia (Willd.) Miers Using Novel g-SSR Markers and Their Comparison with EST-SSR and SCoT Markers for Genetic Diversity Study

In the present study, novel genomic-SSR (g-SSR) markers generated in our laboratory were used to characterize Tinospora cordifolia and related species. The g-SSR marker was also compared with EST-SSR and SCoT markers used earlier in our laboratory to assess the genetic diversity of T. cordifolia. A total of 26 accessions of T. cordifolia and 1 accession each of Tinospora rumphii and Tinospora sinensis were characterized using 65 novel g-SSR markers. A total of 125 alleles were detected with 49 polymorphic g-SSR markers. The number of alleles per locus varied from 1–4 with a mean value of 2.55 alleles per locus. Mean PIC, gene diversity, and heterozygosity were estimated to be 0.33, 0.41, and 0.65, respectively. The two species, namely T. rumphii and T. sinensis, showed cross-species transferability of g-SSRs developed in T. cordifolia. The success rate of cross-species transferability in T. rumphii was 95.3% and 93.8% in T. sinensis, proving the usefulness of this marker in genetic diversity studies of related species. The Tinospora accessions were also used for molecular characterization using SCoT and EST-SSR markers and compared for genetic diversity and cross-species transferability. The PIC, gene diversity, heterozygosity, and principal coordinate analysis showed that g-SSR is the better maker for a genetic diversity study of T. cordifolia. Additionally, high cross-species transferability of g-SSRs was found (95.3% and 93.8%) compared to EST-SSRs (68.8% and 67.7%) in T. rumphii and T. sinensis, respectively.

Phylogenetic Analysis Based on DNA Barcoding and Genetic Diversity Assessment of Morinda officinalis How in Vietnam Inferred by Microsatellites

Morinda officinalis How is well-known as a valuable medicinal plant found in some regions of Vietnam. This species is mainly used for treating male impotence, irregular menstruation, and rheumatoid arthritis. This study aimed to identify the species of and genetic diversity in three M. officinalis populations: one each in Quang Binh (QB), Thua Thien Hue (TTH), and Quang Nam (QN). In this study, four DNA barcoding markers (ITS1, ITS2, matK, and rbcL) were used to identify the species and 22 microsatellite markers were applied for population structure and diversity analyses. The results showed that the sequences of gene regions studied in M. officinalis had a high similarity (>95%) to the ITS1, ITS2, matK, and rbcL sequences of M. officinalis on BLAST. Of the four DNA barcoding markers used, ITS1 and ITS2 showed higher efficiency in DNA amplification of M. officinalis. From this study, 27 GenBank codes were published on BLAST. The results also revealed high levels of genetic diversity in populations. The average observed and expected heterozygosity values were HO = 0.513 and HE = 0.612, respectively. The average FST value was 0.206. Analysis of molecular variance (AMOVA) showed 70% variation within populations and 30% among populations. The population structure of M. officinalis inferred in STRUCTURE revealed that the optimum number of genetic groups for the admixture model was K = 2. These findings provided vital background information for future studies in the conservation of M. officinalis in both ex situ and in situ plans.

De novo transcriptome assembly, gene annotation, and EST-SSR marker development of an important medicinal and edible crop, Amomum tsaoko (Zingiberaceae)

BACKGROUND

Amomum tsaoko is a medicinal and food dual-use crop that belongs to the Zingiberaceae family. However, the lack of transcriptomic and genomic information has limited the understanding of the genetic basis of this species. Here, we performed transcriptome sequencing of samples from different A. tsaoko tissues, and identified and characterized the expressed sequence tag-simple sequence repeat (EST-SSR) markers.

RESULTS

A total of 58,278,226 high-quality clean reads were obtained and de novo assembled to generate 146,911 unigenes with an N50 length of 2002 bp. A total of 128,174 unigenes were successfully annotated by searching seven protein databases, and 496 unigenes were identified as annotated as putative terpenoid biosynthesis-related genes. Furthermore, a total of 55,590 EST-SSR loci were detected, and 42,333 primer pairs were successfully designed. We randomly selected 80 primer pairs to validate their polymorphism in A. tsaoko; 18 of these primer pairs produced distinct, clear, and reproducible polymorphisms. A total of 98 bands and 96 polymorphic bands were amplified by 18 pairs of EST-SSR primers for the 72 A. tsaoko accessions. The Shannon's information index (I) ranged from 0.477 (AM208) to 1.701 (AM242) with an average of 1.183, and the polymorphism information content (PIC) ranged from 0.223 (AM208) to 0.779 (AM247) with an average of 0.580, indicating that these markers had a high level of polymorphism. Analysis of molecular variance (AMOVA) indicated relatively low genetic differentiation among the six A. tsaoko populations. Cross-species amplification showed that 14 of the 18 EST-SSR primer pairs have transferability between 11 Zingiberaceae species.

CONCLUSIONS

Our study is the first to provide transcriptome data of this important medicinal and edible crop, and these newly developed EST-SSR markers are a very efficient tool for germplasm evaluation, genetic diversity, and molecular marker-assisted selection in A. tsaoko.

Target Capture Reveals the Complex Origin of Vietnamese Ginseng

The global market of the medicinal plant ginseng is worth billions of dollars. Many ginseng species are threatened in the wild and effective sustainable development initiatives are necessary to preserve biodiversity at species and genetic level whilst meeting the demand for medicinal produce. This is also the case of Panax vietnamensis Ha & Grushv., an endemic and threatened ginseng species in Vietnam that is locally cultivated at different scales and has been the object of national breeding programs. To investigate the genetic diversity within cultivated and wild populations of P. vietnamensis we captured 353 nuclear markers using the Angiosperm-353 probe set. Genetic diversity and population structure were evaluated for 319 individuals of Vietnamese ginseng across its area of distribution and from wild and a varying range of cultivated areas. In total, 319 individuals were sampled. After filtering, 1,181 SNPs were recovered. From the population statistics, we observe high genetic diversity and high genetic flow between populations. This is also supported by the STRUCTURE analysis. The intense gene flow between populations and very low genetic differentiation is observed regardless of the populations' wild or cultivated status. High levels of admixture from two ancestral populations exist in both wild and cultivated samples. The high gene flow between populations can be attributed to ancient and on-going practices of cultivation, which exist in a continuum from understorey, untended breeding to irrigated farm cultivation and to trade and exchange activities. These results highlight the importance of partnering with indigenous peoples and local communities and taking their knowledge into account for biodiversity conservation and sustainable development of plants of high cultural value.

De novo whole-genome assembly and discovery of genes involved in triterpenoid saponin biosynthesis of Vietnamese ginseng (Panax vietnamensis Ha et Grushv.)

Vietnamese ginseng (Panax vietnamensis Ha et Grushv.), also known as Ngoc Linh ginseng, is a high-value herb in Vietnam. Vietnamese ginseng has been proven to be effective in enhancing the immune system, human memory, anti-stress, anti-inflammatory, anti-cancer, and prevent aging. The present study reports the first draft whole-genome of Vietnamese ginseng and the identification of potential genes involved in the triterpenoid metabolic pathway. De novo whole-genome assembly was performed successfully from a data of approximately 139 Gbps of 394,802,120 high quality reads to generate 9815 scaffolds with an N50 value of 572,722 bp from the leaf of Vietnamese ginseng. The assembled genome of Vietnamese ginseng is 3,001,967,204 bp long containing 79,374 gene models. Among them, there are 55,012 genes (69.30%) were annotated by various public molecular biology databases. The potential genes involved in triterpenoid saponin biosynthesis in Vietnamese ginseng and their metabolic pathway were also predicted." Three genes encoding squalene monooxygenase isozymes in Vietnamese ginseng were cloned, sequenced and characterized. Moreover, expression levels of several key genes involved in terpenoid biosynthesis in different parts of Vietnamese ginseng were also analyzed. The SSR markers were detected by various programs from both of assembly full dataset of Vietnamese ginseng genome and predicted genes. The present work provided important data of the draft whole-genome of Vietnamese ginseng for further studies to understand the role of genes involved in ginsenoside biosynthesis and their metabolic pathway at the molecular level of this rare medicinal species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01076-1.

Research Progress and Future Development Trends in Medicinal Plant Transcriptomics

Transcriptomics is one of the most popular topics in biology in recent times. Transcriptome sequencing (RNA-Seq) is a high-throughput, high-sensitivity, and high-resolution technique that can be used to study model and non-model organisms. Transcriptome sequencing is also an important method for studying the genomes of medicinal plants, a topic on which limited information is available. The study of medicinal plants through transcriptomics can help researchers analyze functional genes and regulatory mechanisms of medicinal plants and improve breeding selection and cultivation techniques. This article analyzes and compares the applications of transcriptome sequencing in medicinal plants over the past decade and briefly introduces the methods of transcriptome sequencing and analysis, their applications in medicinal plant research, and potential development trends. We will focus on the research and application progress of transcriptome sequencing in the following four areas: the mining of functional genes in medicinal plants, development of molecular markers, biosynthetic pathways of secondary metabolites, and developmental mechanisms of medicinal plants. Our review will provide ideas for the mining of functional genes of medicinal plants and breeding new varieties.

A fast likelihood approach for estimation of large phylogenies from continuous trait data

Despite the recent availability of large-scale genomic data for many individuals, few methods for phylogenetic inference are both computationally efficient and highly accurate for trees with hundreds of taxa. Model-based methods such as those developed in the maximum likelihood and Bayesian frameworks are especially time-consuming, as they involve both computationally intensive calculations on fixed phylogenies and searches through the space of possible phylogenies, and they are known to scale poorly with the addition of taxa. Here, we propose a fast approximation to the maximum likelihood estimator that directly uses continuous trait data, such as allele frequency data. The approximation works by first computing the maximum likelihood estimates of some internal branch lengths, and then inferring the tree-topology using these estimates. Our approach is more computationally efficient than existing methods for such data while still achieving comparable accuracy. This method is innovative in its use of the mathematical properties of tree-topologies for inference, and thus serves as a useful addition to the collection of methods available for estimating phylogenies from continuous trait data.