Double migration of the endangered Tricyrtis formosana (Liliaceae) in Japan.

The Ryukyu Islands of Japan are a biodiversity hotspot due to geographical and historical factors. Tricyrtis formosana is a perennial herbaceous plant that commonly found in Taiwan. But only a few populations have been identified in a limited habitat on Iriomote Island, while populations of unknown origin occur near human settlements in an area on the main island of Okinawa. To better understand these populations of the phylogenetic uniqueness and intrinsic vulnerability, we conducted comparative analyses including (1) phylogeny and population structure with MIG-seq data, (2) photosynthesis-related traits of plants grown under common conditions and (3) transcriptome analysis to detect deleterious variations. Results revealed that T. formosana was split into two clades by the congeners and that Iriomote and Okinawa populations independently derived from ancestral Taiwanese populations in each clade. Photosynthetic efficiency was lowest in the Iriomote population, followed by Okinawa and Taiwan. Transcriptome analysis showed that the Iriomote population accumulated more deleterious variations, suggesting intrinsic vulnerability. These results indicate that each T. formosana population in Japan is phylogenetically unique and has been independently dispersed from Taiwan, and that the Iriomote population presents a high conservation difficulty with a unique photosynthesis-related characteristic and a larger amount of deleterious variations.

Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments.

Paris polyphylla var. yunnanensis is a well-known medicinal plant that is mainly distributed in Southwest China; however, its genetic diversity and biodiversity processes are poorly understood. In this study, the sequences of cpDNA trnL-trnF fragments of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis were amplified, sequenced, and aligned to study the population genetics of this species. Genetic diversity was analyzed based on nucleotide diversity, haplotype diversity, Watterson diversity, population-level diversity, and species-level genetic diversity. Genetic structure and genetic differentiation were explored using haplotype distribution maps and genetic distance matrices. A total of 15 haplotypes were identified in the 32 populations of P. polyphylla var. yunnanensis. Five unique haplotypes were identified from the fourteen haplotypes of the cultivated populations, while only one unique haplotype was identified from the ten haplotypes of the wild populations. The haplotype richness and genetic diversity of the cultivated populations were higher than those of the wild populations (HT = 0.900 vs. 0.861). In addition, there were no statistically significant correlations between geographic distance and genetic distance in the cultivated populations (r = 0.16, p > 0.05), whereas there was a significant correlation between geographical distance and genetic structure in the wild populations (r = 0.32, p > 0.05), indicating that there was a geographical and genetic connection between the wild populations. There was only 2.5% genetic variation between the wild populations and cultivated populations, indicating no obvious genetic differentiation between the wild and cultivated populations. Overall, the genetic background of the cultivated populations was complex, and it was hypothesized that the unique haplotypes and higher diversity of the cultivated populations were caused by the mixed provenance of the cultivated populations.

Dynamic variation of Paris polyphylla root-associated microbiome assembly with planting years.

MAIN CONCLUSION: P. polyphylla selectively enriches beneficial microorganisms to help their growth. Paris polyphylla (P. polyphylla) is an important perennial plant for Chinese traditional medicine. Uncovering the interaction between P. polyphylla and the related microorganisms would help to utilize and cultivate P. polyphylla. However, studies focusing on P. polyphylla and related microbes are scarce, especially on the assembly mechanisms and dynamics of the P. polyphylla microbiome. High-throughput sequencing of the 16S rRNA genes was implemented to investigate the diversity, community assembly process and molecular ecological network of the bacterial communities in three root compartments (bulk soil, rhizosphere, and root endosphere) across three years. Our results demonstrated that the composition and assembly process of the microbial community in different compartments varied greatly and were strongly affected by planting years. Bacterial diversity was reduced from bulk soils to rhizosphere soils to root endosphere and varied over time. Microorganisms benefit to plants was selectively enriched in P. polyphylla roots as was its core microbiome, including Pseudomonas, Rhizobium, Steroidobacter, Sphingobium and Agrobacterium. The network's complexity and the proportion of stochasticity in the community assembly process increased. Besides, nitrogen metabolism, carbon metabolism, phosphonate and phosphinate metabolism genes in bulk soils increased over time. These findings suggest that P. polyphylla exerts a selective effect to enrich the beneficial microorganisms and proves the sequential increasing selection pressure with P. polyphylla growth. Our work adds to the understanding of the dynamic processes of plant-associated microbial community assembly, guides the selection and application timing of P. polyphylla-associated microbial inoculants and is vital for sustainable agriculture.

