Phylogeography of the Sino-Himalayan fern Lepisorus clathratus on "the roof of the world"

Uncovering the determinants of biodiversity hotspots in China: Evidence from the drivers of multiple diversity metrics on insect assemblages and implications for conservation

Understanding large-scale biodiversity patterns and underlying mechanisms during the formation process is essential for guiding conservation efforts. However, previous studies on the identification and formation mechanism of diversity hotspots in China were often limited to a single (alpha) diversity metric, while multiple (beta or zeta) diversity has rarely been used for exploring drivers and conservation actions. Here, a comprehensive species distribution dataset consisting of representative families of three insect orders was compiled to explore biodiversity hotspots based on different algorithms. Furthermore, to assess the effects of environmental factors on hotspots, we fitted generalized additive mixed-effects models (GAMMs) for species richness, generalized dissimilarity models (GDMs) and multi-site generalized dissimilarity modeling (MS-GDM) for the total beta and zeta diversity. Our results showed that biodiversity hotspots were mainly concentrated in central and southern China, especially in mountainous areas with complex topography, which indicated the insects' affinity to montane environments. Further analyses based on multiple models showed that water-energy factors exerted the strongest explanatory power for the insect assemblage diversity in hotspots of both alpha and beta (or zeta) levels. Additionally, anthropogenic factors also exerted a significant effect on hotspots, and this effect was higher for beta diversity than for alpha diversity. Overall, our study elucidates a comprehensive analysis of the identification and underlying mechanism of biodiversity hotspots in China. Despite several limitations, we still believe that our findings can provide some new insights for conservation efforts in Chinese hotspots.

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

A comparative plastomic analysis of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow and implication of the origin of Chinese jujube

Comparative plastomics can be used to explicitly dissect various types of plastome variation. In the present study, the plastome variation pattern of Ziziphus jujuba var. spinosa (also called sour jujube) and its phylogenomic relationship with Chinese jujube were investigated. Plastomes of 21 sour jujube individuals were sequenced and assembled. The length of the sour jujube plastomes ranged between 159399 and 161279 bp. The plastomes exhibited collinearity of structure, gene order and content. The most divergent regions were located in the intergenic spacers, such as trnR-UCU-atpA and psbZ-trnG-UCC. Sliding window analysis demonstrated that the sequence variation among the sour jujube plastomes was relatively low. Sixty-two to 76 SSRs with 4 motif types were identified in the sour jujube plastomes with a predominant motif type of A/T. Three protein-coding genes exhibited higher nonsynonymous/synonymous substitution ratios, indicating that these genes may undergo positive selection. A total of 80 SNPs were detected and 1266 potential RNA editing sites of 23 protein-coding genes were predicted. In the phylogenomic tree constructed, sour jujube has a sister relationship to Chinese jujube, which indicates that Chinese jujube may have originated or been domesticated from sour jujube. The present study explicitly investigated the individual-level plastome variation of sour jujube and provides potential valuable molecular markers for future genetic-related study of this lineage.

Phylogeographical Pattern and Population Evolution History of Indigenous Elymus sibiricus L. on Qinghai-Tibetan Plateau

Elymus sibiricus L. is a perennial allotetraploid belonging to Triticeae of Poaceae, Elymus L., as the type species of genus Elymus L. The existing geographical distribution pattern and genetic spatial structure of E. sibiricus on Qinghai-Tibetan Plateau (QTP) are not yet clear. In this study, population genetic structure and demography history of 216 individuals from 44 E. sibiricus populations on QTP were studied used specific-locus amplified fragment sequencing (SLAF-seq). The result of genetic diversity showed that there was no single genetic diversity center was observed across all E. sibiricus populations. The results of genetic variation showed that 44 populations were clearly divided into the following three groups: Qinghai Plateau (Group I), South Tibet (Group II), and Hengduan Mountains (Group III). From the three analyses of AMOVA, Mantel test and Treemix, strong genetic differentiation across all populations and low genetic differentiation among populations within three groups. Molecular dating indicated that E. sibiricus diverged at 16.08 Ma (during the early Miocene) can be linked to the Himalayan Motion stage of QTP uplift. It is speculated that the reasons affecting the current phylogeographical pattern are as follows: (1) The environmental changes due to the uplift of the QTP; (2) The geographic distance between the populations (Groups I and III are close in geographic distance, and gene flow are frequent); (3) Geographical barriers (the Tanggula and Bayangela Mountains between Groups I and II). This study provides new evidence and historical perspective to the future exploration of the evolution and geographic distribution pattern of Elymus L.

Study on the morpho-anatomy of Lepisorus species through light microscopy and scanning electron microscopy and its systematic implications

The genus Lepisorus is often associated with complex lineages in Polypodiaceae, which leads to difficulty in understanding taxonomic relationships among the species. The taxa is mostly epiphytic with some distinct features like rhizomes covered with clathrate scales, simple fronds, sori in a single row, intermixed with peltate paraphyses. In this study, we precisely focus on the Lepisorus species and seek a broader understanding of the taxonomic relationships that prevail in this genus. The morpho-anatomical traits of rhizomes, stipe, lamina, and sporangia were explored in nine species with the aid of light microscopy and scanning electron microscopy techniques for gathering the qualitative and quantitative data sets. These character traits were resolved into groups based on degree of similarity and principal component analysis to comprehend the covariance among the variables. The data were then employed to prepare an artificial dichotomous key and hierarchical cluster analysis was established that revealed five related clades with eight members whereas one member emerged as outgroup. The results correspond with other recent phylogenetic studies of the concerned genus hence confirming immense reliability and thrust of light microscopy, scanning electron microscopy, and morphology-based studies that are being less utilized in fern taxonomy.

