Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub.

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

The complete chloroplast genome sequences of six Hylotelephium species: Comparative genomic analysis and phylogenetic relationships.

To evaluate the phylogenetic relationships between Hylotelephium and Orostachys, and to provide important information for further studies, we analyzed the complete chloroplast genomes of six Hylotelephium species and compared the sequences to those of published chloroplast genomes of congeneric species and species of the closely related genus, Orostachys. The total chloroplast genome length of nineteen species, including the six Hylotelephium species analyzed in this study and the thirteen Hylotelephium and Orostachys species analyzed in previous studies, ranged from 150,369 bp (O. minuta) to 151,739 bp (H. spectabile). Their overall GC contents were almost identical (37.7-37.8%). The chloroplast genomes of the nineteen species contained 113 unique genes comprising 79 protein-coding genes (PCGs), 30 transfer RNA genes (tRNAs), and four ribosomal RNA genes (rRNAs). Among the annotated genes, fourteen genes contained one intron, and two genes contained two introns. The chloroplast genomes of the nineteen Hylotelephium and Orostachys species had identical structures. Additionally, the large single copy (LSC), inverted repeat (IR), and small single copy (SSC) junction regions were conserved in the Hylotelephium and Orostachys species. The nucleotide diversity between the Hylotelephium chloroplast genomes was extremely low in all regions, and only one region showed a high Pi value (>0.03). In all nineteen chloroplast genomes, six regions had a high Pi value (>0.03). The phylogenetic analysis showed that the genus delimitation could not be clearly observed even in this study because Hylotelephium formed a paraphyly with subsect. Orostachys of the genus Orostachys. Additionally, the data supported the taxonomic position of Sedum taqeutii, which was treated as a synonym for H. viridescens in previous studies, as an independent taxon.

Isolation and Characterization of the Genes Involved in the Berberine Synthesis Pathway in Asian Blue Cohosh, Caulophyllum robustum.

Caulophyllum robustum, commonly named Asian blue cohosh, is a perennial herb in the family Berberidaceae. It has traditionally been used for folk medicine in China. We isolated berberine from the leaves, stem, roots, and fruits of C. robustum, and this is the first report on berberine in this species. Transcriptome analysis was conducted for the characterization of berberine biosynthesis genes in C. robustum, in which, all the genes for berberine biosynthesis were identified. From 40,094 transcripts, using gene ontology (GO) analysis, 26,750 transcripts were assigned their functions in the categories of biological process, molecular function, and cellular component. In the analysis of genes expressed in different tissues, the numbers of genes in the categories of intrinsic component of membrane and transferase activity were up-regulated in leaves versus stem. The berberine synthesis genes in C. robustum were characterized by phylogenetic analysis with corresponding genes from other berberine-producing species. The co-existence of genes from different plant families in the deepest branch subclade implies that the differentiation of berberine synthesis genes occurred early in the evolution of berberine-producing plants. Furthermore, the copy number increment of the berberine synthesis genes was detected at the species level.

Genetic and demographic signatures accompanying the evolution of the selfing syndrome in Daphne kiusiana, an evergreen shrub.

BACKGROUND AND AIMS: The evolution of mating systems from outcrossing to self-fertilization is a common transition in flowering plants. This shift is often associated with the 'selfing syndrome', which is characterized by less visible flowers with functional changes to control outcrossing. In most cases, the evolutionary history and demographic dynamics underlying the evolution of the selfing syndrome remain poorly understood. METHODS: Here, we characterize differences in the demographic genetic consequences and associated floral-specific traits between two distinct geographical groups of a wild shrub, Daphne kiusiana, endemic to East Asia; plants in the eastern region (southeastern Korea and Kyushu, Japan) exhibit smaller and fewer flowers compared to those of plants in the western region (southwestern Korea). Genetic analyses were conducted using nuclear microsatellites and chloroplast DNA (multiplexed phylogenetic marker sequencing) datasets. KEY RESULTS: A high selfing rate with significantly increased homozygosity characterized the eastern lineage, associated with lower levels of visibility and herkogamy in the floral traits. The two lineages harboured independent phylogeographical histories. In contrast to the western lineage, the eastern lineage showed a gradual reduction in the effective population size with no signs of a severe bottleneck despite its extreme range contraction during the last glacial period. CONCLUSIONS: Our results suggest that the selfing-associated morphological changes in D. kiusiana are of relatively old origin (at least 100 000 years ago) and were driven by directional selection for efficient self-pollination. We provide evidence that the evolution of the selfing syndrome in D. kiusiana is not strongly associated with a severe population bottleneck.

