Major repeat components covering one-third of the ginseng (Panax ginseng C.A. Meyer) genome and evidence for allotetraploidy

Telomere-to-telomere reference genome for Panax ginseng highlights the evolution of saponin biosynthesis

Ginseng (Panax ginseng) is a representative of Chinese traditional medicine, also used worldwide, while the triterpene saponin ginsenoside is the most important effective compound within it. Ginseng is an allotetraploid, with complex genetic background, making the study of its metabolic evolution challenging. In this study, we assembled a telomere-to-telomere ginseng reference genome, constructed of 3.45 Gb with 24 chromosomes and 77 266 protein-coding genes. Additionally, the reference genome was divided into two subgenomes, designated as subgenome A and B. Subgenome A contains a larger number of genes, whereas subgenome B has a general expression advantage, suggesting that ginseng subgenomes experienced asymmetric gene loss with biased gene expression. The two subgenomes separated approximately 6.07 million years ago, and subgenome B shows the closest relation to Panax vietnamensis var. fuscidiscus. Comparative genomics revealed an expansion of gene families associated with ginsenoside biosynthesis in both ginseng subgenomes. Furthermore, both tandem duplications and proximal duplications play crucial roles in ginsenoside biosynthesis. We also screened functional genes identified in previous research and found that some of these genes located in colinear regions between subgenomes have divergence functions, revealing an unbalanced evolution in both subgenomes and the saponin biosynthesis pathway in ginseng. Our work provides important resources for future genetic studies and breeding programs of ginseng, as well as the biosynthesis of ginsenosides.

Beyond genome: Advanced omics progress of Panax ginseng

Ginseng is a perennial herb of the genus Panax in the family Araliaceae as one of the most important traditional medicine. Genomic studies of ginseng assist in the systematic discovery of genes related to bioactive ginsenosides biosynthesis and resistance to stress, which are of great significance in the conservation of genetic resources and variety improvement. The transcriptome reflects the difference and consistency of gene expression, and transcriptomics studies of ginseng assist in screening ginseng differentially expressed genes to further explore the powerful gene source of ginseng. Protein is the ultimate bearer of ginseng life activities, and proteomic studies of ginseng assist in exploring the biosynthesis and regulation of secondary metabolites like ginsenosides and the molecular mechanism of ginseng adversity adaptation at the overall level. In this review, we summarize the current status of ginseng research in genomics, transcriptomics and proteomics, respectively. We also discuss and look forward to the development of ginseng genome allele mapping, ginseng spatiotemporal, single-cell transcriptome, as well as ginseng post-translational modification proteome. We hope that this review will contribute to the in-depth study of ginseng and provide a reference for future analysis of ginseng from a systems biology perspective.

Cytokinin signaling promotes root secondary growth and bud formation in Panax ginseng

Background

Panax ginseng, one of the valuable perennial medicinal plants, stores numerous pharmacological substrates in its storage roots. Given its perennial growth habit, organ regeneration occurs each year, and cambium stem cell activity is necessary for secondary growth and storage root formation. Cytokinin (CK) is a phytohormone involved in the maintenance of meristematic cells for the development of storage organs; however, its physiological role in storage-root secondary growth remains unknown.

Methods

Exogenous CK was repeatedly applied to P. ginseng, and morphological and histological changes were observed. RNA-seq analysis was used to elucidate the transcriptional network of CK that regulates P. ginseng growth and development. The HISTIDINE KINASE 3 (PgHK3) and RESPONSE REGULATOR 2 (PgRR2) genes were cloned in P. ginseng and functionally analyzed in Arabidopsis as a two-component system involved in CK signaling.

Results

Phenotypic and histological analyses showed that CK increased cambium activity and dormant axillary bud formation in P. ginseng, thus promoting storage-root secondary growth and bud formation. The evolutionarily conserved two-component signaling pathways in P. ginseng were sufficient to restore CK signaling in the Arabidopsis ahk2/3 double mutant and rescue its growth defects. Finally, RNA-seq analysis of CK-treated P. ginseng roots revealed that plant-type cell wall biogenesis-related genes are tightly connected with mitotic cell division, cytokinesis, and auxin signaling to regulate CK-mediated P. ginseng development.

Conclusion

Overall, we identified the CK signaling-related two-component systems and their physiological role in P. ginseng. This scientific information has the potential to significantly improve the field-cultivation and biotechnology-based breeding of ginseng.

High-yield production of protopanaxadiol from sugarcane molasses by metabolically engineered Saccharomyces cerevisiae

Background

Ginsenosides are Panax plant-derived triterpenoid with wide applications in cardiovascular protection and immunity-boosting. However, the saponins content of Panax plants is fairly low, making it time-consuming and unsustainable by direct extraction. Protopanaxadiol (PPD) is a common precursor of dammarane-type saponins, and its sufficient supply is necessary for the efficient synthesis of ginsenoside.

Results

In this study, a combinational strategy was used for the construction of an efficient yeast cell factory for PPD production. Firstly, a PPD-producing strain was successfully constructed by modular engineering in Saccharomyces cerevisiae BY4742 at the multi-copy sites. Then, the INO2 gene, encoding a transcriptional activator of the phospholipid biosynthesis, was fine-tuned to promote the endoplasmic reticulum (ER) proliferation and improve the catalytic efficiency of ER-localized enzymes. To increase the metabolic flux of PPD, dynamic control, based on a carbon-source regulated promoter P_HXT1, was introduced to repress the competition of sterols. Furthermore, the global transcription factor UPC2-1 was introduced to sterol homeostasis and up-regulate the MVA pathway, and the resulting strain BY-V achieved a PPD production of 78.13 ± 0.38 mg/g DCW (563.60 ± 1.65 mg/L). Finally, sugarcane molasses was used as an inexpensive substrate for the first time in PPD synthesis. The PPD titers reached 1.55 ± 0.02 and 15.88 ± 0.65 g/L in shake flasks and a 5-L bioreactor, respectively. To the best of our knowledge, these results were new records on PPD production.

