Late Miocene lineage divergence and ecological differentiation of rare endemic Juniperus blancoi: clues for the diversification of North American conifers

Potential Geographic Range of the Endangered Reed Parrotbill Paradoxornis heudei under Climate Change

The phenomenon of global climate change can impact the geographic range and biodiversity, thereby heightening the vulnerability of rare species to extinction. The reed parrotbill (Paradoxornis heudei David, 1872) is endemic to central and eastern China, it is mainly distributed in the middle and lower reaches of the Yangtze River Plain and the Northeast Plain. In this study, eight of ten algorithms of the species distribution model (SDM) were used to evaluate the impact of climate change on the potential distribution of P. heudei under current and future climate scenarios and to analyze the possible related climate factors. After checking the collected data, 97 occurrence records of P. heudei were used. The relative contribution rate shows that among the selected climatic variables, temperature annual range (bio7), annual precipitation (bio12), and isothermality (bio3) were the principal climatic factors to limit the habitat suitability of P. heudei. The suitable habitat for P. heudei is primarily concentrated in the central–eastern and northeast plains of China, particularly in the eastern coastal region, spanning a mere area of 57,841 km2. The habitat suitability of P. heudei under different representative concentration pathway (RCP) scenarios was predicted to be different under future climatic conditions, but all of them had a larger range than the current one. The species distribution range could expand by more than 100% on average compared with the current range under the four scenarios in 2050, while it could contract by approximately 30% on average relative to the 2050 range in 2070 under different climate change scenarios. In the future, northeastern China may serve as a potential suitable habitat for P. heudei. The changes in the spatial and temporal distributions of P. heudei’s range are of utmost importance in identifying high-priority conservation regions and devising effective management strategies for its preservation.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Incipient speciation, high genetic diversity, and ecological divergence in the alligator bark juniper suggest complex demographic changes during the Pleistocene

The most recent glacial cycles of the Pleistocene affected the distribution, population sizes, and levels of genetic structure of temperate-forest species in the main Mexican mountain systems. Our objective was to investigate the effects these cycles had on the genetic structure and distribution of a dominant species of the "mexical" vegetation across North and Central America. We studied the genetic diversity of Juniperus deppeana, a conifer distributed from the Southwestern United States to the highlands of Central America. We combined information of one plastid marker and two nuclear markers to infer phylogeographic structure, genetic diversity and demographic changes. We also characterized the climatic niche for each variety to infer the plausible area of suitability during past climatic conditions and to evaluate climatic niche discontinuities along with the species distribution. We found a marked phylogeographic structure separating the populations North and South of the Isthmus of Tehuantepec, with populations to the South of this barrier forming a distinct genetic cluster corresponding to Juniperus deppeana var. gamboana. We also found signals of population expansion in the Northern genetic cluster. Ecological niche modeling results confirmed climatic niche differences and discontinuities among J. deppeana varieties and heterogeneous responses to climatic oscillations. Overall, J. deppeana's genetic diversity has been marked by distribution shifts, population growth and secondary contact the North, and in situ permanence in the South since the last interglacial to the present. High genetic variation suggests a wide and climatically diverse distribution during climatic oscillations. We detected the existence of two main genetic clusters, supporting previous proposals that Juniperus deppeana and Juniperus gamboana may be considered two separate species.

Climatic effects on niche evolution in a passerine bird clade depend on paleoclimate reconstruction method

Climatic niches describe the climatic conditions in which species can persist. Shifts in climatic niches have been observed to coincide with major climatic change, suggesting that species adapt to new conditions. We test the relationship between rates of climatic niche evolution and paleoclimatic conditions through time for 65 Old-World flycatcher species (Aves: Muscicapidae). We combine niche quantification for all species with dated phylogenies to infer past changes in the rates of niche evolution for temperature and precipitation niches. Paleoclimatic conditions were inferred independently using two datasets: a paleoelevation reconstruction and the mammal fossil record. We find changes in climatic niches through time, but no or weak support for a relationship between niche evolution rates and rates of paleoclimatic change for both temperature and precipitation niche and for both reconstruction methods. In contrast, the inferred relationship between climatic conditions and niche evolution rates depends on paleoclimatic reconstruction method: rates of temperature niche evolution are significantly negatively related to absolute temperatures inferred using the paleoelevation model but not those reconstructed from the fossil record. We suggest that paleoclimatic change might be a weak driver of climatic niche evolution in birds and highlight the need for greater integration of different paleoclimate reconstructions.

