Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Nearctic-Neotropical transition

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

Forest cover and geographical distance influence fine-scale genetic structure of leaf-toed geckos in the tropical dry forests of western Mexico

The biodiversity within tropical dry forests (TDFs) is astounding and yet poorly catalogued due to inadequate sampling and the presence of cryptic species. In the Mexican TDF, endemic species are common, and the landscape has been continually altered by geological and anthropogenic changes. To understand how landscape and environmental variables have shaped the population structure of endemic species, we studied the recently described species of leaf-toed gecko, Phyllodactylus benedettii, in coastal western Mexico. Using double-digest restriction site-associated DNA sequencing data, we first explore population structure and estimate the number of ancestral populations. The results indicate a high degree of genetic structure with little admixture, and patterns corresponding to both latitudinal and altitudinal gradients. We find that genetic structure cannot be explained purely by geographical distance, and that ecological corridors may facilitate dispersal and gene flow. We then model the spatial distribution of P. benedettii in the TDF through time and find that the coastline has been climatically suitable for the species since the Last Glacial Maximum. Landscape genetic analyses suggest that the combination of isolation by distance (IBD) and isolation by resistance (IBR; forest cover) has influenced the spatial genetic structure of the species. Overall, our genomic data demonstrate fine-scale population structure in TDF habitat, a complex colonization history, and spatial patterns consistent with both IBD and other ecological factors. These results further highlight the Mexican TDF as a diversity hotspot and suggest that continued anthropogenic changes are likely to affect native fauna.

A tale about vipers’ tails: phylogeography of black-tailed rattlesnakes

The phylogenetic relationships among black-tailed rattlesnakes remain poorly understood and some authors indicated that the diversity of this group has been underestimated and additional analyses are required to clarify the biogeographic patterns throughout its distribution in Mexico. Therefore, the aim of this study was to elucidate the phylogenetic relationships among black-tailed rattlesnakes across their range, identifying relative divergence times among the main clades and reconstructing the biogeographical history of the group. Three partial mitochondrial genes (ND4, cytb and ATPase6) and one nuclear gene (RAG1) were sequenced to infer the phylogenetic relationships, through the maximum likelihood and Bayesian inference-based methods; demographic history reconstruction was investigated through Bayesian Skyline plot analysis and the ancestral area reconstruction was carried out considering a Bayesian framework. We found strong evidence that the black-tailed rattlesnakes’ group is composed of six clades, which is in agreement with subspecies previously reported. Divergence time estimation indicated that the origin of the C. molossus group could be traced to the middle of the Miocene (~7.71 Mya). Ancestral area reconstruction indicated that early divergence events occurred in Central Mexico, probably related to the geological dynamics of the Trans-Mexican Volcanic Belt. The lineage C. m. oaxacus is the basal member of the C. molossus group. Furthermore, the combination of geological events and changes in Quaternary vegetation may have contributed to the divergence of C. molossus clades. Our results suggest several clades within C. molossus complex could be potentially recognized as separate species.

Phylogenomics of alligator lizards elucidate diversification patterns across the Mexican Transition Zone and support the recognition of a new genus

Genomic data continue to advance our understanding of species limits and biogeographic patterns. However, there is still no consensus regarding appropriate methods of phylogenomic analysis that make the best use of these heterogeneous data sets. In this study, we used thousands of ultraconserved element (UCE) loci from alligator lizards in the genus Gerrhonotus to compare and contrast species trees inferred using multiple contemporary methods and provide a time frame for biological diversification across the Mexican Transition Zone (MTZ). Concatenated maximum likelihood (ML) and Bayesian analyses provided highly congruent results, with differences limited to poorly supported nodes. Similar topologies were inferred from coalescent analyses in Bayesian Phylogenetics and Phylogeography and SVDquartets, albeit with lower support for some nodes. All divergence times fell within the Miocene, linking speciation to local Neogene vicariance and/or global cooling trends following the mid-Miocene Climatic Optimum. We detected a high level of genomic divergence for a morphologically distinct species restricted to the arid mountains of north-eastern Mexico, and erected a new genus to better reflect evolutionary history. In summary, our results further advocate leveraging the strengths and weaknesses of concatenation and coalescent methods, provide evidence for old divergences for alligator lizards, and indicate that the MTZ continues to harbour substantial unrecognized diversity.