Reticulate evolution in the Pteris fauriei group (Pteridaceae)

The Pteris fauriei group (Pteridaceae) has a wide distribution in Eastern Asia and includes 18 species with similar but varied morphology. We collected more than 300 specimens of the P. fauriei group and determined ploidy by flow cytometry and inferred phylogenies by molecular analyses of chloroplast and nuclear DNA markers. Our results reveal a complicated reticulate evolution, consisting of seven parental taxa and 58 hybrids. The large number of hybrid taxa have added significant morphological complexity to the group leading to difficult taxonomic issues. The hybrids generally had broader ranges and more populations than their parental taxa. Genetic combination of different pairs of parental species created divergent phenotypes of hybrids, exhibited by both morphological characteristics and ecological fidelities. Niche novelty could facilitate hybrid speciation. Apogamy is common in this group and potentially contributes to the sustainability of the whole group. We propose that frequent hybridizations among members of the P. fauriei group generate and maintain genetic diversity, via novel genetic combinations, niche differentiation, and apogamy.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Comparative Analyses of 35 Complete Chloroplast Genomes from the Genus Dalbergia (Fabaceae) and the Identification of DNA Barcodes for Tracking Illegal Logging and Counterfeit Rosewood

The genus Dalbergia contains more than 200 species, several of which are trees that produce traditional medicines and extremely high-value timber commonly referred to as rosewood. Due to the rarity of these species in the wild, the high value of the timber, and a growing international illicit trade, CITES (Convention on International Trade in Endangered Species of Wild Fauna and Flora) has listed the entire genus in appendix II and the species Dalbergia nigra in appendix I because species in this genus are considered at risk of extinction. Given this, and the fact that species or even genus level determination is nearly impossible from cut timber morphology, alternative molecular methods are needed to identify and track intercepted rosewood. To better identify rosewood using molecular methods, we sequenced and assembled eight chloroplast genomes including D. nigra as well as conducted comparative analyses with all other available chloroplast genomes in Dalbergia and closely related lineages. From these analyses, numerous repeats including simple sequence repeats (SSR) and conserved nucleotide polymorphisms unique to subclades within the genus were detected. From phylogenetic analysis based on the CDS from 77 chloroplast genes, the groups Siam rosewood and scented rosewood resolved as monophyletic, supporting the morphological traits used to delimit these species. In addition, several instances of paraphyly and polyphyly resulting from mismatches between taxonomic determinations and phylogenetic tree topology were identified. Ultimately, the highly variable regions in the chloroplast genomes will provide useful plastid markers for further studies regarding the identification, phylogeny, and population genetics of Dalbergia species, including those frequently intercepted in illegal trade.

Phylogeographical Analyses of a Relict Fern of Palaeotropical Flora (Vandenboschia speciosa): Distribution and Diversity Model in Relation to the Geological and Climate Events of the Late Miocene and Early Pliocene

Fern phylogeographic studies have mostly focused on the influence of the Pleistocene climate on fern distributions and the prevalence of long-distance dispersal. The effect of pre-Pleistocene events on the distributions of fern species is largely unexplored. Here, we elucidate a hypothetical scenario for the evolutionary history of Vandenboschia speciosa, hypothesised to be of Tertiary palaeotropical flora with a peculiar perennial gametophyte. We sequenced 40 populations across the species range in one plastid region and two variants of the nuclear gapCp gene and conducted time-calibrated phylogenetic, phylogeographical, and species distribution modelling analyses. Vandenboschia speciosa is an allopolyploid and had a Tertiary origin. Late Miocene aridification possibly caused the long persistence in independent refugia on the Eurosiberian Atlantic and Mediterranean coasts, with the independent evolution of gene pools resulting in two evolutionary units. The Cantabrian Cornice, a major refugium, could also be a secondary contact zone during Quaternary glacial cycles. Central European populations resulted from multiple post-glacial, long-distance dispersals. Vandenboschia speciosa reached Macaronesia during the Pliocene–Pleistocene, with a phylogeographical link between the Canary Islands, Madeira, and southern Iberia, and between the Azores and northwestern Europe. Our results support the idea that the geological and climate events of the Late Miocene/Early Pliocene shifted Tertiary fern distribution patterns in Europe.

Macaronesia Acts as a Museum of Genetic Diversity of Relict Ferns: The Case of Diplazium caudatum (Athyriaceae)

Macaronesia has been considered a refuge region of the formerly widespread subtropical lauroid flora that lived in Southern Europe during the Tertiary. The study of relict angiosperms has shown that Macaronesian relict taxa preserve genetic variation and revealed general patterns of colonization and dispersal. However, information on the conservation of genetic diversity and range dynamics rapidly diminishes when referring to pteridophytes, despite their dominance of the herbaceous stratum in the European tropical palaeoflora. Here we aim to elucidate the pattern of genetic diversity and phylogeography of Diplazium caudatum, a hypothesized species of the Tertiary Palaeotropical flora and currently with its populations restricted across Macaronesia and disjunctly in the Sierras de Algeciras (Andalusia, southern Iberian Peninsula). We analysed 12 populations across the species range using eight microsatellite loci, sequences of a region of plastid DNA, and carry out species-distribution modelling analyses. Our dating results confirm the Tertiary origin of this species. The Macaronesian archipelagos served as a refuge during at least the Quaternary glacial cycles, where populations of D. caudatum preserved higher levels of genetic variation than mainland populations. Our data suggest the disappearance of the species in the continent and the subsequent recolonization from Macaronesia. The results of the AMOVA analysis and the indices of clonal diversity and linkage disequilibrium suggest that D. caudatum is a species in which inter-gametophytic outcrossing predominates, and that in the Andalusian populations there was a shift in mating system toward increased inbreeding and/or clonality. The model that best explains the genetic diversity distribution pattern observed in Macaronesia is, the initial and recurrent colonization between islands and archipelagos and the relatively recent diversification of restricted area lineages, probably due to the decrease of favorable habitats and competition with lineages previously established. This study extends to ferns the concept of Macaronesia archipelagos as refugia for genetic variation.

