Evolution of the Araliaceae family inferred from complete chloroplast genomes and 45S nrDNAs of 10 Panax-related species

Phylogenetic exploration, codon usage bias, and genomic divergence in Hydrocotyle: a comparative plastome study across different geographical locations

Hydrocotyle himalaica from Bhutan, a perennial herb that thrives from 1500 to 2600 m, possesses both ecological importance and medicinal properties. The plastome analysis revealed a length of 153,383 bp, showing variation from conspecific taxa in China. Its standard structure comprises two IR regions (18,336 bp IRa and 18,336 bp IRb), an LSC region of 97,944 bp, and an SSC region of 18,767 bp, with a GC content of 37.63%. Non-coding regions showed higher mutation susceptibility, with Pi values from 0.006 to 0.107. An AT-rich codon bias was consistent across all 18 Hydrocotyle species. Nucleotide composition and GC% in coding sequences differed among the species. The codon preference in Hydrocotyle is shaped by multiple factors, with natural selection being the primary influence, as indicated by the ENC-plot, PR2-plot, and Neutrality-plot. Codon usage patterns varied, with RSCU values from 0 to 2.23. Codons ending in A or U had RSCU > 1, while those ending in C or G had RSCU < 1. GC2 content surpassed GC3 and GC1 in most genes. The phylogenetic analysis placed H. himalaica, sourced from Kanglung, Bhutan, within the monophyly of the Hydrocotyloideae subfamily. However, the species showed weaker bootstrap support (BS < 50) with H. javanica and H. hookeri subsp., a deviation from a prior report on the same species from Jiangkou, Guizhou, China. This analysis highlighted the genomic characteristics and evolutionary relationships of H. himalaica from Bhutan, underscoring the need for a comprehensive phylogenetic, ecological, and botanical characterization to confirm intra-specific variation within Hydrocotyle species.

Manniosides G-J, New Ursane- and Lupane-Type Saponins from Schefflera mannii (Hook.f.) Harms

Four previously unreported triterpenoid saponins named 3β-hydroxy-23-oxours-12-en-28-oic acid 28-O-β-D-glucopyranosyl ester (mannioside G) (1), 23-O-acetyl-3β-hydroxyurs-12-en-28-oic acid 28-O-β-D-glucopyranosyl ester (mannioside H) (2), ursolic acid 28-O-[α-L-rhamnopyranosyl-(1→4)-β-D-glucopyranosyl-(1→6)-β-D-glucopyranosyl] ester (mannioside I) (3), and 3β-hydroxy-23-oxolup-20(29)-en-28-oic acid 28-O-β-D-glucopyranosyl ester (mannioside J) (4) were isolated as minor constituents from the EtOAc soluble fraction of the MeOH extract of the leaves of Schefflera mannii along with the known compounds 23-hydroxyursolic acid 28-O-β-D-glucopyranosyl ester (5), ursolic acid 28-O-β-D-glucopyranosyl ester (6), pulsatimmoside B (7) betulinic acid 28-O-[α-L-rhamnopyranosyl-(1→4)-β-D-glucopyranosyl-(1→6)-β-D-glucopyranosyl] ester (8), 23-hydroxy-3-oxo-urs-12-en-28-oic acid (9), hederagenin (10), ursolic acid (11), betulinic acid (12), and lupeol (13). Their structures were elucidated by a combination of 1D and 2D NMR analysis and mass spectrometry. The MeOH extract, the EtOAc and n-BuOH fractions, and some of the isolated compounds were evaluated for their antibacterial activity against four bacteria: Staphylococcus aureus ATCC1026, Staphylococcus epidermidis ATCC 35984, Escherichia coli ATCC10536, and Klepsiella pnemoniae ATCC13882. They were also screened for their antioxidant properties, but no significant results were obtained.

The mitochondrial genomes of Panax notoginseng reveal recombination mediated by repeats associated with DNA replication

Panax notoginseng is one of the most valuable medicinal species. However, its mitochondrial genome has not been reported yet. We aimed to determine the mitogenome sequence of P. notoginseng. We de novo assembled the mitogenome with Illumina short reads and Nanopore long reads. The mitochondrial genome of P. notoginseng has a multipartite structure consisting of interconversion between a "master circle" and numerous "subgenomic circles" through recombinations mediated by 64 pairs of repetitive sequences. Among the multipartite structure, seven subgenomic circles were best supported. Six of the seven subgenomic circles shared an 852 bp conserved fragment. The complete mitogenome of P. notoginseng was 662,479 bp long including 34 mitochondrial protein-coding genes (PCGs), three rRNA, and 19 tRNA genes. We identified 166 microsatellite repeats and 26 long-tandem repeats. Phylogenetic analysis resolved a tree that was mostly congruent with the phylogeny of Apiales species described in the APG IV system and the tree built with the chloroplast genome sequences. A total of 12 mitochondrial plastid DNA fragments were identified. Lastly, we predicted 591C-to-U RNA editing sites in the coding regions of mitochondrial PCGs. The mitochondrial genome will lay the foundation for understanding the evolution of Panax species.

New Insights into Phylogenetic Relationship of Hydrocotyle (Araliaceae) Based on Plastid Genomes

Hydrocotyle, belonging to the Hydrocotyloideae of Araliaceae, consists of 95 perennial and 35 annual species. Due to the lack of stable diagnostic morphological characteristics and high-resolution molecular markers, the phylogenetic relationships of Hydrocotyle need to be further investigated. In this study, we newly sequenced and assembled 13 whole plastid genomes of Hydrocotyle and performed comparative plastid genomic analyses with four previously published Hydrocotyle plastomes and phylogenomic analyses within Araliaceae. The plastid genomes of Hydrocotyle exhibited typical quadripartite structures with lengths from 152,659 bp to 153,669 bp, comprising a large single-copy (LSC) region (83,958-84,792 bp), a small single-copy (SSC) region (18,585-18,768 bp), and a pair of inverted repeats (IRs) (25,058-25,145 bp). Each plastome encoded 113 unique genes, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Comparative analyses showed that the IR boundaries of Hydrocotyle plastomes were highly similar, and the coding and IR regions exhibited more conserved than non-coding and single-copy (SC) regions. A total of 2932 simple sequence repeats and 520 long sequence repeats were identified, with specificity in the number and distribution of repeat sequences. Six hypervariable regions were screened from the SC region, including four intergenic spacers (IGS) (ycf3-trnS, trnS-rps4, petA-psbJ, and ndhF-rpl32) and two coding genes (rpl16 and ycf1). Three protein-coding genes (atpE, rpl16, and ycf2) were subjected to positive selection only in a few species, implying that most protein-coding genes were relatively conserved during the plastid evolutionary process. Plastid phylogenomic analyses supported the treatment of Hydrocotyle from Apiaceae to Araliaceae, and topologies with a high resolution indicated that plastome data can be further used in the comprehensive phylogenetic research of Hydrocotyle. The diagnostic characteristics currently used in Hydrocotyle may not accurately reflect the phylogenetic relationships of this genus, and new taxonomic characteristics may need to be evaluated and selected in combination with more comprehensive molecular phylogenetic results.