Conclusion

The high-level of PPD production in this study and the successful comprehensive utilization of low-cost carbon source -sugarcane molassesindicate that the constructed yeast cell factory is an excellent candidate strain for the production of high-value-added PPD and its derivativeswith great industrial potential.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01949-4.

Prolonged Exposure to High Temperature Inhibits Shoot Primary and Root Secondary Growth in Panax ginseng

High temperature is one of the most significant abiotic stresses reducing crop yield and quality by inhibiting plant growth and development. Global warming has recently increased the frequency of heat waves, which negatively impacts agricultural fields. Despite numerous studies on heat stress responses and signal transduction in model plant species, the molecular mechanism underlying thermomorphogenesis in Panax ginseng remains largely unknown. Here, we investigated the high temperature response of ginseng at the phenotypic and molecular levels. Both the primary shoot growth and secondary root growth of ginseng plants were significantly reduced at high temperature. Histological analysis revealed that these decreases in shoot and root growth were caused by decreases in cell elongation and cambium stem cell activity, respectively. Analysis of P. ginseng RNA-seq data revealed that heat-stress-repressed stem and root growth is closely related to changes in photosynthesis, cell wall organization, cell wall loosening, and abscisic acid (ABA) and jasmonic acid (JA) signaling. Reduction in both the light and dark reactions of photosynthesis resulted in defects in starch granule development in the storage parenchymal cells of the main tap root. Thus, by combining bioinformatics and histological analyses, we show that high temperature signaling pathways are integrated with crucial biological processes that repress stem and root growth in ginseng, providing novel insight into the heat stress response mechanism of P. ginseng.

Gibberellin Signaling Promotes the Secondary Growth of Storage Roots in Panax ginseng

Gibberellins (GAs) are an important group of phytohormones associated with diverse growth and developmental processes, including cell elongation, seed germination, and secondary growth. Recent genomic and genetic analyses have advanced our knowledge of GA signaling pathways and related genes in model plant species. However, functional genomics analyses of GA signaling pathways in Panax ginseng, a perennial herb, have rarely been carried out, despite its well-known economical and medicinal importance. Here, we conducted functional characterization of GA receptors and investigated their physiological roles in the secondary growth of P. ginseng storage roots. We found that the physiological and genetic functions of P. ginseng gibberellin-insensitive dwarf1s (PgGID1s) have been evolutionarily conserved. Additionally, the essential domains and residues in the primary protein structure for interaction with active GAs and DELLA proteins are well-conserved. Overexpression of PgGID1s in Arabidopsis completely restored the GA deficient phenotype of the Arabidopsis gid1a gid1c (atgid1a/c) double mutant. Exogenous GA treatment greatly enhanced the secondary growth of tap roots; however, paclobutrazol (PCZ), a GA biosynthetic inhibitor, reduced root growth in P. ginseng. Transcriptome profiling of P. ginseng roots revealed that GA-induced root secondary growth is closely associated with cell wall biogenesis, the cell cycle, the jasmonic acid (JA) response, and nitrate assimilation, suggesting that a transcriptional network regulate root secondary growth in P. ginseng. These results provide novel insights into the mechanism controlling secondary root growth in P. ginseng.

Modification of ginsenoside saponin composition via the CRISPR/Cas9-mediated knockout of protopanaxadiol 6-hydroxylase gene in Panax ginseng

Background

The roots of Panax ginseng contain two types of tetracyclic triterpenoid saponins, namely, protopanaxadiol (PPD)-type saponins and protopanaxatiol (PPT)-type saponins. In P. ginseng, the protopanaxadiol 6-hydroxylase (PPT synthase) enzyme catalyses protopanaxatriol (PPT) production from protopanaxadiol (PPD). In this study, we constructed homozygous mutant lines of ginseng by CRISPR/Cas9-mediated mutagenesis of the PPT synthase gene and obtained the mutant ginseng root lines having complete depletion of the PPT-type ginsenosides.

Methods

Two sgRNAs (single guide RNAs) were designed for target mutations in the exon sequences of the two PPT synthase genes (both PPTa and PPTg sequences) with the CRISPR/Cas9 system. Transgenic ginseng roots were generated through Agrobacterium-mediated transformation. The mutant lines were screened by ginsenoside analysis and DNA sequencing.

Result

Ginsenoside analysis revealed the complete depletion of PPT-type ginsenosides in three putative mutant lines (Cr4, Cr7, and Cr14). The reduction of PPT-type ginsenosides in mutant lines led to increased accumulation of PPD-type ginsenosides. The gene editing in the selected mutant lines was confirmed by targeted deep sequencing.

Conclusion

We have established the genome editing protocol by CRISPR/Cas9 system in P. ginseng and demonstrated the mutated roots producing only PPD-type ginsenosides by depleting PPT-type ginsenosides. Because the pharmacological activity of PPD-group ginsenosides is significantly different from that of PPT-group ginsenosides, the new type of ginseng mutant producing only PPD-group ginsenosides may have new pharmacological characteristics compared to wild-type ginseng. This is the first report to generate target-induced mutations for the modification of saponin biosynthesis in Panax species using CRISPR–Cas9 system.

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.

