Combining Biogeographic and Phylogenetic Data to Examine Primate Speciation: An Example Using Cercopithecin Monkeys

We combined phylogenetic and biogeographic data to examine the mode of speciation in a group of African monkeys, the Cercopithecini. If allopatric speciation is the major force producing species, then there should be a positive relationship between the relative divergence time of taxa and their degree of geographic range overlap. Alternatively, an opposite relationship between divergence time and geographic range overlap is consistent with sympatric speciation as the main mechanism underlying the cercopithecin radiation. We collected biogeographic and phylogenetic data for 19 guenon species from the literature. We digitized geographic range maps and utilized three different phylogenetic hypotheses based on Y chromosome, X chromosome, and mitochondrial (mtDNA) data. We used regressions with Monte Carlo simulation to examine the relationship between the relative time since divergence of taxa and their degree of geographic range overlap. We found that there was a positive relationship between relative divergence time and the proportion of geographic range overlap between taxa using all three molecular data sets. Our findings provide evidence for allopatric speciation being the common mode of diversification in the cercopithecin clade. Because most of these primates are forest adapted mammals, the cyclical contraction and expansion of African forests from the late Miocene to the present has likely been an important factor driving allopatric speciation. In addition, geographic barriers such as the Congo and Sanaga rivers have probably played a complementary role in producing new species within the clade.     

HISTORICAL DIVERGENCE AND GENE FLOW: COALESCENT ANALYSES OF MITOCHONDRIAL,AUTOSOMAL AND SEX‐LINKED LOCI IN PASSERINA BUNTINGS

Quantifying the role of gene flow during the divergence of closely related species is crucial to understanding the process of speciation. We collected DNA sequence data from 20 loci (one mitochondrial, 13 autosomal, and six sex‐linked) for population samples of Lazuli Buntings (Passerina amoena) and Indigo Buntings (Passerina cyanea) (Aves: Cardinalidae) to test explicitly between a strict allopatric speciation model and a model in which divergence occurred despite postdivergence gene flow. Likelihood ratio tests of coalescent‐based population genetic parameter estimates indicated a strong signal of postdivergence gene flow and a strict allopatric speciation model was rejected. Analyses of partitioned datasets (mitochondrial, autosomal, and sex‐linked) suggest the overall gene flow patterns are driven primarily by autosomal gene flow, as there is no evidence of mitochondrial gene flow and we were unable to reject an allopatric speciation model for the sex‐linked data. This pattern is consistent with either a parapatric divergence model or repeated periods of allopatry with gene flow occurring via secondary contact. These results are consistent with the low fitness of female avian hybrids under Haldane's rule and demonstrate that sex‐linked loci likely are important in the initial generation of reproductive isolation, not just its maintenance.     

Complex histories of repeated gene flow in Cameroon crater lake cichlids cast doubt on one of the clearest examples of sympatric speciation

One of the most celebrated examples of sympatric speciation in nature are monophyletic radiations of cichlid fishes endemic to Cameroon crater lakes. However, phylogenetic inference of monophyly may not detect complex colonization histories involving some allopatric isolation, such as double invasions obscured by genome‐wide gene flow. Population genomic approaches are better suited to test hypotheses of sympatric speciation in these cases. Here, we use comprehensive sampling from all four sympatric crater lake cichlid radiations in Cameroon and outgroups across Africa combined with next‐generation sequencing to genotype tens of thousands of SNPs. We find considerable evidence of gene flow between all four radiations and neighboring riverine populations after initial colonization. In a few cases, some sympatric species are more closely related to outgroups than others, consistent with secondary gene flow facilitating their speciation. Our results do not rule out sympatric speciation in Cameroon cichlids, but rather reveal a complex history of speciation with gene flow, including allopatric and sympatric phases, resulting in both reproductively isolated species and incipient species complexes. The best remaining non‐cichlid examples of sympatric speciation all involve assortative mating within microhabitats. We speculate that this feature may be necessary to complete the process of sympatric speciation in nature.     

DOES NICHE CONSERVATISM PROMOTE SPECIATION? A CASE STUDY IN NORTH AMERICAN SALAMANDERS

