The role of glacial cycles in promoting genetic diversity in the Neotropics: the case of cloud forests during the Last Glacial Maximum

Climatic niche shifts and ecological sky-island dynamics in Mesoamerican montane birds

An ongoing challenge in evolutionary and ecological research focuses on testing biogeographic hypotheses for the understanding of both species' distributional patterns and of the factors influencing range limits. In this study, we described the climatic niches of Neotropical humid montane forest birds through the analysis of factors driving their evolution at inter- and intraspecific levels; and tested for differences among allopatric lineages within Aulacorhynchus, Chlorospingus, Cardellina, and Eupherusa. We employed ecological niche models (ENMs) along with an ordination approach with kernel smoothing to perform niche overlap analyses and test hypotheses of niche equivalence/similarity among lineages. In addition, we described the potential distributions of each lineage during the Late Pleistocene climate fluctuations, identifying historical range expansions, connectivity, and stability. Overall, we observed differences in environmental variables influencing climatic requirements and distributional patterns for our selected species. We detected the highest values of niche overlap mainly between Eupherusa and some Chlorospingus lineages. At both interspecific and intraspecific levels, sister lineages showed non-identical environmental niches. Our results offer weak support to a moist forest model, in which populations followed the expansion and contraction cycles of montane forests, leading to a lack of niche conservatism among lineages (they tend to occupy not identical climatic environments) throughout Mesoamerica. Therefore, historical climatic conditions may act as ecological barriers determining the distributional ranges of these species.

Montane diversification as a mechanism of speciation in neotropical butterflies

The mountains in the Atlantic Forest domain are environments that harbor a high biodiversity, including species adapted to colder climates that were probably influenced by the climatic variations of the Pleistocene. To understand the phylogeographic pattern and assess the taxonomic boundaries between two sister montane species, a genomic study of the butterflies Actinote mantiqueira and A. alalia (Nymphalidae: Acraeini) was conducted. Analyses based on partial sequences of the mitochondrial gene COI (barcode region) failed to recover any phylogenetic or genetic structure discriminating the two species or sampling localities. However, single nucleotide polymorphisms gathered using Genotyping-by-Sequencing provided a strong isolation pattern in all analyses (genetic distance, phylogenetic hypothesis, clustering analyses, and F ST statistics) which is consistent with morphology, separating all individuals of A. alalia from all populations of A. mantiqueira. The three sampled mountain ranges where A. mantiqueira populations occur-Serra do Mar, Serra da Mantiqueira, and Poços de Caldas Plateau-were identified as three isolated clusters. Paleoclimate simulations indicate that both species' distributions changed according to climatic oscillations in the Pleistocene period, with the two species potentially occurring in areas of lower altitude during glacial periods when compared to the interglacial periods (as the present). Besides, a potential path between their distribution through the Serra do Mar Mountain range was inferred. Therefore, the Pleistocene climatic fluctuation had a significant impact on the speciation process between A. alalia and A. mantiqueira, which was brought on by isolation at different mountain summits during interglacial periods, as shown by the modeled historical distribution and the observed genetic structure.

Recent genetic, phenetic and ecological divergence across the Mesoamerican highlands: a study case with Diglossa baritula (Aves: Thraupidae)

The topographical, geological, climatic and biodiversity complexity of Mesoamerica has made it a primary research focus. The Mesoamerican highlands is a region with particularly high species richness and within-species variation. The Cinnamon-bellied Flowerpiercer, Diglossa baritula (Wagler, 1832), is a species endemic to the Mesoamerican highlands, with three allopatric subspecies currently recognized. To characterize divergence within this species, we integrated genomics, morphology, coloration and ecological niche modeling approaches, obtained from sampling individuals across the entire geographic distribution of the species. Our results revealed a clear genomic divergence between the populations to the east versus the west of the Isthmus of Tehuantepec. In contrast to the genomic results, morphology and coloration analyses showed intermediate levels of differentiation, indicating that population groups within D. baritula have probably been under similar selective pressures. Our morphology results indicated that the only sexually dimorphic morphological variable is the wing chord, with males having a longer wing chord than females. Finally, ecological data indicated that there are differences in ecological niche within D. baritula. Our data suggest that D. baritula could contain two or more incipient species at the intermediate phase of the speciation continuum. These results highlight the importance of the geographical barrier of the Isthmus of Tehuantepec and Pleistocene climatic events in driving isolation and population divergence in D. baritula. The present investigation illustrates the speciation potential of the D. baritula complex and the capacity of Mesoamerican highlands to create cryptic biodiversity and endemism.

An Elevational Phylogeographic Diversity Gradient in Neotropical Birds Is Decoupled from Speciation Rates

AbstractA key question about macroevolutionary speciation rates is whether they are controlled by microevolutionary processes operating at the population level. For example, does spatial variation in population genetic differentiation underlie geographical gradients in speciation rates? Previous work suggests that speciation rates increase with elevation in Neotropical birds, but underlying population-level gradients remain unexplored. Here, we characterize elevational phylogeographic diversity between montane and lowland birds in the megadiverse Andes-Amazonian system and assess its relationship to speciation rates to evaluate the link between population-level differentiation and species-level diversification. We aggregated and georeferenced nearly 7,000 mitochondrial DNA sequences across 103 species or species complexes in the Andes and Amazonia and used these sequences to describe phylogeographic differentiation across both regions. Our results show increased levels of both discrete and continuous metrics of population structure in the Andean mountains compared with the Amazonian lowlands. However, higher levels of population differentiation do not predict higher rates of speciation in our dataset. Multiple potential factors may lead to our observed decoupling of initial population divergence and speciation rates, including the ephemerality of incipient species and the multifaceted nature of the speciation process, as well as methodological challenges associated with estimating rates of population differentiation and speciation.

Past, present and future in the geographical distribution of Mexican Tepezmaite cycads: Genus Ceratozamia

Ceratozamia morettii, C. brevifrons, and C. tenuis are cycads considered endangered in montane forests in the center of Veracruz state. However, the amount of theoretical and empirical information available on the historical distribution of these species and how they could be affected in the future by the effects of climate change still needs to be increased. Our objective was to generate information on the spatial distribution of the species since the last glacial maximum, present, and future. To map the spatial distribution of species, we created a potential distribution model for each species. The spatial data used for the models included 19 bioclimatic data variables in the present, at the last glacial maximum using two models (CCSM4 and MIROC), and in the future (2080) using two models of the RCP 8.5 scenario of climate change (HadGEM2-CC and MIROC5). We found that each species occupies a unique ecoregion and climatic niche. Ceratozamia morettii and C. tenuis have a similar pattern with an expansion of their distribution area since the last glacial maximum with a larger distribution area in the present and a projected reduction in their distribution under future climatic conditions. For C. brevifrons, we also showed an increase in their distributional area since the last glacial maximum. We also showed that this expansion will continue under future climatic conditions when the species reaches its maximum distributional area. Projections about the future of these endemic cycad species show changes in their habitat, highlighting that temperate zone species (C. morettii and C. tenuis) will face imminent extinction if no effort is made to protect them. On the other hand, the tropical climate species (C. brevifrons) will be favored.