The complete plastid genome sequence of Chloranthus fortunei (A. Gray) Solms-Laub. in Chloranthaceae.

Chloranthus fortunei (A. Gray) Solms-Laub. is a perennial herb in a basal angiosperm family Chloranthaceae. Here, we reported the complete plastid genome of C. fortunei using Illumina short-read data. The total genome size was 157,063 bp in length, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The gene content and order were consistent with previously reported Chloranthus plastid genomes. The overall GC content of the C. fortunei plastid genome was 39.0%. In the phylogenetic result, genus Chloranthus was monophyletic and divided into two subclades: C. japonicus+C. angustifolius+C. fortunei, and C. henryi+C. spicatus+C. erectus. Our phylogenetic result was consistent with previous phylogenetic studies, and was supported by a previously proposed infrageneric classification of the genus Chloranthus.

The complete chloroplast genome sequence of Orostachys minuta (Crassulaceae).

The chloroplast (cp) genome sequence is determined and analyzed for Orostachys minuta for the first time. The cp genome was 150,369 bp in length, containing a large single-copy (LSC) of 82,795 bp and a small single-copy (SSC) of 16,854 bp, which were separated by a pair of 25,360 bp inverted repeats (IRs). The overall G + C content of the O. minuta cp genome amounted to 37.7%. In total, 113 unique genes were annotated, consisting of 79 protein-coding genes (PCGs), 30 transfer RNAs (tRNAs), and four ribosomal RNAs (rRNAs). Among these genes, 18 contained one or two introns. A maximum-likelihood (ML) phylogenetic analysis based on 33 taxa showed that O. minuta formed a clade with O. japonica. This study will provide a baseline as well as valuable molecular phylogenomic information for various future studies to determine the taxonomic position and phylogenetic relationships of the genus Orostachys.

The complete chloroplast genome sequences of eight Orostachys species: Comparative analysis and assessment of phylogenetic relationships.

We analyzed the complete chloroplast genomes of eight Orostachys species and compared the sequences to those of published chloroplast genomes of the congeneric and closely related genera, Meterostachys and Hylotelephium. The total chloroplast genome length of thirteen species, including the eight species analyzed in this study and the five species analyzed in previous studies, ranged from 149,860 (M. sikokianus) to 151,707 bp (H. verticillatum). The overall GC contents of the genomes were almost identical (37.6 to 37.8%). The thirteen chloroplast genomes each contained 113 unique genes comprising 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Among the annotated genes, sixteen genes contained one or two introns. Although the genome structures of all Orostachys and Hylotelephium species were identical, Meterostachys differed in structure due to a relatively large gene block (trnS-GCU-trnS-GGA) inversion. The nucleotide diversity among the subsect. Orostachys chloroplast genomes was extremely low in all regions, and among the subsect. Appendiculatae, genus Orostachys, and all thirteen chloroplast genomes showed high values of Pi (>0.03) in one, five, or three regions. The phylogenetic analysis showed that Orostachys formed polyphyly, and subsect. Orostachys and Appendiculatae were clustered with Hylotelephium and Meterostachys, respectively, supporting the conclusion that each subsection should be considered as an independent genus. Furthermore, the data supported the taxonomic position of O. margaritifolia and O. iwarenge f. magnus, which were treated as synonyms for O. iwarenge in a previous study, as independent taxa. Our results suggested that O. ramosa and O. japonica f. polycephala were individual variations of O. malacophylla and O. japonica, respectively. The exact taxonomic position of O. latielliptica and the phylogenetic relationship among the three species, O. chongsunensis, O. malacophylla and O. ramosa, should be a topic of future study.