Complete genome sequence analysis of a new potyvirus isolated from Paris polyphylla var. yunnanensis.

The complete genome sequence of a new potyvirus from Paris polyphylla var. yunnanensis was determined. Its genomic RNA consists of 9571 nucleotides (nt), excluding the 3'-terminal poly(A) tail, containing the typical open reading frame (ORF) of potyviruses and encoding a putative large polyprotein of 3061 amino acids. The virus shares 54.20%-59.60% nt sequence identity and 51.80%-57.90% amino acid sequence identity with other potyviruses. Proteolytic cleavage sites and conserved motifs of potyviruses were identified in the polyprotein and within individual proteins. Phylogenetic analysis indicated that the virus was most closely related to lily yellow mosaic virus. The results suggest that the virus should be classified as a member of a novel species within the genus Potyvirus, and we have tentatively named this virus "Paris yunnanensis mosaic chlorotic virus" (PyMCV).

Comparative transcriptome analysis reveals key genes for polyphyllin difference in five Paris species.

Paris species accumulate a large amount of steroidal saponins, which have numerous pharmacological activities and have become an essential component in many patented drugs. However, only two among all Paris species. Paris are identified as official sources due to high level of bioactive compounds. To clarify the composition of steroidal saponins and the molecular basis behind the differences between species, we investigated transcriptome and metabolic profiles of leaves and rhizomes in Paris polyphylla var. chinensis (PPC), Paris polyphylla var. yunnanensis (PPY), Paris polyphylla var. stenophylla (PPS), Paris fargesii (PF), and Paris mairei (PM). Phytochemical results displayed that the accumulation of steroidal saponins was tissue- and species-specific. PF and PPS contained more steroidal saponins in leaves than rhizomes, while PPY accumulated more steroidal saponins in rhizomes than leaves. PPC and PM contained similar amounts of steroidal saponins in leaves and rhizomes. Transcriptome analysis illustrated that most differentially expressed genes related to the biosynthesis of steroidal saponins were abundantly expressed in rhizomes than leaves. Meanwhile, more biosynthetic genes had significant correlations with steroidal saponins in rhizomes than in leaves. The result of CCA indicated that ACAT, DXS, DWF1, and CYP90 constrained 97.35% of the variance in bioactive compounds in leaves, whereas CYP72, UGT73, ACAT, and GPPS constrained 98.61% of the variance in phytochemicals in rhizomes. This study provided critical information for enhancing the production of steroidal saponins by biotechnological approaches and methodologies.

RNA-Seq analysis reveals the important co-expressed genes associated with polyphyllin biosynthesis during the developmental stages of Paris polyphylla.

BACKGROUND: Plants synthesize metabolites to adapt to a continuously changing environment. Metabolite biosynthesis often occurs in response to the tissue-specific combinatorial developmental cues that are transcriptionally regulated. Polyphyllins are the major bioactive components in Paris species that demonstrate hemostatic, anti-inflammatory and antitumor effects and have considerable market demands. However, the mechanisms underlying polyphyllin biosynthesis and regulation during plant development have not been fully elucidated. RESULTS: Tissue samples of P. polyphylla var. yunnanensis during the four dominant developmental stages were collected and investigated using high-performance liquid chromatography and RNA sequencing. Polyphyllin concentrations in the different tissues were found to be highly dynamic across developmental stages. Specifically, decreasing trends in polyphyllin concentration were observed in the aerial vegetative tissues, whereas an increasing trend was observed in the rhizomes. Consistent with the aforementioned polyphyllin concentration trends, different patterns of spatiotemporal gene expression in the vegetative tissues were found to be closely related with polyphyllin biosynthesis. Additionally, molecular dissection of the pathway components revealed 137 candidate genes involved in the upstream pathway of polyphyllin backbone biosynthesis. Furthermore, gene co-expression network analysis revealed 74 transcription factor genes and one transporter gene associated with polyphyllin biosynthesis and allocation. CONCLUSIONS: Our findings outline the framework for understanding the biosynthesis and accumulation of polyphyllins during plant development and contribute to future research in elucidating the molecular mechanism underlying polyphyllin regulation and accumulation in P. polyphylla.

