Impacts of biogeographic history and marginal population genetics on species range limits: a case study of Liriodendron chinense

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri

Evaluating the contributions of geographic distance and environmental heterogeneity to the genetic divergence can inform the demographic history and responses to environmental change of natural populations. The isolation-by-distance (IBD) reveals that genetic differentiation among populations increases with geographic distance, while the isolation-by-environment (IBE) assumes a linear relationship between genetic variation and environmental differences among populations. Here, we sampled and genotyped 330 individuals from 18 natural populations of Cephalotaxus oliveri throughout the species' distribution. Twenty-eight EST-SSR markers were applied to analyze population genetics, for the investigation of the driving factors that shaped spatial structure. In addition, we identified the outlier loci under positive selection and tested their association with environmental factors. The results showed a moderate genetic diversity in C. oliveri and high genetic differentiation among populations. Population structure analyses indicated that 18 populations were clustered into two major groups. We observed that the genetic diversity of central populations decreased and the genetic differentiation increased towards the marginal populations. Additionally, the signatures of IBD and IBE were detected in C. oliveri, and IBE provided a better contribution to genetic differentiation. Six outlier loci under positive selection were demonstrated to be closely correlated with environmental variables, among which bio8 was associated with the greatest number of loci. Genetic evidence suggests the consistency of the central-marginal hypothesis (CMH) for C. oliveri. Furthermore, our results suggest that temperature-related variables played an important role in shaping genetic differentiation.

Assessing Molecular Diversity in Native and Introduced Populations of Red Wood Ant Formica paralugubris

The Formica rufa group comprises several ant species which are collectively referred to as "red wood ants" and play key roles in boreal forest ecosystems, where they are ecologically dominant and greatly influence habitat dynamics. Owing to their intense predatory activity, some of these species are used as biocontrol agents against several forest insect pests and for this aim in Italy, nearly 6000 ant nests were introduced from their native areas in the Alps to several Appeninic sites during the last century. In this work, we assessed and compared the genetic variability and structure of native and introduced populations of F. paralugubris, thus evaluating the extent of genetic drift that may have occurred since the time of introduction, using amplified fragment length polymorphism (AFLP) markers. PCR amplification with a fam_EcoRI-TAC/MseI-ATG primers combination produced a total of 147 scorable bands, with 17 identified as outlier loci. The genetic variation was higher in the introduced population compared to the native ones that, on the other hand, showed a higher diversity between nests. AMOVA results clearly pointed out that the overall genetic structure was dominated by among-worker variation, considering all populations, the Alpine vs. Apennine groups and the comparison among native and related introduced populations (all ranging between 77.84% and 79.84%). Genetic analyses unveiled the existence of six main different groups that do not entirely mirror their geographic subdivision, pointing towards a wide admixture between populations, but, at the same time, rapid diversification of some Apennine populations. Future studies based on high-throughput genomic methods are needed to obtain a thorough understanding of the effects of environmental pressure on the genetic structure and mating system of these populations.

Genetic Diversity in Marginal Populations of Nitraria schoberi L. from Romania

Nitraria schoberi L. (Nitrariaceae) is a halophytic plant with a continuous range in Central Asia and with only two populations in the westernmost distribution limit of species, in Romania. Currently, there is no documented explanation for the species’ presence in Europe, outside the main distribution area. Considering that marginal populations genetics are important in establishing range limits and species adaptative potential, genetic diversity was assessed using Inter-simple sequence repeat markers (ISSR). Both the Shannon’s Information Index (I) and Expected Heterozygosity (He) suggested a relatively low level of genetic diversity within the two populations. However, the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) dendrogram and Principal Coordinates Analysis clearly distinguished the two populations. Our presumptions, based on current results, are that the marginal westernmost population of N. schoberi was established due to the unique conditions from the “islands of desert” developed in a temperate continental climate. The European establishment of this species was likely accidental and probably due to ornithochory. Genetic relatedness between populations could be a consequence of their common origin, presumably from proximal Asian N. schoberi populations, while the separation can be explained by the lack of genetic material exchange between the two populations.

Alterations in Population Distribution of Liriodendron chinense (Hemsl.) Sarg. and Liriodendron tulipifera Linn. Caused by Climate Change