Host genotype explains rhizospheric microbial community composition: the case of wild cotton metapopulations (Gossypium hirsutum L.) in Mexico

The rhizosphere provides several benefits to the plant host being a strong determinant for its health, growth and productivity. Nonetheless, the factors behind the assembly of the microbial communities associated with the rhizosphere such as the role of plant genotypes are not completely understood. In this study, we tested the role that intraspecific genetic variation has in rhizospheric microbial community assemblages, using genetically distinct wild cotton populations as a model of study. We followed a common garden experiment including five wild cotton populations, controlling for plant genotypes, environmental conditions and soil microbial community inoculum, to test for microbial differences associated with genetic variation of the plant hosts. Microbial communities of the treatments were characterized by culture-independent 16S rRNA gene amplicon sequencing with Illumina MiSeq platform. We analyzed microbial community diversity (alpha and beta), and diversity structure of such communities, determined by co-occurrence networks. Results show that different plant genotypes select for different and specific microbial communities from a common inoculum. Although we found common amplicon sequence variants (ASVs) to all plant populations (235), we also found unique ASVs for different populations that could be related to potential functional role of such ASVs in the rhizosphere.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of π_a/π_s and K_a/K_s with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

Population differentiation and phylogeography in Lycianthes moziniana (Solanaceae: Capsiceae), a perennial herb endemic to the Mexican Transition Zone

Lycianthes moziniana (Solanaceae: Capsiceae) is a perennial herb with edible fruits that is endemic to Mexico. Three varieties are recognized, all known in the Mexican Transition Zone. Lycianthes moziniana var. margaretiana grows in the Sierra Madre Oriental, whereas L. moziniana var. moziniana is common along the Trans-Mexican Volcanic Belt and in the Sierra Madre Occidental. Lycianthes moziniana var. oaxacana is found exclusively in the Sierra Madre del Sur. The Mexican Transition Zone is a complex geological, climatic and biogeographical area, the result of tectonic and volcanic activity that has promoted genetic divergence and speciation. We determined the genetic variation and structure of L. moziniana. Using phylogeographical approaches, we described the demographic history and evolutionary processes leading its divergence. The intergenic spacers rpl32-trnL and ycf1 were sequenced for 133 individuals pertaining to 15 populations. The genealogical relationships were analysed using haplotype networks. Finally, based on ecological niche models, we inferred the palaeodistribution of L. moziniana during the Pleistocene. The genetic differences and the haplogroups matched the three described varieties. Geological and climatic events of the Mexican Transition Zone facilitated these results. The Trans-Mexican Volcanic Belt isolated the populations of the Sierra Madre Oriental and the Sierra Madre del Sur, while allowing the migration to the Sierra Madre Occidental, during the middle Holocene.

Multilocus, phenotypic, behavioral, and ecological niche analyses provide evidence for two species within Euphonia affinis (Aves, Fringillidae)

The integration of genetic, morphological, behavioral, and ecological information in the analysis of species boundaries has increased, allowing integrative systematics that better reflect the evolutionary history of biological groups. In this context, the goal of this study was to recognize independent evolutionary lineages within Euphonia affinis at the genetic, morphological, and ecological levels. Three subspecies have been described: E. affinis godmani, distributed in the Pacific slope from southern Sonora to Guerrero; E. affinis affinis, from Oaxaca, Chiapas and the Yucatan Peninsula to Costa Rica; and E. affinis olmecorum from Tamaulipas and San Luis Potosi east to northern Chiapas (not recognized by some authors). A multilocus analysis was performed using mitochondrial and nuclear genes. These analyses suggest two genetic lineages: E. godmani and E. affinis, which diverged between 1.34 and 4.3 My, a period in which the ice ages and global cooling fragmented the tropical forests throughout the Neotropics. To analyze morphometric variations, six morphometric measurements were taken, and the Wilcoxon Test was applied to look for sexual dimorphism and differences between the lineages. Behavioral information was included, by performing vocalization analysis which showed significant differences in the temporal characteristics of calls. Finally, Ecological Niche Models were estimated with MaxEnt, and then compared using the method of Broennimann. These analyses showed that the lineage distributed in western Mexico (E. godmani) has a more restricted niche than the eastern lineage (E. affinis) and thus we rejected the hypotheses of niche equivalence and similarity. Based on the combined evidence from genetic, morphological, behavioral, and ecological data, it is concluded that E. affinis (with E. olmecorum as its synonym) and E. godmani represent two independent evolutionary lineages.