Complete genome sequence analysis of Paris alphapartitivirus 1: a novel member of the genus Alphapartitivirus infecting Paris polyphylla var. yunnanensis.

A novel double-stranded RNA (dsRNA) virus, tentatively named "Paris alphapartitivirus 1" (ParAPV1, OL960006-OL960007), was detected in Paris polyphylla var. yunnanensis plants exhibiting leaf chlorosis and shrinkage symptoms in Yunnan. Its complete genome sequence was determined using Illumina and Sanger sequencing. ParAPV1 has a bipartite genome that consists of dsRNA1 (1,917 bp) encoding the viral RNA-dependent RNA polymerase (RdRp), and dsRNA2 (1,818 bp) encoding the putative coat protein (CP). Sequence comparisons showed that the RdRp and CP of ParAPV1 are most similar to those of pear alphapartitivirus (PpPV2), with 69.97% and 54.21% amino acid sequence identities respectively. Phylogenetic analysis of the RdRp amino acid sequences of ParAPV1 and other partitiviruses showed that ParAPV1 cluster with viruses in a clade containing alphapartitiviruses, and that its closest known relatives are PpPV2 (BBA66577) and rose partitivirus (RoPV, ANQ45203S). Taken together, these results suggest that ParAPV1 should be regarded as a new member of genus Alphapartitivirus in the family Partitiviridae. This is the first report of a partitivirus infecting P. polyphylla var. yunnanensis.

Four Novel Mycoviruses from the Hypovirulent Botrytis cinerea SZ-2-3y Isolate from Paris polyphylla: Molecular Characterisation and Mitoviral Sequence Transboundary Entry into Plants.

A hypovirulent SZ-2-3y strain isolated from diseased Paris polyphylla was identified as Botrytis cinerea. Interestingly, SZ-2-3y was coinfected with a mitovirus, two botouliviruses, and a 3074 nt fusarivirus, designated Botrytis cinerea fusarivirus 8 (BcFV8); it shares an 87.2% sequence identity with the previously identified Botrytis cinerea fusarivirus 6 (BcFV6). The full-length 2945 nt genome sequence of the mitovirus, termed Botrytis cinerea mitovirus 10 (BcMV10), shares a 54% sequence identity with Fusarium boothii mitovirus 1 (FbMV1), and clusters with fungus mitoviruses, plant mitoviruses and plant mitochondria; hence BcMV10 is a new Mitoviridae member. The full-length 2759 nt and 2812 nt genome sequences of the other two botouliviruses, named Botrytis cinerea botoulivirus 18 and 19 (BcBoV18 and 19), share a 40% amino acid sequence identity with RNA-dependent RNA polymerase protein (RdRp), and these are new members of the Botoulivirus genus of Botourmiaviridae. Horizontal transmission analysis showed that BcBoV18, BcBoV19 and BcFV8 are not related to hypovirulence, suggesting that BcMV10 may induce hypovirulence. Intriguingly, a partial BcMV10 sequence was detected in cucumber plants inoculated with SZ-2-3y mycelium or pXT1/BcMV10 agrobacterium. In conclusion, we identified a hypovirulent SZ-2-3y fungal strain from P. polyphylla, coinfected with four novel mycoviruses that could serve as potential biocontrol agents. Our findings provide evidence of cross-kingdom mycoviral sequence transmission.

Transcriptomic Insight into Underground Floral Differentiation in Erythronium japonicum.

Erythronium japonicum Decne (Liliaceae) flowers in early spring after overwintering. Its sexual reproduction process includes an underground development process of floral organs, but the underlying molecular mechanisms are obscure. The present study is aimed at exploring the transcriptional changes and key genes involved at underground floral developmental stages, including flower primordium differentiation, perianth differentiation, stamen differentiation, and pistil differentiation in E. japonicum. Multistage high-quality transcriptomic data resulted in identifying putative candidate genes for underground floral differentiation in E. japonicum. A total of 174,408 unigenes were identified, 28,508 of which were differentially expressed genes (DEGs) at different floral developmental stages, while only 44 genes were identified with conserved regulation between different stages. Further annotation of DEGs resulted in the identification of 270 DEGs specific to floral differentiation. In addition, ELF3, PHD, cullin 1, SE14, ZSWIM3, GIGNATEA, and SERPIN B were identified as potential candidate genes involved in the regulation of floral differentiation. Besides, we explored transcription factors with differential regulation at different developmental stages and identified bHLH, FAR1, mTERF, MYB-related, NAC, Tify, and WRKY TFs for their potential involvement in the underground floral differentiation process. Together, these results laid the foundation for future molecular works to improve our understanding of the underground floral differentiation process and its genetic regulation in E. japonicum.

