Cryptic Neogene vicariance and Quaternary dispersal of the red-spotted toad (Bufo punctatus): insights on the evolution of North American warm desert biotas

Population genetic structure of a Chihuahuan Desert endemic mammal, the desert pocket gopher, Geomys arenarius

The biogeographic history of the Chihuahuan Desert is complex, driven by numerous physiographic events and climatic changes. This dynamic history would have influenced the flora and fauna of the region including the desert pocket gopher, Geomys arenarius, a subterranean rodent endemic to the northern Chihuahuan Desert. G. arenarius is restricted to sandy soils and are considered to have a disjunct distribution. Two subspecies are recognized: G. a. arenarius and G. a. brevirostris. We used multilocus nuclear (amplified fragment length polymorphisms) and mitochondrial DNA (ND2) sequence data to uncover patterns of genetic diversity within and among populations of G. arenarius. We evaluated correspondence of genetic patterns to traditionally accepted subspecies boundaries, mapped the distribution of potentially suitable soils to identify barriers or corridors to dispersal and to guide future survey efforts, provided evidence that could be used to recognize distinct population segments, and quantified genetic diversity within populations. Both datasets were largely concordant and demonstrated hierarchical patterns of genetic divergence. The greatest divergence was consistent with the two recognized subspecies. Mapping of potentially habitable soils revealed likely barriers to dispersal contributing to the allopatric pattern of geographic distribution and areas, which may be occupied by G. arenarius but not yet documented. Because G. arenarius is restricted to soils with high sand content, and these habitable soils are disjunct within the region occupied by this species, historical factors that impacted soil deposition and deflation likely contributed to the observed patterns of genetic divergence. Genetic diversity was higher within populations of the southern subspecies (G. a. arenarius) compared to G. a. brevirostris. This may be due to a greater availability of continuous suitable soils within the range of G. a. arenarius or higher density due to greater food availability (currently or historically)-both of which could allow for a higher effective population size.

Evolutionary and Ecological Drivers of Local Adaptation and Speciation in a North American Avian Species Complex

Throughout the speciation process, genomic divergence can be differentially impacted by selective pressures, as well as gene flow and genetic drift. Disentangling the effects of these evolutionary mechanisms remains challenging, especially for non-model organisms. Accounting for complex evolutionary histories and contemporary population structure often requires sufficient sample sizes, for which the expense of full genomes remains prohibitive. Here, we demonstrate the utility of partial-genome sequence data for range-wide samples to shed light into the divergence process of two closely related ducks, the Mexican duck (Anas diazi) and mallard (A. platyrhynchos). We determine the role of selective and neutral processes during speciation of Mexican ducks by integrating evolutionary and demographic modelling with genotype-environment and genotype-phenotype association testing. First, evolutionary models and demographic analyses support the hypothesis that Mexican ducks originally diverged ~300,000 years ago in a climate refugia arising during a glacial period in in a southwestern North America, and that subsequent environmental selective pressures played a key role in divergence. Mexican ducks then showed cyclical demographic patterns that likely reflected repeated range expansions and contractions, along with bouts of gene flow with mallards during glacial cycles. Finally, we provide evidence that sexual selection acted on several phenotypic traits as a co-evolutionary process, facilitating the development of reproductive barriers that initially arose due to strong ecological selection. More broadly, this work reveals that the genomic and phenotypic patterns observed across species complexes are the result of myriad factors that contribute in dynamic ways to the evolutionary trajectories of a lineage.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

Diversification and post-glacial range expansion of giant North American camel spiders in genus Eremocosta (Solifugae: Eremobatidae)

Species of camel spiders in the family Eremobatidae are an important component of arthropod communities in arid ecosystems throughout North America. Recently, research demonstrated that the evolutionary history and biogeography of the family are poorly understood. Herein we explore the biogeographic history of this group of arachnids using genome-wide single nucleotide polymorphism (SNP) data, morphology, and distribution modelling to study the eremobatid genus Eremocosta, which contains exceptionally large species distributed throughout North American deserts. Relationships among sampled species were resolved with strong support and they appear to have diversified within distinct desert regions along an east-to-west progression beginning in the Chihuahuan Desert. The unexpected phylogenetic position of some samples suggests that the genus may contain additional, morphologically cryptic species. Geometric morphometric analyses reveal a largely conserved cheliceral morphology among Eremocosta spp. Phylogeographic analyses indicate that the distribution of E. titania was substantially reduced during the last glacial maximum and the species only recently colonized much of the Mojave Desert. Results from this study underscore the power of genome-wide data for unlocking the genetic potential of museum specimens, which is especially promising for organisms like camel spiders that are notoriously difficult to collect.