Spatial phylogenetics of two topographic extremes of the Hengduan Mountains in southwestern China and its implications for biodiversity conservation

Previous attempts to elucidate the drivers of speciation mechanisms and spatial distribution patterns of biodiversity in mountain regions have treated different floras within a single geological region as one flora, ignoring the potential contributions of high habitat/ecosystem heterogeneity. Furthermore, current conservation strategies largely focus on forest ecosystems and/or specific flagship species, ignoring marginal ecosystems, leaving species in these ecosystems at risk. Here, we compared the spatial patterns of biodiversity and the potential drivers of these patterns in the river valley and subnival ecosystems of the Hengduan Mountains region (HDM) in southwestern China. Specifically, we compared spatial patterns of diversity, endemism, and threatened species in these ecosystems based on both traditional measurements and recent phylogenetic approaches. We then examined how those patterns were related to environmental factors and human activity in these same regions. We found that the middle-southern HDM supports the highest diversity and endemism for the river valley and subnival ecosystems; however, the distribution patterns of neo- and paleo-endemism in these two ecosystems differ. Regression models indicate that habitat diversity and paleo-climatic fluctuation are important drivers of diversity and endemism for these two ecosystems. Temperature and precipitation, however, showed different influences on the spatial patterns in different ecosystems. Categorical analysis of neo- and paleo-endemism (CANAPE) indicated that most endemism centers are not covered by current nature reserves. Moreover, the intensity of human activity is highest in the southern and southeastern HDM, which coincides with the distribution patterns of diversity, mixed-endemism and high-priority (and threatened) species. These findings suggest that different floras within a single geographic/floristic region respond differently to environmental factors and show different spatial phylogenetic patterns. We, therefore, recommend that future research into the drivers of biodiversity consider the contributions of various ecosystem types within a single geological region. This study also provides a theoretical basis for protecting habitat diversity. Our work confirms that current conservation efforts are insufficient to protect ecosystem diversity in the river valley and subnival ecosystems of the Hengduan Mountains. Therefore, we recommend the establishment of nature reserves in the regions identified in this study; furthermore, we strongly recommend improving current and establishing new management policies for biodiversity conservation in this region.

Chloroplast genomes in Populus (Salicaceae): comparisons from an intensively sampled genus reveal dynamic patterns of evolution

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

A global phylogeny of Stegnogramma ferns (Thelypteridaceae): generic and sectional revision, historical biogeography and evolution of leaf architecture

The thelypteroid fern genus Stegnogramma s.l. contains around 18-35 species and has a global, cross-continental distribution ranging from tropical to temperate regions. Several genera and infrageneric sections have been recognized previously in Stegnogramma s.l., but their phylogenetic relationships are still unclear. In this study, we present a global phylogeny of Stegnogramma s.l. with the most comprehensive sampling to date and aim to pinpoint the phylogenetic positions of biogeographically and taxonomically important taxa. Based on the reconstructed historical biogeography and character evolution, we propose a new (infra)generic classification and discuss the diversification of Stegnogramma s.l. in a biogeographical context. New names or combinations are made for 12 (infra)species, including transferring the monotypic species of Craspedosorus to Leptogramma. Finally, we discuss a possible link between leaf architecture and ecological adaptation, and hypothesize that the increase in leaf dissection and free-vein proportion is an adaptive feature to cool climates in Stegnogramma s.l.

The rise of the Andes promoted rapid diversification in Neotropical Phlegmariurus (Lycopodiaceae)

Tropical mountains are disproportionately biodiverse relative to their surface area, but the processes underlying their exceptional diversity require further study. Here, we use comparative phylogenetic methods to examine the impact of the Andean orogeny on the diversification of Neotropical Phlegmariurus, a species-rich lycophyte clade. We generated a time-calibrated phylogeny of 105 species of Neotropical Phlegmariurus and estimated lineage diversification rates. We tested for correlations between lineage diversification rates and species range size, niche breadth, elevational range amplitude, and mean elevation of occurrence. A recently developed macroevolutionary model was used to incorporate geological data and test for an association between diversification rates and the Andean uplift. Diversification rates of Neotropical Phlegmariurus are negatively correlated with species range size and positively correlated with mean elevation of species occurrence. The rise of the Andes is strongly associated with increased rates of diversification in Neotropical Phlegmariurus during the last 10 Myr. Our study demonstrates the importance of mountain-building events and geographical isolation of alpine populations as drivers of rapid diversification, even in spore-dispersed plants. This work also highlights the usefulness of combined phylogenetic, geological and ecological datasets, and the promise of comparative environment-dependent diversification models in better understanding the evolutionary origins of biodiversity.

The Complete Chloroplast Genome of Heimia myrtifolia and Comparative Analysis within Myrtales