Climate change has a significant impact on species population size and distribution, global biodiversity, and ecological status. The Liriodendron genus contains two species: Liriodendron chinense and Liriodendron tulipifera, both playing important roles in timber, medicinal, and landscape purposes. However, little is known about their population distribution characteristics and important climatic factors shaping their suitability. In this research, we used the geological record data, 19 climate components, MaxEnt, and ArcGIS to recreate and analyze the potential population distribution and their alterations of Liriodendron within the world beneath the current and future scenarios of RCP 2.6, RCP 4.5, and RCP 8.5 in 2050 and 2070. Our results showed that: Liriodendron is suitable to grow in subtropical monsoon climate areas, and that the climatic factor of precipitation of warmest quarter exerts the greatest impact on L. chinense, with a contribution rate of 57.6%. Additionally, we showed that the climatic factor of precipitation of the driest month exerts the greatest impact on L. tulipifera, with a contribution rate of 60.5%. Further analysis exhibited that low temperature and temperature fluctuations are major temperature factors affecting L. chinense and L. tulipifera, respectively. Therefore, we predicted that by the 2050s and 2070s, the areas of Liriodendron suitable habitats would increase first and then decrease in three scenarios; except the area of L. tulipifera suitable habitats under RCP8.5, which shows a slight increase. We then conclude that the Liriodendron suitable areas would shift to high latitudes due to global climate warming. The information gained from this study will provide a reference for developing forest cultivation, management, and conservation strategies for these two important tree species, and also a basis for subsequent biogeographic research.

Genome-wide identification of the Liriodendron chinense WRKY gene family and its diverse roles in response to multiple abiotic stress

BACKGROUND

Liriodendron chinense (Lchi) is a tree species within the Magnoliaceae family and is considered a basal angiosperm. The too low or high temperature or soil drought will restrict its growth as the adverse environmental conditions, thus improving L. chinense abiotic tolerance was the key issues to study. WRKYs are a major family of plant transcription factors known to often be involved in biotic and abiotic stress responses. So far, it is still largely unknown if and how the LchiWRKY gene family is tied to regulating L. chinense stress responses. Therefore, studying the involvement of the WRKY gene family in abiotic stress regulation in L. chinense could be very informative in showing how this tree deals with such stressful conditions.

RESULTS

In this research, we performed a genome-wide analysis of the Liriodendron chinense (Lchi) WRKY gene family, studying their classification relationships, gene structure, chromosomal locations, gene duplication, cis-element, and response to abiotic stress. The 44 members of the LchiWRKY gene family contain a significant amount of sequence diversity, with their lengths ranging from 525 bp to 40,981 bp. Using classification analysis, we divided the 44 LchiWRKY genes into three phylogenetic groups (I, II, II), with group II then being further divided into five subgroups (IIa, IIb, IIc, IId, IIe). Comparative phylogenetic analysis including the WRKY families from 17 plant species suggested that LchiWRKYs are closely related to the Magnolia Cinnamomum kanehirae WRKY family, and has fewer family members than higher plants. We found the LchiWRKYs to be evenly distributed across 15 chromosomes, with their duplication events suggesting that tandem duplication may have played a major role in LchiWRKY gene expansion model. A Ka/Ks analysis indicated that they mainly underwent purifying selection and distributed in the group IId. Motif analysis showed that LchiWRKYs contained 20 motifs, and different phylogenetic groups contained conserved motif. Gene ontology (GO) analysis showed that LchiWRKYs were mainly enriched in two categories, i.e., biological process and molecular function. Two group IIc members (LchiWRKY10 and LchiWRKY37) contain unique WRKY element sequence variants (WRKYGKK and WRKYGKS). Gene structure analysis showed that most LchiWRKYs possess 3 exons and two different types of introns: the R- and V-type which are both contained within the WRKY domain (WD). Additional promoter cis-element analysis indicated that 12 cis-elements that play different functions in environmental adaptability occur across all LchiWRKY groups. Heat, cold, and drought stress mainly induced the expression of group II and I LchiWRKYs, some of which had undergone gene duplication during evolution, and more than half of which had three exons. LchiWRKY33 mainly responded to cold stress and LchiWRKY25 mainly responded to heat stress, and LchiWRKY18 mainly responded to drought stress, which was almost 4-fold highly expressed, while 5 LchiWRKYs (LchiWRKY5, LchiWRKY23, LchiWRKY14, LchiWRKY27, and LchiWRKY36) responded equally three stresses with more than 6-fold expression. Subcellular localization analysis showed that all LchiWRKYs were localized in the nucleus, and subcellular localization experiments of LchiWRKY18 and 36 also showed that these two transcription factors were expressed in the nucleus.

CONCLUSIONS

This study shows that in Liriodendron chinense, several WRKY genes like LchiWRKY33, LchiWRKY25, and LchiWRKY18, respond to cold or heat or drought stress, suggesting that they may indeed play a role in regulating the tree's response to such conditions. This information will prove a pivotal role in directing further studies on the function of the LchiWRKY gene family in abiotic stress response and provides a theoretical basis for popularizing afforestation in different regions of China.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

