Diversity and evolution of major Panax species revealed by scanning the entire chloroplast intergenic spacer sequences

Ginseng is an important medicinal plant, but almost no genomic information is known for it. For the primary step to understand the Panax genome, we inspected the chloroplast genome sequence diversity and used that to infer the evolution of Panax species using them. We inspected a total of 101 intergenic spacers (IGS) covering 44,563 bp (96.8 % of the total IGS) from four Panax species. Diversity was inspected by three steps: gel electrophoresis, high resolution melting (HRM) analysis, and further confirmation by sequencing. Even though low levels of InDel polymorphism were detected by gel electrophoresis, high levels of reproducible polymorphisms were identified by HRM analysis. Sequencing of the HRM-polymorphic spaces and pair-wise multiple sequence alignments revealed up to 247 sequence variations between species. We have identified 62 IGS showing polymorphism between species. Among them, the trnE-trnT, trnT-psbD, ndhF-rpl32, and rpl14-rpl16 spaces are more informative for studying the diversity of Panax relatives. Phylogenetic analysis and molecular evolution studies revealed that P. notoginseng is most diverged from the other Panax species, with a nucleotide substitution rate of 0.0039. P. quinquefolius shows a close relationship with P. ginseng and P. japonicus, with a nucleotide substitution rate of 0.0009. Meanwhile, no sequence variation was detected between P. ginseng and P. japonicus. Calculation of molecular clocks revealed that P. notoginseng diverged more than 1.30 million years ago from the other Panax species and then, P. quinquefolius diverged from P. ginseng by migration in the American continent more than 0.29 million years ago.