Growth and bioactive phytochemicals of Panax ginseng sprouts grown in an aeroponic system using plasma-treated water as the nitrogen source

Ginseng (Panax ginseng Meyer) sprouts are grown to whole plants in 20 to 25 days in a soil-less cultivation system and then used as a medicinal vegetable. As a nitrogen (N) source, plasma-treated water (PTW) has been used to enhance the seed germination and seedling growth of many crops but has not been investigated for its effects on ginseng sprouts. This study established an in-situ system for N-containing water production using plasma technology and evaluated the effects of the PTW on ginseng growth and its bioactive phytochemicals compared with those of an untreated control. The PTW became weakly acidic 30 min after the air discharge at the electrodes because of the formation of nitrate (NO₃^‒) and nitrite (NO₂^‒) in the water. The NO₃^‒ and NO₂^‒ in the PTW, together with potassium ions (K⁺), enhanced the shoot biomass of the ginseng sprout by 26.5% compared to the untreated control. The ginseng sprout grown in the PTW had accumulated more free amino acids and ginsenosides in the sprout at 25 days after planting. Therefore, PTW can be used as a liquid N fertilizer for P. ginseng growth and phytochemical accumulation during sprouting under aeroponic conditions.

Identification of 5S and 45S rDNA sites in Chrysanthemum species by using oligonucleotide fluorescence in situ hybridization (Oligo-FISH)

Fluorescence in situ hybridization (FISH) is a conventional method used to visualize the distribution of DNA elements within a genome. To examine the relationships within the Chrysanthemum genus, ribosomal DNA (rDNA), a popular cytogenetic marker, was utilized as a probe for FISH within this genus. Based on the genome data of Chrysanthemum nankingense, C. seticuspe and its allied genera in the Compositae(Asteraceae), we explored rDNA sequences to design oligonucleotide probes and perform oligonucleotide fluorescence in situ hybridization (Oligo-FISH) in eight Chrysanthemum accessions. The results showed that the majority of 5S rDNA signals were located in subterminal chromosome regions and that the number of 5S rDNA sites might be tightly associated with ploidy. For 45S rDNA sites, the number and intensity of signals differed from those of previously investigated Chrysanthemum resources. These findings may provide an optimally reliable method of examining the chromosome composition and structural variation of Chrysanthemum and its related species and allow researchers to understand the evolutionary history and phylogenetic relationships of Chrysanthemum.

Age-Dependent and Tissue-Specific Alterations in the rDNA Clusters of the Panax ginseng C. A. Meyer Cultivated Cell Lines

Long-term cultivation of Panax ginseng cell lines leads to a decreasing synthesis of the biologically active substances used in traditional medicine. To gain insight into the cellular mechanisms which may influence this process, we analyzed variations within the rDNA cluster of the Oriental ginseng cell lines. The cell lines were cultivated for 6 and 24 years; the number of nucleoli and chromosomes was analyzed. The complete 18S rDNA sequences were cloned and sequenced. The nucleotide polymorphism and phylogenetic relations of the sequences were analyzed, and the secondary structures for separate 18S rRNA regions were modeled. The 18S rDNA accumulated mutations during cell cultivation that correlate well with an increase in the number of chromosomes and nucleoli. The patterns of nucleotide diversity are culture-specific and the increasing polymorphism associates with cytosine methylation sites. The secondary structures of some 18S rRNA regions and their interaction can alter during cultivation. The phylogenetic tree topologies are particular for each cell line.The observed alterations in rDNA clusters are associated with a somaclonal variation, leading to changes in the pattern of intracellular synthesis during cell cultivation. The identified divergent rRNAs could provide additional gene expression regulation in P. ginseng cells by forming heterogeneous ribosomes.

Characterization of B-Genome Specific High Copy hAT MITE Families in Brassica nigra Genome

Miniature inverted-repeat transposable elements (MITEs) are non-autonomous class II transposons which have been shown to influence genome evolution. Brassica nigra L. (B-genome) is one of three Brassica diploids cultivated primarily as an oil crop, which harbors novel alleles important for breeding. Two new high copy hAT MITE families (BniHAT-1 and BniHAT-2) from the B-genome were characterized and their prevalence assessed in the genomes of the related diploids, rapa L. (A) and Brassica oleracea L. (C). Both novel MITE families were present at high copy numbers in the B-genome with 434 and 331 copies of BniHAT-1 and BniHAT-2, respectively. Yet less than 20 elements were identified in the genome assemblies of the A, and C -genomes, supporting B-genome specific proliferation of these MITE families. Although apparently randomly distributed across the genome, 68 and 70% of the B-genome MITEs were present within 2 kb flanking regions of annotated genes suggesting they might influence gene expression and/or function. In addition, MITE derived microRNAs and transcription factor binding sites suggested a putative role in gene regulation. Age of insertion analysis revealed that the major proliferation of these elements occurred during 2-3 million years ago. Additionally, site-specific polymorphism analyses showed that 44% MITEs were undergoing active amplification into the B-genome. Overall, this study provides a comprehensive analysis of two high copy MITE families, which were specifically amplified in the B-genome, suggesting a potential role in shaping the Brassica B-genome.

Genetic diversity among cultivated and wild Panax ginseng populations revealed by high-resolution microsatellite markers

Background

Ginseng (Panax ginseng Meyer) is one of the world's most valuable medicinal plants with numerous pharmacological effects. Ginseng has been cultivated from wild mountain ginseng collections for a few hundred years. However, the genetic diversity of cultivated and wild ginseng populations is not fully understood.

Methods