Phylogeny and Comparative Analysis for the Plastid Genomes of Five Tulipa (Liliaceae).

Species of Tulipa (Liliaceae) are of great horticultural importance and are distributed across Europe, North Africa, and Asia. The Tien Shan Mountain is one of the primary diversity centres of Tulipa, but the molecular studies of Tulipa species from this location are lacking. In our study, we assembled four Tulipa plastid genomes from the Tien Shan Mountains, T. altaica, T. iliensis, T. patens, and T. thianschanica, combined with the plastid genome of T. sylvestris to compare against other Liliaceae plastid genomes. We focussed on the species diversity and evolution of their plastid genomes. The five Tulipa plastid genomes proved highly similar in overall size (151,691-152,088 bp), structure, gene order, and content. With comparative analysis, we chose 7 mononucleotide SSRs from the Tulipa species that could be used in further population studies. Phylogenetic analyses based on 24 plastid genomes robustly supported the monophyly of Tulipa and the sister relationship between Tulipa and Amana, Erythronium. T. iliensis, T. thianschanica, and T. altaica were clustered together, and T. patens was clustered with T. sylvestris, with our results clearly demonstrating the relationships between these five Tulipa species. Our results provide a more comprehensive understanding of the phylogenomics and comparative genomics of Tulipa.

[Phylogeography of Paris poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences].

Phylogeography is a research hotspot in the field of the genetic diversity and core germplasm construction of endangered rare plants. Paris polyphylla var. yunnanensis is a rare plant species mainly distributed in China. Wild individuals have been overexploited for the last few decades because of increasing demand for such medicines. Therefore, it is great significance to study the phylogeography of P. poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences. In this study, chloroplast genes trnL-trnF were used in the phylogeography analysis of 15 wild and 17 cultivated populations of P. polyphylla var. yunnanensis. This study revealed that based on the results of neutrality tests and mismatch analysis, the rapid expansion of wild population has not been detected in P. polyphylla var. yunnanensis. After aligning and sorting the obtained cpDNA sequences, a total of 15 haplotypes were detected in all 32 populations. One haplotype was unique to the wild population, and 5 haplotypes were unique to the cultivated population. It can be seen that the haplotype richness of cultivated population was higher than that of wild population. The wild populations of P. polyphylla var. yunnanensis were divided into two groups according to evolutionary relationship of haplotypes and distribution map of haplotypes. The haplotype of branch Ⅰ was mainly distributed in Guizhou, and the haplotype of branch Ⅱ was located in Yunnan and Huidong, Sichuan. Therefore, it's speculated that Guizhou and the west Yunnan region may be glacial refuge in the evolutionary history of wild populations of P. polyphylla var. yunnanensis, and in order to protect the wild resources more effectively, wild populations of P. polyphylla var. yunnanensis in these two areas should be included in the protection zone.

A cycloartenol synthase from the steroidal saponin biosynthesis pathway of Paris polyphylla.

Steroidal saponins named polyphyllin are the major active components of Paris polyphylla. Cycloartenol synthase (CAS) is a key enzyme that catalyzes the formation of the sterol scaffold. In this study, we cloned a putative CAS gene from Paris polyphylla. Heterologous expression in yeast indicated that PpCAS can convert 2,3-oxidosqualene into cycloartenol. qRT-PCR analysis showed that the expression of PpCAS was highest in leaves and lowest in roots. To our best knowledge, this is the first report of the functional characterization of cycloartenol synthase from Paris polyphylla, which lays the foundation for further analysis of the biosynthesis pathway of polyphyllins.[Formula: see text].

Phylogeography of Paris poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences / 中国中药杂志