Nuclear and chloroplast DNA phylogeography reveals high genetic diversity and postglacial range expansion in Quercus mexicana

PREMISE

Phylogeographical studies are fundamental for understanding factors that influence the spatial distribution of genetic lineages within species. Population expansions and contractions, distribution shifts, and climate changes are among the most important factors shaping the genetic compositions of populations.

METHODS

We investigated the phylogeography of an endemic oak, Quercus mexicana (Fagaceae), which has a restricted distribution in northeastern Mexico along the Sierra Madre Oriental and adjacent areas. Nuclear and chloroplast DNA microsatellite markers were used to describe the genetic diversity and structure of 39 populations of Q. mexicana along its entire distribution area. We tested whether population expansion or contraction events influenced the genetic diversity and structure of the species. We also modeled the historical distributional range of Q. mexicana (for the Mid Holocene, the Last Glacial Maximum, and the Last Interglacial) to estimate the extent to which climate fluctuations have impacted the distribution of this oak species.

RESULTS

Our results revealed high genetic diversity and low genetic structure in Q. mexicana populations. Ecological niche models suggested historical fluctuations in the distributional range of Q. mexicana. Historical range changes, gene flow, and physical barriers seem to have played an important role in shaping the phylogeographic structure of Q. mexicana.

CONCLUSIONS

Our study indicates that the genetic structure of Q. mexicana may have been the result of responses of oak trees not only to heterogeneous environments present in the Sierra Madre Oriental and adjacent areas, but also to elevational and latitudinal shifts in response to climate changes in the past.

Determinants of genetic diversity in Neotropical salamanders (Plethodontidae: Bolitoglossini)

Genetic diversity is the raw material of evolution, yet the reasons why it varies among species remain poorly understood. While studies at deeper phylogenetic scales point to the influence of life history traits on genetic diversity, it appears to be more affected by population size but less predictable at shallower scales. We used proxies for population size, mutation rate, direct selection, and linked selection to test factors affecting genetic diversity within a diverse assemblage of Neotropical salamanders, which vary widely for these traits. We estimated genetic diversity of noncoding loci using ddRADseq and coding loci using RNAseq for an assemblage of Neotropical salamanders distributed from northern Mexico to Costa Rica. Using ddRADseq loci, we found no significant association with genetic diversity, while for RNAseq data we found that environmental heterogeneity and proxies of population size predict a substantial portion of the variance in genetic diversity across species. Our results indicate that diversity of coding loci may be more predictable than that of noncoding loci, which appears to be mostly unpredictable at shallower phylogenetic scales. Our results suggest that coding loci may be more appropriate for genetic diversity estimates used in conservation planning because of the lack of any association between the variables we used and genetic diversity of noncoding loci.

Plastome phylogenomics reveals an early Pliocene North- and Central America colonization by long-distance dispersal from South America of a highly diverse bromeliad lineage

Understanding the spatial and temporal frameworks of species diversification is fundamental in evolutionary biology. Assessing the geographic origin and dispersal history of highly diverse lineages of rapid diversification can be hindered by the lack of appropriately sampled, resolved, and strongly supported phylogenetic contexts. The use of currently available cost-efficient sequencing strategies allows for the generation of a substantial amount of sequence data for dense taxonomic samplings, which together with well-curated geographic information and biogeographic models allow us to formally test the mode and tempo of dispersal events occurring in quick succession. Here, we assess the spatial and temporal frameworks for the origin and dispersal history of the expanded clade K, a highly diverse Tillandsia subgenus Tillandsia (Bromeliaceae, Poales) lineage hypothesized to have undergone a rapid radiation across the Neotropics. We assembled full plastomes from Hyb-Seq data for a dense taxon sampling of the expanded clade K plus a careful selection of outgroup species and used them to estimate a time- calibrated phylogenetic framework. This dated phylogenetic hypothesis was then used to perform biogeographic model tests and ancestral area reconstructions based on a comprehensive compilation of geographic information. The expanded clade K colonized North and Central America, specifically the Mexican transition zone and the Mesoamerican dominion, by long-distance dispersal from South America at least 4.86 Mya, when most of the Mexican highlands were already formed. Several dispersal events occurred subsequently northward to the southern Nearctic region, eastward to the Caribbean, and southward to the Pacific dominion during the last 2.8 Mya, a period characterized by pronounced climate fluctuations, derived from glacial-interglacial climate oscillations, and substantial volcanic activity, mainly in the Trans-Mexican Volcanic Belt. Our taxon sampling design allowed us to calibrate for the first time several nodes, not only within the expanded clade K focal group but also in other Tillandsioideae lineages. We expect that this dated phylogenetic framework will facilitate future macroevolutionary studies and provide reference age estimates to perform secondary calibrations for other Tillandsioideae lineages.

Late Pleistocene Altitudinal Segregation and Demography Define Future Climate Change Distribution of the Peromyscus mexicanus Species Group: Conservation Implications

Mountains harbor a significant number of the World's biodiversity, both on tropical and temperate regions. Notably, one crucial gap in conservation is the consideration of historical and contemporary patterns influencing differential distribution in small mammal mountain species and how climate change will affect their distribution and survival. The mice Peromyscus mexicanus species group is distributed across mountains in Guatemala-Chiapas and Central America, which experienced significant effects of glacial and interglacial cycles. We determined phylogeographic and demographic patterns of lowlands and highlands mountain lineages, revealing that the radiation of modern P. mexicanus lineages occurred during the Pleistocene (ca. 2.6 mya) along Nuclear Central America. In concert with climatic cycles and the distribution of habitats, lowland and highland lineages showed recent population size increase and decrease, respectively. We also estimated the current and future distribution ranges for six lineages, finding marked area size increase for two lineages for which vegetation type and distribution would facilitate migrating towards higher elevations. Contrastingly, three lineages showed range size decrease; their ecological requirements make them highly susceptible to future habitat loss. Our findings are clear evidence of the negative impacts of future climate change, while our ability to manage and conserve these vulnerable ecosystems and mountain species is contingent on our understanding of the implications of climate change on the distribution, ecology, and genetics of wildlife populations.

Predicting geographic distribution and habitat suitability of Opuntia streptacantha in paleoclimatic, current, and future scenarios in Mexico

Mexican territory is one of the centers of origin and dispersion of the genus Opuntia, where several of its species have been an important plant resource for people in arid and semiarid zones. Opuntia streptacantha is widely distributed in Mexico; however, precise aspects of its geographic distribution and ecological status are still unknown. Here, we modeled its potential distribution under paleoclimatic, current, and future conditions through maximum entropy and predictions from 824 records and seven environmental variables. Potential distribution of O. streptacantha in the interglacial period was contracted and slightly north than current distribution, with 44,773 km2 of optimal habitat. In other past periods, the central location of potential distribution coincides with the actual current distribution, but the period of the last glacial maximum was characterized by 201 km2 of very suitable habitat, absent in interglacial, current, and future periods. The future model suggests that potential distribution will move toward the south of the Mexican territory. Synthesis and applications. The potential distribution of O. streptacantha can be applied for the conservation and management of the species, and also in the selection of areas with crassicaule scrubs for protection, conservation, and reproduction of species resistant to the hostile conditions of arid and semiarid zones of Mexican territory, where the structure and composition of the vegetation will be affected in the next 100 years.