Impact of the Mid-Pleistocene Revolution and Anthropogenic Factors on the Dispersion of Asian Black-Spined Toads (Duttaphrynus melanostictus)

Divergence-time estimation critically improves the understanding of biogeography processes underlying the distribution of species, especially when fossil data is not available. We hypothesise that the Asian black-spined toad, Duttaphrynus melanostictus, expanded into the Eastern Indomalaya following the Quaternary glaciations with the subsequent colonisation of new landscapes during the Last Glacial Maximum. Divergence dating inferred from 364 sequences of mitochondrial tRNAGly ND3 supported the emergence of a common ancestor to the three D. melanostictus clades around 1.85 (±0.77) Ma, matching with the Lower to Mid-Pleistocene transition. Duttaphrynus melanostictus then dispersed into Southeast Asia from the central Indo-Pacific and became isolated in the Southern Sundaic and Wallacea regions 1.43 (±0.10) Ma through vicariance as a result of sea level oscillations. The clade on the Southeast Asian mainland then colonised the peninsula from Myanmar to Vietnam and expanded towards Southeastern China at the end of the Mid-Pleistocene Revolution 0.84 (±0.32) Ma. Population dynamics further highlight an expansion of the Southeast Asian mainland population towards Taiwan, the Northeastern edge of the species' range after the last interglacial, and during the emergence of the Holocene human settlements around 7000 BP. Thus, the current divergence of D. melanostictus into three segregated clades was mostly shaped by Quaternary glaciations, followed by natural dispersion events over land bridges and accelerated by anthropogenic activities.

Genomic divergence in allopatric Northern Cardinals of the North American warm deserts is linked to behavioral differentiation

Biogeographic barriers are considered important in initiating speciation through geographic isolation, but they rarely indiscriminately and completely reduce gene flow across entire communities. Explicitly demonstrating which factors are associated with gene-flow levels across barriers would help elucidate how speciation is initiated and isolation maintained. Here, we investigated the association of behavioral isolation on population differentiation in Northern Cardinals (Cardinalis cardinalis) distributed across the Cochise Filter Barrier, a region of transitional habitat which separates the Sonoran and Chihuahuan deserts of North America. Using genomewide markers, we modeled demographic history by fitting the data to isolation and isolation-with-migration models. The best-fit model indicated that desert populations diverged in the Pleistocene with low, historic, and asymmetric gene flow across the barrier. We then tested behavioral isolation using reciprocal call-broadcast experiments to compare song recognition between deserts, controlling for song dialect changes within deserts. We found that male Northern Cardinals in both deserts were most aggressive to local songs and failed to recognize across-barrier songs. A correlation of genomic differentiation and strong song discrimination is consistent with a model where speciation is initiated across a barrier and maintained by behavioral isolation.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

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.

Phylogeography of a widespread small carnivore, the western spotted skunk (Spilogale gracilis) reveals temporally variable signatures of isolation across western North America

We analyzed phylogeographic patterns in the western spotted skunk, Spilogale gracilis Merriam, 1890 (Carnivora: Mephitidae) in relation to historical events associated with Pre-Pleistocene Divergence (PPD) and Quaternary climate change (QCC) using mitochondrial DNA from 97 individuals distributed across Western North America. Divergence times were generated using BEAST to estimate when isolation in putative refugia occurred. Patterns and timing of demographic expansion was performed using Bayesian skyline plot. Putative climatic refugia resulting from Quaternary climate change were identified using paleoecological niche modeling and divergence dates compared to major vicariant events associated with Pre-Pleistocene conditions. We recovered three major mitochondrial clades corresponding to western North America (California, Baja, and across the Great Basin), east-central North America (Texas, central Mexico, New Mexico), and southwestern Arizona/northwestern Mexico. Time to most recent common ancestor for S. gracilis occurred ~1.36 Ma. Divergence times for each major clade occurred between 0.25 and 0.12 Ma, with signature of population expansion occurring 0.15 and 0.10 Ma. Ecological niche models identified three potential climatic refugia during the Last Interglacial, (1) west coast of California and Oregon, (2) northwestern Mexico, and (3) southern Texas/northeastern Mexico as well as two refugia during the Last Glacial Maximum, (1) western USA and (2) southern Texas/northeastern Mexico. This study supports PPD in shaping species-level diversity compared to QCC-driven changes at the intraspecific level for Spilogale, similar to the patterns reported for other small mammals (e.g., rodents and bats). Phylogeographic patterns also appear to have been shaped by both habitat and river vicariance, especially across the desert southwest. Further, continuing climate change during the Holocene coupled with anthropogenic modifications during the Anthropocene appears to be removing both of these barriers to current dispersal of western spotted skunks.

