Population structure of Nouelia insignis (Asteraceae), an endangered species in southwestern China, based on chloroplast DNA sequences: recent demographic shrinking

Characterization of the complete chloroplast genome sequence of Chinese endemic species of Nouelia insignis (Hyalideae, Asteraceae) and its phylogenetic implications

This study was the first report complete chloroplast genome of Nouelia insignis (Asteraceae, Hyalideae), the large shrubs to small trees endemic to China. The circular whole cp genome of N. insignis was 151,524 bp in length, containing a large single-copy (LSC) region of 83,145 bp and a small single-copy (SSC) region of 18,261 bp. These two regions were separated by a pair of inverted repeat regions (IRa and IRb), each of them 25,060 bp in length. A total of 135 functional genes were encoded, consisting of 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The overall GC content of the chloroplast genome sequence was 37.8%, and the GC contents of the LSC, SSC, and IR regions were 35.9, 31.5, and 43.2%, respectively. The phylogenetic analysis by the Bayesian analysis showed that the species of N.insignis was sister group with Gerbera jamesonii by strong support values, and thus was closely related to members of subfamilies of Cichorioideae and Pertyoideae. These results will be useful for the future studies of Asteraceae in the worldwide.

The impacts of flowering phenology on the reproductive success of the narrow endemic Nouelia insignis Franch. (Asteraceae)

Nouelia insignis Franch. (Asteraceae) is a short, narrow endemic and endangered tree, growing with a natural population in the dry and hot valley of the Jinsha River in the southwest area of China. In this work, flowering phenology (time and duration), floral biology, visit frequency and behavior of pollinators, and pollination characteristics were studied based on investigation in the field and analysis in the laboratory with the help of a stereomicroscope, and the relationship between seed setting rate and reproductive traits, as well as the relationship between flowering time and rainfall before flowering, was tested using the method of general linear regression model. The results showed that natural population of N. insignis exhibited high flowering synchrony with relatively stable flowering duration, and the flowering time fluctuated greatly depending on the rainfall 5 months before flowering. The pollination of N. insignis required pollinators, and insect activities played a very important role in the pollination process. However, lack of the pollinators was not a limitation for reproductive fitness in N. insignis, although the number of pollinators was small and the frequency of visits was low. In addition, no pollen limitation was found during pollination. The average seed setting rate of N. insignis in the natural condition was only 1.52%-3.73%, and it was generally affected by changes in flowering phenology between years and had a higher seed set in early flowering year. The annual variation of seed set might be related to the annual variations of stamen and pistil functions, such as changes of pollen viability and stigma receptivity, which were closely related to flowering time. The results of this study are of value for further conservation actions on natural population of this threatened endemic plant.

Genetic characterization of the entire range of Cycas panzhihuaensis (Cycadaceae)

Cycas panzhihuaensis L. Zhou & S. Y. Yang (Cycadaceae) is an endangered gymnosperm species endemic to the dry-hot valley of the Jinsha River basin in southwest China. Although the wild C. panzhihuaensis population from Panzhihua Cycad Natural Reserve is well protected and its genetic diversity has been well assessed, the genetic characteristics of populations outside the nature reserve, which face larger risks of extinction, remain unknown. Furthermore, the population genetics and historical dynamics of this endemic and endangered species have not been examined across its entire range. In this study, to analyze the genetic diversity, phylogeographical structure and demographic history of C. panzhihuaensis from all its seven known locations, we sequenced and compared molecular data from chloroplastic DNA (psbA-trnH, psbM-trnD, and trnS-trnG), single-copy nuclear genes (PHYP, AC5, HSP70, and AAT) from 61 individuals, as well as 11 nuclear microsatellite loci (SSR) from 102 individuals. We found relatively high genetic diversity within populations and high genetic differentiation among populations of C. panzhihuaensis, which is consistent with the patterns of other Asian inland cycads. Although no significant phylogeographical structure was detected, we found that small and unprotected populations possess higher genetic diversity and more unique haplotypes, which revises our understanding of diversity within this species and deserves due attention. Analysis of demographic dynamics suggest that human activity might be the key threat to C. panzhihuaensis. Based on the genetic characterization of C. panzhihuaensis, we propose several practical guidelines for the conservation of this species, especially for the populations with small sizes.

Ecological and Genetic Divergences with Gene Flow of Two Sister Species (Leucomeris decora and Nouelia insignis) Driving by Climatic Transition in Southwest China

Understanding of the processes of divergence and speciation is a major task for biodiversity researches and may offer clearer insight into mechanisms generating biological diversity. Here, we employ an integrative approach to explore genetic and ecological differentiation of Leucomeris decora and Nouelia insignis distributed allopatrically along the two sides of the biogeographic boundary 'Tanaka Line' in Southwest China. We addressed these questions using ten low-copy nuclear genes and nine plastid DNA regions sequenced among individuals sampled from 28 populations across their geographic ranges in China. Phylogenetic, coalescent-based population genetic analyses, approximate Bayesian computation (ABC) framework and ecological niche models (ENMs) were conducted. We identified a closer phylogenetic relationship in maternal lineage of L. decora with N. insignis than that between L. decora and congeneric Leucomeris spectabilis. A deep divergence between the two species was observed and occurred at the boundary between later Pliocene and early Pleistocene. However, the evidence of significant chloroplast DNA gene flow was also detected between the marginal populations of L. decora and N. insignis. Niche models and statistical analyses showed significant ecological differentiation, and two nuclear loci among the ten nuclear genes may be under divergent selection. These integrative results imply that the role of climatic shift from Pliocene to Pleistocene may be the prominent factor for the divergence of L. decora and N. insignis, and population expansion after divergence may have given rise to chloroplast DNA introgression. The divergence was maintained by differential selection despite in the face of gene flow.

