Evolutionary unraveling: new insights into the Persicaria amphibia complex.

The Persicaria amphibia complex exhibits significant morphological variation depending on its habitat, existing in either aquatic or terrestrial forms. Traditionally, four distinct elements have been recognized based on morphological features along with their distinct geographical distributions. Recent studies suggest that the Asian element may be genetically distinct from the European and American elements. However, a comprehensive study on the genetic differentiation among all four elements remains lacking. This study aimed to leverage whole plastid genome sequences and ITS2 haplotypes to comprehensively assess the genomic diversity within the P. amphibia complex. Notably, we included multiple individuals from New York State to resolve the ongoing debate regarding the taxonomic status of two American elements - whether they represent a single species or distinct entities. Our analysis revealed a well-supported monophyletic clade encompassing all four elements, endorsing their own section, Amphibia. Notably, the terrestrial form of the American element is sister to all other elements, suggesting it deserves its own species status. This reinstates its historical name, P. coccinea, separating it from the broader P. amphibia. Furthermore, distinct compositions of the ITS2 haplotypes differentiated the four elements, although the European element should be further investigated with more sampling. The most intriguing discovery is the identification of putative hybrids between the two American elements. In one population out of four putative hybrid populations, all three entities - the two parent species and their hybrid offspring - thrive together, showcasing a fascinating microcosm of ongoing evolutionary processes. Unraveling the intricate genetic tapestry within each American species and their hybrid populations remains a compelling next step. By delving deeper into their genetic makeup, we can gain a richer understanding of their evolutionary trajectories and the intricacies of their interactions. Finally, it is estimated that the two species of sect. Amphibia diverged approximately 4.02 million years ago during the Pliocene epoch, when there was a significant global cooling and drying trend.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data.

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Sophora genomes provide insight into the evolution of alkaloid metabolites along with small-scale gene duplication.

The genus Sophora (Fabaceae) includes medicinal plants that have been used in East Asian countries since antiquity. Sophora flavescens is a perennial herb indigenous to China, India, Japan, Korea, and Russia. Its dried roots have antioxidant, anti-inflammatory, antibacterial, apoptosis-modulating, and antitumor efficacy. The congeneric S. koreensis is endemic to Korea and its genome is less than half the size of that of S. flavescens. Nevertheless, this discrepancy can be used to assemble and validate the S. flavescens genome. A comparative genomic study of the two genomes can disclose the recent evolutionary divergence of the polymorphic phenotypic profiles of these species. Here, we used the PacBio sequencing platform to sequence and assemble the S. koreensis and S. flavescens genomes. We inferred that it was mainly small-scale duplication that occurred in S. flavescens. A KEGG analysis revealed pathways that might regulate the pharmacologically important secondary metabolites in S. flavescens and S. koreensis. The genome assemblies of Sophora spp. could be used in comparative genomics and data mining for various plant natural products.

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere.

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.

Structure and antiviral activity of a pectic polysaccharide from the root of Sanguisorba officinalis against enterovirus 71 in vitro/vivo.

The increasing prevalence and pandemic risk of viral diseases warrant the development of safe and effective treatments. In this study, we aimed to elucidate the structure and anti-enterovirus 71 (EV71) effects of polysaccharides isolated from the roots of Sanguisorba officinalis (SO), traditionally used for infectious diseases. The purified polysaccharide (S-a3) was a homogenous macromolecule (260.4 kDa) with a concave and porous surface. Linkage and NMR analyses confirmed that S-a3 is a polysaccharide interlinked with homogalacturonan, rhamnogalacturonan-I, 1,4-α-glucan, and arabinogalactan. S-a3 significantly inhibited cell death and viral gene expression in EV71-infected Vero cells, and alleviated EV71-induced body weight loss, death, and paralysis in the hSCARB2-transgenic mouse model. The effective dose of S-a3 was non-toxic to cells and mice. The antiviral mechanism of S-a3 was associated with the disruption of EV71 attachment to host cells. Our findings demonstrate that polysaccharides from SO can be a safe and effective treatment for EV71 infection.

A chromosome-scale genome assembly and annotation of the spring orchid (Cymbidium goeringii).