The complete chloroplast genome of Ulmus mianzhuensis with insights into structural variations, adaptive evolution, and phylogenetic relationships of Ulmus (Ulmaceae)

BACKGROUND

Ulmus mianzhuensis is an endemic tree species in China with high ornamental and economic value. Currently, little is known regarding its genomic architecture, phylogenetic position, or adaptive evolution. Here, we sequenced the complete chloroplast genome (cp genome) of U. mianzhuensis and further compared the variations in gene organization and structure within Ulmus species to define their genomic evolution, then reconstructed the phylogenomic relationship of 31 related Ulmus species to explore the systematic position of U. mianzhuensis and the utility of cp genome for resolving phylogenetics among Ulmus species.

RESULTS

Our results revealed that all the Ulmus species exhibited a typical quadripartite structure, with a large single copy (LSC) region of 87,170 - 88,408 bp, a small single copy (SSC) region of 18,650 - 19,038 bp and an inverted repeat (IR) region of 26,288 - 26,546 bp. Within Ulmus species, gene structure and content of cp genomes were highly conserved, although slight variations were found in the boundary of SC/IR regions. Moreover, genome-wide sliding window analysis uncovered the variability of ndhC-trnV-UAC, ndhF-rpl32, and psbI-trnS-GCU were higher among 31 Ulmus that may be useful for the population genetics and potential DNA barcodes. Two genes (rps15 and atpF) were further detected under a positive selection of Ulmus species. Comparative phylogenetic analysis based on the cp genome and protein-coding genes revealed consistent topology that U. mianzhuensis is a sister group to U. parvifolia (sect. Microptelea) with a relatively low-level nucleotide variation of the cp genome. Additionally, our analyses also found that the traditional taxonomic system of five sections in Ulmus is not supported by the current phylogenomic topology with a nested evolutionary relationship between sections.

CONCLUSIONS

Features of the cp genome length, GC content, organization, and gene order were highly conserved within Ulmus. Furthermore, molecular evidence from the low variation of the cp genome suggested that U. mianzhuensis should be merged into U. parvifolia and regarded as a subspecies of U. parvifolia. Overall, we demonstrated that the cp genome provides valuable information for understanding the genetic variation and phylogenetic relationship in Ulmus.

The Complete Chloroplast Genomes of Gynostemma Reveal the Phylogenetic Relationships of Species within the Genus

Gynostemma is an important medicinal and food plant of the Cucurbitaceae family. The phylogenetic position of the genus Gynostemma in the Cucurbitaceae family has been determined by morphology and phylogenetics, but the evolutionary relationships within the genus Gynostemma remain to be explored. The chloroplast genomes of seven species of the genus Gynostemma were sequenced and annotated, of which the genomes of Gynostemma simplicifolium, Gynostemma guangxiense and Gynostemma laxum were sequenced and annotated for the first time. The chloroplast genomes ranged from 157,419 bp (Gynostemma compressum) to 157,840 bp (G. simplicifolium) in length, including 133 identical genes: 87 protein-coding genes, 37 tRNA genes, eight rRNA genes and one pseudogene. Phylogenetic analysis showed that the genus Gynostemma is divided into three primary taxonomic clusters, which differs from the traditional morphological classification of the genus Gynostemma into the subgenus Gynostemma and Trirostellum. The highly variable regions of atpH-atpL, rpl32-trnL, and ccsA-ndhD, the repeat unilts of AAG/CTT and ATC/ATG in simple sequence repeats (SSRs) and the length of overlapping regions between rps19 and inverted repeats(IRb) and between ycf1 and small single-copy (SSC) were found to be consistent with the phylogeny. Observations of fruit morphology of the genus Gynostemma revealed that transitional state species have independent morphological characteristics, such as oblate fruit and inferior ovaries. In conclusion, both molecular and morphological results showed consistency with those of phylogenetic analysis.

The first report describes features of the chloroplast genome of Withania frutescens

Genomic studies not only help researcher not only to identify genomic features in organisms, but also facilitate understanding of evolutionary relationships. Species in the Withania genus have medicinal benefits, and one of them is Withania frutescens, which is used to treat various diseases. This report investigates the nucleotides and genic features of chloroplast genome of Withania frutescens and trying to clarify the evolutionary relationship with Withania sp and family Solanaceae. We found that the total size of Withania frutescens chloroplast genome was 153.771 kb (the smallest chloroplast genome in genus Withania). A large single-copy region (91.285 kb), a small single-copy region (18.373 kb) form the genomic region, and are distinct from each other by a large inverted repeat (22.056 kb). 137 chloroplast genes are found including 4 rRNAs, 38 tRNAs and 83 protein-coding genes. The Withania frutescens chloroplast genome as well as four closest relatives was compared for features such as structure, nucleotide composition, simple sequence repeats (SSRs) and codon bias. Compared to other Withania species, Withania frutescens has unique characteristics. It has the smallest chloroplast genome of any Withania species, isoleucine is the major amino acid, and tryptophan is the minor, In addition, there are no ycf3 and ycf4 genes, fourth, there are only fifteen replicative genes, while in most other species there are more. Using fast minimum evolution and neighbor joining, we have reconstructed the trees to confirm the relationship with other Solanacaea species. The Withania frutescens chloroplast genome is submitted under accession no. ON153173.

High-throughput discovery of plastid genes causing albino phenotypes in ornamental chimeric plants