Heimia myrtifolia is an important medicinal plant with several pharmacologically active alkaloids and is also used as an ornamental landscape plant. The purpose of this study is to complete and characterize the chloroplast (cp) genome of H. myrtifolia and compare genomic features to other Myrtales species' cp genomes. The analysis showed that H. myrtifolia has a total length of 159,219 bp with a typical quadripartite structure containing two identical inverted repeats (IRs) of 25,643 bp isolated by one large single copy (LSC) of 88,571 bp and one small single copy (SSC) of 18,822 bp. The H. myrtifolia cp genome contains 129 genes with eight ribosomal RNAs, 30 transfer RNAs, and 78 protein coding genes, in which 17 genes are duplicated in two IR regions. The genome organization including gene type and number and guanine-cytosine (GC) content is analyzed among the 12 cp genomes in this study. Approximately 255 simple sequence repeats (SSRs) and 16 forward, two reverses, and two palindromic repeats were identified in the H. myrtifolia cp genome. By comparing the whole H. myrtifolia cp genome with 11 other Myrtales species, the results showed that the sequence similarity was high between coding regions while sequence divergence was high between intergenic regions. By employing the full cp genomes for phylogenetic analysis, structural and sequence differences were characterized between H. myrtifolia and 11 Myrtales species illustrating what patterns are common in the evolution of cp genomes within the Myrtales. The first entire cp genome in the genus Heimia provides a valuable resource for further studies in these medicinally and ornamentally important taxa.

Population genetic dynamics of Himalayan-Hengduan tree peonies, Paeonia subsect. Delavayanae

According to the present taxonomical treatment, Paeonia subsect. Delavayanae consists of only two species (P. delavayi and P. ludlowii) endemic to the Himalayan-Hengduan Mountains. Although P. ludlowii can be distinguished from P. delavayi on the basis of a series of morphological characters, the species delimitation remains controversial because the more widespread one, P. delavayi, exhibits considerable morphological diversity. Both chloroplast DNA markers and nuclear microsatellites or simple sequence repeats (nSSR) are used herein to reveal genetic diversity and relationships of the two taxa included in this subsection, and ecological niche modeling (ENM) is employed to get insights into their paleodistribution. Our results show that genetic boundaries between the two currently recognized species are unclear, probably due to recent divergence. Paeonia ludlowii is budding from P. delavayi, probably by genetic isolation but also by shifting its niche to the harsher upland Tibetan conditions. Paeonia delavayi itself would be, however, under active speciation, showing significant genetic differentiation and morphological diversity. Whereas P. ludlowii would have endured the Pleistocene glacial periods by in situ persistence in local, small refugia, a 'dual' model seems to apply for P. delavayi (in situ persistence and retreat to refugia). The rarity of P. ludlowii and high evolutionary potential of P. delavayi imply high priority for in situ conservation of both taxa. The Himalayan-Hengduan Mountains are an ideal arena for differentiation within subsect. Delavayanae of Paeonia, by means of expansions/contractions/displacements, vertical migrations, and local survival/extinctions in response to the Neogene climate fluctuations and geological changes.

Characterization of the whole chloroplast genome of Chikusichloa mutica and its comparison with other rice tribe (Oryzeae) species

Chloroplast genomes are a significant genomic resource in plant species and have been used in many research areas. The complete genomic information from wild crop species could supply a valuable genetic reservoir for breeding. Chikusichloa mutica is one of the most important wild distant relatives of cultivated rice. In this study, we sequenced and characterized its complete chloroplast (cp) genome and compared it with other species in the same tribe. The whole cp genome sequence is 136,603 bp in size and exhibits a typical quadripartite structure with large and small single-copy regions (LSC, 82,327 bp; SSC, 12,598 bp) separated by a pair of 20,839-bp inverted repeats (IRA, B). A total of 110 unique genes are annotated, including 76 protein-coding genes, 4 ribosomal RNA genes and 30 tRNA genes. The genome structure, gene order, GC content, and other features are similar to those of other angiosperm cp genomes. When comparing the cp genomes between Oryzinae and Zizaniinae subtribes, the main differences were found between the junction regions and distribution of simple sequence repeats (SSRs). In comparing the two Chikusichloa species, the genomes were only 40 bp different in length and 108 polymorphic sites, including 83 single nucleotide substitutions (SNPs) and 25 insertion-deletions (Indels), were found between the whole cp genomes. The complete cp genome of C. mutica will be an important genetic tool for future breeding programs and understanding the evolution of wild rice relatives.

Studies on ontogeny and reproductive behaviour of Lepisorus nudus (Hook.) Ching (Polypodiaceae)

In-vitro studies of the ontogeny and mating system of the gametophytes of Lepisorus nudus were carried out through multispore and isolate cultures lasting 23 weeks. Spore germination begins early, on day 5-6. Spore germination pattern was Vittaria type and the germination percentage reached 82.69% (± 3.20%). Filamentous gametophyte did not branch and never produce separate prothalli. Occasionally the branching and separate prothalli were produced from mature and cordate gametophytes. Prothallial development was Drynaria type (cordate gametophytes with notched apex) contrary to other known species of Lepisorus, where gametophyte development was Kaulinia type (strap gametophytes without apical notch). Gametophyte production in multispore cultures reached up-to 75.6% (± 18.85%). All isolates initially produced archegonia and antheridia only after a prolonged cessation of production in archegonia. In contrast, only 37.2% (±12.63%) of individuals in multispore culture exhibited the same pattern with 29.8% (±7.56%) developing as males that did not produce archegonia by the end of the study. Only 37.2% (±12.63%) of archegoniate gametophytes developed antheridia by the end of the study and only once archegonia had degenerated; i.e., a temporal gap existed in expression of female and male gametangia. In multispore culture, only 26.21% (±5.70%) sporophytes developed on 160th day by fusion of female and male gametes that were derived from matings between sib gametophytes. In contrast, isolated gametophytes did not produce sporophytes. In isolate gametophytes, mature archegonia could not take delivery of male gametangia because antheridia were produced sequentially. This study suggests that the sequential expression of gametangia and absence of the intragametophytic selfing may also be a possible cause of reproductive barriers. Lepisorus nudus promotes inter-gametophytic selfing as an adaptive mechanism for reproductive success in multispore culture. This study presents a detailed account on reproductive biology of the taxa whose population is decreasing at distressing rate.

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai-Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.