Cymbidium goeringii, commonly known as the spring orchid, has long been favoured for horticultural purposes in Asian countries. It is a popular orchid with much demand for improvement and development for its valuable varieties. Until now, its reference genome has not been published despite its popularity and conservation efforts. Here, we report the de novo assembly of the C. goeringii genome, which is the largest among the orchids published to date, using a strategy that combines short- and long-read sequencing and chromosome conformation capture (Hi-C) information. The total length of all scaffolds is 3.99 Gb, with an N50 scaffold size of 178.2 Mb. A total of 29,556 protein-coding genes were annotated and 3.55 Gb (88.87% of genome) repetitive sequences were identified. We constructed pseudomolecular chromosomes using Hi-C, incorporating 89.4% of the scaffolds in 20 chromosomes. We identified 220 expanded and 106 contracted genes families in C. goeringii after divergence from its close relative. We also identified new gene families, resistance gene analogues and changes within the MADS-box genes, which control a diverse set of developmental processes during orchid evolution. Our high quality chromosomal-level assembly of C. goeringii can provide a platform for elucidating the genomic evolution of orchids, mining functional genes for agronomic traits and for developing molecular markers for accelerated breeding as well as accelerating conservation efforts.

The complete chloroplast genome sequence of Viola verecunda (Violaceae).

Viola verecunda is a perennial plant native to the mountainous areas of Northeast Asia. Here the complete chloroplast genome of V. verecunda and its phylogenetic relationships to other species within the genus Viola are reported. The complete chloroplast genome of V. verecunda is 157,843 bp in length and circular in structure with four regions: a large single-copy region (86,345 bp), a small single-copy region (17,292 bp), and a pair of inverted repeat regions (27,103 bp each). The chloroplast genome contains 111 unique genes comprising 77 protein-coding, 30 unique tRNA, and 4 unique rRNA genes. Based on the protein-coding gene sequences from eight Viola chloroplast genomes, with Balanops balansae designated as the outgroup, maximum likelihood tree analysis indicates that V. verecunda is more closely related to V. raddeana than to other Viola species. The complete chloroplast genome of V. verecunda contributes to a better understanding of the phylogenetic relationships among Viola species.

The complete plastid genome sequence of Quercus acuta (Fagaceae), an evergreen broad-leaved oak endemic to East Asia.

We are reporting the complete plastid genome (plastome) of Quercus acuta, an evergreen broad-leaved oak endemic to East Asia. This species is important for maintaining the warm-temperate evergreen forest biome in East Asia. The Q. acuta plastome is 160,522 base pairs (bp) long, with two inverted repeat (IR) regions (25,839 bp each) that separate a large single copy (LSC) region (90,199 bp) and a small single copy (SSC) region (18,645 bp). The phylogenetic tree shows that Quercus acuta is closely related to Quercus sichourensis with strong bootstrap support.

Three complete plastome sequences from the families of Lamiaceae, Mazaceae, and Phrymaceae (Lamiales).

In this study, we report the first complete plastome sequence of Vitex rotundifolia (Lamiaceae) (MT937186). In addition, the plastome sequences of Phryma leptostachya subsp. asiatica (Phrymaceae) (153,324 bp; MT948145) and Mazus pumilus (Mazaceae) (152,847 bp; MT937187) are also included. The gene orders and structures of the three plastomes are collinear with those of the typical plastome of angiosperm. The plastome size of V. rotundifolia is 154,370 bp in length and consists of a large single-copy region of 85,079 bp and a small single-copy region of 17,917 bp, which are separated by a pair of 25,687 bp-long inverted repeat regions. In addition, the plastome sizes of P. leptostachya subsp. asiatica and M. pumilus are 153,324 bp and 152,847 bp, respectively. The three plastomes contain 113 genes, including 79 protein-coding, 30 tRNA, and four rRNA genes. Sixteen genes contain one intron and two genes have two introns. A total of 41 simple sequence repeat loci was identified in the V. rotundifolia plastome. Phylogenetic analysis shows that Viticoideae is a sister group of the last of Lamiaceae except Nepetoideae. The Mazaceae are a sister group of Lamiaceae, while Phrymaceae form a sister group to the Paulowniaceae-Orobanchaceae clade.

The complete chloroplast genome of a fern genus Thelypteris interrupta.