Colonization in North American Arid Lands: The Journey of Agarito (Berberis trifoliolata) Revealed by Multilocus Molecular Data and Packrat Midden Fossil Remains

Here we conduct research to understand the evolutionary history of a shrubby species known as Agarito (Berberis trifoliolata), an endemic species to the Chihuahuan Desert. We identify genetic signatures based on plastid DNA and AFLP markers and perform niche modelling and spatial connectivity analyses as well as niche modelling based on records in packrats to elucidate whether orogenic events such as mountain range uplift in the Miocene or the contraction/expansion dynamics of vegetation in response to climate oscillations in the Pliocene/Pleistocene had an effect on evolutionary processes in Agarito. Our results of current niche modelling and palaeomodelling showed that the area currently occupied by Berberis trifoliolata is substantially larger than it was during the Last Interglacial period and the Last Glacial Maximum. Agarito was probably confined to small areas in the Northeastern and gradually expanded its distribution just after the Last Glacial Maximum when the weather in the Chihuahuan Desert and adjacent regions became progressively warmer and drier. The most contracted range was predicted for the Interglacial period. Populations remained in stable areas during the Last Glacial Maximum and expanded at the beginning of the Holocene. Most genetic variation occured in populations from the Sierra Madre Oriental. Two groups of haplotypes were identified: the Mexican Plateau populations and certain Northeastern populations. Haplogroups were spatially connected during the Last Glacial Maximum and separated during interglacial periods. The most important prediction of packrat middens palaeomodelling lies in the Mexican Plateau, a finding congruent with current and past niche modelling predictions for agarito and genetic results. Our results corroborate that these climate changes in the Pliocene/Pleistocene affected the evolutionary history of agarito. The journey of agarito in the Chihuahuan Desert has been dynamic, expanding and contracting its distribution range and currently occupying the largest area in its history.

Evolutionary History of a Desert Shrub Ephedra przewalskii (Ephedraceae): Allopatric Divergence and Range Shifts in Northwestern China

Based on two chloroplast DNA sequences, psbA-trnH and trnT-trnF, phylogeographical patterns of a desert shrub, Ephedra przewalskii, were examined across most of its geographic range in northwestern China. A total of sixteen haplotypes were detected. There was a common haplotype in each basin, that was haplotype A in Tarim Basin, haplotype G in Junggar Basin, and haplotype M in Qaidam Basin. Genetic variance mainly occurred among populations, geographic regions, and eleven geographic groups subdivided by SAMOVA analysis. E. przewalskii likely had a smaller and more fragmented geographic range during the Last Glacial Maximum, which was determined based on ecological niche modelling. Three groups of E. przewalskii populations were detected to have experience range expansion, and this was based on significant values of Fu’s F_S, Tajima’s D, and unimodel mismatch distributions. The cold and dry climate during the glacial period of the Quaternary is postulated to have been a driver for significant genetic isolation and divergence among populations or groups in E. przewalskii, whereas the warmer and wetter climate during the interglacial period is speculated to have provided favourable conditions for range expansion of the species.

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.

A new phylogeographic pattern of endemic Bufo bankorensis in Taiwan Island is attributed to the genetic variation of populations

Aim

To comprehend the phylogeographic patterns of genetic variation in anurans at Taiwan Island, this study attempted to examine (1) the existence of various geological barriers (Central Mountain Ranges, CMRs); and (2) the genetic variation of Bufo bankorensis using mtDNA sequences among populations located in different regions of Taiwan, characterized by different climates and existing under extreme conditions when compared available sequences of related species B. gargarizans of mainland China.

Methodology/Principal Findings