Recent speciation research has generally focused on how lineages that originate in allopatry evolve intrinsic reproductive isolation, or how ecological divergence promotes nonallopatric speciation. However, the ecological basis of allopatric isolation, which underlies the most common geographic mode of speciation, remains poorly understood and largely unstudied. Here, we explore the ecological and evolutionary factors that promote speciation in Desmognathus and Plethodon salamanders from temperate eastern North America. Based on published molecular phylogenetic estimates and the degree of geographic range overlap among extant species, we find strong evidence for a role for geographic isolation in speciation. We then examine the relationship between climatic variation and speciation in 16 sister-taxon pairs using geographic information system maps of climatic variables, new methods for modeling species' potential geographic distributions, and data on geographic patterns of genetic variation. In contrast to recent studies in tropical montane regions, we found no evidence for parapatric speciation along climatic gradients. Instead, many montane sister taxa in the Appalachian Highlands inhabit similar climatic niches and seemingly are allopatric because they are unable to tolerate the climatic conditions in the intervening lowlands. This temporal and spatial-ecological pattern suggests that niche conservatism, rather than niche divergence, plays the primary role in promoting allopatric speciation and montane endemism in this species-rich group of vertebrates. Our results demonstrate that even the relatively subtle climatic differences between montane and lowland habitats in eastern North America may play a key role in the origin of new species.     

Latitude, elevational climatic zonation and speciation in New World vertebrates

Many biodiversity hotspots are located in montane regions, especially in the tropics. A possible explanation for this pattern is that the narrow thermal tolerances of tropical species and greater climatic stratification of tropical mountains create more opportunities for climate-associated parapatric or allopatric speciation in the tropics relative to the temperate zone. However, it is unclear whether a general relationship exists among latitude, climatic zonation and the ecology of speciation. Recent taxon-specific studies obtained different results regarding the role of climate in speciation in tropical versus temperate areas. Here, we quantify overlap in the climatic distributions of 93 pairs of sister species of mammals, birds, amphibians and reptiles restricted to either the New World tropics or to the Northern temperate zone. We show that elevational ranges of tropical- and temperate-zone species do not differ from one another, yet the temperature range experienced by species in the temperate zone is greater than for those in the tropics. Moreover, tropical sister species tend to exhibit greater similarity in their climatic distributions than temperate sister species. This pattern suggests that evolutionary conservatism in the thermal niches of tropical taxa, coupled with the greater thermal zonation of tropical mountains, may result in increased opportunities for allopatric isolation, speciation and the accumulation of species in tropical montane regions. Our study exemplifies the power of combining phylogenetic and spatial datasets of global climatic variation to explore evolutionary (rather than purely ecological) explanations for the high biodiversity of tropical montane regions.     

The geography of mammalian speciation: mixed signals from phylogenies and range maps

The importance of geographic isolation in speciation has been debated since the 19th century. Since the beginning of the 20th century, the consensus has been that most speciation involves divergence in allopatry. This consensus was based largely on decades of observations by naturalists and verbal arguments against speciation without isolation. Recent attempts to quantify the importance of allopatric versus sympatric speciation using comparative methods called "age-range correlation" (ARC) suggest that allopatric speciation is more common than sympatric speciation. However, very few taxa have been studied and there are concerns about the adequacy of the methods. We propose methodological improvements including changes in the way overlap between clades is quantified and Monte Carlo methods to test the null hypothesis of no relationship between phylogenetic relatedness and geographic range overlap. We analyze 14 clades of mammals, chosen because of the availability of data and the consensus among mammalogists that speciation is routinely allopatric. Although data from a few clades clearly indicate allopatric speciation, divergence with gene flow is plausible in others and many results are inconclusive. The relative rarity of significant correlations between phylogenetic distance and range overlap may have three distinct causes: (1) post-speciation range changes, (2) relative rarity of range overlap, and (3) a mixture of geographic modes of speciation. Our results support skepticism about ARC's power for inferring the biogeography of speciation. Yet, even if few clades provide clear signals, meta-analytic approaches such as ARC may set bounds on the prevalence of alternative modes of speciation.     

The role of Quaternary environmental change in plant macroevolution: the exception or the rule?