We developed 92 polymorphic microsatellite markers based on whole-genome sequence data. We selected five markers that represent clear allele diversity for each of their corresponding loci to elucidate genetic diversity. These markers were applied to 147 individual plants, including cultivars, breeding lines, and wild populations in Korea and neighboring countries.

Results

Most of the 92 markers displayed multiple-band patterns, resulting from genome duplication, which causes confusion in interpretation of their target locus. The five high-resolution markers revealed 3 to 8 alleles from each single locus. The proportion of heterozygosity (H_e) ranged from 0.027 to 0.190, with an average of 0.132, which is notably lower than that of previous studies. Polymorphism information content of the markers ranged from 0.199 to 0.701, with an average of 0.454. There was no statistically significant difference in genetic diversity between cultivated and wild ginseng groups, and they showed intermingled positioning in the phylogenetic relationship.

Conclusion

Ginseng has a relatively high level of genetic diversity, and cultivated and wild groups have similar levels of genetic diversity. Collectively, our data demonstrate that current breeding populations have abundant genetic diversity for breeding of elite ginseng cultivars.

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine protein-coding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.

Discovery and annotation of a novel transposable element family in Gossypium

BACKGROUND

Fluorescence in situ hybridization (FISH) is an efficient cytogenetic technology to study chromosome structure. Transposable element (TE) is an important component in eukaryotic genomes and can provide insights in the structure and evolution of eukaryotic genomes.

RESULTS

A FISH probe derived from bacterial artificial chromosome (BAC) clone 299N22 generated striking signals on all 26 chromosomes of the cotton diploid A genome (AA, 2x=26) but very few on the diploid D genome (DD, 2x=26). All 26 chromosomes of the A sub genome (At) of tetraploid cotton (AADD, 2n=4x=52) also gave positive signals with this FISH probe, whereas very few signals were observed on the D sub genome (Dt). Sequencing and annotation of BAC clone 299N22, revealed a novel Ty3/gypsy transposon family, which was named as 'CICR'. This family is a significant contributor to size expansion in the A (sub) genome but not in the D (sub) genome. Further FISH analysis with the LTR of CICR as a probe revealed that CICR is lineage-specific, since massive repeats were found in A and B genomic groups, but not in C-G genomic groups within the Gossypium genus. Molecular evolutionary analysis of CICR suggested that tetraploid cottons evolved after silence of the transposon family 1-1.5 million years ago (Mya). Furthermore, A genomes are more homologous with B genomes, and the C, E, F, and G genomes likely diverged from a common ancestor prior to 3.5-4 Mya, the time when CICR appeared. The genomic variation caused by the insertion of CICR in the A (sub) genome may have played an important role in the speciation of organisms with A genomes.

CONCLUSIONS

The CICR family is highly repetitive in A and B genomes of Gossypium, but not amplified in the C-G genomes. The differential amount of CICR family in At and Dt will aid in partitioning sub genome sequences for chromosome assemblies during tetraploid genome sequencing and will act as a method for assessing the accuracy of tetraploid genomes by looking at the proportion of CICR elements in resulting pseudochromosome sequences. The timeline of the expansion of CICR family provides a new reference for cotton evolutionary analysis, while the impact on gene function caused by the insertion of CICR elements will be a target for further analysis of investigating phenotypic differences between A genome and D genome species.

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.

Rerouting of NADPH synthetic pathways for increased protopanaxadiol production in Saccharomyces cerevisiae

Ginseng (Panax ginseng) and its bioactive components, ginsenosides, are popular medicinal herbal products, exhibiting various pharmacological effects. Despite their advocated use for medication, the long cultivation periods of ginseng roots and their low ginsenoside content prevent mass production of this compound. Yeast Saccharomyces cerevisiae was engineered for production of protopanaxadiol (PPD), a type of aglycone characterizing ginsenoside. PPD-producing yeast cell factory was further engineered by obtaining a balance between enzyme expressions and altering cofactor availability. Different combinations of promoters (PGPD, PCCW12, and PADH2) were utilized to construct the PPD biosynthetic pathway. Rerouting the redox metabolism to improve NADPH availability in the engineered S. cerevisiae also increased PPD production. Combining these approaches resulted in more than an 11-fold increase in PPD titer over the initially constructed strain. The series of metabolic engineering strategies of this study provides a feasible approach for the microbial production of PPD and development of microbial platforms producing other industrially-relevant terpenoids.

Panax ginseng genome examination for ginsenoside biosynthesis

Ginseng, which contains ginsenosides as bioactive compounds, has been regarded as an important traditional medicine for several millennia. However, the genetic background of ginseng remains poorly understood, partly because of the plant's large and complex genome composition. We report the entire genome sequence of Panax ginseng using next-generation sequencing. The 3.5-Gb nucleotide sequence contains more than 60% repeats and encodes 42 006 predicted genes. Twenty-two transcriptome datasets and mass spectrometry images of ginseng roots were adopted to precisely quantify the functional genes. Thirty-one genes were identified to be involved in the mevalonic acid pathway. Eight of these genes were annotated as 3-hydroxy-3-methylglutaryl-CoA reductases, which displayed diverse structures and expression characteristics. A total of 225 UDP-glycosyltransferases (UGTs) were identified, and these UGTs accounted for one of the largest gene families of ginseng. Tandem repeats contributed to the duplication and divergence of UGTs. Molecular modeling of UGTs in the 71st, 74th, and 94th families revealed a regiospecific conserved motif located at the N-terminus. Molecular docking predicted that this motif captures ginsenoside precursors. The ginseng genome represents a valuable resource for understanding and improving the breeding, cultivation, and synthesis biology of this key herb.

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 /.