Long-Distance Dispersal after the Last Glacial Maximum (LGM) Led to the Disjunctive Distribution of Pedicularis kansuensis (Orobanchaceae) between the Qinghai-Tibetan Plateau and Tianshan Region

Quaternary climate fluctuations have profoundly affected the current distribution patterns and genetic structures of many plant and animal species in the Qinghai-Tibetan Plateau (QTP) and adjacent mountain ranges, e.g. Tianshan (TSR), Altay, etc. In this greater area disjunct distributions are prominent but have nevertheless received little attention with respect to the historical processes involved. Here, we focus on Pedicularis kansuensis to test whether the current QTP and TSR disjunction is the result of a recent Holocene range expansion involving dispersal across arid land bridge(s) or a Pleistocene range fragmentation involving persistence in refugia. Two chloroplast DNA spacers were sequenced for 319 individuals from 34 populations covering the entire distribution range of this species in China. We found a total of 17 haplotypes of which all occurred in the QTP, and only five in the TSR. Overall genetic diversity was high (H_T = 0.882, H_S = 0.559) and higher in the QTP than in the TSR. Genetic differentiation among regions and populations was relatively low (G_ST = 0.366) and little evidence for a phylogeographic pattern emerged. The divergence times for the four main lineages could be dated to the early Pleistocene. Surprisingly, the two ubiquitous haplotypes diverged just before or around the Last Glacial Maximum (LGM) and were found in different phylogenetic lineages. The Species Distribution Model suggested a disappearance of P. kansuensis from the TSR during the LGM in contrast to a relatively constant potential distribution in the QTP. We conclude that P. kansuensis colonized the TSR after the LGM. The improbable long-distance dispersal by wind or water across arid land seed flow may well have had birds or men as vector.

Isolation and characterization of microsatellite markers in the Lepisorus clathratus complex (Polypodiaceae)

PREMISE OF THE STUDY

Microsatellites were designed and characterized in the Sino-Himalayan fern Lepisorus clathratus complex (Polypodiaceae) to further study the phylogeography and reproductive ecology of this species.

METHODS AND RESULTS

From a genomic library obtained by next-generation sequencing, 10 polymorphic and six monomorphic microsatellite loci were developed. In one population of L. clathratus from Taibaishan in central China, the number of alleles observed for these microsatellites ranged from seven to 29, and observed and expected heterozygosity ranged from 0.463 to 0.919 and from 0.797 to 0.947, respectively. Cross-amplification in other taxa within this complex was successful, but cross-amplification was poor for other congeneric species.

CONCLUSIONS

This set of newly developed microsatellite markers will be useful for assessing genetic diversity, population structure, and mating system, and to infer polyploid origin in the L. clathratus complex.

The complete chloroplast genome of queen's crape-myrtle(Lagerstroemia macrocarpa)

The whole complete chloroplast genome of Lagerstroemia macrocarpa was assembled in this study. Total genome is 152,472 bp in length consisting of two inverted repeats of 17,562 bp separated by a large single-copy region and a small single-copy region of 84,050 bp and 33,295 bp, respectively. This genome contains 112 unique genes including 78 protein-coding genes, 4 ribosomal RNA genes and 30 transfer RNA genes. In 78 protein-coding genes, 8 genes (atpF, ndhA, ndhB, petB, petD, rpl16, rpoC1, rps16) contain one intron and three genes with two introns each (clpP, rps12 and ycf3). This newly sequenced chloroplast genome supply highly variable information of polymorphisms within Lagerstroemia species.

The complete chloroplast genome of Arabidopsis lyrata

We report the complete chloroplast DNA (cpDNA) of Arabidopsis lyrata (Brassicaceae), a less studied relative of A. thaliana, by employing next-generation sequencing reads and de novo assembly. The length of the closed circular cpDNA is 154,604 bp with a typical quadripartite structure. The genome is composed of one large single copy and one small single copy regions of 84,209 bp and 17,871 bp, respectively, and separated by a pair of inverted repeats of 26,262 bp in length. The overall GC content is 36.35% and the GC content of the LSC, IRs and SSC regions are 34.12%, 42.30% and 29.38%, separately. The gene content and the number for A. lyrata are the same as other published species in Brassicaceae with 112 annotated known unique genes including 78 protein-coding genes, 30 tRNA genes and four rRNA genes. The complete cpDNA of A. lyrata will provide valuable molecular resources for further phylogenetic and evolutionary analysis in the model Arabidopsis genus.

The complete plastid genome of the wild rice species Oryza brachyantha (Poaceae)

The whole plastid genome of wild rice (Oryza brachyantha) is characterized in this study. The genome is 134 604 bp in length and is arranged in a typical circular structure, including a pair of inverted repeats (IRs) of 20 832 bp in size separated by a large single-copy region (LSC) of 80 411 bp in length and a small single-copy region (SSC) of 12 529 bp in length. The overall GC content of the genome is 38.98%. One hundred and ten unique genes were annotated from the chloroplast genome, including 76 protein-coding genes, 4 ribosomal RNA genes and 30 tRNA genes. A total of 20 of these genes are duplicated in the IR regions, 13 genes contain 1 intron and 2 genes (rps12 and ycf3) have 2 introns.