Population Genetic Structure and Diversity of Cryptic Species of the Plant Genus Macrocarpaea (Gentianaceae) from the Tropical Andes

The Pleistocene climatic oscillations (PCO) that provoked several cycles of glacial-interglacial periods are thought to have profoundly affected species distribution, richness and diversity around the world. While the effect of the PCO on population dynamics at temperate latitudes is well known, considerable questions remain about its impact on the biodiversity of neotropical mountains. Here, we use amplified fragment length polymorphism molecular markers (AFLPs) to investigate the phylogeography and genetic structure of 13 plant species belonging to the gentian genus Macrocarpaea (Gentianaceae) in the tropical Andes. These woody herbs, shrubs or small trees show complex and potentially reticulated relationships, including cryptic species. We show that populations of M. xerantifulva in the dry system of the Rio Marañón in northern Peru have lower levels of genetic diversity compared to other sampled species. We suggest that this is due to a recent demographic bottleneck resulting from the contraction of the montane wet forests into refugia because of the expansion of the dry system into the valley during the glacial cycles of the PCO. This may imply that the ecosystems of different valleys of the Andes might have responded differently to the PCO.

Evolutionary origins of the eastern North American-Mesoamerican floristic disjunction: Current status and future prospects

Biogeographic disjunction patterns, where multiple taxa are shared between isolated geographic areas, represent excellent systems for investigating the historical assembly of modern biotas and fundamental biological processes such as speciation, diversification, niche evolution, and evolutionary responses to climate change. Studies on plant genera disjunct across the northern hemisphere, particularly between eastern North America (ENA) and eastern Asia (EAS), have yielded tremendous insight on the geologic history and assembly of rich temperate floras. However, one of the most prevalent disjunction patterns involving ENA forests has been largely overlooked: that of taxa disjunct between ENA and cloud forests of Mesoamerica (MAM), with examples including Acer saccharum, Liquidambar styraciflua, Cercis canadensis, Fagus grandifolia, and Epifagus virginiana. Despite the remarkable nature of this disjunction pattern, which has been recognized for over 75 years, there have been few recent efforts to empirically examine its evolutionary and ecological origins. Here I synthesize previous systematic, paleobotanical, phylogenetic, and phylogeographic studies to establish what is known about this disjunction pattern to provide a roadmap for future research. I argue that this disjunction pattern, and the evolution and fossil record of the Mexican flora more broadly, represents a key missing piece in the broader puzzle of northern hemisphere biogeography. I also suggest that the ENA-MAM disjunction represents an excellent system for examining fundamental questions about how traits and life history strategies mediate plant evolutionary responses to climate change and for predicting how broadleaf temperate forests will respond to the ongoing climatic pressures of the Anthropocene.

Genomic-based microsatellite development for Ternstroemia (Pentaphylacaceae) and transferability to other Ericales

BACKGROUND

The genus Ternstroemia is associated with the vulnerable tropical montane cloud forest in Mexico and with other relevant vegetation types worldwide. It contains threatened and pharmacologically important species and has taxonomic issues regarding its species limits. This study describes 38 microsatellite markers generated using a genomic-based approach.

METHODS AND RESULTS

We tested 23 of these markers in a natural population of Ternstroemia lineata. These markers are highly polymorphic (all loci polymorphic with 3-14 alleles per locus and expected heterozygosity between 0.202 and 0.908), most of them (19 out of 23) are in Hardy-Weinberg Equilibrium and free of null alleles (18 out of 23). Also we found no evidence of linkage among them. Finally, we tested the transferability to six other American species of Ternstroemia, two other Pentaphylacaceae species, and four species from different families within the order Ericales.

CONCLUSIONS

These molecular resources are promising tools to investigate genetic diversity loss and as barcodes for ethnopharmacological applications and species delimitation in the family Pentaphylacaceae and some Ericales, among other applications.

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.

Phylogeography supports lineage divergence for an endemic rattlesnake (Crotalus ravus) of the Neotropical montane forest in the Trans-Mexican Volcanic Belt

The formation of the Trans-Mexican Volcanic Belt (TMVB) and Pleistocene climatic fluctuations have been shown to influence the diversification of lineages and species distributed throughout central Mexico. In some taxa, however, evidence of lineage diversification is not easily recognized, as often is the case in reptiles. Here we present a phylogeographic study on a Mexican endemic rattlesnake species (Crotalus ravus), with the aim of understanding how distinct lineages are distributed across the TMVB. Genetic (mtDNA) and genomic (ddRADseq) data were generated from samples across the species’ range to evaluate phylogeographic structure, estimate phylogenetic relationships and divergence times, and perform environmental niche modeling (ENM). Both datasets recover strong phylogeographic structuring of two distinct lineages on an east-west axis, with an estimated Pleistocene divergence (~1.47 Myr). The ENM suggest that the distribution of the two lineages experienced expansion and reduction events throughout recent evolutionary time. We attribute the diversification of C. ravus lineages to geological events associated with the formation of the TMVB, as well as Quaternary climate changes, both of which have been previously recognized in co-distributed taxa in the TMVB. This work emphasizes the existence of cryptic diversification processes in a morphologically conserved species distributed in a region of complex climatic and orogenic heterogeneity.

Himalayan orogeny and monsoon intensification explain species diversification in an endemic ginger (Hedychium: Zingiberaceae) from the Indo-Malayan Realm

The Indo-Malayan Realm is a biogeographic realm that extends from the Indian Subcontinent to the islands of Southeast Asia (Malay Archipelago). Despite being megadiverse, evolutionary hypotheses explaining taxonomic diversity in this region have been rare. Here, we investigate the role of geoclimatic events such as Himalayan orogeny and monsoon intensification in the diversification of the ginger-lilies (Hedychium J.Koenig: Zingiberaceae). We first built a comprehensive, time-calibrated phylogeny of Hedychium with 75% taxonomic and geographic sampling. We found that Hedychium is a very young lineage that originated in Northern Indo-Burma, in the Late Miocene (c. 10.6 Ma). This was followed by a late Neogene and early Quaternary diversification, with multiple dispersal events to Southern Indo-Burma, Himalayas, Peninsular India, and the Malay Archipelago. The most speciose clade IV i.e., the predominantly Indo-Burmese clade also showed a higher diversification rate, suggesting its recent rapid radiation. Our divergence dating and GeoHiSSE results demonstrate that the diversification of Hedychium was shaped by both the intensifications in the Himalayan uplift as well as the Asian monsoon. Ancestral character-state reconstructions identified the occurrence of vegetative dormancy in both clades I and II, whereas the strictly epiphytic growth behavior, island dwarfism, lack of dormancy, and a distinct environmental niche were observed only in the predominantly island clade i.e., clade III. Finally, we show that the occurrence of epiphytism in clade III corresponds with submergence due to sea-level changes, suggesting it to be an adaptive trait. Our study highlights the role of recent geoclimatic events and environmental factors in the diversification of plants within the Indo-Malayan Realm and the need for collaborative work to understand biogeographic patterns within this understudied region. This study opens new perspectives for future biogeographic studies in this region and provides a framework to explain the taxonomic hyperdiversity of the Indo-Malayan Realm.