Chimeric plants composed of green and albino tissues have great ornamental value. To unveil the functional genes responsible for albino phenotypes in chimeric plants, we inspected the complete plastid genomes (plastomes) in green and albino leaf tissues from 23 ornamental chimeric plants belonging to 20 species, including monocots, dicots, and gymnosperms. In nine chimeric plants, plastomes were identical between green and albino tissues. Meanwhile, another 14 chimeric plants were heteroplasmic, showing a mutation between green and albino tissues. We identified 14 different point mutations in eight functional plastid genes related to plastid-encoded RNA polymerase (rpo) or photosystems which caused albinism in the chimeric plants. Among them, 12 were deleterious mutations in the target genes, in which early termination appeared due to small deletion-mediated frameshift or single nucleotide substitution. Another was single nucleotide substitution in an intron of the ycf3 and the other was a missense mutation in coding region of the rpoC2 gene. We inspected chlorophyll structure, protein functional model of the rpoC2, and expression levels of the related genes in green and albino tissues of Reynoutria japonica. A single amino acid change, histidine-to-proline substitution, in the rpoC2 protein may destabilize the peripheral helix of plastid-encoded RNA polymerase, impairing the biosynthesis of the photosynthesis system in the albino tissue of R. japonica chimera plant.

Analysis of the complete plastomes and nuclear ribosomal DNAs from Euonymus hamiltonianus and its relatives sheds light on their diversity and evolution

Euonymus hamiltonianus and its relatives (Celastraceae family) are used for ornamental and medicinal purposes. However, species identification in Euonymus is difficult due to their morphological diversity. Using plastid genome (plastome) data, we attempt to reveal phylogenetic relationship among Euonymus species and develop useful markers for molecular identification. We assembled the plastome and nuclear ribosomal DNA (nrDNA) sequences from five Euonymus lines collected from South Korea: three Euonymus hamiltonianus accessions, E. europaeus, and E. japonicus. We conducted an in-depth comparative analysis using ten plastomes, including other publicly available plastome data for this genus. The genome structures, gene contents, and gene orders were similar in all Euonymus plastomes in this study. Analysis of nucleotide diversity revealed six divergence hotspots in their plastomes. We identified 339 single nucleotide polymorphisms and 293 insertion or deletions among the four E. hamiltonianus plastomes, pointing to abundant diversity even within the same species. Among 77 commonly shared genes, 9 and 33 were identified as conserved genes in the genus Euonymus and E. hamiltonianus, respectively. Phylogenetic analysis based on plastome and nrDNA sequences revealed the overall consensus and relationships between plastomes and nrDNAs. Finally, we developed six barcoding markers and successfully applied them to 31 E. hamiltonianus lines collected from South Korea. Our findings provide the molecular basis for the classification and molecular taxonomic criteria for the genus Euonymus (at least in Korea), which should aid in more objective classification within this genus. Moreover, the newly developed markers will be useful for understanding the species delimitation of E. hamiltonianus and closely related species.

Evidence for saponin diversity-mycobiome links and conservatism of plant-fungi interaction patterns across Holarctic disjunct Panax species

Although interplays between plant and coevolved microorganisms are believed to drive landscape formation and ecosystem services, the relationships between the mycobiome and phytochemical evolution and the evolutionary characteristics of plant-mycobiome interaction patterns are still unclear. The present study explored fungal communities from 405 multiniche samples of three Holarctic disjunct Panax species. The overall mycobiomes showed compartment-dominated variations and dynamic universality. Neutral models were fitted for each compartment at the Panax genus (I) and species (II) levels to infer the community assembly mechanism and identify fungal subgroups potentially representing different plant-fungi interaction results, i.e., the potentially selected, opposed, and neutral taxa. Selection contributed more to the endosphere than to external compartments. The nonneutral taxa showed significant phylogenetic clustering. In Model I, the opposed subgroups could best reflect Panax saponin diversities (r = 0.69), and genera with highly positive correlations to specific saponins were identified using machine learning. Although mycobiomes in the three species differed significantly, subgroups in Model II were phylogenetically clustered based on potential interaction type rather than plant species, indicating potentially conservative plant-fungi interactions. In summary, the finding of strong links between invaders and saponin diversity can help explore the underlying mechanisms of saponin biosynthesis evolution from microbial insights, which is important to understanding the formation of the current landscape. The potential conservatism of plant-fungi interaction patterns suggests that the related genetic modules and selection pressures were convergent across Panax species, advancing our understanding of plant interplay with biotic environments.

Chloroplast Genomes of Genus Tilia: Comparative Genomics and Molecular Evolution

Tilia is a complex genus in the family Malvaceae that has high ecological and economical values. Owing to the lack of sufficient distinguishable morphological and molecular characteristics, interspecific relationships in this genus are not clear. Chloroplast (cp) genomes are small, meanwhile most angiosperms usually undergo matrilineal inheritance. Consequently, they can be used in molecular systematics and phylogenetic analyses. Here, we sequenced and assembled cp genomes from T. endochrysea, T. tomentosa, T. miqueliana, T. americana and T. cordata, and compared them with those of seven previously reported Tilia species. Similar gene contents, gene orders and GC contents existed among the 12 cp genomes, which ranged from 162,564 to 162,855 bp and encoded 113 unique genes. Abundant simple sequence repeats (119–127) and dispersed repeats (97–135) were detected in Tilia cp genomes. In total, 11 hypervariable regions were identified that could be suitable for species identification and phylogenetic studies. A phylogenetic analysis of Malvaceae based on 5 hypervariable genes (matK + ndhF + rpoB + rpoC2+ycf1) revealed that all eight subfamilies were monophyletic groups. Additionally, the genus Tilia was divided into three groups on the basis of all 521 molecular variation loci. The current study provides valuable insights into the genomic evolution of the genus Tilia.

The complete plastid genome of the endangered shrub Brassaiopsis angustifolia (Araliaceae): Comparative genetic and phylogenetic analysis