Positive Selection Driving Cytoplasmic Genome Evolution of the Medicinally Important Ginseng Plant Genus Panax

Panax L. (the ginseng genus) is a shade-demanding group within the family Araliaceae and all of its species are of crucial significance in traditional Chinese medicine. Phylogenetic and biogeographic analyses demonstrated that two rounds of whole genome duplications accompanying with geographic and ecological isolations promoted the diversification of Panax species. However, contributions of the cytoplasmic genomes to the adaptive evolution of Panax species remained largely uninvestigated. In this study, we sequenced the chloroplast and mitochondrial genomes of 11 accessions belonging to seven Panax species. Our results show that heterogeneity in nucleotide substitution rate is abundant in both of the two cytoplasmic genomes, with the mitochondrial genome possessing more variants at the total level but the chloroplast showing higher sequence polymorphisms at the genic regions. Genome-wide scanning of positive selection identified five and 12 genes from the chloroplast and mitochondrial genomes, respectively. Functional analyses further revealed that these selected genes play important roles in plant development, cellular metabolism and adaptation. We therefore conclude that positive selection might be one of the potential evolutionary forces that shaped nucleotide variation pattern of these Panax species. In particular, the mitochondrial genes evolved under stronger selective pressure compared to the chloroplast genes.

The mRNA and miRNA transcriptomic landscape of Panax ginseng under the high ambient temperature

BACKGROUND

Ginseng is a popular traditional herbal medicine in north-eastern Asia. It has been used for human health for over thousands of years. With the rise in global temperature, the production of Korean ginseng (Panax ginseng C.A.Meyer) in Korea have migrated from mid to northern parts of the Korean peninsula to escape from the various higher temperature related stresses. Under the high ambient temperature, vegetative growth was accelerated, which resulted in early flowering. This precocious phase change led to yield loss. Despite of its importance as a traditional medicine, biological mechanisms of ginseng has not been well studied and even the genome sequence of ginseng is yet to be determined due to its complex genome structure. Thus, it is challenging to investigate the molecular biology mechanisms at the transcript level.

RESULTS

To investigate how ginseng responds to the high ambient temperature environment, we performed high throughput RNA sequencing and implemented a bioinformatics pipeline for the integrated analysis of small-RNA and mRNA-seq data without a reference genome. By performing reverse transcriptase (RT) PCR and sanger sequencing of transcripts that were assembled using our pipeline, we validated that their sequences were expressed in our samples. Furthermore, to investigate the interaction between genes and non-coding small RNAs and their regulation status under the high ambient temperature, we identified potential gene regulatory miRNAs. As a result, 100,672 contigs with significant expression level were identified and 6 known, 214 conserved and 60 potential novel miRNAs were predicted to be expressed under the high ambient temperature.

CONCLUSION

Collectively, we have found that development, flowering and temperature responsive genes were induced under high ambient temperature, whereas photosynthesis related genes were repressed. Functional miRNAs were down-regulated under the high ambient temperature. Among them are miR156 and miR396 that target flowering (SPL6/9) and growth regulating genes (GRF) respectively.

Genome-Wide Variation Patterns Uncover the Origin and Selection in Cultivated Ginseng (Panax ginseng Meyer)

Chinese ginseng (Panax ginseng Meyer) is a medicinally important herb and plays crucial roles in traditional Chinese medicine. Pharmacological analyses identified diverse bioactive components from Chinese ginseng. However, basic biological attributes including domestication and selection of the ginseng plant remain under-investigated. Here, we presented a genome-wide view of the domestication and selection of cultivated ginseng based on the whole genome data. A total of 8,660 protein-coding genes were selected for genome-wide scanning of the 30 wild and cultivated ginseng accessions. In complement, the 45s rDNA, chloroplast and mitochondrial genomes were included to perform phylogenetic and population genetic analyses. The observed spatial genetic structure between northern cultivated ginseng (NCG) and southern cultivated ginseng (SCG) accessions suggested multiple independent origins of cultivated ginseng. Genome-wide scanning further demonstrated that NCG and SCG have undergone distinct selection pressures during the domestication process, with more genes identified in the NCG (97 genes) than in the SCG group (5 genes). Functional analyses revealed that these genes are involved in diverse pathways, including DNA methylation, lignin biosynthesis, and cell differentiation. These findings suggested that the SCG and NCG groups have distinct demographic histories. Candidate genes identified are useful for future molecular breeding of cultivated ginseng.

Elucidating the major hidden genomic components of the A, C, and AC genomes and their influence on Brassica evolution

Decoding complete genome sequences is prerequisite for comprehensive genomics studies. However, the currently available reference genome sequences of Brassica rapa (A genome), B. oleracea (C) and B. napus (AC) cover 391, 540, and 850 Mbp and represent 80.6, 85.7, and 75.2% of the estimated genome size, respectively, while remained are hidden or unassembled due to highly repetitive nature of these genome components. Here, we performed the first comprehensive genome-wide analysis using low-coverage whole-genome sequences to explore the hidden genome components based on characterization of major repeat families in the B. rapa and B. oleracea genomes. Our analysis revealed 10 major repeats (MRs) including a new family comprising about 18.8, 10.8, and 11.5% of the A, C and AC genomes, respectively. Nevertheless, these 10 MRs represented less than 0.7% of each assembled reference genome. Genomic survey and molecular cytogenetic analyses validates our insilico analysis and also pointed to diversity, differential distribution, and evolutionary dynamics in the three Brassica species. Overall, our work elucidates hidden portions of three Brassica genomes, thus providing a resource for understanding the complete genome structures. Furthermore, we observed that asymmetrical accumulation of the major repeats might be a cause of diversification between the A and C genomes.