Historical biogeography of the genus Rhadinaea (Squamata: Dipsadinae)

Multiple geological and climatic events have created geographical or ecological barriers associated with speciation events, playing a role in biological diversification in North and Central America. Here, we evaluate the influence of the Neogene and Quaternary geological events, as well as the climatic changes in the diversification of the colubrid snake genus Rhadinaea using molecular dating and ancestral area reconstruction. A multilocus sequence dataset was generated for 37 individuals of Rhadinaea from most of the biogeographical provinces where the genus is distributed, representing 19 of the 21 currently recognized species, and two undescribed species. Our analyses show that the majority of the Rhadinaea species nest in two main clades, herein identified as "Eastern" and "Southern". These clades probably diverged from each other in the early Miocene, and their divergence was followed by 11 divergences during the middle to late Miocene, three divergences during the Pliocene, and six divergences in the Pleistocene. The ancestral distribution of Rhadinaea was reconstructed across the Sierra Madre del Sur. Our phylogenetic analyses do not support the monophyly of Rhadinaea. The Miocene and Pliocene geomorphology, perhaps in conjunction with climate change, appears to have triggered the diversification of the genus, while the climatic changes during the Miocene probably induced the diversification of Rhadinaea in the Sierra Madre del Sur. Our analysis suggests that the uplifting of the Trans-Mexican Volcanic Belt and Chiapan-Guatemalan highlands in this same period resulted in northward and southward colonization events. This was followed by more recent, independent colonization events in the Pliocene and Pleistocene involving the Balsas Basin, Chihuahuan Desert, Pacific Coast, Sierra Madre Occidental, Sierra Madre Oriental, Sierra Madre del Sur, Trans-Mexican Volcanic Belt, and Veracruz provinces, probably driven by the climatic fluctuations of the time.

Insights into the importance of areas of climatic stability in the evolution and maintenance of avian diversity in the Mesoamerican dry forests

We analysed the phylogeographic structure of five resident bird lineages distributed in the seasonally dry tropical forests (SDTF) of Mesoamerica to test whether they show patterns of synchronous and geographically coincident genetic divergence during the Quaternary. We generated phylogenetic trees, estimated divergence times and analysed the genetic structure of populations (based on sequences of mitochondrial genes), as well as estimating historical distributions (range extension and areas of long-term climate stability) during the Late Pleistocene. We tested and selected the phylogeographic divergence scenarios that best explain the current divergence patterns of taxa using the Approximate Bayesian Computation (ABC) approach. For most species, phylogenetic trees and haplotype networks showed a clear genetic structure associated with geographical distribution. Overall, the divergence times ranged from 0.29–2.0 Mya, suggesting that diversification of populations occurred at different times during the Pleistocene. The palaeodistribution models predicted at least two areas of climatic stability within the current SDTF that probably allowed glacial-interglacial persistence of isolated bird populations along the Mexican Pacific, thus promoting their genetic divergence. The results provide information relevant to the identification of diversification hotspots for the Mesoamerican SDTF avifauna.

Introgression obscures lineage boundaries and phylogeographic history in the western banded gecko, Coleonyx variegatus (Squamata: Eublepharidae)

The geomorphological formation of the Baja California peninsula and the Gulf of California is a principal driver of diversification for the reptiles of North America’s warm deserts. The western banded gecko, Coleonyx variegatus, is distributed throughout the Mojave, Sonoran and Peninsular deserts. In this study we use multilocus sequence data to address deep phylogeographic structure within C. variegatus. Analyses of mtDNA data recover six divergent clades throughout the range of C. variegatus. Topology of the mtDNA gene tree suggests separate origins of peninsular populations with an older lineage in the south and a younger one in the north. In contrast, analyses of multilocus nuclear data provide support for four lineages, corresponding to the subspecies C. v. abbotti, C. v. peninsularis, C. v. sonoriensis and C. v. variegatus. Phylogenetic analyses of the nuclear data recover C. v. abbotti and C. v. peninsularis as a clade, indicating a single origin of the peninsular populations. Discordance between the nuclear and mtDNA data is largely the result of repeated episodes of mtDNA introgression that have obscured both lineage boundaries and biogeographic history. Dating analyses of the combined nuclear and mtDNA data suggest that the peninsular clade diverged from the continental group in the Late Miocene.

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.