Population structure in Neotropical plants: Integrating pollination biology, topography and climatic niches

Animal pollinators mediate gene flow among plant populations, but in contrast to well-studied topographic and (Pleistocene) environmental isolating barriers, their impact on population genetic differentiation remains largely unexplored. Comparing how these multifarious factors drive microevolutionary histories is, however, crucial for better resolving macroevolutionary patterns of plant diversification. Here we combined genomic analyses with landscape genetics and niche modelling across six related Neotropical plant species (424 individuals across 33 localities) differing in pollination strategy to test the hypothesis that highly mobile (vertebrate) pollinators more effectively link isolated localities than less mobile (bee) pollinators. We found consistently higher genetic differentiation (FST ) among localities of bee- than vertebrate-pollinated species with increasing geographical distance, topographic barriers and historical climatic instability. High admixture among montane populations further suggested relative climatic stability of Neotropical montane forests during the Pleistocene. Overall, our results indicate that pollinators may differentially impact the potential for allopatric speciation, thereby critically influencing diversification histories at macroevolutionary scales.

Comparing the population history of Neotropical annual species: The role of climate change and hybridization between Tithonia tubaeformis and T. rotundifolia (Asteraceae)

Historical climate changes and interspecific gene flow have played an important role in shaping the distribution and genetic diversity of the biota in the Neotropics. In this study, we explored the role of both the Pleistocene climate changes and hybridization on the contemporary geographic structure of two Neotropical weedy species: Tithonia tubaeformis and T. rotundifolia. Distribution shifts under past and current climate conditions were explored through ecological niche modelling (ENM). We then tested these hypotheses using chloroplast microsatellite (cpSSR) data in T. tubaeformis and compared the patterns deduced with those previously reported for T. rotundifolia using the same cpSSR loci. Lastly, we searched for shared haplotypes between species. Both species exhibited significant downwards altitudinal shifts during the last interglacial (LIG) and the Last Glacial Maximum (LGM). For instance, T. rotundifolia showed large suitable habitat areas since the LIG at xeric conditions in western Mesoamerica. Tithonia tubaeformis showed a ~five-fold range contraction during the LIG compared to current climate conditions, followed by a range expansion in the LGM. Despite the large shared refugial areas predicted through ENM, we found a low number of shared haplotypes, suggesting a minor role of hybridization in shaping the geographic structure of these species. Our results provide additional patterns of the population history of the northern Neotropics during the Quaternary, and we suggest that weedy widespread species are a well-suited group for the study of the effects of historic climatic changes on the biota of this region.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Evolutionary conservatism will limit responses to climate change in the tropics

Rapid species turnover in tropical mountains has fascinated biologists for centuries. A popular explanation for this heightened beta diversity is that climatic stability at low latitudes promotes the evolution of narrow thermal tolerance ranges, leading to local adaptation, evolutionary divergence and parapatric speciation along elevational gradients. However, an emerging consensus from research spanning phylogenetics, biogeography and behavioural ecology is that this process rarely, if ever, occurs. Instead, closely related species typically occupy a similar elevational niche, while species with divergent elevational niches tend to be more distantly related. These results suggest populations have responded to past environmental change not by adapting and diverging in place, but instead by shifting their distributions to tightly track climate over time. We argue that tropical species are likely to respond similarly to ongoing and future climate warming, an inference supported by evidence from recent range shifts. In the absence of widespread in situ adaptation to new climate regimes by tropical taxa, conservation planning should prioritize protecting large swaths of habitat to facilitate movement.

Discordance in maternal and paternal genetic markers in lesser long-nosed bat Leptonycteris yerbabuenae, a migratory bat: recent expansion to the North and male phylopatry

Leptonycteris yerbabuenae, the lesser long-nosed bat is an abundant migratory nectar-feeding bat found in most of Mexico, and in some areas of northern Central America and small sections of southwestern USA. We analyzed the distribution of the maternal and paternal lineages of this species with phylogeographic methods based on two mitochondrial markers, Cyt-b and D-loop, and a marker located in the Y chromosome, DBY. We obtained tissue samples from 220 individuals from 23 localities. Levels of genetic diversity (haplotype diversity, Hd ) were high (Cyt-b = 0.757; D-loop = 0.8082; DBY = 0.9137). No clear patterns of population genetic structure were found for mitochondrial markers, while male genetic differentiation suggested the presence of two lineages: one from Mexican Pacific coast states and another from central-southern Mexico; in accordance to strong male philopatry and higher female migration. We used genealogical reconstructions based on Bayesian tools to calculate divergence times, and to test coalescent models to explain changes in L. yerbabuenae historical demography. Our results show that recent demographic changes were consistent with global climatic changes (∼130,000 kyr ago for Cyt-b and ∼160,000 kyr for D-loop) and divergence times dated from molecular genealogies exhibited older divergence times, Cyt-b (4.03 mya), D-loop (10.26 mya) and DBY (12.23 mya). Accordingly, the female lineage underwent demographic expansion associated to Pleistocene climate change, whereas the male lineage remained constant.

Influence of Pleistocene climatic oscillations on the phylogeography and demographic history of endemic vulnerable trees (section Magnolia) of the Tropical Montane Cloud Forest in Mexico

The Tropical Montane Cloud Forest (TMCF) is a highly dynamic ecosystem that has undergone frequent spatial changes in response to the interglacial-glacial cycles of the Pleistocene. These climatic fluctuations between cold and warm cycles have led to species range shifts and contractions-expansions, resulting in complex patterns of genetic structure and lineage divergence in forest tree species. In this study, we sequenced four regions of the chloroplast DNA (trnT-trnL, trnK5-matk, rpl32-trnL, trnS-trnG) for 20 populations and 96 individuals to evaluate the phylogeography, historical demography, and paleodistributions of vulnerable endemic TMCF trees in Mexico: Magnolia pedrazae (north-region), M. schiedeana (central-region), and M. schiedeana population Oaxaca (south-region). Our data recovered 49 haplotypes that showed a significant phylogeographic structure in three regions: north, central, and south. Bayesian Phylogeographic and Ecological Clustering (BPEC) analysis also supported the divergence in three lineages and highlighted the role of environmental factors (temperature and precipitation) in genetic differentiation. Our historical demography analyses revealed demographic expansions predating the Last Interglacial (LIG, ~125,000 years ago), while Approximate Bayesian Computation (ABC) simulations equally supported two contrasting demographic scenarios. The BPEC and haplotype network analyses suggested that ancestral haplotypes were geographically found in central Veracruz. Our paleodistributions modeling showed evidence of range shifts and expansions-contractions from the LIG to the present, which suggested the complex evolutionary dynamics associated to the climatic oscillations of the Pleistocene. Habitat management of remnant forest fragments where large and genetically diverse populations occur in the three TMCF regions analyzed would be key for the conservation of these magnolia populations.