Comparative Transcriptomics for Genes Related to Berberine and Berbamine Biosynthesis in Berberidaceae.

Berberine and berbamine are bioactive compounds of benzylisoquinoline alkaloids (BIAs) present in Berberis species. The contents of berbamine are 20 times higher than berberine in leaf tissues in three closely related species: Berberis koreana, B. thunbergii and B. amurensis. This is the first report on the quantification of berberine compared to the berbamine in the Berberis species. Comparative transcriptome analyses were carried out with mRNAs from the leaf tissues of the three-species. The comparison of the transcriptomes of B. thunbergii and B. amurensis to those of B. koreana, B. thunbergii showed a consistently higher number of differentially expressed genes than B. amurensis in KEGG and DEG analyses. All genes encoding enzymes involved in berberine synthesis were identified and their expressions were variable among the three species. There was a single copy of CYP80A/berbamunine synthase in B. koreana. Methyltransferases and cytochrome P450 mono-oxidases (CYPs) are key enzymes for BIA biosynthesis. The current report contains the copy numbers and other genomic characteristics of the methyltransferases and CYPs in Berberis species. Thus, the contents of the current research are valuable for molecular characterization for the medicinal utilization of the Berberis species.

Characterization of the chloroplast genome of Meterostachys sikokianus (Makino) Nakai (Crassulaceae) and its phylogenetic analysis.

Meterostachys is a monotypic genus of Crassulaceae, though its phylogenetic position remains unclear. Here, we report the complete chloroplast (cp) genome sequence of Meterostachys sikokianus using the Illumina high-throughput sequencing approach. The cp genome was 149,860 bp in length, containing a large single copy (LSC) of 82,293 bp and a small single copy (SSC) of 16,879 bp, which were separated by a pair of 25,344 bp inverted repeats (IRs). The overall GC content of the M. sikokianus cp genome was 37.6%. A total of 113 unique genes were annotated, consisting of 79 protein coding genes (PCGs), 30 transfer RNAs (tRNAs), and four ribosomal RNAs (rRNAs). Among these genes, eighteen contained one or two introns. A maximum-likelihood (ML) phylogenetic analysis based on 30 accessions of Crassulaceae showed that M. sikokianus was most closely related to Orostachys japonica and Orostachys fimbriata.

Gene Expression and Isoform Identification of PacBio Full-Length cDNA Sequences for Berberine Biosynthesis in Berberis koreana.

Berberis koreana is a medicinal plant containing berberine, which is a bioactive compound of the benzylisoquinoline alkaloid (BIA) class. BIA is widely used in the food and drug industry for its health benefits. To investigate the berberine biosynthesis pathway, gene expression analysis was performed in leaves, flowers, and fruits at different stages of growth. This was followed by full-length cDNA sequencing analysis using the PacBio sequencer platform to determine the number of isoforms of those expressed genes. We identified 23,246 full-length unigenes, among which 8479 had more than one isoform. The number of isoforms ranged between two to thirty-one among all genes. Complete isoform analysis was carried out on the unigenes encoding BIA synthesis. Thirteen of the sixteen genes encoding enzymes for berberine synthesis were present in more than one copy. This demonstrates that gene duplication and translation into isoforms may contribute to the functional specificity of the duplicated genes and isoforms in plant alkaloid synthesis. Our study also demonstrated the streamlining of berberine biosynthesis via the absence of genes for enzymes of other BIAs, but the presence of all the genes for berberine biosynthesize in B. koreana. In addition to genes encoding enzymes for the berberine biosynthesis pathway, the genes encoding enzymes for other BIAs were not present in our dataset except for those encoding corytuberine synthase (CTS) and berbamunine synthase (BS). Therefore, this explains how B. koreana produces berberine by blocking the pathways leading to other BIAs, effectively only allowing the pathway to lead to berberine synthesis.

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.

The complete chloroplast genome of Campanula takesimana Nakai from Dokdo Island in Korea (Campanulaceae).