Brassaiopsis angustifolia K.M. Feng belongs to the family Araliaceae, and is an endangered shrub species in southwest China. Despite the importance of this species, the plastid genome has not been sequenced and analyzed. In this study, the complete plastid genome of B. angustifolia was sequenced, analyzed, and compared to the eight species in the Araliaceae family. Our study reveals that the complete plastid genome of B. angustifolia is 156,534 bp long, with an overall GC content of 37.9%. The chloroplast genome (cp) encodes 133 genes, including 88 protein-coding genes, 37 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. All protein-coding genes consisted of 21,582 codons. Among the nine species of Araliaceae, simple sequence repeats (SSRs) and five large repeat sequences were identified with total numbers ranging from 37 to 46 and 66 to 78, respectively. Five highly divergent regions were successfully identified that could be used as potential genetic markers of Brassaiopsis and Asian Palmate group. Phylogenetic analysis of 47 plastomes, representing 19 genera of Araliaceae and two related families, was performed to reconstruct highly supported relationships for the Araliaceae, which highlight four well-supported clades of the Hydrocotyle group, Greater Raukaua group, Aralia-Panax group, and Asian Palmate group. The genus Brassaiopsis can be divided into four groups using internal transcribed spacer (ITS) data. The results indicate that plastome and ITS data can contribute to investigations of the taxonomy, and phylogeny of B. angustifolia. This study provides a theoretical basis for species identification and future biological research on resources of the genus Brassaiopsis.

Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata

Araliaceae species produce various classes of triterpene and triterpenoid saponins, such as the oleanane-type triterpenoids in Aralia species and dammarane-type saponins in Panax, valued for their medicinal properties. The lack of genome sequences of Panax relatives has hindered mechanistic insight into the divergence of triterpene saponins in Araliaceae. Here, we report a chromosome-level genome of Aralia elata with a total length of 1.05 Gb. The loss of 12 exons in the dammarenediol synthase (DDS)-encoding gene in A. elata after divergence from Panax might have caused the lack of dammarane-type saponin production, and a complementation assay shows that overexpression of the PgDDS gene from Panax ginseng in callus of A. elata recovers the accumulation of dammarane-type saponins. Tandem duplication events of triterpene biosynthetic genes are common in the A. elata genome, especially for AeCYP72As, AeCSLMs, and AeUGT73s, which function as tailoring enzymes of oleanane-type saponins and aralosides. More than 13 aralosides are de novo synthesized in Saccharomyces cerevisiae by overexpression of these genes in combination. This study sheds light on the diversity of saponins biosynthetic pathway in Araliaceae and will facilitate heterologous bioproduction of aralosides.

Phytochemical constituents and protective efficacy of Schefflera arboricola L. leaves extract against thioacetamide-induced hepatic encephalopathy in rats

Context: Hepatic encephalopathy (HE) is a severe neuropsychiatric syndrome resulting from liver failure. Objective: To evaluate the protective effect of Schefflera arboricola L. leaves methanol extract against thioacetamide (TAA) induced HE in rats. Materials and methods: GC/MS, LC-ESI-MS and the total phenolic and flavonoid contents were determined. The methanol extract was orally administrated (100 and 200 mg/kg) for 21 days. TAA (200 mg/kg) was given intraperitoneally on day 19 and continued for three days. The evaluation was done by measuring alanine aminotransferases (ALT), alkaline phosphatase (ALP), ammonia, reduced glutathione (GSH), malondialdehyde (MDA), nitric oxide (NO) alpha tumor necrotic factor (TNFα), toll like receptor (TLR4), interleukin 1 beta (IL-1β), interlukin 6 (IL-6), cyclooxygenase 2(COX2), B cell lymphoma (BCL2), alpha smooth muscle actin (α-SMA) and cluster of differentiation 163 (CD163). The histological features of liver and brain were conducted. Results: Forty five compounds were identified from the n-hexane fraction, while twenty nine phenolic compounds were determined from the methanol extract. Pretreatment with the plant extract returned most of the measurements under investigation to nearly normal. Conclusion: Due to its richness with bioactive compounds, Schefflera arboricola L. leaves extract succeeded to exert anti-fibrotic, anti-inflammatory and antioxidants properties in TAA-induced HE in rats with more efficacy to its high protective dose.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Characterization of the chloroplast genome of Lonicera ruprechtiana Regel and comparison with other selected species of Caprifoliaceae

Lonicera ruprechtiana Regel is widely used as a greening tree in China and also displays excellent pharmacological activities. The phylogenetic relationship between L. ruprechtiana and other members of Caprifoliaceae remains unclear. In this study, the complete cp genome of L. ruprechtiana was identified using high-throughput Illumina pair-end sequencing data. The circular cp genome was 154,611 bp long and has a large single-copy region of 88,182 bp and a small single-copy region of 18,713 bp, with the two parts separated by two inverted repeat (IR) regions (23,858 bp each). A total of 131 genes were annotated, including 8 ribosomal RNAs, 39 transfer RNAs, and 84 protein-coding genes (PCGs). In addition, 49 repeat sequences and 55 simple sequence repeat loci of 18 types were also detected. Codon usage analysis demonstrated that the Leu codon is preferential for the A/U ending. Maximum-likelihood phylogenetic analysis using 22 Caprifoliaceae species revealed that L. ruprechtiana was closely related to Lonicera insularis. Comparison of IR regions revealed that the cp genome of L. ruprechtiana was largely conserved with that of congeneric species. Moreover, synonymous (Ks) and non-synonymous (Ka) substitution rate analysis showed that most genes were under purifying selection pressure; ycf3, and some genes associated with subunits of NADH dehydrogenase, subunits of the cytochrome b/f complex, and subunits of the photosystem had been subjected to strong purifying selection pressure (Ka/Ks < 0.1). This study provides useful genetic information for future study of L. ruprechtiana evolution.

Genome structure and diversity among Cynanchum wilfordii accessions

BACKGROUND

Cynanchum wilfordii (Cw) and Cynanchum auriculatum (Ca) have long been used in traditional medicine and as functional food in Korea and China, respectively. They have diverse medicinal functions, and many studies have been conducted, including pharmaceutical efficiency and metabolites. Especially, Cw is regarded as the most famous medicinal herb in Korea due to its menopausal symptoms relieving effect. Despite the high demand for Cw in the market, both species are cultivated using wild resources with rare genomic information.

RESULTS