Population dynamics linked to glacial cycles in Cercis chuniana F. P. Metcalf (Fabaceae) endemic to the montane regions of subtropical China

The mountains of subtropical China are an excellent system for investigating the processes driving the geographical distribution of biodiversity and radiation of plant populations in response to Pleistocene climate fluctuations. How the major mountain ranges in subtropical China have affected the evolution of plant species in the subtropical evergreen broadleaved forest is an issue with long-term concern. Here, we focused on Cercis chuniana, a woody species endemic to the southern mountain ranges in subtropical China, to elucidate its population dynamics. We used genotyping by sequencing (GBS) to investigate the spatial pattern of genetic variation among 11 populations. Geographical isolation was detected between the populations located in adjacent mountain ranges, thought to function as geographical barriers due to their complex physiography. Bayesian time estimation revealed that population divergence occurred in the middle Pleistocene, when populations in the Nanling Mts. separated from those to the east. The orientation and physiography of the mountain ranges of subtropical China appear to have contributed to the geographical pattern of genetic variation between the eastern and western populations of C. chuniana. Complex physiography plus long-term stable ecological conditions across glacial cycles facilitated the demographic expansion in the Nanling Mts., from which contemporary migration began. The Nanling Mts. are thus considered as a suitable area for preserving population diversity and large population sizes of C. chuniana compared with other regions. As inferred by ecological niche modeling and coalescent simulations, secondary contact occurred during the warm Lushan-Tali Interglacial period, with intensified East Asia summer monsoon and continuous habitat available for occupation. Our data support the strong influence of both climatic history and topographic characteristics on the high regional phytodiversity of the subtropical evergreen broadleaved forest in subtropical China.

The different fates of two Asian horseshoe crab species with different dispersal abilities

Impending anthropogenic climate change will severely impact coastal organisms at unprecedented speed. Knowledge on organisms' evolutionary responses to past sea-level fluctuations and estimation of their evolutionary potential is therefore indispensable in efforts to mitigate the effects of future climate change. We sampled tens of thousands of genomic markers of ~300 individuals in two of the four extant horseshoe crab species across the complex archipelagic Singapore Straits. Carcinoscorpius rotundicauda Latreille, a less mobile mangrove species, has finer population structure and lower genetic diversity compared with the dispersive deep-sea Tachypleus gigas Müller. Even though the source populations of both species during the last glacial maximum exhibited comparable effective population sizes, the less dispersive C. rotundicauda seems to lose genetic diversity much more quickly because of population fragmentation. Contra previous studies' results, we predict that the more commonly sighted C. rotundicauda faces a more uncertain conservation plight, with a continuing loss in evolutionary potential and higher vulnerability to future climate change. Our study provides important genomic baseline data for the redirection of conservation measures in the face of climate change and can be used as a blueprint for assessment and mitigation of the adverse effects of impending sea-level rise in other systems.

Characterization of the Liriodendron Chinense MYB Gene Family and Its Role in Abiotic Stress Response

Liriodendron chinense (Lchi) is a Magnoliaceae plant, which is a basic angiosperm left behind by the Pleistocene and mainly distributed in the south of the Yangtze River. Liriodendron hybrids has good wood properties and is widely used in furniture and in other fields. It is not clear if they can adapt to different environmental conditions, such as drought and high and low temperatures, and the molecular mechanisms for this adaptation are unknown. Among plant transcription factors (TFs), the MYB gene family is one of the largest and is often involved in stress or adversity response signaling, growth, and development. Therefore, studying the role of MYBTFs in regulating abiotic stress signaling, growth, and development in Lchi is helpful to promote afforestation in different environments. In our research, a genome-wide analysis of the LchiMYB gene family was performed, including the phylogenetic relationship tree, gene exon-intron structure, collinearity, and chromosomal position. According to the evolutionary tree, 190 LchiMYBs were divided into three main branches. LchiMYBs were evenly distributed across 19 chromosomes, with their collinearity, suggesting that segment duplication events may have contributed to LchiMYB gene expansion. Transcriptomes from eight tissues, 11 stages of somatic embryogenesis, and leaves after cold, heat, and drought stress were used to analyze the function of the MYB gene family. The results of tissue expression analysis showed that most LchiMYB genes regulated bark, leaf, bud, sepal, stigma, and stamen development, as well as the four important stages (ES3, ES4, ES9, and PL) of somatic embryogenesis. More than 60 LchiMYBs responded to heat, cold, and drought stress; some of which underwent gene duplication during evolution. LchiMYB3 was highly expressed under all three forms of stress, while LchiMYB121 was strongly induced by both cold and heat stress. Eight genes with different expression patterns were selected and verified by quantitative real-time PCR (qRT-PCR) experiments. The results suggested that these LchiMYBs may regulate Lchi growth development and resistance to abiotic stress. This study shows the cross-regulatory function of LchiMYBs in the growth and development, asexual reproduction, and abiotic resistance of Lchi. This information will prove pivotal to directing further studies on the biological function of Lchi MYBTFs in genetic improvement and abiotic stress response.