Genetic imprints of Brosimum alicastrum Sw. in Mexico

PREMISE

The mechanisms generating the geographical distributions of genetic diversity are a central theme in evolutionary biology. The amount of genetic diversity and its distribution are controlled by several factors, including dispersal abilities, physical barriers, and environmental and climatic changes. We investigated the patterns of genetic diversity and differentiation among populations of the widespread species Brosimum alicastrum in Mexico.

METHODS

Using nuclear DNA microsatellite data, we tested whether the genetic structure of B. alicastrum was associated with the roles of the Trans-Mexican Volcanic Belt and the Isthmus of Tehuantepec as geographical barriers to gene flow and to infer the role of past events in the genetic diversity patterns. We further used a maximum-likelihood population-effects mixed model (MLPE) to identify the main factor affecting population differentiation in B. alicastrum.

RESULTS

Our results suggested that Mexican B. alicastrum is well differentiated into three main lineages. Patterns of the genetic structure at a finer scale did not fully correspond to the current geographical barriers to gene flow. According to the MLPE mixed model, isolation by distance is the best model for explaining the genetic differentiation of B. alicastrum in Mexico.

CONCLUSIONS

We propose that the differentiation patterns might reflect (1) an ancient differentiation that occurred in Central and South America, (2) the effects of past climatic changes, and (3) the functions of some physical barriers to gene flow. This study provides insights into the possible mechanisms underlying the geographic genetic variation of B. alicastrum along a moisture gradient in tropical lowland forests.

Potential distribution of Amblyomma mixtum (Koch, 1844) in climate change scenarios in the Americas

Amblyomma mixtum is a Neotropical generalist tick of medical and veterinary importance which is widely distributed from United States of America to Ecuador. The aim of this study was to evaluate changes in the geographic projections of the ecological niche models of A. mixtum in climate change scenarios in America. We constructed a database of published scientific publications, personal collections, personal communications, and online databases. Ecological niche modelling was performed with 15 Bioclimatic variables using kuenm in R and was projected to three time periods (Last Glacial Maximum, Current and 2050) for America. Our model indicated a wide distribution for A. mixtum, with higher probability of occurrence along the Gulf of Mexico and occurring in a lesser proportion in the Pacific states, Central America, and the northern part of South America. The areas of new invasion are located mainly on the border of Mexico with Guatemala and Belize, some regions of Central America and Colombia. We conclude that the ecological niche modelling are effective tools to infer the potential distribution of A. mixtum in America, in addition to helping to propose future measures of epidemiological control and surveillance in the new potential areas of invasion.

Rugged relief and climate promote isolation and divergence between two neotropical cold-associated birds

The role of historical factors in establishing patterns of diversity in tropical mountains is of interest to understand the buildup of megadiverse biotas. In these regions, the historical processes of range fragmentation and contraction followed by dispersal are thought to be mediated by the interplay between rugged relief (complex topography) and climate fluctuations and likely explain most of the dynamics of diversification in plants and animals. Although empirical studies addressing the interaction between climate and topography have provided invaluable insights into population divergence and speciation patterns in tropical montane organisms, a more detailed and robust test of such processes in an explicit spatio-temporal framework is still lacking. Consequently, our ability to gain insights into historical range shifts over time and the genomic footprint left by them is limited. Here, we used niche modeling and subgenomic population-level datasets to explore the evolution of two species of warbling finches (genus Microspingus) disjunctly distributed across the Montane Atlantic Forest, a Neotropical region with complex geological and environmental histories. Population structure inferences suggest a scenario of three genetically differentiated populations, which are congruent with both geography and phenotypic variation. Demographic simulations support asynchronous isolation of these populations as recently as ∼40,000 years ago, relatively stable population sizes over recent time, and past gene flow subsequent to divergence. Throughout the last 800,000 years, niche models predicted extensive expansion into lowland areas with increasing overlap of species distributions during glacial periods, with prominent retractions and isolation into higher altitudes during interglacials, which are in line with signs of introgression of currently isolated populations. These results support a dual role of cyclical climatic changes: population divergence and persistence in mountain tops during warm periods followed by periods of expansion and admixture in lower elevations during cold periods. Our results underscore the role of the interplay between landscape and climate as an important mechanism in the evolution of the Neotropical montane biota.

The DNA history of a lonely oak: Quercus humboldtii phylogeography in the Colombian Andes

The climatic and geological changes that occurred during the Quaternary, particularly the fluctuations during the glacial and interglacial periods of the Pleistocene, shaped the population demography and geographic distribution of many species. These processes have been studied in several groups of organisms in the Northern Hemisphere, but their influence on the evolution of Neotropical montane species and ecosystems remains unclear. This study contributes to the understanding of the effect of climatic fluctuations during the late Pleistocene on the evolution of Andean mountain forests. First, we describe the nuclear and plastidic DNA patterns of genetic diversity, structure, historical demography, and landscape connectivity of Quercus humboldtii, which is a typical species in northern Andean montane forests. Then, these patterns were compared with the palynological and evolutionary hypotheses postulated for montane forests of the Colombian Andes under climatic fluctuation scenarios during the Quaternary. Our results indicated that populations of Q. humboldtii have high genetic diversity and a lack of genetic structure and that they have experienced a historical increase in connectivity from the last glacial maximum (LGM) to the present. Furthermore, our results showed a dramatic reduction in the effective population size followed by an expansion before the LGM, which is consistent with the results found by palynological studies, suggesting a change in dominance in Andean forests that may be related to ecological factors rather than climate change.

Morphological and genetic diversification of Russula floriformis, sp. nov., along the Isthmus of Panama

Species of Russula are ubiquitous members of ectomycorrhizal fungal communities in tropical ecosystems. However, an important part of the total tropical diversity of this genus and its biogeographic patterns is unknown due to the lack of studies on Russula in tropical ecosystems. We combined molecular, morphological, ecological, and biogeographic data to elaborate concepts for two new subspecies of R. floriformis (subsection Substriatinae). Russula floriformis subsp. floriformis and R. floriformis subsp. symphoniae are described as new from montane forest dominated by Quercus and/or Oreomunnea (Fagales) from Colombia and Panama, respectively. Phylogenies were constructed using nuc rDNA internal transcribed spacer region ITS1-5.8S-ITS2 (ITS), D1-D2 domains of nuc 28S rDNA (28S), and partial regions of the second largest subunit of RNA polymerase II (rpb2) and translation elongation factor 1-alpha (tef1). Similar environmental conditions, similar morphology, and an ITS sequence similarity higher than 99% with only three different positions indicate that these two subspecies are closely related. Detailed observations of microscopic structures and analyses of further DNA loci, however, revealed morphological and molecular characteristics that allow distinguishing the two subspecies of R. floriformis. Spatial distribution and phylogenetic proximity of the two Russula subspecies and their ectomycorrhizal hosts, i.e., species of Quercus, suggest that their diversification is a result of comigration, adaptation, and geographic isolation along the Isthmus of Panama during the Pliocene and Pleistocene.

Idiosyncratic responses to drivers of genetic differentiation in the complex landscapes of Isthmian Central America