In this study, we report the complete chloroplast (cp) genome of Thelypteris interrupta, a fern member, and comparative analysis with its related family members. The cp genome was 155,983 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (25,614 bp) separated by a large (82,769 bp) and small (21,986 bp) single-copy (SC) region. The genome encodes a total of 88 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Additionally, we identified 86 RNA editing sites in 52 genes; most of the substitution was U to C (52 sites), while C to U conversion occurred in 34 positions. The phylogenetic analysis strongly supported the relationship of T. interrupta with Ampelopteris prolifera and Christella appendiculata of Thelypteridoideae family.

Identification of plastid genomic regions inferring species identity from de novo plastid genome assembly of 14 Korean-native Iris species (Iridaceae).

Iris is one of the largest genera in the family Iridaceae, comprising hundreds of species, including numerous economically important horticultural plants used in landscape gardening and herbal medicine. Improved taxonomic classification of Iris species, particularly the endangered Korean-native Iris, is needed for correct species delineation. To this end, identification of diverse genetic markers from Iris genomes would facilitate molecular identification and resolve ambiguous classifications from molecular analyses; however, only two Iris plastid genomes, from Iris gatesii and Iris sanguinea, have been sequenced. Here, we used high-throughput next-generation sequencing, combined with Sanger sequencing, to construct the plastid genomes of 14 Korean-native Iris species with one outgroup and predict their gene content. Using these data, combined with previously published plastid genomes from Iris and one outgroup (Sisyrinchium angustifolium), we constructed a Bayesian phylogenetic tree showing clear speciation among the samples. We further identified sub-genomic regions that have undergone neutral evolution and accurately recapitulate Bayesian-inferred speciation. These contain key markers that could be used to identify and classify Iris samples into taxonomic clades. Our results confirm previously reported speciation patterns and resolve questionable relationships within the Iris genus. These data also provide a valuable resource for studying genetic diversity and refining phylogenetic relationships between Iris species.

A Disjunctive Marginal Edge of Evergreen Broad-Leaved Oak (Quercus gilva) in East Asia: The High Genetic Distinctiveness and Unusual Diversity of Jeju Island Populations and Insight into a Massive, Independent Postglacial Colonization.

Jeju Island is located at a marginal edge of the distributional range of East Asian evergreen broad-leaved forests. The low genetic diversity of such edge populations is predicted to have resulted from genetic drift and reduced gene flow when compared to core populations. To test this hypothesis, we examined the levels of genetic diversity of marginal-edge populations of Quercus gilva, restricted to a few habitats on Jeju Island, and compared them with the southern Kyushu populations. We also evaluated their evolutionary potential and conservation value. The genetic diversity and structure were analyzed using 40 polymorphic microsatellite markers developed in this study. Ecological Niche Modeling (ENM) has been employed to develop our insights, which can be inferred from historical distribution changes. Contrary to our expectations, we detected a similar level of genetic diversity in the Jeju populations, comparable to that of the southern Kyushu populations, which have been regarded as long-term glacial refugia with a high genetic variability of East Asian evergreen trees. We found no signatures of recent bottlenecks in the Jeju populations. The results of STRUCTURE, neighbor-joining phylogeny, and Principal Coordinate Analysis (PCoA) with a significant barrier clearly demonstrated that the Jeju and Kyushu regions are genetically distinct. However, ENM showed that the probability value for the distribution of the trees on Jeju Island during the Last Glacial Maximum (LGM) converge was zero. In consideration of these results, we hypothesize that independent massive postglacial colonization from a separate large genetic source, other than Kyushu, could have led to the current genetic diversity of Jeju Island. Therefore, we suggest that the Jeju populations deserve to be separately managed and designated as a level of management unit (MU). These findings improve our understanding of the paleovegetation of East Asian evergreen forests, and the microevolution of oaks.

Corrigendum: Plastome Evolution and Phylogeny of Orchidaceae, With 24 New Sequences.

[This corrects the article DOI: 10.3389/fpls.2020.00022.].

Plastome Evolution and Phylogeny of Orchidaceae, With 24 New Sequences.