Population Genetics of Odontarrhena (Brassicaceae) from Albania: The Effects of Anthropic Habitat Disturbance, Soil, and Altitude on a Ni-Hyperaccumulator Plant Group from a Major Serpentine Hotspot

Albanian taxa and populations of the genus Odontarrhena are most promising candidates for research on metal tolerance and Ni-agromining, but their genetic structure remains unknown. We investigated phylogenetic relationships and genetic differentiation in relation to distribution and ploidy of the taxa, anthropic site disturbance, elevation, soil type, and trace metals at each population site. After performing DNA sequencing of selected accessions, we applied DNA-fingerprinting to analyze the genetic structure of 32 populations from ultramafic and non-ultramafic outcrops across Albania. Low sequence divergence resulted in poorly resolved phylograms, but supported affinity between the two diploid serpentine endemics O. moravensis and O. rigida. Analysis of molecular variance (AMOVA) revealed significant population differentiation, but no isolation by distance. Among-population variation was higher in polyploids than in diploids, in which genetic distances were lower. Genetic admixing at population and individual level occurred especially in the polyploids O. chalcidica, O. decipiens, and O. smolikana. Admixing increased with site disturbance. Outlier loci were higher in serpentine populations but decreased along altitude with lower drought and heat stress. Genetic variability gained by gene flow and hybridization at contact zones with "resident" species of primary ultramafic habitats promoted expansion of the tetraploid O. chalcidica across anthropogenic sites.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Ethanol extract of Liriodendron chinense (Hemsl.) Sarg barks attenuates hyperuricemic nephropathy by inhibiting renal fibrosis and inflammation in mice

ETHNOPHARMACOLOGICAL RELEVANCE

Liriodendron chinense (Hemsl.) Sarg, known as the Chinese tulip tree, has a long history of cultivation and utilization in many Asia countries, especially in China to use in traditional Chinese medicine for expelling "wind and dampness", a term corresponding to rheumatic fever and rheumatoid arthritis. Interestingly, the barks of Liriodendron chinense (Hemsl.) Sarg was also found in folk to treat gout. However, further experimental studies remained to confirm its uric acid-lowering effects.

AIM OF THE STUDY

The aim of the study was to evaluate the protective effect of ethanol extract of the barks of Liriodendron chinense (Hemsl.) Sarg (EELC) in a mouse model of hyperuricemic nephropathy (HN) and the involved mechanisms.

MATERIALS AND METHODS

EELC at a respective dose of 250 mg/kg/d or 500 mg/kg/d were orally administered to HN mice induced by a mixture of adenine (160 mg/kg/d)/potassium oxonate (2.4 g/kg/d) for 21 days. At the end of the treatment, serum uric acid, kidney functions (serum creatinine, blood urea nitrogen and urine microalbumin), 24-h urine uric acid excretion, as well as kidney pathological changes were investigated by biochemical assay, histopathological score, immunofluorescence and histochemistry, RT-qPCR, and western blotting analysis.

RESULTS AND DISCUSSION

Oral administration of EELC significantly lowered serum uric acid level at 500 mg/kg (185.75 ± 15.49 μmol/L of EELC vs. 238.28 ± 20.97 μmol/L of HN model, p < 0.01) in HN mice. EELC at 500 mg/kg also remarkably reduced the levels of serum creatinine (82.92 ± 7.86 μmol/L of EELC vs. 92.08 ± 6.13 μmol/L of HN model, p < 0.0001), blood urea nitrogen (21.50 ± 1.87 mmol/L of EELC vs. 29.40 ± 3.95 mmol/L of HN model, p < 0.001) and urine microalbumin (4.25 ± 0.40 mg/L of EELC vs. 5.95 ± 0.33 mg/L of HN model, p < 0.001) to improve renal function. It also attenuated renal fibrosis, especially the high-dose of EELC. Furthermore, EELC could inhibit the activation of NF-κB, ASK1/JNK/c-Jun, JAK2/STAT3 signaling pathways and reduce the release of pro-inflammatory cytokine TNF-α in the kidneys of HN mice. Additionally, EELC remarkably increased urine uric acid excretion of HN mice, which may be achieved by the upregulation of organic anion transporter 1 (OAT1), OAT3 and ATP-binding cassette subfamily G member 2 (ABCG2) proteins.