Phylogeography is a research hotspot in the field of the genetic diversity and core germplasm construction of endangered rare plants. Paris polyphylla var. yunnanensis is a rare plant species mainly distributed in China. Wild individuals have been overexploited for the last few decades because of increasing demand for such medicines. Therefore, it is great significance to study the phylogeography of P. poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences. In this study, chloroplast genes trnL-trnF were used in the phylogeography analysis of 15 wild and 17 cultivated populations of P. polyphylla var. yunnanensis. This study revealed that based on the results of neutrality tests and mismatch analysis, the rapid expansion of wild population has not been detected in P. polyphylla var. yunnanensis. After aligning and sorting the obtained cpDNA sequences, a total of 15 haplotypes were detected in all 32 populations. One haplotype was unique to the wild population, and 5 haplotypes were unique to the cultivated population. It can be seen that the haplotype richness of cultivated population was higher than that of wild population. The wild populations of P. polyphylla var. yunnanensis were divided into two groups according to evolutionary relationship of haplotypes and distribution map of haplotypes. The haplotype of branch Ⅰ was mainly distributed in Guizhou, and the haplotype of branch Ⅱ was located in Yunnan and Huidong, Sichuan. Therefore, it's speculated that Guizhou and the west Yunnan region may be glacial refuge in the evolutionary history of wild populations of P. polyphylla var. yunnanensis, and in order to protect the wild resources more effectively, wild populations of P. polyphylla var. yunnanensis in these two areas should be included in the protection zone.

Pre-quaternary diversification and glacial demographic expansions of Cardiocrinum (Liliaceae) in temperate forest biomes of Sino-Japanese Floristic Region.

The Sino-Japanese Floristic Region (SJFR) in East Asia is one of the most diverse temperate floras in the world. However, the relative influence of Neogene palaeogeographical changes and Quaternary climatic fluctuations as causal mechanisms on species diversification remains largely controversial, because most divergence time estimates were inferred from single-locus data and have limited geographic or taxonomic sampling. To evaluate these influences, we use SNP markers from restriction site-associated DNA sequencing (RAD-Seq) loci and expressed sequence tags-simple sequence repeat (EST-SSR) markers to investigate the levels of genetic variation, speciation and demographic history of the temperate-deciduous forest (TDF) endemic Cardiocrinum (Endlicher) Lindley (Liliaceae), a genus comprising three species in China (C. giganteum, C. cathayanum) and Japan (C. cordatum). Phylogenomic and population genomic coalescent-based analyses demonstrated that Late Neogene tectonic/climatic events triggered speciation of Cardiocrinum, and Pleistocene climatic fluctuations had limited influence on its divergence history. Population demographic inference using Approximate Bayesian Computation from EST-SSRs and palaeoclimatic niche models both indicated that all three Cardiocrinum species experienced population expansions during the transition from the LIG to the LGM. We also discussed the implications of these results on the conservation of montane TDF species in the SJFR under ongoing environmental change. Our results improve our understanding of how the constituents of montane TDF across the SJFR responded to previous periods of rapid climate and environmental change in terms of speciation and population demographic processes.

Transcriptome analyses of Paris polyphylla var. chinensis, Ypsilandra thibetica, and Polygonatum kingianum characterize their steroidal saponin biosynthesis pathway.

Steroidal saponins, one of the most diverse groups of plant-derived natural products, elicit biological and pharmacological activities; however, the genes involved in their biosynthesis and the corresponding biosynthetic pathway in monocotyledon plants remain unclear. This study aimed to identify genes involved in the biosynthesis of steroidal saponins by performing a comparative analysis among transcriptomes of Paris polyphylla var. chinensis (PPC), Ypsilandra thibetica (YT), and Polygonatum kingianum (PK). De novo transcriptome assemblies generated 57,537, 140,420, and 151,773 unigenes from PPC, YT, and PK, respectively, of which 56.54, 47.81, and 44.30% were successfully annotated, respectively. Among the transcriptomes for PPC, YT, and PK, we identified 194, 169, and 131; 17, 14, and 26; and, 80, 122, and 113 unigenes corresponding to terpenoid backbone biosynthesis; sesquiterpenoid and triterpenoid biosynthesis; and, steroid biosynthesis pathways, respectively. These genes are putatively involved in the biosynthesis of cholesterol that is the primary precursor of steroidal saponins. Phylogenetic analyses indicated that lanosterol synthase may be exclusive to dicotyledon plant species, and the cytochrome P450 unigenes were closely related to clusters CYP90B1 and CYP734A1, which are UDP-glycosyltransferases unigenes homologous with the UGT73 family. Thus, unigenes of ß-glucosidase may be candidate genes for catalysis of later period modifications of the steroidal saponin skeleton. Our data provide evidence to support the hypothesis that monocotyledons biosynthesize steroidal saponins from cholesterol via the cycloartenol pathway.