Two complete chloroplast genomes of white campion (Silene latifolia) from male and female individuals

In this study, we assembled two individuals' complete chloroplast genomes with one male and one female from the dioecious plant species white campion (Silene latifolia). The two chloroplast genomes have an identical composition with each other as a circular molecule of 150 931 bp in length, with an overall GC content of 36.4%. The genomes consist of a pair of inverted repeats (IRs) of 25 503 bp, separated by a large single-copy (LSC) region and a small single-copy (SSC) region of 82 708 and 17 217 bp, respectively. The genomes contain 111 single copy genes, including 77 protein-coding genes, 4 ribosomal RNA genes and 30 transfer RNA genes. In these protein-coding genes, eight genes (rpl16, rpoC1, rps16, petD, petB, ndhB, ndhA and atpF) contain a single intron and three genes (rps12, clpP and ycf3) contain two introns. The two newly sequenced chloroplast genomes provide valuable information for detecting polymorphisms within species and between sexes.

Limited variation across two chloroplast genomes with finishing chloroplast genome of Capsella grandiflora

The complete chloroplast genome of Capsella grandiflora is finished in this study, which consists of 154 638 base pairs (bp) in size containing a pair of inverted repeats (IRa and IRb) of 26 462 bp each and a small and large single-copies (SSC and LSC) of 17 835 and 83 879 bp, respectively. The overall GC content is 36.54% and the GC contents of LSC, IRs, and SSC are 34.34%, 42.38%, and 29.61% separately. The gene contents and numbers are the same with other two published species in genus Capsella with 112 annotated unique genes including 78 protein-coding genes, 30 tRNA genes, and four rRNA genes. Among these, 16 are duplicated in the inverted repeat regions, 15 genes contained one intron, and three genes (rps12, clpP, and ycf3) comprising two introns. Based on the whole genome comparison, only 28 SNPs and 16 Indels (insertion and deletion) are detected between two closed relatives' chloroplast genomes.

Phylogeography of Thlaspi arvense (Brassicaceae) in China Inferred from Chloroplast and Nuclear DNA Sequences and Ecological Niche Modeling

Thlaspi arvense is a well-known annual farmland weed with worldwide distribution, which can be found from sea level to above 4000 m high on the Qinghai-Tibetan Plateau (QTP). In this paper, a phylogeographic history of T. arvense including 19 populations from China was inferred by using three chloroplast (cp) DNA segments (trnL-trnF, rpl32-trnL and rps16) and one nuclear (n) DNA segment (Fe-regulated transporter-like protein, ZIP). A total of 11 chloroplast haplotypes and six nuclear alleles were identified, and haplotypes unique to the QTP were recognized (C4, C5, C7 and N4). On the basis of molecular dating, haplotypes C4, C5 and C7 have separated from others around 1.58 Ma for cpDNA, which corresponds to the QTP uplift. In addition, this article suggests that the T. arvense populations in China are a mixture of diverged subpopulations as inferred by hT/vT test (hT ≤ vT, cpDNA) and positive Tajima's D values (1.87, 0.05 < p < 0.10 for cpDNA and 3.37, p < 0.01 for nDNA). Multimodality mismatch distribution curves and a relatively large shared area of suitable environmental conditions between the Last Glacial Maximum (LGM) as well as the present time recognized by MaxEnt software reject the sudden expansion population model.

The complete chloroplast genome of Capsella rubella

The whole nucleotide sequence of the chloroplast genome from Capsella rubella is determined in this study using short Illumina sequence data from public database. The circular double-stranded DNA, which consists of 154,601 base pairs (bp) in size, contains a pair of inverted repeats (IRa and IRb) of 26,462 bp each, which are separated by a small and large single-copies (SSC and LSC) of 17,855 and 83,822 bp, respectively. The overall GC content of the chloroplast genome is 36.54% and the GC contents of LSC, IRs and SSC are 34.33%, 42.38% and 29.61% separately. One hundred and twelve unique genes were annotated, including 78 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Among these, 16 are duplicated in the inverted repeat regions, 15 genes contained 1 intron, and 3 genes (rps12, clpP and ycf3) comprised of 2 introns. Two protein genes (rps19 and ycf1) span in the boundaries of LSC-IR and IR-SSC to produce two partial pseudogenes.

The whole chloroplast genome of shrub willows (Salix suchowensis)

The whole chloroplast genome of Salix suchowensis was determined in this study using next generation sequencing data. The total genome size was 155,214 bp in length, containing a pair of inverted repeats (IRs) of 27,459 bp, which were separated by large single copy (LSC) and small single copy (SSC) of 84,077 bp and 16,219 bp, respectively. The overall GC contents of the chloroplast genome were 36.73%. One hundred and ten unique genes were annotated, including 76 protein coding genes, 30 tRNA genes and 4 rRNA genes. Among these, 18 are duplicated in the inverted repeat regions, 14 genes contained 1 intron, and 3 genes (rps12, clpP and ycf3) comprised of 2 introns.

DNA barcodes and species distribution models evaluate threats of global climate changes to genetic diversity: a case study from Nanorana parkeri (Anura: Dicroglossidae)