CONCLUSIONS

EELC alleviated the progression of HN by suppressing the activation of NF-κB, ASK1/JNK/c-Jun and JAK2/STAT3 signaling pathway, reducing the infiltration of inflammatory factors and uric acid accumulation in the kidney.

Genetic evidence for allopatric speciation of the Siberian ibex Capra sibirica in India

Montane systems, formed by a series of climatic oscillations and temporal topographic metamorphoses, have broken up the contiguous distribution of widespread species and accelerated allopatric speciation. We used a partial fragment of the mitochondrial cytochrome b (cytb) gene to address speciation across the entire range of the Siberian ibex Capra sibirica. We demonstrated that the Siberian ibex is a polytypic species, plausibly formed by a combination of at least 2 species and/or 3 to 4 sub-species. Bayesian phylogeny showed that the Indian-Tajikistan (I-T) clade is adequately diverged from the other clades based on the mean intra-specific distance criterion, and warrants recognition as a distinct species. We provide pragmatic evidence for the endorsement of the I-T clade as a distinct species of Siberian ibex and urge prioritization of the conservation of this species at global and regional scales.

Chloroplast population genetics reveals low levels of genetic variation and conformation to the central-marginal hypothesis in Taxus wallichiana var. mairei, an endangered conifer endemic to China

The central-marginal hypothesis predicts that geographically peripheral populations should exhibit reduced genetic diversity and increased genetic differentiation than central populations due to smaller effective population size and stronger geographical isolation. We evaluated these predictions in the endangered conifer Taxus wallichiana var. mairei. Eight plastid simple sequence repeats (cpSSRs) were used to investigate plastid genetic variation in 22 populations of Taxus wallichiana var. mairei, encompassing nearly its entire distribution range. Low levels of plastid genetic variation and differentiation were detected in the populations, and the findings were attributed to low mutation rates, small population sizes, habitat fragmentation and isolation, and effective pollen or seed dispersal. Hunan and Hubei were identified as major refugia based on the number of private haplotypes and species distribution modeling. Trends in plastid genetic diversity and genetic differentiation from central to peripheral populations supported the predictions of the central-marginal hypothesis. In scenarios wherein the future climate becomes warmer, we predict that some peripheral populations will disappear and southern and southeastern regions will become significantly less habitable. Factors that include the levels of precipitation during the driest month, annual precipitation level, and annual temperature range will be decisive in shaping the future distribution of these populations. This study provides a theoretical basis for the conservation of T. wallichiana var. mairei.

Integrating Phylogeographic Analysis and Geospatial Methods to Infer Historical Dispersal Routes and Glacial Refugia of Liriodendron chinense

Liriodendron chinense (Hemsl.), a Tertiary relic tree, is mainly distributed in subtropical China. The causes of the geographical distribution pattern of this species are poorly understood. In this study, we inferred historical dispersal routes and glacial refugia of this species by combining genetic data (chloroplast DNA (cpDNA), nuclear ribosomal DNA (nrDNA), and nuclear DNA (nDNA)) and geospatial data (climate and geology) with the methods of landscape genetics. Additionally, based on sequence variation at multiple loci, we employed GenGIS and Barrier software to analyze L. chinense population genetic structure. Dispersal corridors and historical gene flow between the eastern and western populations were detected, and they were located in mountainous regions. Based on species distribution model (SDMs), the distribution patterns in paleoclimatic periods were consistent with the current pattern, suggesting the presence of multiple refuges in multiple mountainous regions in China. The genetic structure analysis clustered most eastern populations into a clade separated from the western populations. Additionally, a genetic barrier was detected between the eastern and western populations. The dispersal corridors and historical gene flow detected here suggested that the mountains acted as a bridge, facilitating gene flow between the eastern and western populations. Due to Quaternary climatic fluctuations, the habitats and dispersal corridors were frequently inhabited by warm-temperate evergreen forests, which may have fragmented L. chinense habitats and exacerbated the differentiation of eastern and western populations. Ultimately, populations retreated to multiple isolated mountainous refugia, shaping the current geographical distribution pattern. These dispersal corridors and montane refugia suggested that the mountains in subtropical China play a crucial role in the conservation of genetic resources and migration of subspecies or related species in this region.

New insight into the phylogeographic pattern of Liriodendron chinense (Magnoliaceae) revealed by chloroplast DNA: east-west lineage split and genetic mixture within western subtropical China