Phylogenetic analyses of the dataset with mitochondrial DNA (mtDNA) D-loop gene (348 bp) recovered a close relationship between B. bankorensis and B. gargarizans, identified three distinct lineages. Furthermore, the network of mtDNA D-loop gene (564 bp) amplified (279 individuals, 27 localities) from Taiwan Island indicated three divergent clades within B. bankorensis (Clade W, E and S), corresponding to the geography, thereby verifying the importance of the CMRs and Kaoping River drainage as major biogeographic barriers. Mismatch distribution analysis, neutrality tests and Bayesian skyline plots revealed that a significant population expansion occurred for the total population and Clade W, with horizons dated to approximately 0.08 and 0.07 Mya, respectively. These results suggest that the population expansion of Taiwan Island species B. bankorensis might have resulted from the release of available habitat in post-glacial periods, the genetic variation on mtDNA showing habitat selection, subsequent population dispersal, and co-distribution among clades.

Conclusions

The multiple origins (different clades) of B. bankorensis mtDNA sequences were first evident in this study. The divergent genetic clades found within B. bankorensis could be independent colonization by previously diverged lineages; inferring B. bankorensis originated from B. gargarizans of mainland China, then dispersal followed by isolation within Taiwan Island. Highly divergent clades between W and E of B. bankorensis, implies that the CMRs serve as a genetic barrier and separated the whole island into the western and eastern phylogroups.

Cryptic genetic diversity and complex phylogeography of the boreal North American scorpion, Paruroctonus boreus (Vaejovidae)

Diverse studies in western North America have revealed the role of topography for dynamically shaping genetic diversity within species though vicariance, dispersal and range expansion. We examined patterns of phylogeographical diversity in the widespread but poorly studied North American vaejovid scorpion, Paruroctonus boreus Girard 1854. We used mitochondrial sequence data and parsimony, likelihood, and Bayesian inference to reconstruct phylogenetic relationships across the distributional range of P. boreus, focusing on intermontane western North America. Additionally, we developed a species distribution model to predict its present and historical distributions during the Last Glacial Maximum and the Last Interglacial Maximum. Our results documented complex phylogeographic relationships within P. boreus, with multiple, well-supported crown clades that are either geographically-circumscribed or widespread and separated by short, poorly supported internodes. We also observed subtle variation in predicted habitat suitability, especially at the northern, eastern and southern edges of the predicted distributional range under past climatic conditions. The complex phylogenetic relationships of P. boreus suggests that historical isolation and expansion of populations may have occurred. Variation in the predicted distributional range over time may implicate past climatic fluctuations in generating the patterns of genetic diversity observed in P. boreus. These findings highlight both the potential for cryptic biodiversity in widespread North American scorpion species and the importance of phylogeographical studies for understanding the factors responsible for generating the biodiversity of western North America.

Low phylogeographic structure of Rhinella arunco (Anura: Bufonidae), an endemic amphibian from the Chilean Mediterranean hotspot

Background

The Mediterranean zone of central Chile (30° to 38°S) is one of the 25 diversity hotspots in the world. However, there are few phylogeographic studies which identify the factors that have influenced population diversification in the fauna of this area. In this study, we investigated the phylogeographic structure of Rhinella arunco, an anuran endemic to Mediterranean Chile, using 160 individuals from 23 localities representative of its entire distribution (32° to 38°S).

Results

The haplotype network revealed four haplogroups, three of which overlap geographically and only one of which has an exclusive geographic distribution. An analysis of molecular variance indicated that neither watershed limits nor the main rivers in the current distribution of this species have been important geographic barriers against the dispersion of individuals. Finally, the Geneland analysis showed three population units, one of which concurs with one of the haplogroups found in the haplotype network. Together, these analyses indicated a low level of phylogeographic structure for this species. On the other hand, the highest levels of intrapopulational genetic variation were concentrated in the central part of the distribution (33° to 34°S), which may indicate an effect of Pleistocene glaciations on the genetic diversity of the populations in the extreme south of its range.

Conclusions

The low phylogeographic structure observed in R. arunco is a rarely documented pattern for amphibians and contrasts with the phylogeographic studies of other vertebrates which inhabit the same zone. This result may be attributed to a series of attributes of bufonids, related mainly to water retention and their reproductive biology, which have allowed them to disperse and colonize an enormous variety of environments.

Diversity and niche evolution along aridity gradients in north american lizards (phrynosomatidae)