The Quaternary has been described as an important time for genetic diversification and speciation. This is based on the premise that Quaternary climatic conditions fostered the isolation of populations and, in some instances, allopatric speciation. However, the 'Quaternary Ice-Age speciation model' rests on two key assumptions: (i) that biotic responses to climate change during the Quaternary were significantly different from those of other periods in Earth's history; and (ii) that the mechanisms of isolation during the Quaternary were sufficient in time and space for genetic diversification to foster speciation. These assumptions are addressed by examining the plant fossil record for the Quaternary (in detail) and for the past 410 Myr, which encompasses previous intervals of icehouse Earth. Our examination of the Quaternary record indicates that floristic responses to climate changes during the past 1.8 Myr were complex and that a distinction has to be made between those plants that were able to withstand the extremes of glacial conditions and those that could not. Generation times are also important as are different growth forms (e.g. herbaceous annuals and arborescent perennials), resulting in different responses in terms of genetic divergence rates during isolation. Because of these variations in the duration of isolation of populations and genomic diversification rates, no canonical statement about the predominant floristic response to climatic changes during the Quaternary (i.e. elevated rates of speciation or extinction, or stasis) is currently possible. This is especially true because of a sampling bias in terms of the fossil record of tree species over that of species with non-arborescent growth forms. Nevertheless, based on the available information, it appears that the dominant response of arborescent species during the Quaternary was extinction rather than speciation or stasis. By contrast, our examination of the fossil record of vascular plants for the past 410 Myr indicates that speciation rates often increased during long intervals of icehouse Earth (spanning up to 50 Myr). Therefore, longer periods of icehouse Earth than those occurring during the Quaternary may have isolated plant populations for sufficiently long periods of time to foster genomic diversification and allopatric speciation. Our results highlight the need for more detailed study of the fossil record in terms of finer temporal and spatial resolution than is currently available to examine the significance of intervals of icehouse Earth. It is equally clear that additional and detailed molecular studies of extant populations of Quaternary species are required in order to determine the extent to which these 'relic' species have genomically diversified across their current populations.     

LATITUDINAL VARIATION IN SPECIATION MECHANISMS IN FROGS

Speciation often has a strong geographical and environmental component, but the ecological factors that potentially underlie allopatric and parapatric speciation remain understudied. Two ecological mechanisms by which speciation may occur on geographic scales are allopatric speciation through niche conservatism and parapatric or allopatric speciation through niche divergence. A previous study on salamanders found a strong latitudinal pattern in the prevalence of these mechanisms, with niche conservatism dominating in temperate regions and niche divergence dominating in the tropics, and related this pattern to Janzen's hypothesis of greater climatic zonation between different elevations in the tropics. Here, we test for latitudinal patterns in speciation in a related but more diverse group of amphibians, the anurans. Using data from up to 79 sister-species pairs, we test for latitudinal variation in elevational and climatic overlap between sister species, and evaluate the frequency of speciation via niche conservatism versus niche divergence in relation to latitude. In contrast to salamanders, we find no tendency for greater niche divergence in the tropics or for greater niche conservatism in temperate regions. Although our results support the idea of greater climatic zonation in tropical regions, they show that this climatic pattern does not lead to straightforward relationships between speciation, latitude, and niche evolution.     

Phylogeography of the Rhabditis (Pellioditis) marina species complex: evidence for long-distance dispersal, and for range expansions and restricted gene flow in the northeast Atlantic

Pinpointing processes that structure the geographical distribution of genetic diversity of marine species and lead to speciation is challenging because of the lack of obvious dispersal barriers and the likelihood of substantial (passive) dispersal in oceans. In addition, cryptic radiations with sympatric distributions abound in marine species, challenging the allopatric speciation mechanism. Here, we present a phylogeographical study of the marine nematode species complex Rhabditis ( Pellioditis ) marina to investigate processes shaping genetic structure and speciation. Rhabditis ( P .) marina lives on decaying macroalgae in the intertidal, and may therefore disperse over considerable distances. Rhabditis ( P .) marina consists of several cryptic species sympatrically distributed at a local scale. Genetic variation in the COI gene was screened in 1362 specimens from 45 locations around the world. Two nuclear DNA genes (ITS and D2D3) were sequenced to infer phylogenetic species. We found evidence for ten sympatrically distributed cryptic species, seven of which show a strong genetic structuring. A historical signature showed evidence for restricted gene flow with occasional long-distance dispersal and range expansions pre-dating the last glacial maximum. Our data also point to a genetic break around the British Isles and a contact zone in the Southern Bight of the North Sea. We provide evidence for the transoceanic distribution of at least one cryptic species (PmIII) and discuss the dispersal capacity of marine nematodes. The allopatric distribution of some intraspecific phylogroups and of closely related cryptic species points to the potential for allopatric speciation in R. ( P .) marina .     