Isthmian Central America (ICA) is one of the most biodiverse regions in the world, hosting an exceptionally high number of species per unit area. ICA was formed <25 million years ago and, consequently, its biotic assemblage is relatively young and derived from both colonization and in situ diversification. Despite intensive taxonomic work on the local fauna, the potential forces driving genetic divergences and ultimately speciation in ICA remain poorly studied. Here, we used a landscape genetics approach to test whether isolation by distance, topography, habitat suitability, or environment drive the genetic diversity of the regional frog assemblage. To this end, we combined data on landscape features and mitochondrial DNA sequence variation for nine codistributed amphibian species with disparate life histories. In five species, we found that at least one of the factors tested explained patterns of genetic divergence. However, rather than finding a general pattern, our results revealed idiosyncratic responses to historical and ecological processes, indicating that intrinsic life-history characteristics may determine the effect of different drivers of isolation on genetic divergence in ICA. Our work also suggests that the convergence of several factors promoting isolation among populations over a heterogeneous landscape might maximize genetic differentiation, despite short geographical distances. In conclusion, abiotic factors and geographical features have differentially affected the genetic diversity across the regional frog assemblage. Much more complex models (i.e., considering multiple drivers), beyond simple vicariance of Caribbean and Pacific lineages, are needed to better understand the evolutionary history of ICA's diverse biotas.

Niche divergence and paleo-distributions of Lutzomyia longipalpis mitochondrial haplogroups (Diptera: Psychodidae)

Lutzomyia longipalpis is a complex of species which has a wide but discontinuous distribution from southeastern Mexico to northern Argentina and Uruguay. To date, eight mitochondrial haplogroups have been identified along its distribution although key environmental tolerances and ecological niche models have been analyzed only at the complex level. The aim of the present study was to analyze whether genetic diversification using three mitochondrial genes of the Lu. longipalpis complex is associated with niche divergence and to explore evolution of distributional projections of all haplogroups between the Last Glacial Maximum (LGM; 21,000 yrs ago) and the present. Current occurrence of all haplogroups was used to develop ecological niche models (ENM) and these were projected in both periods to quantify and identify geographic area shifts. Environmental space was used to estimate niche similarity between major clades and pairwise between individual haplogroups. The two major Lu. longipalpis clades (Mex, CA, Col and Ven vs Arg and Bra) had significantly different environmental space, indicating niche divergence. Environmental space overlap of southern haplogroups was variable, with divergent niche, except between Arg and ArgBra. The most suitable regions for the ArgBra haplogroup were northeastern and southeastern Brazil, and the Gran Chaco region. In contrast, ENM of haplogroups within the northern major clade have significantly similar niche, with highest geographic ENM suitability along both the Caribbean and Pacific coasts. The intensity and coverage of high suitability areas in the LGM decreased for most haplogroups in the present. Integrating ENM and phylogenetic analyses has allowed us to test hypotheses of niche similarity between Lu. longipalpis haplogroups and major clades, and to identify conserved distributional areas of haplogroups since the LGM, with the exception of Arg. Evidence for distributional shifts and overlap of haplogroups is important to analyze Leishmaniasis´ eco-epidemiology and to successfully monitor and control transmission.

The role of environment, local adaptation, and past climate fluctuation on the amount and distribution of genetic diversity in two subspecies of Mexican wild Zea mays

PREMISE

Past climate fluctuations during the Holocene and Pleistocene shaped the distribution of several plant species in temperate areas over the world. Wild maize, commonly known as teosinte, is a good system to evaluate the effects of historical climate fluctuations on genetic diversity due to its wide distribution in Mexico with contrasting environmental conditions. We explored the influence of contemporary factors and historical environmental shifts on genetic diversity, including present and three historical periods using neutral markers.

METHODS

We used 22 nuclear microsatellite loci to examine the genetic diversity of 14 populations of Zea mays subsp. parviglumis and 15 populations of Zea mays subsp. mexicana (527 individuals total). We implemented genetic structure analyses to evaluate genetic differentiation between and within subspecies. We applied coalescent-based demographic analysis and species distribution modeling to evaluate the effects of historical environmental shifts.

RESULTS

We found 355 alleles in total for the two subspecies and variable levels of diversity in each (Z. mays subsp. parviglumis expected heterozygosity H_E = 0.3646-0.7699; Z. mays subsp. mexicana H_E = 0.5885-0.7671). We detected significant genetic structure among populations (D_EST = 0.4332) with significant heterozygote deficiency (F_IS = 0.1796), and variable selfing rates (sg₂ = 0.0-0.3090). The Bayesian assignment analysis differentiated four genetic groups. Demographic and species distribution modeling analysis suggested that environmental shifts were influential in the amount of genetic diversity.

CONCLUSIONS

Our analyses suggest that the current genetic diversity in teosinte is shaped by factors such as local adaptation and genetic isolation, along with historical environmental fluctuations.

Long-term cloud forest response to climate warming revealed by insect speciation history

Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change.

Altitudinal, latitudinal and longitudinal responses of cloud forest species to Quaternary glaciations in the northern Neotropics

The glaciations of the Quaternary caused changes in the geographical distributions of species associated with tropical montane cloud forests. The most obvious effect of the glacial conditions was the downward displacement of cloud forest species, thus giving opportunities for population connectivity in the lowlands. Considerable attention has been paid to these altitudinal changes, but latitudinal and longitudinal movements remain poorly understood in the northern Neotropics. Here, I use ecological niche modelling to generate palaeodistributions of small-eared shrews (Mammalia: Soricidae) closely associated with cloud forests in the mountain systems of Mexico and then retrodict their range shifts during the Last Glacial Maximum (LGM), one of the coldest periods of the Quaternary. The results suggest that cloud forest species not only migrated downwards in response to global cooling and dryness but also migrated latitudinally and longitudinally onto those slopes that maintained moist conditions (other slopes remained unsuitable during the LGM), thus revealing a hitherto unknown route for postglacial colonization of cloud forest species. This scenario of past distributional change probably had genetic and demographic implications and has repercussions for the identification of areas of refugia and postglacial colonization routes of cloud forest species.

High genetic diversity and stable Pleistocene distributional ranges in the widespread Mexican red oak Quercus castanea Née (1801) (Fagaceae)