Isoform Sequencing Provides a More Comprehensive View of the Panax ginseng Transcriptome

Korean ginseng (Panax ginseng C.A. Meyer) has been widely used for medicinal purposes and contains potent plant secondary metabolites, including ginsenosides. To obtain transcriptomic data that offers a more comprehensive view of functional genomics in P. ginseng, we generated genome-wide transcriptome data from four different P. ginseng tissues using PacBio isoform sequencing (Iso-Seq) technology. A total of 135,317 assembled transcripts were generated with an average length of 3.2 kb and high assembly completeness. Of those unigenes, 67.5% were predicted to be complete full-length (FL) open reading frames (ORFs) and exhibited a high gene annotation rate. Furthermore, we successfully identified unique full-length genes involved in triterpenoid saponin synthesis and plant hormonal signaling pathways, including auxin and cytokinin. Studies on the functional genomics of P. ginseng seedlings have confirmed the rapid upregulation of negative feed-back loops by auxin and cytokinin signaling cues. The conserved evolutionary mechanisms in the auxin and cytokinin canonical signaling pathways of P. ginseng are more complex than those in Arabidopsis thaliana. Our analysis also revealed a more detailed view of transcriptome-wide alternative isoforms for 88 genes. Finally, transposable elements (TEs) were also identified, suggesting transcriptional activity of TEs in P. ginseng. In conclusion, our results suggest that long-read, full-length or partial-unigene data with high-quality assemblies are invaluable resources as transcriptomic references in P. ginseng and can be used for comparative analyses in closely related medicinal plants.

Analysis of the Panax ginseng stem/leaf transcriptome and gene expression during the leaf expansion period

Ginseng (Panax ginseng C.A Meyer) is a widely used herbal remedy, however, the majority of studies have focused on the roots, with less known about the aerial regions of the plant. As the stems and leaves are the primary aerial tissues, the present study characterized their transcriptional profiles using Illumina next‑generation sequencing technology. The gene expression profiles and the functional genes of ginseng stems (GS) and leaves (GL) were analyzed during the leaf‑expansion period. cDNA libraries of the GS and GL of 5‑year‑old ginseng plants were separately constructed. In the GS library, 38,000,000 sequencing reads were produced. These reads were assembled into 99,809 unique sequences with a mean size of 572 bp, and 57,371 sequences were identified based on similarity searches against known proteins. The assembled sequences were annotated using Gene Ontology terms, Clusters of Orthologous Groups classifications and Kyoto Encyclopedia of Genes and Genomes pathways. For GL, >118,000,000 sequencing reads were produced, which were assembled into 73,163 unique sequences, from which 50,523 sequences were identified. Additionally, several genes involved in the regulation of growth‑related, stress‑related, pathogenesis‑related, and chlorophyll metabolism‑associated proteins were found and expressed at high levels, with low expression levels of ginsenoside biosynthesis enzymes also found. The results of the present study provide a valuable useful sequence resource for ginseng in general, and specifically for further investigations of the functional genomics and molecular genetics of GS and GL during early growth.

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.

Intercontinental and intracontinental biogeography of the eastern Asian - Eastern North American disjunct Panax (the ginseng genus, Araliaceae), emphasizing its diversification processes in eastern Asia

The intercontinental biogeography between eastern Asia and eastern North America has attracted much attention from evolutionary biologists. Further insights into understanding the evolution of the intercontinental disjunctions have been hampered by the lack of studies on the intracontinental biogeography in eastern Asia, a region with complex geology, geography, climates and habitats. Herein we studied the biogeographic history of the eastern Asian-eastern North American disjunct genus Panax with special emphasis on the investigation of its uneven diversification in Asia. This study reconstructs the diversification history of Panax and also emphasizes a large clade of Panax taxa, which has a wide distribution in eastern Asia, but was unresolved in previous studies. We examined the noncoding plastid DNA fragments of trnH-psbA, rps16, and psbM-trnD, the mitochondrial b/c intron of NAD1, and the nuclear ribosomal internal transcribed spacer (ITS) region of 356 samples from 47 populations. The results revealed the subtropical Northern Hemisphere origin (Asia or Asia and North America) of Panax in the Paleocene. Intercontinental disjunctions between eastern Asia and eastern North America formed twice in Panax, once estimated in early Eocene for the split of P. trifolius and another in mid-Miocene for the divergence of P. quinquefolius. Intercontinental diversifications in Panax showed temporal correlation with the increase of global temperature. The evolutionary radiation of the P. bipinnatifidus species complex occurred around the boundary of Oligocene and Miocene. Strong genetic structure among populations of the species complex was detected and the populations may be isolated by distance. The backbone network and the Bayesian clustering analysis revealed a major evolutionary radiation centered in the Hengduan Mountains of western China. Our results suggested that the evolutionary radiation of Panax was promoted by geographic barriers, including mountain ranges (Hengduan Mountains, Nanling Mountains and Wuyishan Mountains), oceans and altitudinal shifts, which further contribute to the knowledge of the uneven species diversification between eastern Asia and North America.

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