What drives genetic and phenotypic divergence in the Red-crowned Ant tanager (Habia rubica, Aves: Cardinalidae), a polytypic species?

AIM

The effects of geographic and environmental variables on patterns of genetic and phenotypic differentiation have been thoroughly studied. Ecological speciation involves reproductive isolation due to divergent natural selection that can result in a positive correlation between genetic divergence and adaptive phenotypic divergence (isolation by adaptation, IBA). If the phenotypic target of selection is unknown or not easily measured, environmental variation can be used as a proxy, expecting positive correlation between genetic and environmental distances, independent of geographic distances (isolation by environment, IBE). The null model is that the amount of gene flow between populations decreases as the geographic distance between them increases, and genetic divergence is due simply to the neutral effects of genetic drift (isolation by distance, IBD). However, since phenotypic differentiation in natural populations may be autocorrelated with geographic distance, it is often difficult to distinguish IBA from the neutral expectation of IBD. In this work, we test hypotheses of IBA, IBE, and IBD in the Red-crowned Ant tanager (Habia rubica).

LOCATION

Mesoamerica (Mexico-Central America) and South America.

TAXON

Habia rubica (Aves: Cardinalidae).

METHODS

We compiled genetic data, coloration, and morphometric data from specimens from collections in Mexico and the United States. We used the Multiple Matrix Regression with Randomization (MMRR) approach to evaluate the influence of geographic and environmental distances on genetic and phenotypic differentiation of H. rubica at both phylogroup and population levels.

RESULTS

Our results provide strong evidence that geographic distance is the main driver of genetic variation in H. rubica. We did not find evidence that climate variation is driving population differentiation in this species across a widespread geographic region.

MAIN CONCLUSIONS

Our data point to geographic isolation as the main factor structuring genetic variation within populations of H. rubica and suggest that climate is not playing a major role in genetic differentiation within this species.

Genetic differentiation in the southern population of the Fathead Minnow Pimephales promelas Rafinesque (Actinopterygii: Cyprinidae)

The North American cyprinid Pimephales promelas is a species with a wide distribution range, occurring in distinct hydrographic basins in Mexico, Canada, and the United States. Previous morphological and meristic analyses of P. promelas concluded that at least three subspecies exist in the midwestern and northeast region of the United States. No studies have been carried out on the Mexican population of P. promelas, but the findings of cryptic diversity in United States populations of this species, as well as in other codistributed fish species in Mexico could be an indication that Mexican populations of P. promelas consist of cryptic species. Using the mitochondrial gene cyt b and the first intron of the S7 ribosomal protein-coding nuclear gene we carried out phylogenetic and phylogeographic analyses of populations of P. promelas across its distribution range in northwestern Mexico. Using this information were analyzed the structure and differentiation level between populations of P. promelas from distinct river basins in the region in identifying cryptic diversity. Twenty-four sequences were obtained for cyt b, and 30 for S7, which included the two heterozygous alleles. The results revealed the existence of four well-differentiated lineages: (1) Yaqui in the Pacific slope; (2) Santa Maria, and (3) Casas Grandes in the Guzman Basin; and (4) Nazas+Conchos in Chihuahua state. This challenges the current taxonomy of P. promelas. Differences in the relationships between markers and the small sample size for the Santa Maria population (n = 1), indicate that our results must be corroborated with more data and morphological analyses. Biogeographic analysis of these findings suggest that the evolutionary history of P. promelas is associated with the fragmentation of the ancestral Rio Grande river system since Miocene in northwestern Mexico consistent with findings for codistributed fish species.

Introgressive hybridization in a Spiny-Tailed Iguana, Ctenosaura pectinata, and its implications for taxonomy and conservation