Calochortus gunnisonii furthers evidence for the complex genetic legacy of historical climate change in the southern Rocky Mountains.

PREMISE OF THE STUDY: Climate cycles of the Quaternary have impacted plants at a global scale, leaving behind a complex genetic legacy. Species of the northern Rocky Mountains of North America were exposed to more uniform glacial patterns than the central and southern ranges, where synergistic relationships between temperature and precipitation caused differences in the timing and extent of glacier onset. We examined the genetic impacts of climate oscillations on Calochortus gunnisonii (Liliaceae) in the central and southern Rocky Mountains. METHODS: Populations were sampled from disjunct mountain ranges across the basins of Wyoming and northern and central Colorado. Allelic data from nuclear microsatellites and plastid sequences (trnV-ndhC, petA-psbJ, and rpl16) were used to examine patterns of genetic structure between and among populations along the southern Rocky Mountain corridor. KEY RESULTS: We infer considerable population structure concordant with mountain range of origin. Clustering analysis supports separate north and south genetic clusters on either side of major basins in Wyoming, suggesting that populations were maintained in two distinct refugia. Additionally, populations within the Sierra Madre Range of southern Wyoming show localized, divergent genetic signal indicative of a third potential glacial refugium. By contrast, recent genetic admixture is observed in the Laramie, Medicine Bow, and Front ranges, where population expansion from glacial refugia has likely occurred. CONCLUSIONS: We conclude that during climate cycles of the Quaternary, C. gunnisonii experienced periods of population expansion and reduction, habitat fragmentation, isolation in three or more refugia, and admixture mirroring genetic impacts of other southern Rocky Mountains organisms.

Molecular phylogenetics and historical biogeography of the tribe Lilieae (Liliaceae): bi-directional dispersal between biodiversity hotspots in Eurasia.

Background and Aims: The role played by the Qinghai-Tibet Plateau (QTP) in the organismal diversification and biogeography of plants in the Northern Hemisphere has attracted much attention from evolutionary biologists. Here we use tribe Lilieae (Liliaceae), including primarily temperate and alpine lineages with disjunct distributions in the North Temperate Zone, as a case study to shed light upon these processes. Methods: Using 191 taxa (five outgroup taxa) comprising more than 60 % of extant Lilieae species across the entire geographical range, we analyse phylogenetic relationships based on three plastid markers (matK, rbcL, rpl16) and nuclear ITS. Divergence time estimation and ancestral range reconstruction were further inferred. Key Results: The results support a monophyletic Lilieae divided into four clades. Lilium is nested within Fritillaria, which is paraphyletic and partitioned into two clades, New World and Old World, in the chloroplast DNA (cpDNA) analysis. Incongruences between the ITS and cpDNA trees may be explained by divergent ITS paralogues and hybridization. Lilieae originated around 40-49 (28-67) Mya and probably diversified in the QTP region with four major clades that were established during the Oligocene and the Early Miocene. Uplift of the QTP and climatic changes probably drove early diversification of Lilieae in the QTP region. A rapid radiation occurred during the Late Miocene and the Pleistocene, coinciding temporally with recent orogenic process in the QTP region and climatic oscillations. Several lineages dispersed out of the QTP. Conclusions: Lineage persistence and explosive radiation were important processes for establishing high species diversity of Lilieae in the QTP region. Both long-distance dispersal and migration across Beringia probably contributed to the modern distribution range of Lilieae. Our study shows that biotic interchanges between the QTP region and Irano-Turanian region and the Mediterranean Basin were bi-directional, suggesting the latter was a secondary centre of diversity.

Comparison of genetic variation between northern and southern populations of Lilium cernuum (Liliaceae): Implications for Pleistocene refugia.