Karyotypic diversity and speciation in Agrodiaetus butterflies

That chromosomal rearrangements may play an important role in maintaining postzygotic isolation between well-established species is part of the standard theory of speciation. However, little evidence exists on the role of karyotypic change in speciation itself--in the establishment of reproductive barriers between previously interbreeding populations. The large genus Agrodiaetus (Lepidoptera: Lycaenidae) provides a model system to study this question. Agrodiaetus butterflies exhibit unusual interspecific diversity in chromosome number, from n= 10 to n= 134; in contrast, the majority of lycaenid butterflies have n= 23/24. We analyzed the evolution of karyotypic diversity by mapping chromosome numbers on a thoroughly sampled mitochondrial phylogeny of the genus. Karyotypic differences accumulate gradually between allopatric sister taxa, but more rapidly between sympatric sister taxa. Overall, sympatric sister taxa have a higher average karyotypic diversity than allopatric sister taxa. Differential fusion of diverged populations may account for this pattern because the degree of karyotypic difference acquired between allopatric populations may determine whether they will persist as nascent biological species in secondary sympatry. This study therefore finds evidence of a direct role for chromosomal rearrangements in the final stages of animal speciation. Rapid karyotypic diversification is likely to have contributed to the explosive speciation rate observed in Agrodiaetus, 1.6 species per million years.     

The population genetics of the origin and divergence of the Drosophila simulans complex species

The origins and divergence of Drosophila simulans and close relatives D. mauritiana and D. sechellia were examined using the patterns of DNA sequence variation found within and between species at 14 different genes. D. sechellia consistently revealed low levels of polymorphism, and genes from D. sechellia have accumulated mutations at a rate that is approximately 50% higher than the same genes from D. simulans. At synonymous sites, D. sechellia has experienced a significant excess of unpreferred codon substitutions. Together these observations suggest that D. sechellia has had a reduced effective population size for some time, and that it is accumulating slightly deleterious mutations as a result. D. simulans and D. mauritiana are both highly polymorphic and the two species share many polymorphisms, probably since the time of common ancestry. A simple isolation speciation model, with zero gene flow following incipient species separation, was fitted to both the simulans/mauritiana divergence and the simulans/sechellia divergence. In both cases the model fit the data quite well, and the analyses revealed little evidence of gene flow between the species. The exception is one gene copy at one locus in D. sechellia, which closely resembled other D. simulans sequences. The overall picture is of two allopatric speciation events that occurred quite near one another in time.     

THE RATE TEST OF SPECIATION: ESTIMATING THE LIKELIHOOD OF NON‐ALLOPATRIC SPECIATION FROM REPRODUCTIVE ISOLATION RATES IN DROSOPHILA

Among the most debated subjects in speciation is the question of its mode. Although allopatric (geographical) speciation is assumed the null model, the importance of parapatric and sympatric speciation is extremely difficult to assess and remains controversial. Here I develop a novel approach to distinguish these modes of speciation by studying the evolution of reproductive isolation (RI) among taxa. I focus on the Drosophila genus, for which measures of RI are known. First, I incorporate RI into age‐range correlations. Plots show that almost all cases of weak RI are between allopatric taxa whereas sympatric taxa have strong RI. This either implies that most reproductive isolation (RI) was initiated in allopatry or that RI evolves too rapidly in sympatry to be captured at incipient stages. To distinguish between these explanations, I develop a new “rate test of speciation” that estimates the likelihood of non‐allopatric speciation given the distribution of RI rates in allopatry versus sympatry. Most sympatric taxa were found to have likely initiated RI in allopatry. However, two putative candidate species pairs for non‐allopatric speciation were identified (5% of known Drosophila). In total, this study shows how using RI measures can greatly inform us about the geographical mode of speciation in nature.     

Radiating on oceanic islands: patterns and processes of speciation in the land snail genus Theba (Risso 1826)