In order to understand the evolution of the orchid plastome, we annotated and compared 124 complete plastomes of Orchidaceae representing all the major lineages in their structures, gene contents, gene rearrangements, and IR contractions/expansions. Forty-two of these plastomes were generated from the corresponding author's laboratory, and 24 plastomes-including nine genera (Amitostigma, Bulbophyllum, Dactylorhiza, Dipodium, Galearis, Gymnadenia, Hetaeria, Oreorchis, and Sedirea)-are new in this study. All orchid plastomes, except Aphyllorchis montana, Epipogium aphyllum, and Gastrodia elata, have a quadripartite structure consisting of a large single copy (LSC), two inverted repeats (IRs), and a small single copy (SSC) region. The IR region was completely lost in the A. montana and G. elata plastomes. The SSC is lost in the E. aphyllum plastome. The smallest plastome size was 19,047 bp, in E. roseum, and the largest plastome size was 178,131 bp, in Cypripedium formosanum. The small plastome sizes are primarily the result of gene losses associated with mycoheterotrophic habitats, while the large plastome sizes are due to the expansion of noncoding regions. The minimal number of common genes among orchid plastomes to maintain minimal plastome activity was 15, including the three subunits of rpl (14, 16, and 36), seven subunits of rps (2, 3, 4, 7, 8, 11, and 14), three subunits of rrn (5, 16, and 23), trnC-GCA, and clpP genes. Three stages of gene loss were observed among the orchid plastomes. The first was ndh gene loss, which is widespread in Apostasioideae, Vanilloideae, Cypripedioideae, and Epidendroideae, but rare in the Orchidoideae. The second stage was the loss of photosynthetic genes (atp, pet, psa, and psb) and rpo gene subunits, which are restricted to Aphyllorchis, Hetaeria, Hexalectris, and some species of Corallorhiza and Neottia. The third stage was gene loss related to prokaryotic gene expression (rpl, rps, trn, and others), which was observed in Epipogium, Gastrodia, Lecanorchis, and Rhizanthella. In addition, an intermediate stage between the second and third stage was observed in Cyrtosia (Vanilloideae). The majority of intron losses are associated with the loss of their corresponding genes. In some orchid taxa, however, introns have been lost in rpl16, rps16, and clpP(2) without their corresponding gene being lost. A total of 104 gene rearrangements were counted when comparing 116 orchid plastomes. Among them, many were concentrated near the IRa/b-SSC junction area. The plastome phylogeny of 124 orchid species confirmed the relationship of {Apostasioideae [Vanilloideae (Cypripedioideae (Orchidoideae, Epidendroideae))]} at the subfamily level and the phylogenetic relationships of 17 tribes were also established. Molecular clock analysis based on the whole plastome sequences suggested that Orchidaceae diverged from its sister family 99.2 mya, and the estimated divergence times of five subfamilies are as follows: Apostasioideae (79.91 mya), Vanilloideae (69.84 mya), Cypripedioideae (64.97 mya), Orchidoideae (59.16 mya), and Epidendroideae (59.16 mya). We also released the first nuclear ribosomal (nr) DNA unit (18S-ITS1-5.8S-ITS2-28S-NTS-ETS) sequences for the 42 species of Orchidaceae. Finally, the phylogenetic tree based on the nrDNA unit sequences is compared to the tree based on the 42 identical plastome sequences, and the differences between the two datasets are discussed in this paper.

The complete chloroplast genome of a Korean endemic species Sophora koreensis, Nakai.

In this study, we report the complete chloroplast (cp) genome of Sophora koreensis and its relation with other species within the Fabaceae family. The cp genome was 154,870 bp long, with a typical quadripartite structure including a pair of inverted repeat regions (25,866 bp) separated by a large (85,037 bp) and small (18,101 bp) single-copy (SC) region. The genome encodes a total of 84 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Phylogenetic analysis suggested that S. koreensis is closely related to genus Sophora alopecuroides var. alopecuroides within Fabaceae.

Neuroprotective effects of Dendropanax morbifera leaves on glutamate-induced oxidative cell death in HT22 mouse hippocampal neuronal cells.