Anthropogenic global climate changes are one of the greatest threats to biodiversity. Distribution modeling can predict the effects of climate changes and potentially their effects on genetic diversity. DNA barcoding quickly identifies patterns of genetic diversity. As a case study, we use DNA barcodes and distribution models to predict threats under climate changes in the frog Nanorana parkeri, which is endemic to the Qinghai-Tibetan Plateau. Barcoding identifies major lineages W and E. Lineage W has a single origin in a refugium and Lineage E derives from three refugia. All refugia locate in river valleys and each greatly contributes to the current level of intraspecific genetic diversity. Species distribution models suggest that global climate changes will greatly influence N. parkeri, especially in the level of genetic diversity, because two former refugia will fail to provide suitable habitat. Our pipeline provides a novel application of DNA barcoding and has important implications for the conservation of biodiversity in southern areas of the Qinghai-Tibetan Plateau.

Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) is the highest and one of the most extensive plateaus in the world. Phylogenetic, phylogeographic, and ecological studies support plant diversifications on the QTP through multiple mechanisms such as allopatric speciation via geographic isolation, climatic oscillations and divergences, pollinator-mediated isolation, diploid hybridization and introgression, and allopolyploidy. These mechanisms have driven spectacular radiations and/or species diversifications in various groups of plants such as Pedicularis L., Saussurea DC., Rhododendron L., Primula L., Meconopsis Vig., Rhodiola L., and many lineages of gymnosperms. Nevertheless, much work is needed toward understanding the evolutionary mechanisms of plant diversifications on the QTP. Well-sampled biogeographic analyses of the QTP plants in the broad framework of the Northern Hemisphere as well as the Southern Hemisphere are still relatively few and should be encouraged in the next decade. This paper reviews recent evidence from phylogenetic and biogeographic studies in plants, in the context of rapid radiations, mechanisms of species diversifications on the QTP, and the biogeographic significance of the QTP in the broader context of both the Northern and Southern Hemisphere biogeography. Integrative multidimensional analyses of phylogeny, morphological innovations, geography, ecology, development, species interactions and diversifications, and geology are needed and should shed insights into the patterns of evolutionary assembly and radiations in this fascinating region.

Present, past and future of the European rock fern Asplenium fontanum: combining distribution modelling and population genetics to study the effect of climate change on geographic range and genetic diversity

BACKGROUND AND AIMS

Climate change is expected to alter the geographic range of many plant species dramatically. Predicting this response will be critical to managing the conservation of plant resources and the effects of invasive species. The aim of this study was to predict the response of temperate homosporous ferns to climate change.

METHODS

Genetic diversity and changes in distribution range were inferred for the diploid rock fern Asplenium fontanum along a South-North transect, extending from its putative last glacial maximum (LGM) refugia in southern France towards southern Germany and eastern-central France. This study reconciles observations from distribution models and phylogeographic analyses derived from plastid and nuclear diversity.

KEY RESULTS

Genetic diversity distribution and niche modelling propose that genetic diversity accumulates in the LGM climate refugium in southern France with the formation of a diversity gradient reflecting a slow, post-LGM range expansion towards the current distribution range. Evidence supports the fern's preference for outcrossing, contradicting the expectation that homosporous ferns would populate new sites by single-spore colonization. Prediction of climate and distribution range change suggests that a dramatic loss of range and genetic diversity in this fern is possible. The observed migration is best described by the phalanx expansion model.

CONCLUSIONS

The results suggest that homosporous ferns reproducing preferentially by outcrossing accumulate genetic diversity primarily in LGM climate refugia and may be threatened if these areas disappear due to global climate change.

Chloroplast DNA phylogeography reveals repeated range expansion in a widespread aquatic herb Hippuris vulgaris in the Qinghai-Tibetan Plateau and adjacent areas

Background

The Qinghai-Tibetan Plateau (QTP) is one of the most extensive habitats for alpine plants in the world. Climatic oscillations during the Quaternary ice age had a dramatic effect on species ranges on the QTP and the adjacent areas. However, how the distribution ranges of aquatic plant species shifted on the QTP in response to Quaternary climatic changes remains almost unknown.

Methodology and Principal Findings

We studied the phylogeography and demographic history of the widespread aquatic herb Hippuris vulgaris from the QTP and adjacent areas. Our sampling included 385 individuals from 47 natural populations of H. vulgaris. Using sequences from four chloroplast DNA (cpDNA) non-coding regions, we distinguished eight different cpDNA haplotypes. From the cpDNA variation in H. vulgaris, we found a very high level of population differentiation (G_ST = 0.819) but the phylogeographical structure remained obscure (N_ST = 0.853>G_ST = 0.819, P>0.05). Phylogenetic analyses revealed two main cpDNA haplotype lineages. The split between these two haplotype groups can be dated back to the mid-to-late Pleistocene (ca. 0.480 Myr). Mismatch distribution analyses showed that each of these had experienced a recent range expansion. These two expansions (ca. 0.12 and 0.17 Myr) might have begun from the different refugees before the Last Glacial Maximum (LGM).

Conclusions/Significance

This study initiates a research on the phylogeography of aquatic herbs in the QTP and for the first time sheds light on the response of an alpine aquatic seed plant species in the QTP to Quaternary climate oscillations.

The rise of the Himalaya enforced the diversification of SE Asian ferns by altering the monsoon regimes

BACKGROUND

The rise of high mountain chains is widely seen as one of the factors driving rapid diversification of land plants and the formation of biodiversity hotspots. Supporting evidence was reported for the impact of the rapid rise of the Andean mountains but this hypothesis has so far been less explored for the impact of the "roof of the world". The formation of the Himalaya, and especially the rise of the Qinghai-Tibetan Plateau in the recent 20 million years, altered the monsoon regimes that dominate the current climates of South East Asia. Here, we infer the hypothesis that the rise of Himalaya had a strong impact on the plant diversity in the biodiversity hotspot of the Southwest Chinese Mountains.