Using demographic model selection to untangle allopatric divergence and diversification mechanisms in the Rheum palmatum complex in the Eastern Asiatic Region

Allopatric divergence is often initiated by geological uplift and restriction to sky-islands, climate oscillations, or river capture. However, it can be difficult to establish which mechanism was the most likely to generate the current phylogeographical structure of a species. Recently, genomic data in conjunction with a model testing framework have been applied to address this issue in animals. To test whether such an approach is also likely to be successful in plants, we used population genomic data of the Rheum palmatum complex from the Eastern Asiatic Region, in conjunction with biogeographical reconstruction and demographic model selection, to identify the potential mechanism(s) which have led to the current level of divergence. Our results indicate that the R. palmatum complex originated in the central Hengduan Mts and possibly in regions further east, and then dispersed westward and eastward resulting in genetically distinct lineages. Populations are likely to have diverged in refugia during climate oscillations followed by subsequent expansion and secondary contact. However, model simulations within the western lineage of the R. palmatum complex cannot reject a restriction to sky-islands as a possible mechanism of diversification due to the genetically ambiguous position of one population. This highlights that genetically mixed populations might introduce ambiguity regarding the best diversification model in some cases. Although it might be possible to resolve this ambiguity using other data, sometimes this could prove to be difficult in complex biogeographical areas.

Multilocus Analyses Reveal Postglacial Demographic Shrinkage of Juniperus morrisonicola (Cupressaceae), a Dominant Alpine Species in Taiwan

Postglacial climate changes alter geographical distributions and diversity of species. Such ongoing changes often force species to migrate along the latitude/altitude. Altitudinal gradients represent assemblage of environmental, especially climatic, variable factors that influence the plant distributions. Global warming that triggered upward migrations has therefore impacted the alpine plants on an island. In this study, we examined the genetic structure of Juniperus morrisonicola, a dominant alpine species in Taiwan, and inferred historical, demographic dynamics based on multilocus analyses. Lower levels of genetic diversity in north indicated that populations at higher latitudes were vulnerable to climate change, possibly related to historical alpine glaciers. Neither organellar DNA nor nuclear genes displayed geographical subdivisions, indicating that populations were likely interconnected before migrating upward to isolated mountain peaks, providing low possibilities of seed/pollen dispersal across mountain ranges. Bayesian skyline plots suggested steady population growth of J. morrisonicola followed by recent demographic contraction. In contrast, most lower-elevation plants experienced recent demographic expansion as a result of global warming. The endemic alpine conifer may have experienced dramatic climate changes over the alternation of glacial and interglacial periods, as indicated by a trend showing decreasing genetic diversity with the altitudinal gradient, plus a fact of upward migration.

Evolutionary origin and demographic history of an ancient conifer (Juniperus microsperma) in the Qinghai-Tibetan Plateau

All Qinghai-Tibetan Plateau (QTP) endemic species are assumed to have originated recently, although very rare species most likely diverged early. These ancient species provide an excellent model to examine the origin and evolution of QTP endemic plants in response to the QTP uplifts and the climate changes that followed in this high altitude region. In this study, we examined these hypotheses by employing sequence variation from multiple nuclear and chloroplast DNA of 239 individuals of Juniperus microsperma and its five congeners. Both phylogenetic and population genetic analyses revealed that J. microsperma diverged from its sister clade comprising two species with long isolation around the Early Miocene, which corresponds to early QTP uplift. Demographic modeling and coalescent tests suggest that J. microsperma experienced an obvious bottleneck event during the Quaternary when the global climate greatly oscillated. The results presented here support the hypotheses that the QTP uplifts and Quaternary climate changes played important roles in shaping the evolutionary history of this rare juniper.