Introgression, the transmission of genetic material of one taxon into another through hybridization, can have various evolutionary outcomes. Previous studies have detected signs of introgression between western populations of the Mexican endemic and threatened spiny-tailed iguana, Ctenosaura pectinata. However, the extent of this phenomenon along the geographic distribution of the species is unknown. Here, we use multilocus data together with detailed geographic sampling to (1) define genotypic clusters within C. pectinata; (2) evaluate geographic concordance between maternally and biparentally inherited markers; (3) examine levels of introgression between genotypic clusters, and (4) suggest taxonomic modifications in light of this information. Applying clustering methods to genotypes of 341 individuals from 49 localities of C. pectinata and the closely related C. acanthura, we inferred the existence of five genotypic clusters. Contact zones between genotypic clusters with signatures of interbreeding were detected, showing different levels of geographic discordance with mtDNA lineages. In northern localities, mtDNA and microsatellites exhibit concordant distributions, supporting the resurrection of C. brachylopha. Similar concordance is observed along the distribution of C. acanthura, confirming its unique taxonomic identity. Genetic and geographic concordance is also observed for populations within southwestern Mexico, where the recognition of a new species awaits in depth taxonomic revision. In contrast, in western localities a striking pattern of discordance was detected where up to six mtDNA lineages co-occur with only two genotypic clusters. Given that the type specimen originated from this area, we suggest that individuals from western Mexico keep the name C. pectinata. Our results have profound implications for conservation, management, and forensics of Mexican iguanas.

Tracking Host Trees: The Phylogeography of Endemic Psittacanthus sonorae (Loranthaceae) Mistletoe in the Sonoran Desert

The host dependence of mistletoes suggests that they track the distributions of their hosts. However, the factors that determine the geographic distribution of mistletoes are not well understood. In this study, the phylogeography of Psittacanthus sonorae was reconstructed by sequencing one nuclear (ITS) and two plastid (trnL-F and atpB-rbcL) regions of 148 plants from populations separated by the Sea of Cortez. Divergence time and gene flow were estimated to gain insight into the historical demography and geographic structuring of genetic variation. We also described and mapped the spatial distribution of suitable habitat occupied by P. sonorae and its most common host Bursera microphylla in the Sonoran Desert, along with their responses to Quaternary climate fluctuations using environmental data and ecological niche modeling (ENM). We detected environmental and genetic differentiation between the peninsular and continental P. sonorae populations. Population divergence occurred during the Pleistocene, around the time of the Last Glacial Maximum. No signals of population growth were detected, with net gene flow moving from the continent to the peninsula. ENM models indicate decoupled responses by the mistletoe and its main host to past climate changes. For the Last Interglacial to the present, most models produce only partial areas of overlap on both the peninsula and the continent. Our results support a scenario of Late-Pleistocene isolation and divergence with asymmetrical gene flow between peninsular and continental P. sonorae populations. Continental populations migrated to the peninsula and the spatial isolation probably produced genetic differentiation under different environmental conditions.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Molecular systematics and phylogeography of the endemic Osgood's deermouse Osgoodomys banderanus (Rodentia: Cricetidae) in the lowlands of western Mexico

Osgoodomys banderanus is a recognized and endemic rodent species of western Mexico, an area known for its high biodiversity and number of endemisms. Phylogeographical relationships within this taxon were analyzed based on mitochondrial (ND3, tRNA-Arginine, ND4L and partial ND4) and nuclear (GHR) nucleotide sequences. We obtained a total of 112 samples from 22 localities, covering the complete distribution of the species. Phylogenetic analyses using Maximum Likelihood and Bayesian inference confirmed that Osgoodomys is a monophyletic group. In addition, phylogenetic and phylogeographic analyses detected three major clades, which do not coincide with the recognized subspecies of O. banderanus. The genetic lineages detected are the western clade (Nayarit, Jalisco and northern Colima), the central clade (Colima, Michoacán, and northern Guerrero) and the eastern clade (central and southern Guerrero). Genetic distances among clades (5-9%) and nucleotide substitutions (30-88) among haplogroups were high, especially in the southern group. Mountain ranges such as the Transmexican Volcanic Belt and the Sierra Madre del Sur, as well as the Balsas River act as geographical barriers for Osgoodomys. Our results suggest the presence of three independent species, which need to be characterized morphologically to confirm our hypothesis.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Phylogeographic assessment of the northern pygmy mouse, Baiomys taylori