Island radiations have played a major role in shaping our current understanding of allopatric, sympatric and parapatric speciation. However, the fact that species divergence correlates with island size emphasizes the importance of geographic isolation (allopatry) in speciation. Based on molecular and morphological data, we investigated the diversification of the land snail genus Theba on the two Canary Islands of Lanzarote and Fuerteventura. Due to the geological history of both islands, this study system provides ideal conditions to investigate the interplay of biogeography, dispersal ability and differentiation in generating species diversity. Our analyses demonstrated extensive cryptic diversification of Theba on these islands, probably driven mainly by non-adaptive allopatric differentiation and secondary gene flow. In a few cases, we observed a complete absence of gene flow among sympatrically distributed forms suggesting an advanced stage of speciation. On the Jandía peninsula genome scans suggested genotype-environment associations and potentially adaptive diversification of two closely related Theba species to different ecological environments. We found support for the idea that genetic differentiation was enhanced by divergent selection in different environments. The diversification of Theba on both islands is therefore best explained by a mixture of non-adaptive and adaptive speciation, promoted by ecological and geomorphological factors.     

Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach

Background  

Marine allopatric speciation is an enigma because pelagic larval dispersal can potentially connect disjunct populations thereby preventing reproductive and morphological divergence. Here we present a new hierarchical approximate Bayesian computation model (HABC) that tests two hypotheses of marine allopatric speciation: 1.) "soft vicariance", where a speciation involves fragmentation of a large widespread ancestral species range that was previously connected by long distance gene flow; and 2.) peripatric colonization, where speciations in peripheral archipelagos emerge from sweepstakes colonizations from central source regions. The HABC approach analyzes all the phylogeographic datasets at once in order to make across taxon-pair inferences about biogeographic processes while explicitly allowing for uncertainty in the demographic differences within each taxon-pair. Our method uses comparative phylogeographic data that consists of single locus mtDNA sequences from multiple co-distributed taxa containing pairs of central and peripheral populations. We use the method on two comparative phylogeographic data sets consisting of cowrie gastropod endemics co-distributed in the Hawaiian (11 taxon-pairs) and Marquesan archipelagos (7 taxon-pairs).     

Reproductive isolation among allopatric Drosophila montana populations

An outstanding goal in speciation research is to trace the mode and tempo of the evolution of barriers to gene flow. Such research benefits from studying incipient speciation, in which speciation between populations has not yet occurred, but where multiple potential mechanisms of reproductive isolation (RI: i.e., premating, postmating‐prezygotic (PMPZ), and postzygotic barriers) may act. We used such a system to investigate these barriers among allopatric populations of Drosophila montana. In all heteropopulation crosses we found premating (sexual) isolation, which was either symmetric or asymmetric depending on the population pair compared. Postmating isolation was particularly strong in crosses involving males from one of the study populations, and while sperm were successfully transferred, stored, and motile, we experimentally demonstrated that the majority of eggs produced were unfertilized. Thus, we identified the nature of a PMPZ incompatibility. There was no evidence of intrinsic postzygotic effects. Measures of absolute and relative strengths of pre‐ and postmating barriers showed that populations differed in the mode and magnitude of RI barriers. Our results indicate that incipient RI among populations can be driven by different contributions of both premating and PMPZ barriers occurring between different population pairs and without the evolution of postzygotic barriers.     

Glacial survival east and west of the 'Mekong-Salween Divide' in the Himalaya-Hengduan Mountains region as revealed by AFLPs and cpDNA sequence variation in Sinopodophyllum hexandrum (Berberidaceae)