RESULTS

Our analyses of the diversification pattern of the derived fern genus Lepisorus recovered evidence for changes in plant diversity that correlated with the strengthening of South East Asian monsoon. Southwest China or Southwest China and Japan was recovered as the putative area of origin of Lepisorus and enhancing monsoon regime were found to shape the early diversification of the genus as well as subsequent radiations during the late Miocene and Pliocene.

CONCLUSIONS

We report new evidence for a coincidence of plant diversification and changes of the climate caused by the uplift of the Himalaya. These results are discussed in the context of the impact of incomplete taxon sampling, uncertainty of divergence time estimates, and limitations of current methods used to assess diversification rates.

Population expanding with the phalanx model and lineages split by environmental heterogeneity: a case study of Primula obconica in subtropical China

BACKGROUND

Current and historical events have both affected the current distribution patterns and intraspecific divergence of plants. While numerous studies have focused on the Qinghai-Tibetan Plateau (QTP), the impacts of such events on the flora of subtropical China remain poorly understood. Subtropical China is famous for its highly complex topography and the limited impact from glaciation during the Pleistocene; this may have resulted in a different genetic legacy for species in this region compared to fully glaciated areas.

METHODOLOGY/PRINCIPAL FINDINGS

We used plastid and nuclear DNA sequence data and distribution modeling to analyze the divergence patterns and demographic history of Primula obconica Hance, a widespread herbaceous montane species in subtropical China. The phylogenetic analysis revealed two major lineages (lineage A and lineage B), representing a west-east split into the Yunnan and Eastern groups, and the Sichuan and Central groups, respectively. The Eastern and Central groups comprised relatively new derived haplotypes. Nested Clade Analysis and Bayesian Skyline Plot analyses both indicated that P. obconica mainly experienced a gradual expansion of populations. In addition, the simulated distribution of P. obconica during the Last Glacial Maximum was slightly larger than its present-day distribution.

CONCLUSION/SIGNIFICANCE

Our results are the first to identify a west-east migration of P. obconica. The gradual expansion pattern and a larger potential distribution range in cold periods detected for P. obconica indicate that the population expansion of this species is consistent with the phalanx model. In addition, the current patterns of genetic differentiation have persisted as a result of the extensive environmental heterogeneity that exists in subtropical China.

Indehiscent sporangia enable the accumulation of local fern diversity at the Qinghai-Tibetan Plateau

BACKGROUND

Indehiscent sporangia are reported for only a few of derived leptosporangiate ferns. Their evolution has been likely caused by conditions in which promotion of self-fertilization is an evolutionary advantageous strategy such as the colonization of isolated regions and responds to stressful habitat conditions. The Lepisorus clathratus complex provides the opportunity to test this hypothesis because these derived ferns include specimens with regular dehiscent and irregular indehiscent sporangia. The latter occurs preferably in well-defined regions in the Himalaya. Previous studies have shown evidence for multiple origins of indehiscent sporangia and the persistence of populations with indehiscent sporangia at extreme altitudinal ranges of the Qinghai-Tibetan Plateau (QTP).

RESULTS

Independent phylogenetic relationships reconstructed using DNA sequences of the uniparentally inherited chloroplast genome and two low-copy nuclear genes confirmed the hypothesis of multiple origins of indehiscent sporangia and the restriction of particular haplotypes to indehiscent sporangia populations in the Lhasa and Nyingchi regions of the QTP. In contrast, the Hengduan Mountains were characterized by high haplotype diversity and the occurrence of accessions with and without indehiscent sporangia. Evidence was found for polyploidy and reticulate evolution in this complex. The putative case of chloroplast capture in the Nyingchi populations provided further evidence for the promotion of isolated but persistent populations by indehiscent sporangia.

CONCLUSIONS

The presented results confirmed the hypothesis that indehiscent sporangia promote the establishment of persistent population in different regions of the QTP. These results are consistent with the expectations of reproductive reassurance by promotion of self-fertilization that played a critical role in the assembly of populations in isolated locations and/or extreme habitats.

Great genetic differentiation among populations of Meconopsis integrifolia and its implication for plant speciation in the Qinghai-Tibetan Plateau

The complex tectonic events and climatic oscillations in the Qinghai-Tibetan Plateau (QTP), the largest and highest plateau in the world, are thought to have had great effects on the evolutionary history of the native plants. Of great interest is to investigate plant population genetic divergence in the QTP and its correlation with the geologic and climatic changes. We conducted a range-wide phylogeographical analysis of M. integrifolia based on the chloroplast DNA (cpDNA) trnL-trnF and trnfM-trnS regions, and defined 26 haplotypes that were phylogenetically divided into six clades dated to the late Tertiary. The six clades correspond, respectively, to highly differentiated population groups that do not overlap in geographic distribution, implying that the mountain ranges acting as corridors or barriers greatly affected the evolutionary history of the QTP plants. The older clade of M. integrifolia only occurs in the southwest of the species' range, whereas the distributions of younger clades extend northeastward in the eastern QTP, suggesting that climatic divergence resulting from the uplift of the QTP triggered the initial divergence of M. integrifolia native to the plateau. Also, the nrDNA ITS region was used to clarify the unexpected phylogenetic relationships of cpDNA haplotypes between M. integrifolia and M. betonicifolia. The topological incongruence between the two phylogenies suggests an ancestral hybridization between the two species. Our study indicates that geographic isolation and hybridization are two important mechanisms responsible for the population differentiation and speciation of Meconopsis, a species-rich genus with complex polyploids.