ETHNOPHARMACOLOGICAL RELEVANCE: Dendropanax morbifera (DM) has long been used as a traditional herbal medicine for migraines. Glutamate toxicity and oxidative stress have emerged as the possible triggers implicated in migraine pathogenesis. AIM OF THE STUDY: We aimed to examine the neuroprotective effects of DM leaves (DML) on glutamate-induced oxidative cell death in HT22 mouse hippocampal neuronal cells. MATERIALS AND METHODS: Molecular authentication of DML was assessed using DNA barcoding analysis. Four different solvent extracts of DML were prepared and subjected to antioxidant activity and phytochemical assays. Neuroprotective effects of DML extracts were evaluated using relevant biochemical and imaging assays that measure cell viability/death, ROS generation, Ca2+ levels, mitochondrial dysfunction, and AIF nuclear translocation. RESULTS: The sequences of matK, rbcL, atpF-H, and psbK-I in DML were identical with those in voucher specimens, confirming that DML was indeed D. morbifera. The ethyl acetate extract of DML (DMLE) showed the highest flavonoid and phenolic content, and prominent DPPH/superoxide radical scavenging and reducing power activities. In the HT22 cell model, glutamate was shown to be the causative agent for apoptotic cell death via elevation of intracellular ROS and Ca2+ levels, induction of mitochondrial depolarization and membrane permeabilization, and translocation of AIF to the nucleus. Of note, N-acetyl-L-cysteine and necrostatin-1, but not z-VAD-fmk, completely prevented glutamate-induced cell death, implying that oxidative stress and AIF translocation were pivotal in glutamate cytotoxicity. DMLE significantly recovered glutamate-induced apoptotic cell death in a concentration-dependent manner. It completely inhibited intracellular/mitochondrial ROS generation, the elevation of Ca2+ levels, and mitochondrial dysfunction induced by glutamate during early exposure within 8 h. It significantly reversed subsequent AIF nuclear translocation after 12 h of treatment. Antioxidant activities of DMLE may be the protective mechanism that regulates homeostatic balance of ROS and Ca2+ as well as maintains mitochondrial function. CONCLUSIONS: DMLE shows significant neuroprotective effects against glutamate-induced oxidative neuronal cell death. Therefore, DM could be a potential therapeutic candidate for neurological disorders propagated by glutamate toxicity.

Plastome evolution and phylogeny of subtribe Aeridinae (Vandeae, Orchidaceae).

Subtribe Aeridinae (Vandeae, Epidendroideae, Orchidaceae) consists of 83 genera and 2,345 species. The present study completely decoded the plastomes and nuclear ribosomal (nr) RNA gene clusters of seven species of Aeridinae belonging to Gastrochilus, Neofinetia, Pelatantheria, and Thrixspermum and compared them with existing data to investigate their genome evolution and phylogeny. Although no large structural variations were observed among the Aeridinae plastomes, 14 small inversions (SI) were found in Orchidaceae for the first time. Therefore, the evolutionary trends and usefulness of SI as molecular identification markers were evaluated. Since all 11 ndh genes in the Aeridinae plastome were lost or pseudogenized, the evolutionary trends of ndh genes are discussed at the tribe and family levels. In the maximum likelihood tree reconstructed from 83 plastome genes, the five Orchidaceae subfamilies were shown to have diverged in the following order: Apostasioideae, Vanilloideae, Cypripedioideae, Orchioideae, Epidendroideaeae. Divergence times for major lineages were found to be more recent, 5-10 Mya, than previous studies, which only used two or three genes. Vandeae, which includes Aeridinae, formed a sister group with Cymbidieae and Epidendreae. The Vandeae, Cymbidieae, and Epidendreae lineages were inferred to have diverged at 25.31 Mya; thus, numerous speciation events within Aeridineae occurred since then. Furthermore, the present study reconstructed a phylogenetic tree from 422 nrITS sequences belonging to Aerdinae and allied taxa and uses it to discuss the phylogenetic positions and species identities of five endangered species.

Dynamic Chloroplast Genome Rearrangement and DNA Barcoding for Three Apiaceae Species Known as the Medicinal Herb "Bang-Poong".