We collected 160 Cw germplasm from local areas of Korea and analyzed their morphological diversity. Five Cw and one Ca of them, which were morphologically diverse, were sequenced, and nuclear ribosomal DNA (nrDNA) and complete plastid genome (plastome) sequences were assembled and annotated. We investigated the genomic characteristics of Cw as well as the genetic diversity of plastomes and nrDNA of Cw and Ca. The Cw haploid nuclear genome was approximately 178 Mbp. Karyotyping revealed the juxtaposition of 45S and 5S nrDNA on one of 11 chromosomes. Plastome sequences revealed 1226 interspecies polymorphisms and 11 Cw intraspecies polymorphisms. The 160 Cw accessions were grouped into 21 haplotypes based on seven plastome markers and into 108 haplotypes based on seven nuclear markers. Nuclear genotypes did not coincide with plastome haplotypes that reflect the frequent natural outcrossing events.

CONCLUSIONS

Cw germplasm had a huge morphological diversity, and their wide range of genetic diversity was revealed through the investigation with 14 molecular markers. The morphological and genomic diversity, chromosome structure, and genome size provide fundamental genomic information for breeding of undomesticated Cw plants.

Bacterial endophytes from ginseng and their biotechnological application

Ginseng has been well-known as a medicinal plant for thousands of years. Bacterial endophytes ubiquitously colonize the inside tissues of ginseng without any disease symptoms. The identification of bacterial endophytes is conducted through either the internal transcribed spacer region combined with ribosomal sequences or metagenomics. Bacterial endophyte communities differ in their diversity and composition profile, depending on the geographical location, cultivation condition, and tissue, age, and species of ginseng. Bacterial endophytes have a significant effect on the growth of ginseng through indole-3-acetic acid (IAA) and siderophore production, phosphate solubilization, and nitrogen fixation. Moreover, bacterial endophytes can protect ginseng by acting as biocontrol agents. Interestingly, bacterial endophytes isolated from Panax species have the potential to produce ginsenosides and bioactive metabolites, which can be used in the production of food and medicine. The ability of bacterial endophytes to transform major ginsenosides into minor ginsenosides using β-glucosidase is gaining increasing attention as a promising biotechnology. Recently, metabolic engineering has accelerated the possibilities for potential applications of bacterial endophytes in producing beneficial secondary metabolites.

Chromosomal-scale genome assembly of Eleutherococcus senticosus provides insights into chromosome evolution in Araliaceae

Siberian ginseng (Eleutherococcus senticosus, also known as ciwujia) belongs to the Araliaceae family, which contains more than 1,500 species in 41 genera with diverse chromosome numbers and genome sizes. General consensus posits that ancient whole-genome duplication events and rapid evolutionary radiation are the driving forces for this variation in genome properties. In an attempt to generate more genomic information for the Araliaceae family, we report a 1.30 Gb high-quality draft genome assembly (contig N50 = 309.43 kb) of E. senticosus via PacBio long reads and Hi-C chromatin interaction maps. We found that transposable elements accounted for 72.82% of the genome and a total of 36,372 protein-coding genes were predicted. Comparative analyses of the E. senticosus, Panax notoginseng and Daucus carota genomes revealed a burst expansion of Tekay chromoviral elements in Araliaceae after its divergence with Apiaceae. We also found that E. senticosus underwent a lineage-specific whole-genome duplication event Es-α and a whole-genome duplication event Araliaceae-β that was probably shared by all Araliaceae species. Even though the rediploidization of the E. senticosus genome is evident, pathway analyses show that these two whole-genome duplication events may have contributed to the adaptation of E. senticosus to a cold environment. Taken together, the high-quality genome assembly of E. senticosus provides a valuable genomic resource for future research on the evolution of Araliaceae.

In Vivo Evaluation of Dendropanax morbifera Leaf Extract for Anti-Obesity and Cholesterol-Lowering Activity in Mice

Metabolic syndrome is a worldwide health problem, and obesity is closely related to type 2 diabetes, cardiovascular disease, hypertension, and cancer. According to WHO in 2018, the prevalence of obesity in 2016 tripled compared to 1975. D. morbifera reduces bad cholesterol and triglycerides levels in the blood and provides various antioxidant nutrients and germicidal sub-stances, as well as selenium, which helps to remove active oxygen. Moreover, D. morbifera is useful for treating cardiovascular diseases, hypertension, hyperlipidemia, and diabetes. Therefore, we study in vivo efficacy of D. morbifera to investigate the prevention effect of obesity and cholesterol. The weight and body fat were effectively reduced by D. morbifera water (DLW) extract administration to high-fat diet-fed C57BL/6 mice compared to those of control mice. The group treated with DLW 500 mg∙kg-1∙d-1 had significantly lower body weights compared to the control group. In addition, High-density lipoprotein (HDL) cholesterol increased in the group treated with DLW 500 mg∙kg-1∙d-1. The effect of DLW on the serum lipid profile could be helpful to prevent obesity. DLW suppresses lipid formation in adipocytes and decreases body fat. In conclusion, DLW can be applied to develop anti-obesity functional foods and other products to reduce body fat.

Diversity and authentication of Rubus accessions revealed by complete plastid genome and rDNA sequences

Complete plastid genome (plastome) and ribosomal DNA (rDNA) sequences of three Rubus accessions (two Rubus longisepalus and one R. hirsutus) were newly assembled using Illumina whole-genome sequences. Rubus longisepalus Nakai and R. longisepalus var. tozawai, described as different varieties, have identical plastomes and rDNA sequences. The plastomes are 155,957 bp and 156,005 bp and the 45S rDNA transcription unit sizes are 5809 bp and 5811 bp in R. longisepalus and R. hirsutus, respectively. The 5S rDNA transcription unit is an identical 121 bp in three Rubus accessions. We developed three DNA markers to authenticate R. longisepalus and R. hirsutus based on plastome diversity. Phylogenomic analysis revealed that the Rubus species classified as two clades and R. longisepalus, R. hirsutus, and R. chingii are the most closely related species in clade 1.

Dynamic evolution of Panax species

Background

Panax ginseng is one of the most valuable medicinal plants in Korea. However, deciphering its full genome sequence information for crop improvement has been hampered due to its complex genomic, genetic, and growth characteristics. Many efforts have been made in the past decade to overcome these limitations and understand the genome structure and the evolutionary history of P. ginseng.

Methods

This review aims to discuss the current status of genomic studies on P. ginseng and related species, and the experimental clues suggesting phylogenetic classification and evolutionary history of the genus Panax.

Conclusion