Background

Subtropical China is a global center of biodiversity and one of the most important refugia worldwide. Mountains play an important role in conserving the genetic resources of species. Liriodendron chinense is a Tertiary relict tree largely endemic to subtropical China. In this study, we aimed to achieve a better understanding of the phylogeographical pattern of L. chinense and to explore the role of mountains in the conservation of L. chinense genetic resources.

Methods

Three chloroplast regions (psbJ-petA, rpl32-ndhF, and trnK5’-matK) were sequenced in 40 populations of L. chinense for phylogeographical analyses. Relationships among chloroplast DNA (cpDNA) haplotypes were determined using median-joining networks, and genetic structure was examined by spatial analysis of molecular variance (SAMOVA). The ancestral area of the species was reconstructed using the Bayesian binary Markov Chain Monte Carlo (BBM) method according to its geographic distribution and a maximum parsimony (MP) tree based on Bayesian methods.

Results

Obvious phylogeographic structure was found in L. chinense. SAMOVA revealed seven groups matching the major landscape features of the L. chinense distribution area. The haplotype network showed three clades distributed in the eastern, southwestern, and northwestern regions. Separate northern and southern refugia were found in the Wu Mountains and Yungui Plateau, with genetic admixture in the Dalou Mountains and Wuling Mountains. BBM revealed a more ancient origin of L. chinense in the eastern region, with a west–east split most likely having occurred during the Mindel glacial stage.

Discussion

The clear geographical distributions of haplotypes suggested multiple mountainous refugia of L. chinense. The east–west lineage split was most likely a process of gradual genetic isolation and allopatric lineage divergence when the Nanling corridor was frequently occupied by evergreen or coniferous forest during Late Quaternary oscillations. Hotspots of haplotype diversity in the Dalou Mountains and Wuling Mountains likely benefited from gene flow from the Wu Mountains and Yungui Plateau. Collectively, these results indicate that mountain regions should be the main units for conserving and collecting genetic resources of L. chinense and other similar species in subtropical China.

The influence of a relict distribution on genetic structure and variation in the Mediterranean tree, Platanus orientalis

The distribution of plant species around the Mediterranean basin is a product of the influence of both geographical barriers and of climatic changes experienced during the Tertiary, with the transition from a warm to cool periods. Several species, once largely distributed across the Northern Hemisphere, retracted to refugial areas in southern Europe where they are described as Tertiary relicts. Platanus orientalis is a typical representative of Tertiary flora in southwest Eurasia; its distribution spreads from the Caucasus to the Mediterranean, with its western border in Southern Italy and Sicily. We analysed genetic diversity and differentiation in the centre and western part of its distribution range using nuclear microsatellites and compared genetic parameters between core and western populations. We found an overall decrease in genetic diversity estimates (heterozygosity, private allelic richness) from central towards western populations, with those from Southern Italy and Sicily showing the lowest values. The low level of genetic diversity probably results from historic range fragmentation experienced by P. orientalis in its westernmost distribution as confirmed by high level genetic isolation of these populations. Ornamental hybrids were genetically distinguished from P. orientalis as contained private alleles, indicating that gene flow with natural populations is rare. Population assignment and neighbour-joining (NJ) analysis of populations identified four groups belonging to two main phyletic groups (the Southern Italian-Sicilian-Balkan and Cretan-Bulgarian-Turkish lineages) that seem to have different biogeographic origin and also excluded an artificial origin for southern Italian and Sicilian populations as previously suggested. These results show that quantifying the genetic variation of a Tertiary relict in a geographical context, and the potential effect of hybridization with introduced ornamental hybrids, can provide useful insights on factors influencing population genetic structure. Such information is crucial to predict how such taxa will respond to increasing anthropogenic influence on the environment.

RAD-Seq Data Point to a Distinct Split in Liriodendron (Magnoliaceae) and Obvious East–West Genetic Divergence in L. chinense

Liriodendron is a Tertiary period relic tree genus with a typical East Asian and North American disjunction distribution pattern. As an angiosperm base group of trees, Liriodendron provides a valuable resource for the study of evolution processes. Here, we reconstruct the phylogeny and population genetic structure of Liriodendron based on the restriction site-associated DNA sequencing (RAD-Seq) of a wide collection of individuals from 16 populations. Our results reveal a clear phylogenetic break between L. chinense and L. tulipifera and obvious genetic divergence between the eastern and western populations of L. chinense, which are consistent with the patterns of geographical distributions. The phylogeographic history and long-term geographical isolation of the genus may be responsible for this pattern. Furthermore, a closer relationship was found between L. tulipifera and the eastern populations of L. chinense, indicating the ancient phylogeny of L. chinense in this area. The results of this study will aid in the development of scientific strategies for the conservation and utilization of the Liriodendron germplasm.