The so-called "Baekdudaegan" (BDDG), a mountain range that stretches along the Korean Peninsula, has been recently proposed as a major "southern" glacial refugium for boreal or temperate plant species based on palaeoecological and, especially, genetic data. Genetic studies comparing genetic variation between population occurring on the BDDG and more northern ones (i.e. in NE China and/or in Russian Far East) are, however, still too few to draw firm conclusions on the role of the BDDG as a refugium and a source for possible northward post-glacial recolonizations. In order to fill this gap, we selected a boreal/temperate herb, Lilium cernuum, and compared levels of allozyme-based genetic diversity of five populations from NE China with five populations from South Korea (home of its hypothesized refuge areas). As a complementary tool, we used the maximum entropy algorithm implemented in MaxEnt to infer the species' potential distribution for the present time, which was projected to different past climate scenarios for the Last Glacial Maximum (LGM). Permutation tests revealed that Korean populations harbored significantly higher levels of within-population genetic variation than those from NE China (expected heterozygosity = 0.173 vs. 0.095, respectively). Our results suggest that the lowered levels of genetic diversity in NE Chinese populations might be due to founder effects associated with post-glacial migration from southern regions. Congruent with genetic data, past distribution models showed higher probability of occurrence in southern ranges than in northern ones during the LGM. In addition, a positive correlation was detected between the expected heterozygosity and environmental LGM suitability. From a conservation perspective, our results further suggest that the southern populations in South Korea may be particularly worthy of protection.

Comparative chloroplast genomes of Paris Sect. Marmorata: insights into repeat regions and evolutionary implications.

BACKGROUND: Species of Paris Sect. Marmorata are valuable medicinal plants to synthesize steroidal saponins with effective pharmacological therapy. However, the wild resources of the species are threatened by plundering exploitation before the molecular genetics studies uncover the genomes and evolutionary significance. Thus, the availability of complete chloroplast genome sequences of Sect. Marmorata is necessary and crucial to the understanding the plastome evolution of this section and facilitating future population genetics studies. Here, we determined chloroplast genomes of Sect. Marmorata, and conducted the whole chloroplast genome comparison. RESULTS: This study presented detailed sequences and structural variations of chloroplast genomes of Sect. Marmorata. Over 40 large repeats and approximately 130 simple sequence repeats as well as a group of genomic hotspots were detected. Inverted repeat contraction of this section was inferred via comparing the chloroplast genomes with the one of P. verticillata. Additionally, almost all the plastid protein coding genes were found to prefer ending with A/U. Mutation bias and selection pressure predominately shaped the codon bias of most genes. And most of the genes underwent purifying selection, whereas photosynthetic genes experienced a relatively relaxed purifying selection. CONCLUSIONS: Repeat sequences and hotspot regions can be scanned to detect the intraspecific and interspecific variability, and selected to infer the phylogenetic relationships of Sect. Marmorata and other species in subgenus Daiswa. Mutation and natural selection were the main forces to drive the codon bias pattern of most plastid protein coding genes. Therefore, this study enhances the understanding about evolution of Sect. Marmorata from the chloroplast genome, and provide genomic insights into genetic analyses of Sect. Marmorata.

Transcriptome-Wide Identification and Prediction of miRNAs and Their Targets in Paris polyphylla var. yunnanensis by High-Throughput Sequencing Analysis.

Long dormancy period of seeds limits the large-scale artificial cultivation of the scarce Paris polyphylla var. yunnanensis, an important traditional Chinese medicine. Characterizing miRNAs and their targets is crucial to understanding the role of miRNAs during seed dormancy in this species. Considering the limited genome information of this species, we first sequenced and assembled the transcriptome data of dormant seeds and their seed coats as the reference genome. A total of 146,671 unigenes with an average length of 923 bp were identified and showed functional diversity based on different annotation methods. Two small RNA libraries from respective seeds and seed coats were sequenced and the combining data indicates that 263 conserved miRNAs belonging to at least 83 families and 768 novel miRNAs in 1174 transcripts were found. The annotations of the predicted putative targets of miRNAs suggest that these miRNAs were mainly involved in the cell, metabolism and genetic information processing by direct and indirect regulation patterns in dormant seeds of P. polyphylla var. yunnanensis. Therefore, we provide the first known miRNA profiles and their targets, which will assist with further study of the molecular mechanism of seed dormancy in P. polyphylla var. yunnanensis.