The development of sequencing technologies made genome sequencing of the large P. ginseng genome possible, providing fundamental information to deciphering the evolutionary history of P. ginseng and related species. P. ginseng went through two rounds of whole genome duplication events after diverging from the closest family Apiaceae, which was unveiled from complete whole genome sequences. Further in-depth comparative genome analysis with other related species and genera will uncover the evolutionary history as well as important morphological and ecological characteristics of Panax species.

Inheritance of chloroplast and mitochondrial genomes in cucumber revealed by four reciprocal F1 hybrid combinations

Both genomes in chloroplasts and mitochondria of plant cell are usually inherited from maternal parent, with rare exceptions. To characterize the inheritance patterns of the organelle genomes in cucumber (Cucumis sativus var. sativus), two inbred lines and their reciprocal F₁ hybrids were analyzed using an next generation whole genome sequencing data. Their complete chloroplast genome sequences were de novo assembled, and a single SNP was identified between the parental lines. Two reciprocal F₁ hybrids have the same chloroplast genomes with their maternal parents. Meanwhile, 292 polymorphic sites were identified between mitochondrial genomes of the two parental lines, which showed the same genotypes with their paternal parents in the two reciprocal F₁ hybrids, without any recombination. The inheritance patterns of the chloroplast and mitochondria genomes were also confirmed in four additional cucumber accessions and their six reciprocal F₁ hybrids using molecular markers derived from the identified polymorphic sites. Taken together, our results indicate that the cucumber chloroplast genome is maternally inherited, as is typically observed in other plant species, whereas the large cucumber mitochondrial genome is paternally inherited. The combination of DNA markers derived from the chloroplast and mitochondrial genomes will provide a convenient system for purity test of F₁ hybrid seeds in cucumber breeding.

Assessing the genetic and chemical diversity of Taraxacum species in the Korean Peninsula

The genetic relationship between Taraxacum species, also known as the dandelion, is complicated because of asexual and mixed sexual apomictic reproduction. The usage of Taraxacum species in traditional medicines make their specialized metabolism important, but interspecific chemical difference has rarely been reported for the genus. In this study, we assembled the chloroplast genome and 45S rDNA of six Taraxacum species that occur in Korea (T. campylodes, T. coreanum, T. erythrospermum, T. mongolicum, T. platycarpum, and T. ussuriense), and performed a comparative analysis, which revealed their phylogenetic relationships and possible natural hybridity. We also performed a liquid chromatography-mass spectrometry-based phytochemical analysis to reveal interspecific chemical diversity. The comparative metabolomics analysis revealed that Taraxacum species could be separated into three chemotypes according to their major defensive specialized metabolites, which were the sesquiterpene lactones, the phenolic inositols, and chlorogenic acid derivatives. The CP DNA- and 45S rDNA-based phylogenetic trees showed a tangled relationship, which supports the notion of ongoing hybridization of wild Taraxacum species. The untargeted LC-MS analysis revealed that each Taraxacum plant exhibits species-specific defensive specialized metabolism. Moreover, 45S rDNA-based phylogenetic tree correlated with the hierarchical cluster relied on metabolite compositions. Given the coincidence between these analyses, we represented that 45S rDNA could well reflect overall nuclear genome variation in Taraxacum species.

Dendropanax Morbiferus and Other Species from the Genus Dendropanax: Therapeutic Potential of Its Traditional Uses, Phytochemistry, and Pharmacology

The Dendropanax genus is a kind of flowering plant in the family of Araliaceae that encompasses approximately 91 to 95 species. Several Dendropanax species are used as traditional medicinal plants, extensively used Korea and South America and other parts of the world. Almost every part of the plant, including the leaves, bark, roots, and stems, can be used as traditional medicine for the prevention and management of a broad spectrum of health disorders. This paper sought to summarizes the ethnopharmacological benefits, biological activities, and phytochemical investigations of plants from the genus Dendropanax, and perhaps to subsequently elucidate potential new perspectives for future pharmacological research to consider. Modern scientific literature suggests that plants of the Dendropanax genus, together with active compounds isolated from it, possess a wide range of therapeutic and pharmacological applications, including antifungal, anti-complement, antioxidant, antibacterial, insect antifeedant, cytotoxic, anti-inflammatory, neuroprotective, anti-diabetic, anti-cancer, and anti-hypouricemic properties. The botanical descriptions of approximately six to 10 species are provided by different scientific web sources. However, only six species, namely, D. morbiferus, D. gonatopodus, D. dentiger, D. capillaris, D. chevalieri, and D. arboreus, were included in the present investigation to undergo phytochemical evaluation, due to the unavailability of data for the remaining species. Among these plant species, a high concentration of variable bioactive ingredients was identified. In particular, D. morbifera is a traditional medicinal plant used for the multiple treatment purposes and management of several human diseases or health conditions. Previous experimental evidence supports that the D. morbifera species could be used to treat various inflammatory disorders, diarrhea, diabetes, cancer, and some microbial infections. It has recently been reported, by our group and other researchers, that D. morbifera possesses a neuroprotective and memory-enhancing agent. A total of 259 compounds have been identified among six species, with 78 sourced from five of these species reported to be bioactive. However, there is no up-to-date information concerning the D. morbifera, its different biological properties, or its prospective benefits in the enhancement of human health. In the present study, we set out to conduct a comprehensive analysis of the botany, traditional medicinal history, and medicinal resources of species of the Dendropanax genus. In addition, we explore several phytochemical constituents identified in different species of the Dendropanax genus and their biological properties. Finally, we offer comprehensive analysis findings of the phytochemistry, medicinal uses, pharmacological actions, and a toxicity and safety evaluation of the D. morbifera species and its main bioactive ingredients for future consideration.

Till 2018: a survey of biomolecular sequences in genus Panax

Ginseng is popularly known to be the king of ancient medicines and is used widely in most of the traditional medicinal compositions due to its various pharmaceutical properties. Numerous studies are being focused on this plant's curative effects to discover their potential health benefits in most human diseases, including cancer- the most life-threatening disease worldwide. Modern pharmacological research has focused mainly on ginsenosides, the major bioactive compounds of ginseng, because of their multiple therapeutic applications. Various issues on ginseng plant development, physiological processes, and agricultural issues have also been studied widely through state-of-the-art, high-throughput sequencing technologies. Since the beginning of the 21st century, the number of publications on ginseng has rapidly increased, with a recent count of more than 6,000 articles and reviews focusing notably on ginseng. Owing to the implementation of various technologies and continuous efforts, the ginseng plant genomes have been decoded effectively in recent years. Therefore, this review focuses mainly on the cellular biomolecular sequences in ginseng plants from the perspective of the central molecular dogma, with an emphasis on genomes, transcriptomes, and proteomes, together with a few other related studies.