The Mexican highlands are areas of high biological complexity where taxa of Nearctic and Neotropical origin and different population histories are found. To gain a more detailed view of the evolution of the biota in these regions, it is necessary to evaluate the effects of historical tectonic and climate events on species. Here, we analyzed the phylogeographic structure, historical demographic processes, and the contemporary period, Last Glacial Maximum (LGM) and Last Interglacial (LIG) ecological niche models of Quercus castanea, to infer the historical population dynamics of this oak distributed in the Mexican highlands. A total of 36 populations of Q. castanea were genotyped with seven chloroplast microsatellite loci in four recognized biogeographic provinces of Mexico: the Sierra Madre Occidental (western mountain range), the Central Plateau, the Trans-Mexican Volcanic Belt (TMVB, mountain range crossing central Mexico from west to east) and the Sierra Madre del Sur (SMS, southern mountain range). We obtained standard statistics of genetic diversity and structure and tested for signals of historical demographic expansions. A total of 90 haplotypes were identified, and 29 of these haplotypes were restricted to single populations. The within-population genetic diversity was high (mean h S = 0.72), and among-population genetic differentiation showed a strong phylogeographic structure (N ST = 0.630 > G ST = 0.266; p < .001). Signals of demographic expansion were identified in the TMVB and the SMS. The ecological niche models suggested a considerable percentage of stable distribution area for the species during the LGM and connectivity between the TMVB and the SMS. High genetic diversity, strong phylogeographic structure, and ecological niche models suggest in situ permanence of Q. castanea populations with large effective population sizes. The complex geological and climatic histories of the TMVB help to explain the origin and maintenance of a large proportion of the genetic diversity in this oak species.

Nuclear phylogeography of the temperate tree species Chiranthodendron pentadactylon (Malvaceae): Quaternary relicts in Mesoamerican cloud forests

Background

The Mexican hand tree or Canac (Chiranthodendron pentadactylon) is a temperate tree species of cloud and pine-oak forests of southern Mexico and Guatemala. Its characteristic hand-shaped flower is used in folk medicine and has constituted the iconic symbol of the Sociedad Botánica de México since 1940. Here, the evolutionary history of this species was estimated through phylogeographic analyses of nuclear DNA sequences obtained through restriction site associated DNA sequencing and ecological niche modeling. Total genomic DNA was extracted from leaf samples obtained from a representative number (5 to 10 per sampling site) of individuals distributed along the species geographic range. In Mexico, population is comprised by spatially isolated individuals which may follow the trends of cloud forest fragmentation. By contrast, in Guatemala Chiranthodendron may constitute a canopy dominant species near the Acatenango volcano. The distributional range of this species encompasses geographic provinces separated by the Isthmus of Tehuantepec.

The objectives of the study were to: (i) estimate its genetic structure to define whether the observed range disjunction exerted by the Isthmus of Tehuantepec translates into separate populations, (ii) link population divergence timing and demographic trends to historical climate change, and (iii) test hypotheses related to Pleistocene refugia.

Results

Patterns of genetic diversity indicated high levels of genetic differentiation between populations separated by the Isthmus. The western and eastern population diverged approximately 0.873 Million years ago (Ma). Demographic analyses supported a simultaneous split from an ancestral population and rapid expansion from a small stock approximately 0.2 Ma corresponding to a glacial period. The populations have remained stable since the LIG (130 Kilo years ago (Ka)). Species distribution modelling (SDM) predicted a decrease in potential distribution in the Last Interglacial (LIG) and an increase during the Last Glacial Maximum (LGM) (22 Ka), Mid-Holocene (6 Ka) and present times.

Conclusions

Divergence time estimations support the hypothesis that populations represent Quaternary relict elements of a species with broader and northernmost distribution. Pleistocene climatic shifts exerted major influence on the distribution of populations allowing dispersion during episodes of suitable climatic conditions and structuring during the first interglacial with a time period length of 100 Kilo years (Kyr) and the vicariant influence of the Isthmus. Limited demographic expansion and population connectivity during the LGM supports the moist forest hypothesis model.

Climatic dynamics and topography control genetic variation in Atlantic Forest montane birds

Montane organisms responded to Quaternary climate change by tracking suitable habitat along elevational gradients. However, it is unclear whether these past climatic dynamics generated predictable patterns of genetic diversity in co-occurring montane taxa. To test if the genetic variation is associated with historical changes in the elevational distribution of montane habitats, we integrated paleoclimatic data and a model selection approach for testing the demographic history of five co-distributed bird species occurring in the southern Atlantic Forest sky islands. We found that changes in historical population sizes and current genetic diversity are attributable to habitat dynamics among time periods and the current elevational distribution of populations. Taxa with populations restricted to the more climatically dynamic southern mountain block (SMB) had, on average, a six-fold demographic expansion, whereas the populations from the northern mountain block (NMB) remained constant. In the current configuration of the southern Atlantic Forest montane habitats, populations in the SMB have more widespread elevational distributions, occur at lower elevations, and harbor higher levels of genetic diversity than NMB populations. Despite the apparent coupling of demographic and climatic oscillations, our data rejected simultaneous population structuring due to historical habitat fragmentation. Demographic modeling indicated that the species had different modes of differentiation, and varied in the timing of divergence and the degree of gene flow across mountain blocks. Our results suggest that the heterogeneous distribution of genetic variation in birds of the Atlantic Forest sky islands is associated with the interplay between topography and climate of distinct mountains, leading to predictable patterns of genetic diversity.

Phylogeography and Genetic Diversity in a Southern North American Desert: Agave kerchovei From the Tehuacán-Cuicatlán Valley, Mexico

The Tehuacán-Cuicatlán Valley, located at the southeast of the state of Puebla and the northeast of the state of Oaxaca in Central Mexico, south of the Trans-Mexican Volcanic Belt (TMVB), is of particular interest for understanding the evolutionary dynamics of arid and semi-arid environments, being one of the main reservoirs of biological diversity for the arid zones of North America, including the highest diversity of Agavaceae worldwide and high levels of endemism. Studying in detail the phylogeography, environmental history and population genetics of representative species will hopefully shed light on the evolutionary and ecological dynamics that generated the tremendous biodiversity and endemism of this important region in Mexico. We sequenced three non-coding regions of chloroplast genome of Agave kerchovei, a representative species of the Tehuacán Valley, generating 2,188 bp from 128 individuals sampled from eight populations throughout the species range. We used this data set to (i) characterize the levels of genetic diversity and genetic structure in A. kerchovei; (ii) predict the distribution of A. kerchovei for the present day, and to reconstruct the past geographical history of the species by constructing ecological niche models (ENM); and (iii) compare the levels of diversity in this species with those estimated for the widely distributed Agave lechuguilla. Agave kerchovei has high levels of total chloroplast genetic variation (Hd = 0.718), especially considering that it is a species with a very restricted distribution. However, intrapopulation diversity is low (zero in some populations), and genetic structure is high (F _ST = 0.928, G _ST = 0.824), which can be expected for endemic species with isolated populations. Our data suggest that Pleistocene glacial cycles have played an important role in the distribution of A. kerchovei, where the climatic variability of the region - likely associated with its topographic complexity - had a significant effect on the levels of genetic diversity and population dynamics, while the potential distribution of the species seems to be stable since the middle Holocene (6 kya). We conclude that in A. kerchovei there is a core group of populations in the Tehuacán Valley, and peripheric populations that appear to be evolving independently and thus the species is fundamentally an endemic species from the Tehuacán Valley while the populations outside the Valley appear to be in the process of incipient speciation.

The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico

Background

Mesoamerica is a remarkable region with a high geological and ecological complexity. Within northern Mesoamerica, the biotic province of the Sierra Madre del Sur (SMS) in southwestern Mexico harbors exceptionally high avian endemism and diversity. Herein, we searched for spatially and temporally concordant phylogeographic patterns, in four bird genera from three distinct avian orders co-distributed across Mesoamerica and investigated their causes through hypothesis testing regarding historical processes. Selected species include endemic and differentiated populations across the montane forests of Mesoamerica, and particularly within the SMS.