We produced complete sequences and conducted comparative analysis of the maternally inherited chloroplast (cp) genomes and bi-parentally inherited 45S nuclear ribosomal RNA genes (nrDNA) from ten Araliaceae species to elucidate the genetic diversity and evolution in that family. The cp genomes ranged from 155,993 bp to 156,730 bp with 97.1-99.6% similarity. Complete 45S nrDNA units were about 11 kb including a 5.8-kb 45S cistron. Among 79 cp protein-coding genes, 74 showed nucleotide variations among ten species, of which infA, rpl22, rps19 and ndhE genes showed the highest Ks values and atpF, atpE, ycf2 and rps15 genes showed the highest Ka/Ks values. Four genes, petN, psaJ, psbF, and psbN, related to photosynthesis and one gene, rpl23, related to the ribosomal large subunit remain conserved in all 10 Araliaceae species. Phylogenetic analysis revealed that the ten species could be resolved into two monophyletic lineages, the Panax-Aralia and the Eleutherococcus-Dendropanax groups, which diverged approximately 8.81-10.59 million years ago (MYA). The Panax genus divided into two groups, with diploid species including P. notoginseng, P. vietnamensis, and P. japonicus surviving in Southern Asia and a tetraploid group including P. ginseng and P. quinquefolius Northern Asia and North America 2.89-3.20 MYA.

A refined Panax ginseng karyotype based on an ultra-high copy 167-bp tandem repeat and ribosomal DNAs

Background

Panax ginseng Meyer (Asian ginseng) has a large nuclear genome size of > 3.5 Gbp in haploid genome equivalent of 24 chromosomes. Tandem repeats (TRs) occupy significant portions of the genome in many plants and are often found in specific genomic loci, making them a valuable molecular cytogenetic tool in discriminating chromosomes. In an effort to understand the P. ginseng genome structure, we characterized an ultrahigh copy 167-bp TR (Pg167TR) and explored its chromosomal distribution as well as its utility for chromosome identification.

Methods

Polymerase chain reaction amplicons of Pg167TR were labeled, along with 5S and 45S rDNA amplicons, using a direct nick-translation method. Direct fluorescence in situ hybridization (FISH) was used to analyze the chromosomal distribution of Pg167TR.

Results

Recently, we reported a method of karyotyping the 24 chromosome pairs of P. ginseng using rDNA and DAPI (4′,6-diamidino-2-phenylindole) bands. Here, a unique distribution of Pg167TR in all 24 P. ginseng chromosomes was observed, allowing easy identification of individual homologous chromosomes. Additionally, direct labeling of 5S and 45S rDNA probes allowed the identification of two additional 5S rDNA loci not previously reported, enabling the refinement of the P. ginseng karyotype.

Conclusion

Identification of individual P. ginseng chromosomes was achieved using Pg167TR-FISH. Chromosome identification is important in understanding the P. ginseng genome structure, and our method will be useful for future integration of genetic linkage maps and genome scaffold anchoring. Additionally, it is a good tool for comparative studies with related species in efforts to understand the evolution of P. ginseng.

Cytohistological study of the leaf structures of Panax ginseng Meyer and Panax quinquefolius L

Background

Both Panax ginseng Meyer and Panax quinquefolius are obligate shade-loving plants whose natural habitats are broadleaved forests of Eastern Asia and North America. Panax species are easily damaged by photoinhibition when they are exposed to high temperatures or insufficient shade. In this study, a cytohistological study of the leaf structures of two of the most well-known Panax species was performed to better understand the physiological processes that limit photosynthesis.

Methods

Leaves of ginseng plants grown in soil and hydroponic culture were sectioned for analysis. Leaf structures of both Panax species were observed using a light microscope, scanning electron microscope, and transmission electron microscope.

Results

The mesostructure of both P. ginseng and P. quinquefolius frequently had one layer of noncylindrical palisade cells and three or four layers of spongy parenchymal cells. P. quinquefolius contained a similar number of stomata in the abaxial leaf surface but more tightly appressed enlarged grana stacks than P. ginseng contained. The adaxial surface of the epidermis in P. quinquefolius showed cuticle ridges with a pattern similar to that of P. ginseng.

Conclusion

The anatomical leaf structure of both P. ginseng and P. quinquefolius shows that they are typical shade-loving sciophytes. Slight differences in chloroplast structure suggests that the two different species can be authenticated using transmission electron microscopy images, and light-resistant cultivar breeding can be performed via controlling photosynthesis efficiency.

A Gossypium BAC clone contains key repeat components distinguishing sub-genome of allotetraploidy cottons

Background

Dissecting genome organization is indispensable for further functional and applied studies. As genome sequences data shown, cotton genomes contain more than 60 % repetitive sequences, so study on repetitive sequences composition, structure, and distribution is the key step to dissect cotton genome.

Results

In this study, a bacterial artificial chromosome (BAC) clone enriched in repetitive sequences, was discovered initiatively by fluorescence in situ hybridization (FISH). FISHing with allotetraploidy cotton as target DNA, dispersed signals on most regions of all A sub-genome chromosomes, and only middle regions of all D sub-genome chromosomes were detected. Further FISHing with other cotton species bearing A or D genome as target DNA, specific signals were viewed. After BAC sequencing and bioinformational analysis, 129 repeat elements, size about 57,172 bp were found, accounting for more than 62 % of the BAC sequence (91,238 bp). Among them, a type of long terminal repeat-retrotransposon (LTR-RT), LTR/Gypsy was the key element causing the specific FISH results. Using the fragments of BAC matching with the identified Gypsy-like LTR as probes, the BAC-57I23-like FISH signals were reappeared. Running BLASTN, the fragments had good match with all chromosomes of G. arboreum (A₂) genome and A sub-genome of G. hirsutum (AD₁), and had relatively inferior match with all chromosomes of D sub-genome of AD₁, but had little match with the chromosomes of G. raimondii (D₅) genome, which was consistent with the FISH results.

Conclusion