The complete chloroplast genome of the Lonicera maackii (Caprifoliaceae), an ornamental plant

Lonicera maackii, is scattered in west and northeast China as well as adjacent Korea, Japan and the Soviet union. Here, we assembled and characterized the complete chloroplast (cp) genome of L. maackii using Illumina sequencing data for the first time. The complete cp genome was 155,337 bp in length, consisting of a pair of inverted repeats of 23,718 bp, a large single-copy region of 89,221 bp and a small single-copy region of 18,680 bp. The genome encoded 113 unique genes, including 79 protein-coding genes, 30 tRNA genes and four rRNA genes. Phylogenetic analysis based on 25 complete cp genome sequences indicated that L. maackii is closely related to Lonicera sachalinensis and Lonicera insularis.

Comprehensive comparative analysis of chloroplast genomes from seven Panax species and development of an authentication system based on species-unique single nucleotide polymorphism markers

BACKGROUND

Panax species are important herbal medicinal plants in the Araliaceae family. Recently, we reported the complete chloroplast genomes and 45S nuclear ribosomal DNA sequences from seven Panax species, two (P . quinqu e folius and P . trifolius) from North America and five (P . ginseng, P . notoginseng, P . japonicus, P . vietnamensis, and P . stipuleanatus) from Asia.

METHODS

We conducted phylogenetic analysis of these chloroplast sequences with 12 other Araliaceae species and comprehensive comparative analysis among the seven Panax whole chloroplast genomes.

RESULTS

We identified 1,128 single nucleotide polymorphisms (SNP) in coding gene sequences, distributed among 72 of the 79 protein-coding genes in the chloroplast genomes of the seven Panax species. The other seven genes (including psaJ, psbN, rpl23, psbF, psbL, rps18, and rps7) were identical among the Panax species. We also discovered that 12 large chloroplast genome fragments were transferred into the mitochondrial genome based on sharing of more than 90% sequence similarity. The total size of transferred fragments was 60,331 bp, corresponding to approximately 38.6% of chloroplast genome. We developed 18 SNP markers from the chloroplast genic coding sequence regions that were not similar to regions in the mitochondrial genome. These markers included two or three species-specific markers for each species and can be used to authenticate all the seven Panax species from the others.

CONCLUSION

The comparative analysis of chloroplast genomes from seven Panax species elucidated their genetic diversity and evolutionary relationships, and 18 species-specific markers were able to discriminate among these species, thereby furthering efforts to protect the ginseng industry from economically motivated adulteration.

Three new diterpenoids from Aralia dumetorum

Three new diterpenoids, dumetoranes A (1) and B (2), melanocane B (3), together with four known ones including melanocane A (4), ent-15S,16-dihydroxypimar-8(14)-en-19-oic acid (5), ent-pimara-8(14),15-diene-19-oic acid (6), and ent-pimara-8(14),15-diene-19-ol (7) were obtained from the ethanol extract of the roots of Aralia dumetorum. Their structure elucidation was achieved by the methods of spectroscopic HRMS, IR, NMR, and by comparison with literature. The cytotoxicities of compounds 1-3 and 5 were assayed by in vitro MTT methods.

Pre-labelled oligo probe-FISH karyotype analyses of four Araliaceae species using rDNA and telomeric repeat

BACKGROUND

The family Araliaceae contains many medicinal species including ginseng of which the whole genome sequencing analyses have been going on these days.

OBJECTIVE

To characterize the chromosomal distribution of 5S and 45S rDNAs and telomeric repeat in four ginseng related species of Aralia elata (Miq.) Seem., Dendropanax morbiferus H. Lév., Eleutherococcus sessiliflorus (Rupr. Et Maxim.) Seem., Kalopanax septemlobus (Thunb. ex A.Murr.) Koidz.

METHOD

Pre-labelled oligoprobe (PLOP)-fluorescence in situ hybridization (FISH) was carried out.

RESULTS

The chromosome number of A. elata was 2n = 24, whereas that of the other three species of D. morbiferus, E. sessiliflorus, and K. septemlobus was 2n = 48, corresponding to diploid and tetraploid, respectively, based on the basic chromosome number x = 12 in Araliaceae. In all four species, one pair of 5S signals were detected in the proximal regions of the short arms of chromosome 3, whereas in K. septemlobus, the 5S rDNA signals localized in the subtelomeric region of short arm of chromosome 3, while all the 45S rDNA signals localized at the paracentromeric region of the short arm of chromosome 1. And the telomeric repeat signals were detected at the telomeric region of both short and long arms of most chromosomes.

CONCLUSION

The PLOP-FISH was very effective compared with conventional FISH method. These results provide useful comparative cytogenetic information to better understand the genome structure of each species and will be useful to trace the history of ginseng genomic constitution.

Comparative mitogenomic and phylogenetic analysis of Apalone spinifera and Apalone ferox (Testudines: Trionychidae)

The soft-shell turtles Apalone spinifera (AS) and Apalone ferox (AF) are two important economic species. AF is found in the Yellow River of China, and is a confirmed member of the Trionychidae family. However, the classification of AS was in dispute. Mitochondrial genomes (mitogenomes) have been widely used for species identification, as well as population and phylogenetic analysis. In order to understand the phylogenetic and mitogenomic features of AS and AF, the complete mitogenomes were sequenced, annotated and analyzed in this study. The complete mitogenomes of AS and AF are 16,817 bp and 16,756 bp in length, respectively. Both mitogenomes contain 37 genes, seven short intergenic spacers and two long intergenic spacers. Comparative analysis showed that there are 1,137 variation sites (6.79%) between the two mitogenomes. AS and AF mitogenomes both show a usage preference in terms of nucleotides, codons and amino acids. In addition, the non-synonymous substitution rate/synonymous substitution rate indicates that all protein-coding genes (PCGs) have undergone a strong purifying selection. Phylogenetic trees constructed by 13 PCGs show a clear phylogenetic relationship of the soft-shell turtles and suggest that AS is a sister species to AF of the genus Apalone. The data could be useful for further research of species identification, population analysis and the mitogenomic features of soft-shell turtles.