Species delimitation, genetic diversity and population historical dynamics of Cycas diannanensis (Cycadaceae) occurring sympatrically in the Red River region of China

Delimitating species boundaries could be of critical importance when evaluating the species' evolving process and providing guidelines for conservation genetics. Here, species delimitation was carried out on three endemic and endangered Cycas species with resembling morphology and overlapped distribution range along the Red River (Yuanjiang) in China: Cycas diananensis Z. T. Guan et G. D. Tao, Cycas parvula S. L. Yang and Cycas multiovula D. Y. Wang. A total of 137 individuals from 15 populations were genotyped by using three chloroplastic (psbA-trnH, atpI-atpH, and trnL-rps4) and two single copy nuclear (RPB1 and SmHP) DNA sequences. Basing on the carefully morphological comparison and cladistic haplotype aggregation (CHA) analysis, we propose all the populations as one species, with the rest two incorporated into C. diannanensis. Genetic diversity and structure analysis of the conflated C. diannanensis revealed this species possessed a relative lower genetic diversity than estimates of other Cycas species. The higher genetic diversity among populations and relative lower genetic diversity within populations, as well as obvious genetic differentiation among populations inferred from chloroplastic DNA (cpDNA) suggested a recent genetic loss within this protected species. Additionally, a clear genetic structure of C. diannanensis corresponding with geography was detected based on cpDNA, dividing its population ranges into "Yuanjiang-Nanhun" basin and "Ejia-Jiepai" basin groups. Demographical history analyses based on combined cpDNA and one nuclear DNA (nDNA) SmHP both showed the population size of C. diannanensis began to decrease in Quaternary glaciation with no subsequent expansion, while another nDNA RPB1 revealed a more recent sudden expansion after long-term population size contraction, suggesting its probable bottleneck events in history. Our findings offer grounded views for clarifying species boundaries of C. diannanensis when determining the conservation objectives. For operational guidelines, the downstream populations which occupy high and peculiar haplotypes should be given prior in-situ conservation. In addition, ex-situ conservation and reintroduction measures for decades of generations are supplemented for improving the population size and genetic diversity of the endemic and endangered species.

Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the 'Tanaka Line' in Southwest China

BACKGROUND

Leucomeris decora and Nouelia insignis (Asteraceae) are narrowly and allopatrically distributed species, separated by the important biogeographic boundary Tanaka Line in Southwest China. Previous morphological, cytogenetic and molecular studies suggested that L. decora is sister to N. insignis. However, it is less clear how the two species diverged, whether in full isolation or occurring gene flow across the Tanaka Line. Here, we performed a molecular study at the population level to characterize genetic differentiation and decipher phylogeographic history in two closely related species based on variation examined in plastid and nuclear DNAs using a coalescent-based approach.

RESULTS

These morphologically distinct species share plastid DNA (cpDNA) haplotypes. In contrast, Bayesian analysis of nuclear DNA (nDNA) uncovered two distinct clusters corresponding to L. decora and N. insignis. Based on the IMa analysis, no strong indication of migration was detected based on both cpDNA and nDNA sequences. The molecular data pointed to a major west-east split in nuclear DNA between the two species corresponding with the Tanaka Line. The coalescent time estimate for all cpDNA haplotypes dated to the Mid-Late Pleistocene. The estimated demographic parameters showed that the population size of L. decora was similar to that of N. insignis and both experienced limited demographic fluctuations recently.

CONCLUSIONS

The study revealed comprehensive species divergence and phylogeographic histories of N. insignis and L. decora divided by the Tanaka Line. The phylogeographic pattern inferred from cpDNA reflected ancestrally shared polymorphisms without post-divergence gene flow between species. The marked genealogical lineage divergence in nDNA provided some indication of Tanaka Line for its role as a barrier to plant dispersal, and lent support to its importance in promoting strong population structure and allopatric divergence.

Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau

Habitat heterogeneity, physical barriers, and the uplift of the Yungui Plateau were found to deeply affect the phylogeographic pattern and evolutionary history of Cephalotaxus oliveri, a perennial conifer endemic to China. In this study, we explored the phylogeography using three chloroplast sequences (trnL-trnF, trnT-trnD and atpB-rbcL) in 22 natural populations of C. oliveri distributed throughout its range. The Yungui Plateau populations of C. oliveri were revealed to origin ca. 9.15Ma by molecular clock estimation, which is consistent with rapid uplift of the Qinghai-Tibetan Plateau (QTP) ca. 8-10Ma. Additionally, geological effects of the Yungui Plateau were suggested to promote the rapid intra-specific differentiation of C. oliveri in the Pliocene and Early Pleistocene. The relatively low level of genetic diversity (h=0.719, θ=1.17×10(-3)) and high population differentiation (NST=0.771 and GST=0.642) implied restricted gene flow among populations, which was confirmed by the Nested Clade Analysis (NCA). Mismatch distribution and haplotypes network provided evidences of recent demographic population expansion. Furthermore, the statistical dispersal-vicariance analysis indicated that the center of origin was in Central China. The comparison of haplotype distribution patterns indicated that the regions of HNHPS and HBLD were the potential refugia during the Pleistocene ice ages. Our results highlighted that habitat heterogeneity and physical barriers presenting in a species range can predict genetic patterns.