Interspecies evolutionary divergence in Liriodendron, evidence from the nucleotide variations of LcDHN-like gene

Background

Liriodendron is a genus of Magnoliaceae, which consists of two relict species, Liriodendron chinense and L. tulipifera. Although the morphologies are highly similar, the two species exhibit different adaptive capacity. Dehydrins (DHNs) are abiotic stresses resistant proteins in planta, which are associated with adaptive evolution. To better understand the evolution divergence between L. chinense and L. tulipifera and how DHN genes are associated with adaptation evolution, we firstly investigated the DNA polymorphisms of the LcDHN-like gene in 21 L. chinense and 6 L. tulipifera populations.

Results

A 707 bp LcDHN-like gene was cloned, which included a 477 bp open reading frame (ORF) and coding 158 amino acids. 311 LcDHN-like gDNA sequences were obtained from 70 L. chinense and 35 L. tulipifera individuals. The AMOVA and phylogenetic relationship analysis showed significant differences between the two species. A higher genetic diversity was observed in L. tulipifera compared to L. chinense, in consistent with the higher adaptive capacity of L. tulipifera. Our data also suggested that the LcDHN-like genes’ polymorphisms were under neutral mutation and purifying selection model in the L. chinense and L. tulipifera populations, respectively. The distinct expanding range and rate between the two species, haplotypes shared only in L.chinense’s nearby populations, and wide dispersals in L. tulipifera could contribute to the obscure east-west separation in L. chinense and entirely unordered phylogeny in L. tulipifera. The completely separated nonsynonymous substitution at position 875 and the higher range scope of aliphatic index in L. tulipifera populations may be related with its higher adaptive capacity. Taken together, our study suggests LcDHN-like gene is a potential mark gene responsible for adaptive evolution divergence in Liriodendron.

Conclusions

Significant differences and completely distinct haplogroups between L. chinense and L. tulipifera showed that the two species have evolved into different directions. The more widely distribution, earlier haplogroups divergence events, and richer SNPs variations in L. tulipifera could imply its stronger adaptation in this species. And potential effect of the allelic variations in LcDHN-like gene may reflect the difference of water stress and chill tolerance between L. chinense and L. tulipifera, which could provide some information for further adaption evolution studies of Liriodendron.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1318-7) contains supplementary material, which is available to authorized users.

AFLP Approach Reveals Variability in Phragmites australis: Implications for Its Die-Back and Evidence for Genotoxic Effects

Phragmites australis is a subcosmopolitan species typical of wetlands being studied in Europe for its disappearance from natural stands, a phenomenon called reed die-back syndrome (RDBS). Although it is conjectured that low genetic variability contributes to RDBS, this aspect remains neglected to this day. Using a molecular fingerprinting approach and a sequence analysis of the trnT-trnL/rbcL-psaI regions of cpDNA, this study aimed to compare the genetic structure of stable vs. RDBS-affected P. australis stands from five wetlands of central Italy. Beforehand, in order to characterize the health condition of reed populations, the occurrence of the main macromorphological descriptors for RDBS was considered on 40 reed stands. Soil samples were also collected to examine the total content of heavy metals. The current study analyzed cpDNA in 19 samples and AFLP profiles in 381 samples to investigate the genetic structure of Phragmites populations. Based on the multinomial-Dirichlet model, an analysis of candidate loci under selective pressure was also performed. The relationships among AFLP data, RDBS descriptors and chemicals were evaluated with the use of Linear Mixed Models. The analysis of the cpDNA shows the occurrence of the haplotypes M (the most widespread), and K here recorded for the first time in Italy. Three new haplotypes were also described. The DNA fingerprinting analysis has produced a total of 322 loci (98% polymorphic) and shows the medium-to-high amount of genetic diversity. The significant genetic differentiation among wetlands (Fst = 0.337) suggests either low gene flow or small effective population size. Moreover, the low amount of outlier loci (only 5; l.5% of the total), seems to indicate the scarce occurrence of selective pressure upon the reed's genome. Genetic diversity increased in relationship to the decrease in diameter and of flowering buds of the reed, two of the trends associated with the die-back. The current study rejects the hypothesis that genetic diversity massively contributed to RDBS. Moreover, significant relationships between genetic diversity and the total concentration of some heavy metals (Cr, Cu, and Zn) were highlighted, indicating possible genotoxic effects on P. australis. The current study represents a fact-finding background useful for the conservation of common reed.