The northern pygmy mouse, Baiomys taylori, occurs throughout the Trans-Mexican Volcanic Belt and southern Altiplano of central Mexico and extends northward in 3 projections into northern Mexico and the United States. We used mitochondrial DNA (from the cytochrome-b and NADH dehydrogenase 2 genes) and morphological data to assess diversity within B. taylori across its geographic range in relation to recognized subspecies and putative physiographic filter-barriers. Our results indicate 5 distinct mitochondrial clades despite little morphological variation across the species’ geographic range. The Sierra Madre Oriental and Sierra Madre Occidental separate clades representing the eastern (Gulf coastal lowlands), central (Altiplano), and Pacific coastal lowlands, which appear to be divided into 3 major mitochondrial clades. Based on a preliminary analysis of cranial morphology, we are able to reject the Balcones Escarpment of Texas as an impediment to the well-documented recent northern expansion of the species, while we are unable to reject a causal role of filter-barriers elsewhere in subspecific differentiation. Revision of subspecific taxonomy must await further genetic sampling, particularly along the western and southeastern portions of the species’ distribution.

El ratón pigmeo norteño, Baiomys taylori, se distribuye a través del Eje Volcánico Transmexicano y la parte sur del Altiplano del centro de México y se extiende hacia el norte de México y los Estados Unidos en tres direcciones. Utilizamos ADN mitocondrial (del citocromo-b y del gen NADH deshidrogenasa) y datos morfológicos para evaluar la diversidad dentro de B. taylori a través de su rango geográfico en relación con las subespecies reconocidas y las supuestas barreras-filtro fisiográficas. Nuestros resultados indican 5 clados mitocondriales distintos a pesar de la escasa variación morfológica a lo largo de su rango geográfico. La Sierra Madre Oriental y la Sierra Madre Occidental separan los clados representantes de las tierras bajas del este (tierras bajas del Golfo), del centro (Altiplano), y las tierras bajas del Pacífico, que parecen estar divididas en 3 clados mitocondriales principales. Basados en un análisis preliminar de la morfología craneal, pudimos rechazar el Escarpe de Balcones de Texas como un impedimento a la bien documentada expansión norteña de la especie, mientras que no fue posible rechazar el rol causal de las otras barreras-filtro en la diferenciación subespecífica. La revisión de la taxonomía subespecífica deberá esperar un mayor muestreo genético, en particular a lo largo de las porciones oeste y sureste de la distribución de la especie.

Phylogeography of Pogonomyrmex barbatus and P. rugosus harvester ants with genetic and environmental caste determination

We present a phylogeographic study of at least six reproductively isolated lineages of new world harvester ants within the Pogonomyrmex barbatus and P. rugosus species group. The genetic and geographic relationships within this clade are complex: Four of the identified lineages show genetic caste determination (GCD) and are divided into two pairs. Each pair has evolved under a mutualistic system that necessitates sympatry. These paired lineages are dependent upon one another because their GCD requires interlineage matings for the production of F1 hybrid workers, and intralineage matings are required to produce queens. This GCD system maintains genetic isolation among these interdependent lineages, while simultaneously requiring co-expansion and emigration as their distributions have changed over time. It has also been demonstrated that three of these four GCD lineages have undergone historical hybridization, but the narrower sampling range of previous studies has left questions on the hybrid parentage, breadth, and age of these groups. Thus, reconstructing the phylogenetic and geographic history of this group allows us to evaluate past insights and hypotheses and to plan future inquiries in a more complete historical biogeographic context. Using mitochondrial DNA sequences sampled across most of the morphospecies' ranges in the U.S.A. and Mexico, we conducted a detailed phylogeographic study. Remarkably, our results indicate that one of the GCD lineage pairs has experienced a dramatic range expansion, despite the genetic load and fitness costs of the GCD system. Our analyses also reveal a complex pattern of vicariance and dispersal in Pogonomyrmex harvester ants that is largely concordant with models of late Miocene, Pliocene, and Pleistocene range shifts among various arid-adapted taxa in North America.

Independent Demographic Responses to Climate Change among Temperate and Tropical Milksnakes (Colubridae: Genus Lampropeltis)

The effects of Late Quaternary climate change have been examined for many temperate New World taxa, but the impact of Pleistocene glacial cycles on Neotropical taxa is less well understood, specifically with respect to changes in population demography. Here, we examine historical demographic trends for six species of milksnake with representatives in both the temperate and tropical Americas to determine if species share responses to climate change as a taxon or by area (i.e., temperate versus tropical environments). Using a multilocus dataset, we test for the demographic signature of population expansion and decline using non-genealogical summary statistics, as well as coalescent-based methods. In addition, we determine whether range sizes are correlated with effective population sizes for milksnakes. Results indicate that there are no identifiable trends with respect to demographic response based on location, and that species responded to changing climates independently, with tropical taxa showing greater instability. There is also no correlation between range size and effective population size, with the largest population size belonging to the species with the smallest geographic distribution. Our study highlights the importance of not generalizing the demographic histories of taxa by region and further illustrates that the New World tropics may not have been a stable refuge during the Pleistocene.

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito-nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation-migration analyses, however, revealed near-zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.

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