Completed chloroplast genome of Campanula takesimana Nakai isolated from Dokdo island in Korea is 169,719 bp long (GC ratio is 38.8%) and has four subregions: 102,381 bp of large single-copy (37.8%) and 7,750 bp of small single-copy (32.6%) regions are separated by 29,794 bp of inverted repeat (41.3%) regions including 131 genes (87 protein-coding genes, eight rRNAs, and 36 tRNAs). Phylogenetic analyses suggested that C. takesimana from Dokdo Island form a clade with C. takesimana from Ulleungdo Island and that chloroplast genomes of the two accessions are diverged.

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.

An unexpected genetic diversity pattern and a complex demographic history of a rare medicinal herb, Chinese asparagus (Asparagus cochinchinensis) in Korea.

Range-wide population studies of wide spread species are often associated with complex diversity patterns resulting from genetically divergent evolutionary significant units (ESUs). The compound evolutionary history creating such a pattern of diversity can be inferred through molecular analyses. Asparagus cochinchinensis, a medicinally important perennial herb, is in decline due to overharvesting in Korea. Eight A. cochinchinensis populations in Korea and three populations from neighboring countries (China, Japan and Taiwan) were examined using nine nuclear microsatellite loci and three chloroplast microsatellite loci to characterize molecular diversity patterns. The average within-population diversity was limited likely due to long-term bottlenecks observed in all eight populations. High pairwise FST values indicated that the populations have largely diverged, but the divergences were not correlated with geographic distances. Clustering analyses revealed a highly complex spatial structure pattern associated with two ESUs. Approximate Bayesian Computation (ABC) suggests that the two ESUs split about 21,000 BP were independently introduced to Korea approximately 1,800 years ago, and admixed in secondary contact zones. The two ESUs found in our study may have different habitat preferences and growth conditions, implying that the two genetically divergent groups should be considered not only for conservation and management but also for breeding programs in agricultural areas.

Genetic diagnosis of a rare myrmecochorous species, Plagiorhegma dubium (Berberidaceae): Historical genetic bottlenecks and strong spatial structures among populations.

Distribution of genetic variation over time and space is relevant to demographic histories and tightly linked to ecological disturbances as well as evolutionary potential of an organism. Therefore, understanding the pattern of genetic diversity is a primary step in conservation and management projects for rare and threatened plant species. We used eight microsatellite markers to examine the level of genetic diversity, spatial structure, and demographic history of Plagiorhegma dubium, a rare myrmecochorous herb, populations sampled across northeast Asia and Siberia. We found low within-population genetic variation associated with historical bottlenecks. Although pairwise F ST values were not much higher than the ones found in similar life form species, STRUCTURE and PCoA revealed a clear broadscale spatial pattern of genetic structure. Bayesian clustering (best K = 6) and PCoA identified three populations that are distinctive from neighboring populations in the Korean peninsula, which suggests potential units for conservation and management plans in Korea. MIGRATE-N and BAYESASS showed that both contemporary (0.003-0.045) and historical migration rates (2 × e-5-4.6 × e-4) were low. Our findings provide a good example, where genetic considerations should be integrated for conservation and management plans of rare and threatened species.

Development and characterization of 30 microsatellite loci for Plagiorhegma dubium (Berberidaceae).

PREMISE OF THE STUDY: Plagiorhegma dubium (Berberidaceae) has been listed as an endangered species in Korea due to extensive collection and destruction of natural habitats. In this study, 30 microsatellite loci, including 25 polymorphic loci, were developed for P. dubium for use in population-level genetic analyses. METHODS AND RESULTS: We carried out transcriptome sequencing and isolated a total of 30 expressed sequence tag-simple sequence repeat markers from P. dubium using Illumina HiSeq high-throughput sequencing. To test utility of the developed markers, we genotyped 60 individuals from three populations and estimated the number of alleles and levels of observed and expected heterozygosity. Expected heterozygosity levels ranged from 0.000 to 0.594, 0.000 to 1.000, and 0.000 to 0.744 in the three populations, respectively. CONCLUSIONS: These transcriptome-derived simple sequence repeat markers are highly polymorphic and can be widely used in characterization of the endangered P. dubium. Population genetic studies with these markers will provide valuable insights for conservation by unraveling evolutionary patterns of P. dubium.

Impatiens bokorensis (Balsaminaceae), a new species from Cambodia.