Results

We gathered mitochondrial DNA sequences for at least one locus from 177 individuals across all species. We assessed genetic structure, demographic history, and defined a framework for the coalescent simulations used in biogeographic hypothesis testing temporal and spatial co-variance. Our analyses suggested shared phylogeographic breaks in areas corresponding to the SMS populations, and between the main montane systems in Mesoamerica, with the Central Valley of Oaxaca and the Nicaragua Depression being the most frequently shared breaks among analyzed taxa. Nevertheless, dating analyses and divergence patterns observed were consistent with the hypothesis of broad vicariance across Mesoamerica derived from mechanisms operating at distinct times across taxa in the SMS.

Conclusions

Our study provides a framework for understanding the evolutionary origins and historical factors enhancing speciation in well-defined regions within Mesoamerica, indicating that the evolutionary history of extant biota inhabiting montane forests is complex and often idiosyncratic.

Sufficient versus optimal climatic stability during the Late Quaternary: using environmental quality to guide phylogeographic inferences in a Neotropical montane system

Quaternary climatic oscillations affected species distributions worldwide, creating cycles of connectivity and isolation that impacted population demography and promoted lineage divergence. These effects have been well studied in temperate regions. Taxa inhabiting mesic montane habitats in tropical ecosystems show high levels of endemism and diversification in the distinct mountain ranges they inhabit; such a pattern has commonly been ascribed to past climatic oscillations, but few phylogeographic studies have tested this hypothesis. Here, we combine ecological niche models of species distributions with molecular data to study phylogeographic patterns in two rodents endemic to the highlands of Costa Rica and western Panama (Reithrodontomys creper and Nephelomys devius). In so doing, we apply a novel approach that incorporates a basic ecological principle: the expected positive relationship between environmental suitability and population abundance. Specifically, we use niche models to predict potential patterns of population connectivity and stability of different suitability levels during climatic extremes of the last glacial–interglacial cycle; we then test these predictions with population genetic analyses of a mitochondrial and a nuclear marker. The detailed predictions arising from the different levels of suitability were moderately to highly congruent with the molecular data depending on the species. Overall, results suggest that in these tropical montane ecosystems, cycles of population connectivity and isolation followed a pattern opposite to that typically described for temperate or lowland tropical ecosystems: namely, higher connectivity during the colder glacials, with isolation in montane refugia during the interglacials, including today. Nevertheless, the individualistic patterns for each species indicate a potentially wide gamut of phylogeographic histories reflecting particularities of their niches. Taken together, this study illustrates how phylogeographic inferences may benefit from niche model outputs that provide more detailed predictions of connectivity and finer characterizations of potential refugia through time.

Contrasting Patterns of Population History and Seed-mediated Gene Flow in Two Endemic Costa Rican Oak Species

Lower Central America is an important area to study recent population history and diversification of Neotropical species due to its complex and dynamic geology and climate. Phylogeographic studies in this region are few in comparison with other regions and even less for tree species. The aim of the present study was to characterize the phylogeographic structure in 2 partially co-distributed endemic oak species (Quercus costaricensis and Q. bumelioides) of the Costa Rican mountains using chloroplast short sequence repeats (cpSSRs), and to test for the effect of geological and palaeoclimatic processes on their population history. Genetic diversity and structure, haplotype networks, patterns of seed-mediated gene flow and historical demography were estimated for both species. Results suggested contrasting patterns. Quercus costaricensis exhibited high values of genetic diversity, a marked phylogeographic structure, a north-to-south genetic diversity gradient and evidence of a demographic expansion during the Quaternary. Quercus bumelioides did not show significant genetic structure and the haplotype network and historical demography estimates suggested a recent population expansion probably during the Pleistocene-Holocene transition. The phylogeographic structure of Q. costaricensis seems to be related to Pleistocene altitudinal migration due to its higher altitudinal distribution. Meanwhile, historical seed-mediated gene flow through the lower altitudinal distribution of Q. bumelioides may have promoted the homogenization of genetic variation. Population expansion and stable availability of suitable climatic areas in both species probably indicate that palaeoclimatic changes promoted downwards altitudinal migration and formation of continuous forests allowing oak species to expand their distribution into the Panamanian mountains during glacial stages.

Species distribution modelling and cpSSR reveal population history of the Neotropical annual herb Tithonia rotundifolia (Asteraceae)

The impacts of the historical geologic and climatic events on the diversity and genetic structure of Neotropical taxa have recently become a subject of study. However, annual plants associated with tropical dry forests remain under-studied. The exploration of additional taxa in contrasting environments will improve the current understanding of responses of the Neotropical biota to these events. Here, we explore the species distribution and geographic structure of the annual herb Tithonia rotundifolia. We sampled 175 individuals from 19 populations of T. rotundifolia. Species distribution modelling and six microsatellite chloroplast loci were used to infer its population history. We identified areas of historical climate suitability and then tested if there is genetic structuring among these areas. Haplotypes showed strong phylogeographic structure. Historical climatic suitability areas were found along the Pacific coast; however, a gap was found at the Isthmus of Tehuantepec (IT). Although Bayesian analysis showed population structuring, amova revealed that the IT is not its main driver. Instead, a subdivision into a higher number of regions had higher F_CT values. Also, populations to the east of the IT showed evidence of recent population expansion and migration in a south-north direction. Pleistocene climate fluctuations partially explain the geographic structure of T. rotundifolia. However, life-history characteristics such as limited seed dispersal and the patchy distribution of suitable habitats explain the high haplotype diversity and population sub-structuring and diversity. Lastly, the absence of geographic structure of some haplotypes may indicate long-distance dispersal, or hybridisation with the closely related T. tubaeformis.

Phylogeography indicates incomplete genetic divergence among phenotypically differentiated montane forest populations of Atlapetesalbinucha (Aves, Passerellidae)

The White-naped Brushfinch (Atlapetesalbinucha) comprises up to eight allopatric subspecies mainly identified by the color of the underparts (gray vs. yellow belly). Yellow and gray bellied forms were long considered two different species (A.albinucha and A.gutturalis), but they are presently considered as one polytypic species. Previous studies in the genus Atlapetes have shown that the phylogeny, based on molecular data, is not congruent with characters such as coloration, ecology, or distributional patterns. The phylogeography of A.albinucha was analyzed using two mitochondrial DNA regions from samples including 24 different localities throughout montane areas from eastern Mexico to Colombia. Phylogeographic analyses using Bayesian inference, maximum likelihood and haplotype network revealed incomplete geographic structure. The genetic diversity pattern is congruent with a recent process of expansion, which is also supported by Ecological Niche Models (ENM) constructed for the species and projected into three past scenarios. Overall, the results revealed an incomplete genetic divergence among populations of A.albinucha in spite of the species’ ample range, which contrasts with previous results of phylogeographic patterns in other Neotropical montane forest bird species, suggesting idiosyncratic evolutionary histories for different taxa throughout the region.