Molecular phylogeographic studies have recently begun to elucidate how plant species from the Qinghai-Tibetan Plateau (QTP) and adjacent regions responded to the Quaternary climatic oscillations. In this regard, however, far less attention has been paid to the southern and south-eastern declivities of the QTP, i.e. the Himalaya-Hengduan Mountains (HHM) region. Here, we report a survey of amplified fragment length polymorphisms (AFLPs) and chloroplast DNA (cpDNA) sequence variation in the HHM endemic Sinopodophyllum hexandrum, a highly selfing alpine perennial herb with mainly gravity-dispersed berries (105 individuals, 19 localities). We specifically aimed to test a vicariant evolutionary hypothesis across the 'Mekong-Salween Divide', a known biogeographic and phytogeographic boundary of north-to-south trending river valleys separating the East Himalayas and Hengduan Mts. Both cpDNA and AFLPs identified two divergent phylogroups largely congruent with these mountain ranges. There was no genetic depauperation in the more strongly glaciated East Himalayas (AFLPs: H(E)=0.031; cpDNA: h(S)=0.133) compared to the mainly ice-free Hengduan Mts. (AFLPs: H(E)=0.037; cpDNA: h(S)=0.082), while population differentiation was consistently higher in the former region (AFLPs: Φ(ST)=0.522 vs. 0.312; cpDNA: Φ(ST)=0.785 vs. 0.417). Our results suggest that East Himalayan and Hengduan populations of S. hexandrum were once fragmented, persisted in situ during glacials in both areas, and have not merged again, except for a major instance of inter-lineage chloroplast capture identified at the MSD boundary. Our coalescent time estimate for all cpDNA haplotypes (c. 0.37-0.48 mya), together with paleogeological evidence, strongly rejects paleo-drainage formation as a mechanism underlying allopatric fragmentation, whereas mountain glaciers following the ridges of the MSD during glacials (and possible interglacials) could have been responsible. This study thus indicates an important role for mountain glaciers in driving (incipient) allopatric speciation across the MSD in the HHM region by causing vicariant lineage divergence and acting as barriers to post-divergence gene flow.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Pleistocene speciation with and without gene flow in Euphaea damselflies of subtropical and tropical East Asian islands

Climatic oscillations during the Pleistocene period could have had a profound impact on the origin of tropical species by the alternation of allopatric isolation and interpopulation gene flow cycles. However, whether tropical speciation involves strictly allopatric isolation, or proceeds in the face of homogenizing gene flow, is relatively unclear. Here, we investigated geographical modes of speciation in four closely related Euphaea damselfly species endemic to the subtropical and tropical East Asian islands using coalescent analyses of a multilocus data set. The reconstructed phylogenies demonstrated distinct species status for each of the four species and the existence of two sister species pairs, Euphaea formosa/E. yayeyamana and E. decorata/E. ornata. The species divergence time of the sibling Euphaea damselflies dates back to within the last one Mya of the Middle to Lower Pleistocene. The speciation between the populous E. formosa of Taiwan and the less numerous E. yayeyamana of the Yaeyama islands occurred despite significant bidirectional, asymmetric gene flow, which is strongly inconsistent with a strictly allopatric model. In contrast, speciation of the approximately equal-sized populations of E. decorata of the southeast Asian mainland and E. ornata of Hainan is inferred to have involved allopatric divergence without gene flow. Our findings suggest that differential selection of natural or sexual environments is a prominent driver of species divergence in subtropical E. formosa and E. yayeyamana; whereas for tropical E. decorata and E. ornata at lower latitudes, allopatric isolation may well be a pivotal promoter of species formation.     

The impact of Quaternary climate oscillations on divergence times and historical population sizes in Thylamys opossums from the Andes

Climate oscillations during the Quaternary altered the distributions of terrestrial animals at a global scale. In mountainous regions, temperature fluctuations may have led to shifts in range size and population size as species tracked their shifting habitats upslope or downslope. This creates the potential for both allopatric speciation and population size fluctuations, as species are either constrained to smaller patches of habitat at higher elevations or able to expand into broader areas at higher latitudes. We considered the impact of climate oscillations on three pairs of marsupial species from the Andes (Thylamys opossums) by inferring divergence times and demographic changes. We compare four different divergence dating approaches, using anywhere from one to 26 loci. Each pair comprises a northern (tropical) lineage and a southern (subtropical to temperate) lineage. We predicted that divergences would have occurred during the last interglacial (LIG) period approximately 125 000 years ago and that population sizes for northern and southern lineages would either contract or expand, respectively. Our results suggest that all three north–south pairs diverged in the late Pleistocene during or slightly after the LIG. The three northern lineages showed no signs of population expansion, whereas two southern lineages exhibited dramatic, recent expansions. We attribute the difference in responses between tropical and subtropical lineages to the availability of ‘montane‐like’ habitats at lower elevations in regions at higher latitudes. We conclude that climate oscillations of the late Quaternary had a powerful impact on the evolutionary history of some of these species, both promoting speciation and leading to significant population size shifts.