Impatiens bokorensis, a new species of family Balsaminaceae from Phnum Bokor National Park in southwestern Cambodia, is described and illustrated. The species is similar to Impatiens patula, but is readily distinguished by the orbicular-obovate dorsal petal, shorter pedicels and larger seeds.

Development of microsatellite markers based on expressed sequence tags in Asparagus cochinchinensis (Asparagaceae).

PREMISE OF THE STUDY: Transcriptome-derived simple sequence repeat (SSR) markers were developed in Asparagus cochinchinensis (Asparagaceae). Due to its application in traditional medicine, its wild populations are threatened by over-collection even in protected areas, requiring immediate conservation efforts. METHODS AND RESULTS: Based on transcriptome data of A. cochinchinensis, 96 primer pairs with two to seven alleles per locus were selected for initial validation; of those, 27 primer pairs amplified across all samples, resulting in 15 polymorphic and 12 monomorphic microsatellite markers. The usefulness of these markers was assessed in 60 individuals representing three populations of A. cochinchinensis. Observed and expected heterozygosity values ranged from 0.050 to 0.950 and 0.049 to 0.626, respectively. Cross-species amplification of the 27 markers was tested in the related species A. rigidulus and A. schoberioides. CONCLUSIONS: These polymorphic, transcriptome-derived SSR markers can be used as molecular markers to study population genetics and ecological conservation in A. cochinchinensis and related taxa.

Development of microsatellite markers for Viscum coloratum (Santalaceae) and their application to wild populations.

PREMISE OF THE STUDY: Microsatellite primers were developed for Viscum coloratum (Santalaceae), a semiparasitic medicinal plant that is known for its anticancer properties. Due to excessive human harvesting and loss of suitable habitat of its populations, it has become a potentially threatened species requiring immediate conservation efforts. METHODS AND RESULTS: Based on transcriptome data for V. coloratum, 124 primer pairs were randomly selected for initial validation, of which 19 yielded polymorphic microsatellite loci, with two to six alleles per locus. The usefulness of these markers was assessed for 60 individuals representing three populations of V. coloratum. Observed and expected heterozygosity values ranged from 0.033 to 0.833 and 0.032 to 0.672, respectively. Cross-species amplification for 19 loci in the related species V. album was conducted. CONCLUSIONS: The 19 newly developed loci are expected to be useful for studying the population genetics and ecological conservation of V. coloratum.

Synthesis and Characterization of Water-Soluble Conjugated Oligoelectrolytes for Near-Infrared Fluorescence Biological Imaging.

Near-infrared (NIR) fluorophores attract increasing attention as a molecular marker (or probe) for in vivo and in vitro biological fluorescence imaging. Three types of new NIR fluorescent conjugated oligoelectrolytes (COEs: Q-FlTBTTFl, Q-FlBBTFl, and Q-FlTBBTTFl) are synthesized with quaternized ammonium ionic groups in their side-chains for water solubility. The emission wavelength is modulated in the range 600-1300 nm, by adjusting the intramolecular charge transfer in the molecular backbone based on the electron-rich fluorene (and/or thiophene) and electron-deficient benzo[2,1,3]thiadiazole (or benzo[1,2-c:4,5-c']bis[1,2,5]thiadiazole) moieties. The COEs show a remarkably larger Stokes shift (147-276 nm) compared to commercial rhodamine and cyanine dyes in water, avoiding self-quenching and interference from the excitation backscattered light. The photoluminescence (PL) quantum efficiency is improved substantially by up to 27.8% in water by fabricating a vesicular complex, COE/v, with a block ionomer, poly[(ethylene oxide)-block-(sodium 2-acrylamido-2-methyl-1-propanesulfonate)]. In vitro cellular uptake images with the COEs are obtained with good biocompatibility by confocal single-photon and two-photon microscopy. The ex vivo and in vivo images of a mouse xenograft model treated with the Q-FlBBTFl/v exhibit a substantially stronger fluorescence signal at the tumor site than at the other organs, highlighting the potential of the COE/v as an NIR fluorescent imaging agent for the diagnosis of cancer.