Western North America and Mexico contain a large number of conifer species. This diversity could be the product of orographic and climate changes of the late Tertiary and Quaternary. In this study, we focus on the evolutionary history of Juniperus blancoi, in order to determine the impact of climate change and environmental heterogeneity on population differentiation. We estimated the population structure, phylogenetic relationships and historical demography of J. blancoi populations using nuclear genes. We correlated genetic structure with ecological differentiation, divergence times and changes in population size. Populations of J. blancoi are differentiated into three lineages that correspond to low-, mid- and high-altitude populations. The three groups diversified in the late Miocene, early Pliocene, with only a few events of gene flow since then. Two lineages in the north exhibited a pattern of population growth during the Pleistocene that could be linked to climate changes. Populations of J. blancoi experienced significant ecological differentiation and early divergence events, which correspond to periods of global cooling and mountain uplift during the Miocene. This suggests that mountain ranges in tropical and subtropical latitudes play an important role in the speciation and persistence of conifer taxa in diversity hotspots, by providing diverse environmental conditions.

Genetic differentiation and species cohesion in two widespread Central American Begonia species

Begonia is one of the ten largest plant genera, with over 1500 species. This high species richness may in part be explained by weak species cohesion, which has allowed speciation by divergence in allopatry. In this study, we investigate species cohesion in the widespread Central American Begonia heracleifolia and Begonia nelumbiifolia, by genotyping populations at microsatellite loci. We then test for post-zygotic reproductive barriers using experimental crosses, and assess whether sterility barriers are related to intraspecific changes in genome size, indicating major genome restructuring between isolated populations. Strong population substructure was found for B. heracleifolia (FST=0.364, F'ST=0.506) and B. nelumbiifolia (FST=0.277, F'ST=0.439), and Bayesian admixture analysis supports the division of most populations into discrete genetic clusters. Moderate levels of inferred selfing (B. heracleifolia s=0.40, B. nelumbiifolia s=0.62) and dispersal limitation are likely to have contributed to significant genetic differentiation (B. heracleifolia Jost's D=0.274; B. nelumbiifolia D=0.294). Interpopulation crosses involving a divergent B. heracleifolia population with a genome size ∼10% larger than the species mean had a ∼20% reduction in pollen viability compared with other outcrosses, supporting reproductive isolation being polymorphic within the species. The population genetic data suggest that Begonia populations are only weakly connected by gene flow, allowing reproductive barriers to accumulate between the most isolated populations. This supports allopatric divergence in situ being the precursor of speciation in Begonia, and may also be a common speciation mechanism in other tropical herbaceous plant groups.

Evolutionary diversification and speciation in rodents of the Mexican lowlands: the Peromyscus melanophrys species group

Despite some studies of the species groups within the genus Peromyscus have been performed, both evolutionary relationships among species within groups and group composition have remained controversial. In this study, we address phylogenetic relationships among species in the Peromyscus melanophrys group (P. melanophrys, P. perfulvus, and P. mekisturus), using a molecular phylogenetic analysis. This analysis is the first to include the poorly known P. mekisturus. We conducted maximum likelihood and Bayesian inference analyses with the ND3, tRNA-Arginine, ND4L, and partial ND4 mitochondrial genes, and the GHR nuclear gene. We consistently recovered a P. melanophrys group that is monophyletic with respect to the set of outgroups. Also, we recovered two distinct clades within P. perfulvus and two within P. melanophrys, one of which contain P. mekisturus among other P. melanophrys, all with geographic consistency. According to our divergence time estimates, the P. melanophrys group diverged during the Pliocene and the main diversification events within the group occurred at the end of the Pliocene and through the Pleistocene.

Phylogeographic pattern and extensive mitochondrial DNA divergence disclose a species complex within the Chagas disease vector Triatoma dimidiata

BACKGROUND

Triatoma dimidiata is among the main vectors of Chagas disease in Latin America. However, and despite important advances, there is no consensus about the taxonomic status of phenotypically divergent T. dimidiata populations, which in most recent papers are regarded as subspecies.

METHODOLOGY AND FINDINGS

A total of 126 cyt b sequences (621 bp long) were produced for specimens from across the species range. Forty-seven selected specimens representing the main cyt b clades observed (after a preliminary phylogenetic analysis) were also sequenced for an ND4 fragment (554 bp long) and concatenated with their respective cyt b sequences to produce a combined data set totalling 1175 bp/individual. Bayesian and Maximum-Likelihood phylogenetic analyses of both data sets (cyt b, and cyt b+ND4) disclosed four strongly divergent (all pairwise Kimura 2-parameter distances >0.08), monophyletic groups: Group I occurs from Southern Mexico through Central America into Colombia, with Ecuadorian specimens resembling Nicaraguan material; Group II includes samples from Western-Southwestern Mexico; Group III comprises specimens from the Yucatán peninsula; and Group IV consists of sylvatic samples from Belize. The closely-related, yet formally recognized species T. hegneri from the island of Cozumel falls within the divergence range of the T. dimidiata populations studied.

CONCLUSIONS

We propose that Groups I-IV, as well as T. hegneri, should be regarded as separate species. In the Petén of Guatemala, representatives of Groups I, II, and III occur in sympatry; the absence of haplotypes with intermediate genetic distances, as shown by multimodal mismatch distribution plots, clearly indicates that reproductive barriers actively promote within-group cohesion. Some sylvatic specimens from Belize belong to a different species - likely the basal lineage of the T. dimidiata complex, originated ~8.25 Mya. The evidence presented here strongly supports the proposition that T. dimidiata is a complex of five cryptic species (Groups I-IV plus T. hegneri) that play different roles as vectors of Chagas disease in the region.

Phylogenetic relationships of the enigmatic longtailed rattlesnakes (Crotalus ericsmithi, C. lannomi, and C. stejnegeri)

The longtailed rattlesnakes of western Mexico represent an enigmatic group of poorly known venomous snake species: Crotalus ericsmithi, C. lannomi, and C. stejnegeri. In the 120 years since their discovery, fewer than twenty individuals have been deposited in natural history collections worldwide. These three species share similar morphological traits, including a particularly long tail that has been interpreted as either an ancestral condition among rattlesnakes or as derived within the longtailed group. An understanding of the phylogenetic distinctiveness and relationships among the longtailed rattlesnakes, and their relationships to other rattlesnake groups, has previously been hampered by a dearth of comparative material and tissues for collection of DNA sequence data. Facilitated by the recent availability of tissue samples from multiple individuals of each species, we estimate the phylogenetic relationships among the longtailed rattlesnakes and their placement among other rattlesnake groups, using DNA sequence data from three mitochondrial and three nuclear gene fragments. We explore phylogenetic signal in our data using Bayesian and maximum likelihood methods, species tree analyses and hypothesis testing. Our results strongly support the monophyly of longtailed rattlesnakes and suggest the three species diverged from each other during the mid to late Pliocene or early Pleistocene (~1.5-5.6 mya). Contrary to prevailing hypotheses, we find no evidence for an early or basal divergence of the longtailed clade within the rattlesnake tree, and instead estimate that it diverged relatively recently (~6.8 mya) from its sister lineage, composed of the diamondback rattlesnakes (C. atrox group) and the prairie rattlesnakes (C. viridis group). With our added sampling of lineages and identification of previously used problematic sequences, we provide a revised hypothesis for relationships among Crotalus species, yet underscore the need for future studies and new data to deliver a well-supported robust estimate of rattlesnake relationships.

A new species of Metacyclops Kiefer, 1927 (Copepoda, Cyclopidae, Cyclopinae) from the Chihuahuan desert, northern Mexico

A new species of the freshwater cyclopoid copepod genus Metacyclops Kiefer, 1927 is described from a single pond in northern Mexico, within the binational area known as the Chihuahuan Desert. This species belongs to a group of Metacyclops species with a 3443 spine formula of swimming legs. It is morphologically similar to Metacyclops lusitanus Lindberg, 1961 but differs from this and other congeners by having a unique combination of characters, including a caudal rami length/width proportion of 3.5–3.8, a innermost terminal seta slightly longer than the outermost terminal seta, intercoxal sclerites of legs 1-4 naked, a strong apical spine of the second endopodal segment of leg 1 and one row of 6-8 small spinules at the insertion of this spine. The finding of this species represents also the first record of the genus in Mexico and the third in North America, where only two other species, Metacyclops gracilis (Lilljeborg, 1853)and Metacyclops cushae Reid, 1991 have been hitherto reported. This is also the first continental record of a species of Metacyclops from an arid environment in the Americas. This species appears to be endemic to the Chihuahuan Desert, thus emphasizing the high endemicity of this area.

Comparative phylogeographic analyses illustrate the complex evolutionary history of threatened cloud forests of northern Mesoamerica

Comparative phylogeography can elucidate the influence of historical events on current patterns of biodiversity and can identify patterns of co-vicariance among unrelated taxa that span the same geographic areas. Here we analyze temporal and spatial divergence patterns of cloud forest plant and animal species and relate them to the evolutionary history of naturally fragmented cloud forests–among the most threatened vegetation types in northern Mesoamerica. We used comparative phylogeographic analyses to identify patterns of co-vicariance in taxa that share geographic ranges across cloud forest habitats and to elucidate the influence of historical events on current patterns of biodiversity. We document temporal and spatial genetic divergence of 15 species (including seed plants, birds and rodents), and relate them to the evolutionary history of the naturally fragmented cloud forests. We used fossil-calibrated genealogies, coalescent-based divergence time inference, and estimates of gene flow to assess the permeability of putative barriers to gene flow. We also used the hierarchical Approximate Bayesian Computation (HABC) method implemented in the program msBayes to test simultaneous versus non-simultaneous divergence of the cloud forest lineages. Our results show shared phylogeographic breaks that correspond to the Isthmus of Tehuantepec, Los Tuxtlas, and the Chiapas Central Depression, with the Isthmus representing the most frequently shared break among taxa. However, dating analyses suggest that the phylogeographic breaks corresponding to the Isthmus occurred at different times in different taxa. Current divergence patterns are therefore consistent with the hypothesis of broad vicariance across the Isthmus of Tehuantepec derived from different mechanisms operating at different times. This study, coupled with existing data on divergence cloud forest species, indicates that the evolutionary history of contemporary cloud forest lineages is complex and often lineage-specific, and thus difficult to capture in a simple conservation strategy.

Deep divergences and extensive phylogeographic structure in a clade of lowland tropical salamanders

BACKGROUND

The complex geological history of Mesoamerica provides the opportunity to study the impact of multiple biogeographic barriers on population differentiation. We examine phylogeographic patterns in a clade of lowland salamanders (Bolitoglossa subgenus Nanotriton) using two mitochondrial genes and one nuclear gene. We use several phylogeographic analyses to infer the history of this clade and test hypotheses regarding the geographic origin of species and location of genetic breaks within species. We compare our results to those for other taxa to determine if historical events impacted different species in a similar manner.

RESULTS

Deep genetic divergence between species indicates that they are relatively old, and two of the three widespread species show strong phylogeographic structure. Comparison of mtDNA and nuclear gene trees shows no evidence of hybridization or introgression between species. Isolated populations of Bolitoglossa rufescens from Los Tuxtlas region constitute a separate lineage based on molecular data and morphology, and divergence between Los Tuxtlas and other areas appears to predate the arrival of B. rufescens in other areas west of the Isthmus of Tehuantepec. The Isthmus appears responsible for Pliocene vicariance within B. rufescens, as has been shown for other taxa. The Motagua-Polochic fault system does not appear to have caused population vicariance, unlike in other systems.

CONCLUSIONS

Species of Nanotriton have responded to some major geological events in the same manner as other taxa, particularly in the case of the Isthmus of Tehuantepec. The deep divergence of the Los Tuxtlas populations of B. rufescens from other populations highlights the contribution of this volcanic system to patterns of regional endemism, and morphological differences observed in the Los Tuxtlas populations suggests that they may represent an undescribed species of Bolitoglossa. The absence of phylogeographic structure in B. nympha, in contrast to the other widespread species in the subgenus, may be due to historical forest contraction and more recent range expansion in the region. Phylogeographic data provide substantial insight into the evolutionary history of these morphologically similar species of salamanders, and contribute to our understanding of factors that have generated the high biodiversity of Mesoamerica.