A repeats-enriched cytogenetic marker to identify A and D sub-genomes of Gossypium was discovered by FISH. Combined sequences analysis with FISH verification, the assembly quality of repetitive sequences in the allotetraploidy cotton draft genome was assessed, and better chromosome belonging was verified. We also found the genomic distribution of the identified Gypsy-LTR-RT was similar to the distribution of heterochromatin. The expansion of this type of Gypsy-LTR-RT in heterochromatic regions may be one of the major reasons for the size gap between A and D genome. The findings showed here will help to understand the composition, structure, and evolution of cotton genome, and contribute to the further perfection of the draft genomes of cotton.

The impacts of polyploidy, geographic and ecological isolations on the diversification of Panax (Araliaceae)

BACKGROUND

Panax L. is a medicinally important genus within family Araliaceae, where almost all species are of cultural significance for traditional Chinese medicine. Previous studies suggested two independent origins of the East Asia and North America disjunct distribution of this genus and multiple rounds of whole genome duplications (WGDs) might have occurred during the evolutionary process.

RESULTS

We employed multiple chloroplast and nuclear markers to investigate the evolution and diversification of Panax. Our phylogenetic analyses confirmed previous observations of the independent origins of disjunct distribution and both ancient and recent WGDs have occurred within Panax. The estimations of divergence time implied that the ancient WGD might have occurred before the establishment of Panax. Thereafter, at least two independent recent WGD events have occurred within Panax, one of which has led to the formation of three geographically isolated tetraploid species P. ginseng, P. japonicus and P. quinquefolius. Population genetic analyses showed that the diploid species P. notoginseng harbored significantly lower nucleotide diversity than those of the two tetraploid species P. ginseng and P. quinquefolius and the three species showed distinct nucleotide variation patterns at exon regions.

CONCLUSION

Our findings based on the phylogenetic and population genetic analyses, coupled with the species distribution patterns of Panax, suggested that the two rounds of WGD along with the geographic and ecological isolations might have together contributed to the evolution and diversification of this genus.

Genetic and Epigenetic Diversities Shed Light on Domestication of Cultivated Ginseng (Panax ginseng)

Chinese ginseng (Panax ginseng) is a medically important herb within Panax and has crucial cultural values in East Asia. As the symbol of traditional Chinese medicine, Chinese ginseng has been used as a herbal remedy to restore stamina and capacity in East Asia for thousands of years. To address the evolutionary origin and domestication history of cultivated ginseng, we employed multiple molecular approaches to investigate the genetic structures of cultivated and wild ginseng across their distribution ranges in northeastern Asia. Phylogenetic and population genetic analyses revealed that the four cultivated ginseng landraces, COMMON, BIANTIAO, SHIZHU, and GAOLI (also known as Korean ginseng), were not domesticated independently and Fusong Town is likely one of the primary domestication centers. In addition, our results from population genetic and epigenetic analyses demonstrated that cultivated ginseng maintained high levels of genetic and epigenetic diversity, but showed distinct cytosine methylation patterns compared with wild ginseng. The patterns of genetic and epigenetic variation revealed by this study have shed light on the domestication history of cultivated ginseng, which may serve as a framework for future genetic improvements.

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Comprehensive analysis of Panax ginseng root transcriptomes

BACKGROUND

Korean ginseng (Panax ginseng C.A. Meyer) is a highly effective medicinal plant containing ginsenosides with various pharmacological activities, whose roots are produced commercially for crude drugs.

RESULTS

Here, we used the Illumina platform to generate over 232 million RNA sequencing reads from four root samples, including whole roots from one-year-old plants and three types of root tissue from six-year-old plants (i.e., main root bodies, rhizomes, and lateral roots). Through de novo assembly and reference-assisted selection, we obtained a non-redundant unigene set consisting of 55,949 transcripts with an average length of 1,250 bp. Among transcripts in the unigene set, 94 % were functionally annotated via similarity searches against protein databases. Approximately 28.6 % of the transcripts represent novel gene sequences that have not previously been reported for P. ginseng. Digital expression profiling revealed 364 genes showing differential expression patterns among the four root samples. Additionally, 32 were uniquely expressed in one-year-old roots, while seven were uniquely expressed in six-year-old root tissues. We identified 38 transcripts encoding enzymes involved in ginsenoside biosynthesis pathways and 189 encoding UDP-glycosyltransferases.

CONCLUSION

Our analysis provides new insights into the role of the root transcriptome in development and secondary metabolite biosynthesis in P. ginseng.

Comprehensive Survey of Genetic Diversity in Chloroplast Genomes and 45S nrDNAs within Panax ginseng Species

We report complete sequences of chloroplast (cp) genome and 45S nuclear ribosomal DNA (45S nrDNA) for 11 Panax ginseng cultivars. We have obtained complete sequences of cp and 45S nrDNA, the representative barcoding target sequences for cytoplasm and nuclear genome, respectively, based on low coverage NGS sequence of each cultivar. The cp genomes sizes ranged from 156,241 to 156,425 bp and the major size variation was derived from differences in copy number of tandem repeats in the ycf1 gene and in the intergenic regions of rps16-trnUUG and rpl32-trnUAG. The complete 45S nrDNA unit sequences were 11,091 bp, representing a consensus single transcriptional unit with an intergenic spacer region. Comparative analysis of these sequences as well as those previously reported for three Chinese accessions identified very rare but unique polymorphism in the cp genome within P. ginseng cultivars. There were 12 intra-species polymorphisms (six SNPs and six InDels) among 14 cultivars. We also identified five SNPs from 45S nrDNA of 11 Korean ginseng cultivars. From the 17 unique informative polymorphic sites, we developed six reliable markers for analysis of ginseng diversity and cultivar authentication.