Phylogeography and genetic effects of habitat fragmentation on endemic Urophysa (Ranunculaceae) in Yungui Plateau and adjacent regions

Urophysa is a Chinese endemic genus with only two species (U. rockii and U. henryi) distributed in Yungui Plateau (Guizhou Province) and adjacent regions (i.e., Provinces of Hunan, Hubei, Chongqing and Sichuan). The aim of this study was to determine the genetic diversity and population differentiation within Urophysa and investigate the effect of the Yungui Plateau uplift and climate oscillations on evolution of Urophysa. In this study, micro-morphological characteristics, nine microsatellite loci (SSR), two nuclear loci (ITS and ETS) and two chloroplast fragments (psbA-trnH and trnL-trnF) were used to analyze the phylogenetic relationships and assess genetic and phylogeographical structure of Urophysa. Isolation by distance (IBD) was performed to research the effects of geographical isolation. We detected high genetic diversity at the species level but low genetic diversity within populations. Striking genetic differentiation (AMOVA) among populations and a significant phylogeographical structure (NST > GST, p < 0.01) were detected among U. henryi populations, along with significant effects of isolation by distance (IBD). Molecular clock estimation using calibration strategy and cpDNA substitution rate indicated that the divergence of U. henryi occurred during late Miocene to early Quaternary, when the orogeny of Yungui Plateau was violent. U. rockii originated at the early Quaternary and further differentiated at early Pleistocene. Our results suggested that habitat fragmentation played an important role in the genetic diversity and population differentiation of U. rockii and U. henryi. Heterogenous geomorphological configuration and complicated environment resulted from rapid uplift of the Yungui Plateau were inferred as important incentives for the modern phylogeograhpical pattern and species divergence of Urophysa. The geographical isolation, limited gene flow, specialized morphologies and the Pleistocene climatic oscillation greatly contributed to the allopatric divergence of U. rockii. Significant genetic drift and inbreeding were detected in these two species, in situ measures should be implemented to protect them.

Two disjunct Pleistocene populations and anisotropic postglacial expansion shaped the current genetic structure of the relict plant Amborella trichopoda

Past climate fluctuations shaped the population dynamics of organisms in space and time, and have impacted their present intra-specific genetic structure. Demo-genetic modelling allows inferring the way past demographic and migration dynamics have determined this structure. Amborella trichopoda is an emblematic relict plant endemic to New Caledonia, widely distributed in the understory of non-ultramafic rainforests. We assessed the influence of the last glacial climates on the demographic history and the paleo-distribution of 12 Amborella populations covering the whole current distribution. We performed coalescent genetic modelling of these dynamics, based on both whole-genome resequencing and microsatellite genotyping data. We found that the two main genetic groups of Amborella were shaped by the divergence of two ancestral populations during the last glacial maximum. From 12,800 years BP, the South ancestral population has expanded 6.3-fold while the size of the North population has remained stable. Recent asymmetric gene flow between the groups further contributed to the phylogeographical pattern. Spatially explicit coalescent modelling allowed us to estimate the location of ancestral populations with good accuracy (< 22 km) and provided indications regarding the mid-elevation pathways that facilitated post-glacial expansion.

AFLP Genome Scanning Reveals Divergent Selection in Natural Populations of Liriodendron chinense (Magnoliaceae) along a Latitudinal Transect

Understanding adaptive genetic variation and its relation to environmental factors are important for understanding how plants adapt to climate change and for managing genetic resources. Genome scans for the loci exhibiting either notably high or low levels of population differentiation (outlier loci) provide one means of identifying genomic regions possibly associated with convergent or divergent selection. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) genome scan and environmental association analysis to test for signals of natural selection in natural populations of Liriodendron chinense (Chinese Tulip Tree; Magnoliaceae) along a latitudinal transect. We genotyped 276 individuals from 11 populations of L. chinense using 987 AFLP markers. Both frequency-based (Dfdist and BayeScan) and correlation-based (MLM) methods were applied to detect outlier loci. Our analyses recovered both neutral and potentially adaptive genetic differentiation among populations of L. chinense. We found moderate genetic diversity within populations and high genetic differentiation among populations with reduced genetic diversity toward the periphery of the species ranges. Nine AFLP marker loci showed evidence of being outliers for population differentiation for both detection methods. Of these, six were strongly associated with at least one climate factor. Temperature, precipitation, and radiation were found to be three important factors influencing local adaptation of L. chinense. The outlier AFLP loci are likely not the target of natural selection, but the neighboring genes of these loci might be involved in local adaptation. Hence, these candidates should be validated by further studies.