Three Apiaceae species Ledebouriella seseloides, Peucedanum japonicum, and Glehnia littoralis are used as Asian herbal medicines, with the confusingly similar common name "Bang-poong". We characterized the complete chloroplast (cp) genomes and 45S nuclear ribosomal DNA (45S nrDNA) sequences of two accessions for each species. The complete cp genomes of G. littoralis, L. seseloides, and P. japonicum were 147,467, 147,830, and 164,633 bp, respectively. Compared to the other species, the P. japonicum cp genome had a huge inverted repeat expansion and a segmental inversion. The 45S nrDNA cistron sequences of the three species were almost identical in size and structure. Despite the structural variation in the P. japonicum cp genome, phylogenetic analysis revealed that G. littoralis diverged 5-6 million years ago (Mya), while P. japonicum diverged from L. seseloides only 2-3 Mya. Abundant copy number variations including tandem repeats, insertion/deletions, and single nucleotide polymorphisms, were found at the interspecies level. Intraspecies-level polymorphism was also found for L. seseloides and G. littoralis. We developed nine PCR barcode markers to authenticate all three species. This study characterizes the genomic differences between L. seseloides, P. japonicum, and G. littoralis; provides a method of species identification; and sheds light on the evolutionary history of these three species.

The complete chloroplast genome sequences of four Viola species (Violaceae) and comparative analyses with its congeneric species.

We report the complete chloroplast genomes of four Viola species (V. mirabilis, V. phalacrocarpa, V. raddeana, and V. websteri) and the results of a comparative analysis between these species and the published plastid genome of the congeneric species V. seoulensis. The total genome length of the five Viola species, including the four species analyzed in this study and the species analyzed in the previous study, ranged from 156,507 (V. seoulensis) to 158,162 bp (V. mirabilis). The overall GC contents of the genomes were almost identical (36.2-36.3%). The five Viola plastomes each contained 111 unique genes comprising 77 protein-coding genes, 30 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes. Among the annotated genes, 16 contained one or two introns. Based on the results of a chloroplast genome structure comparison using MAUVE, all five Viola plastomes were almost identical. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved among the Viola species. A total of 259 exon, intron, and intergenic spacer (IGS) fragments were compared to verify the divergence hotspot regions. The nucleotide diversity (Pi) values ranged from 0 to 0.7544. The IR region was relatively more conserved than the LSC and SSC regions. The Pi values in ten noncoding regions were relatively high (>0.03). Among these regions, all but rps19-trnH, petG-trnW, rpl16-rps3, and rpl2-rpl23 represent useful molecular markers for phylogenetic studies and will be helpful to resolve the phylogenetic relationships of Viola. The phylogenetic tree, which used 76 protein-coding genes from 21 Malpighiales species and one outgroup species (Averrhoa carambola), revealed that Malpighiales is divided into five clades at the family level: Erythroxylaceae, Chrysobalanaceae, Euphorbiaceae, Salicaceae, and Violaceae. Additionally, Violaceae was monophyletic, with a bootstrap value of 100% and was divided into two subclades.

Genetic diversity and population structure of Lychnis wilfordii (Caryophyllaceae) with newly developed 17 microsatellite markers.

Lychnis wilfordii (Regel) Maxim. is a perennial plant designated as an endangered species by the Korean government because of rapid reduction in its population size. Thus, a population genetic study of this species is needed to establish the strategy for management and conservation based on scientific evidences. The goals of this study were to develop useful microsatellite markers for L. wilfordii and to understand current genetic status of L. wilfordii in Korean peninsula. Seventeen microsatellite markers were identified using next-generation sequencing and bioinformatic analysis and then analyzed genetic diversity in one hundred forty-five individuals from Korea (KI1, KI2, and KP), China (CX, CF) and Russia (RP). Analysis of molecular variance (AMOVA), principal coordinates analysis (PCoA) and STRUCTURE results consistently showed discontinuity among L. wilfordii populations. AMOVA showed that the percentage of variation among populations was 53%, which was higher than the variation within populations (19%). PCoA showed that the populations were divided into three genetic clusters, (1) Chinese (CX, CF), (2) Russian (RP) populations and Korean populations (KI1, KI2) excluding KP, and (3) the KP population. In particular, KP, the most southern population on the Korean peninsula, showed significantly lower observed and expected heterozygosity, number of effective alleles, and Shannon index (I) than those of KI1 and KI2. L. wilfordii showed high differentiation between populations with low genetic diversity within populations. Among Korean populations, KP is likely to be affected by genetic drift due to small population size, low genetic diversity and limited gene flow.