Deserts occupy approximately 12% of the Earth's land surface, and are thought to have species poor but highly specialized biotas. However, few studies have examined the evolutionary origins of desert biotas and of diversity patterns along aridity gradients. Further, it is unclear if species occurring in more extreme conditions on a given niche axis (i.e., precipitation) are more specialized for those conditions (i.e., have narrower niche breadths). We address these questions here using a time-calibrated phylogeny and climatic data for 117 species of phrynosomatid lizards. Phrynosomatids are the most species-rich family of lizards in North America, and are found from deserts to rainforests. Surprisingly, we find that phrynosomatids have higher richness in more arid environments. This pattern occurs seemingly because they have been present in more arid habitats longer (~55 million years), and lineages in mesic environments are recently derived from more arid-dwelling ancestors. We find little support for the hypothesis that species in more extreme environments are more specialized. Instead, many desert-dwelling species are broadly distributed, and species in the most mesic environments have the broadest niche breadths. In summary, phrynosomatids offer a counterexample to the idea that arid regions are inhabited by a small number of recent and highly specialized lineages.

Neogene-dominated diversification in neotropical montane lichens: dating divergence events in the lichen-forming fungal genus Oropogon (Parmeliaceae)

PREMISE OF THE STUDY

Diversification in neotropical regions has been attributed to both Tertiary geological events and Pleistocene climatic fluctuations. However, the timing and processes driving speciation in these regions remain unexplored in many important groups. Here, we address the timing of diversification in the neotropical lichenized fungal genus Oropogon (Ascomycota) and assess traditional species boundaries.

METHODS

We analyzed sequence data from three loci to assess phenotypically circumscribed Oropogon species from the Oaxacan Highlands, Mexico. We provide a comparison of dated divergence estimates between concatenated gene trees and a calibrated multilocus species-tree using substitution rates for two DNA regions. We also compare estimates from a data set excluding ambiguously aligned regions and a data set including the hyper-variable regions in two ribosomal markers.

KEY RESULTS

Phylogenetic reconstructions were characterized by well-supported monophyletic clades corresponding to traditionally circumscribed species, with the exception of a single taxon. Divergence estimates indicate that most diversification of the sampled Oropogon species occurred throughout the Oligocene and Miocene, although diversification of a single closely related clade appears to have occurred during the late Pliocene and into the Pleistocene. Divergence estimates calculated from a data set with ambiguously aligned regions removed were much more recent than those from the full data set.

CONCLUSIONS

Overall, our analyses place the majority of divergence events of Oropogon species from the Oaxacan Highlands within the Neogene and provide strong evidence that climatic changes during the Pleistocene were not a major factor driving speciation in the lichenized genus Oropogon in neotropical highlands.

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

Background

Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescent-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies.

Results

We found that divergence estimates, from both concatenated gene tree and coalescent-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages.

Conclusions

Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.

Diversification in North American arid lands: niche conservatism, divergence and expansion of habitat explain speciation in the genus Ephedra

A lineage of 12 arid land shrubby species in the gymnosperm genus Ephedra (Gnetales) from North America is used to evaluate the influence of climate on speciation. With a long evolutionary history, and a well documented fossil record this lineage is an ideal model for understanding the process of speciation under a niche conservatism scenario. Using seven DNA molecular markers, Bayesian inference is carried out to uncover sister species and to estimate time of divergence of the lineages. Ecological niche models are generated for four parapatric and sympatric sister species and two analyses of niche evolution are performed, one based on ecological niche models and another using raw data and multivariate analysis. As previous analyses suggest, the diversification of North America Ephedra species may be the result of a recent secondary radiation. Both parapatric and sympatric species diverged mostly in a scenario of climatic niche conservatism. However, we also found strong evidence for niche divergence for one of the sister species pairs (E. californica-E. trifurca). Moreover, the multivariate analysis found environmental differences for some variables between sister species. The estimated divergence time of three pairs of sister species distributed in southwestern North America (E. cutleri-E. aspera, E. californica-E. trifurca and E. torreyana-E. viridis) is inferred to have occurred in the Late Miocene to Pliocene and for the sister species pair E. antisyphilitica-E. coryi distributed in the southern United States and northeastern Mexico, it was inferred from the Pliocene to Pleistocene. The orogenetic and climatic changes documented for these regions related to expansion of arid lands, may have contributed to the diversification in North American Ephedra, rather than adaptations to new climatic conditions.