Chemical and genomic diversity of six Lonicera species occurring in Korea

Lonicera spp. (Caprifoliaceae) are important not only as a common medicinal herb in East Asia but also as one of the most problematic invasive species in North America. In the present study, we performed a systemic analysis of genomic and chemical diversity among six Lonicera species occurring in Korea, L. japonica, L. maackii, L. insularis, L. sachalinensis, L. praeflorens, and L. vesicaria, using chloroplast DNA whole genome shotgun (WGS) sequencing and LC-MS analyses. The phylogenetic and phylochemical relationships did not coincide with each other, but partial consistency could be found among them. InDel-based cDNA marker for authentication was developed based on the genome sequences. Flavonoids, iridoids, and organic acids were identified in the LC-MS analyses, and their inter-species distribution and localization were also revealed.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

Ginseng Genome Database: an open-access platform for genomics of Panax ginseng

BACKGROUND

The ginseng (Panax ginseng C.A. Meyer) is a perennial herbaceous plant that has been used in traditional oriental medicine for thousands of years. Ginsenosides, which have significant pharmacological effects on human health, are the foremost bioactive constituents in this plant. Having realized the importance of this plant to humans, an integrated omics resource becomes indispensable to facilitate genomic research, molecular breeding and pharmacological study of this herb.

DESCRIPTION

The first draft genome sequences of P. ginseng cultivar "Chunpoong" were reported recently. Here, using the draft genome, transcriptome, and functional annotation datasets of P. ginseng, we have constructed the Ginseng Genome Database http://ginsengdb.snu.ac.kr /, the first open-access platform to provide comprehensive genomic resources of P. ginseng. The current version of this database provides the most up-to-date draft genome sequence (of approximately 3000 Mbp of scaffold sequences) along with the structural and functional annotations for 59,352 genes and digital expression of genes based on transcriptome data from different tissues, growth stages and treatments. In addition, tools for visualization and the genomic data from various analyses are provided. All data in the database were manually curated and integrated within a user-friendly query page.

CONCLUSION

This database provides valuable resources for a range of research fields related to P. ginseng and other species belonging to the Apiales order as well as for plant research communities in general. Ginseng genome database can be accessed at http://ginsengdb.snu.ac.kr /.

Phylogenomics and barcoding of Panax: toward the identification of ginseng species

BACKGROUND

The economic value of ginseng in the global medicinal plant trade is estimated to be in excess of US$2.1 billion. At the same time, the evolutionary placement of ginseng (Panax ginseng) and the complex evolutionary history of the genus is poorly understood despite several molecular phylogenetic studies. In this study, we use a full plastome phylogenomic framework to resolve relationships in Panax and to identify molecular markers for species discrimination.

RESULTS

We used high-throughput sequencing of MBD2-Fc fractionated Panax DNA to supplement publicly available plastid genomes to create a phylogeny based on fully assembled and annotated plastid genomes from 60 accessions of 8 species. The plastome phylogeny based on a 163 kbp matrix resolves the sister relationship of Panax ginseng with P. quinquefolius. The closely related species P. vietnamensis is supported as sister of P. japonicus. The plastome matrix also shows that the markers trnC-rps16, trnS-trnG, and trnE-trnM could be used for unambiguous molecular identification of all the represented species in the genus.

CONCLUSIONS

MBD2 depletion reduces the cost of plastome sequencing, which makes it a cost-effective alternative to Sanger sequencing based DNA barcoding for molecular identification. The plastome phylogeny provides a robust framework that can be used to study the evolution of morphological characters and biosynthesis pathways of ginsengosides for phylogenetic bioprospecting. Molecular identification of ginseng species is essential for authenticating ginseng in international trade and it provides an incentive for manufacturers to create authentic products with verified ingredients.

Rapid amplification of four retrotransposon families promoted speciation and genome size expansion in the genus Panax

Genome duplication and repeat multiplication contribute to genome evolution in plants. Our previous work identified a recent allotetraploidization event and five high-copy LTR retrotransposon (LTR-RT) families PgDel, PgTat, PgAthila, PgTork, and PgOryco in Panax ginseng. Here, using whole-genome sequences, we quantified major repeats in five Panax species and investigated their role in genome evolution. The diploids P. japonicus, P. vietnamensis, and P. notoginseng and the tetraploids P. ginseng and P. quinquefolius were analyzed alongside their relative Aralia elata. These species possess 0.8-4.9 Gb haploid genomes. The PgDel, PgTat, PgAthila, and PgTork LTR-RT superfamilies accounted for 39-52% of the Panax species genomes and 17% of the A. elata genome. PgDel included six subfamily members, each with a distinct genome distribution. In particular, the PgDel1 subfamily occupied 23-35% of the Panax genomes and accounted for much of their genome size variation. PgDel1 occupied 22.6% (0.8 Gb of 3.6 Gb) and 34.5% (1.7 Gb of 4.9 Gb) of the P. ginseng and P. quinquefolius genomes, respectively. Our findings indicate that the P. quinquefolius genome may have expanded due to rapid PgDel1 amplification over the last million years as a result of environmental adaptation following migration from Asia to North America.

Authentication Markers for Five Major Panax Species Developed via Comparative Analysis of Complete Chloroplast Genome Sequences

Ginseng represents a set of high-value medicinal plants of different species: Panax ginseng (Asian ginseng), Panax quinquefolius (American ginseng), Panax notoginseng (Chinese ginseng), Panax japonicus (Bamboo ginseng), and Panax vietnamensis (Vietnamese ginseng). Each species is pharmacologically and economically important, with differences in efficacy and price. Accordingly, an authentication system is needed to combat economically motivated adulteration of Panax products. We conducted comparative analysis of the chloroplast genome sequences of these five species, identifying 34-124 InDels and 141-560 SNPs. Fourteen InDel markers were developed to authenticate the Panax species. Among these, eight were species-unique markers that successfully differentiated one species from the others. We generated at least one species-unique marker for each of the five species, and any of the species can be authenticated by selection among these markers. The markers are reliable, easily detectable, and valuable for applications in the ginseng industry as well as in related research.