Biogeographic barriers,Pleistocene refugia,and climatic gradients in the southeastern Nearctic drive diversification in cornsnakes (Pantherophis guttatus complex)

The southeastern Nearctic is a biodiversity hotspot that is also rich in cryptic species. Numerous hypotheses (e.g., vicariance, local adaptation, and Pleistocene speciation in glacial refugia) have been tested in an attempt to explain diversification and the observed pattern of extant biodiversity. However, previous phylogeographic studies have both supported and refuted these hypotheses. Therefore, while data support one or more of these diversification hypotheses, it is likely that taxa are forming within this region in species‐specific ways. Here, we generate a genomic data set for the cornsnakes (Pantherophis guttatus complex), which are widespread across this region, spanning both biogeographic barriers and climatic gradients. We use phylogeographic model selection combined with hindcast ecological niche models to determine regions of habitat stability through time. This combined approach suggests that numerous drivers of population differentiation explain the current diversity of this group of snakes. The Mississippi River caused initial speciation in this species complex, with more recent divergence events linked to adaptations to ecological heterogeneity and allopatric Pleistocene refugia. Lastly, we discuss the taxonomy of this group and suggest there may be additional cryptic species in need of formal recognition.     

Genealogical portraits of speciation in montane grasshoppers (genus Melanoplus) from the sky islands of the Rocky Mountains

Grasshoppers in the genus Melanoplus have undergone a radiation in the 'sky islands' of western North America, with many species originating during the Pleistocene. Despite their recent origins, phylogenetic analyses indicate that all the species exhibit monophyletic or paraphyletic gene trees. The objectives of this study were to determine whether the monophyletic genealogies are the result of a bottleneck at speciation and to investigate the extent to which the different phylogenetic states of eight species (i.e. monophyletic versus paraphyletic gene trees) can be ascribed to the effects of speciation. A coalescent simulation was used to test for a bottleneck at speciation in each species. The effective population sizes and demographic histories of species were compared across taxa to evaluate the possibility that the paraphyly versus monophyly of the species reflects differential rates of lineage loss rather than speciation mode. While coalescent analyses indicate that the monophyly of Melanoplus species might not be indicative of bottlenecks at speciation, the results suggest that the paraphyletic gene trees may reflect the demography of speciation, involving localized divergences in the ancestral species. With respect to different models of Pleistocene divergence, the data do not support a model of founder-effect speciation but are compatible with divergence in allopatric refugia.     

Adaptive radiation and convergence in subdivisions of the butterfly genus Heliconius (Lepidoptera: Nymphalidae)

The process of adaptive radiation and convergence, usually regarded as a feature of macro-evolution, can be seen in the mimetic colour patterns of the butterflies within the confines of the South American genus Heliconius. This can be shown by dividing the genus into subgroups on the basis of adult, pupal and larval morphology: the theory that the mimicry between species results solely from close systematic relationships is thereby refuted, as members of the same morphological group can display widely divergent mimetic patterns, and conversely mutual mimics may belong to several different morphological groups. Various forms of parallel and convergent evolution are thought to account for the present pattern of mimicry, the process is known to start even before full speciation has taken place. A new subgenus (Neruda) is created to contain three atypical members of the genus.     

Pliocene forest dynamics as a primary driver of African bird speciation

Aim Montane tropics are areas of high endemism, and mechanisms driving this endemism have been receiving increasing attention at a global scale. A general trend is that climatic factors do not explain the species richness of species with small to medium‐sized geographic ranges, suggesting that geological and evolutionary processes must be considered. On the African continent, several hypotheses including both refugial and geographic uplift models have been advanced to explain avian speciation and diversity in the lowland forest and montane regions of central and eastern Africa; montane regions in particular are recognized as hotspots of vertebrate endemism. Here, we examine the possible role of these models in driving speciation in a clade of African forest robins. Location Africa. Methods We constructed the first robustly supported molecular phylogenetic hypothesis of forest robins. On this phylogeny, we reconstructed habitat‐based distributions and geographic distributions relative to the Albertine Rift. We also estimated the timing of lineage divergences via a molecular clock. Results Robust estimates of phylogenetic relationships and clock‐based divergences reject Miocene tectonic uplift and Pleistocene forest refugia as primary drivers of speciation in forest robins. Instead, our data suggest that most forest robin speciation took place in the Late Pliocene, from 3.2 to 2.2 Ma. Distributional patterns are complex, with the Albertine Rift region serving as a general east–west break across the group. Montane distributions are inferred to have evolved four times. Main conclusions Phylogenetic divergence dates coincide with a single period of lowland forest retraction in the late Pliocene, suggesting that most montane speciation resulted from the rapid isolation of populations in montane areas, rather than montane areas themselves being drivers of speciation. This conclusion provides additional evidence that Pliocene climate change was a major driver of speciation in broadly distributed African animal lineages. We further show that lowland forest robins are no older than their montane relatives, suggesting that lowland areas are not museums which house ‘ancient’ taxa; rather, for forest robins, montane areas should be viewed as living museums of a late Pliocene diversification event. A forest refugial pattern is operating in Africa, but it is not constrained to the Pleistocene.     

A molecular phylogeny of the neotropical butterfly genus Anartia (Lepidoptera: Nymphalidae)

While Anartia butterflies have served as model organisms for research on the genetics of speciation, no phylogeny has been published to describe interspecific relationships. Here, we present a molecular phylogenetic analysis of Anartia species relationships, using both mitochondrial and nuclear genes. Analyses of both data sets confirm earlier predictions of sister species pairings based primarily on genital morphology. Yet both the mitochondrial and nuclear gene phylogenies demonstrate that Anartia jatrophae is not sister to all other Anartia species, but rather that it is sister to the Anartia fatima-Anartia amathea lineage. Traditional biogeographic explanations for speciation across the genus relied on A. jatrophae being sister to its congeners. These explanations invoked allopatric divergence of sister species pairs and multiple sympatric speciation events to explain why A. jatrophae flies alongside all its congeners. The molecular phylogenies are more consistent with lineage divergence due to vicariance, and range expansion of A. jatrophae to explain its sympatry with congeners. Further interpretations of the tree topologies also suggest how morphological evolution and eco-geographic adaptation may have set species range boundaries.     

Mutualistic Interactions Drive Ecological Niche Convergence in a Diverse Butterfly Community

Ecological communities are structured in part by evolutionary interactions among their members. A number of recent studies incorporating phylogenetics into community ecology have upheld the paradigm that competition drives ecological divergence among species of the same guild. However, the role of other interspecific interactions, in particular positive interactions such as mutualism, remains poorly explored. We characterized the ecological niche and inferred phylogenetic relationships among members of a diverse community of neotropical Müllerian mimetic butterflies. Müllerian mimicry is one of the best studied examples of mutualism, in which unpalatable species converge in wing pattern locally to advertize their toxicity to predators. We provide evidence that mutualistic interactions can drive convergence along multiple ecological axes, outweighing both phylogeny and competition in shaping community structure. Our findings imply that ecological communities are adaptively assembled to a much greater degree than commonly suspected. In addition, our results show that phenotype and ecology are strongly linked and support the idea that mimicry can cause ecological speciation through multiple cascading effects on species' biology.     

Speciation in North American black basses,Micropterus (Actinopterygii: Centrarchidae)

Abstract.— The Pleistocene Epoch has been frequently cited as a period of intense speciation for a significant portion of temperate continental biotas. To critically assess the role of Pleistocene glaciations on the evolution of the freshwater fish clade Micropterus , we use a phylogenetic analysis of complete gene sequences from two mitochondrial genes (cytochrome b and ND2), and a fossil calibration of the molecular clock to estimate ages of speciation events and rates of diversification. The absence of substantial morphological and ecological divergence together with endemism of five of the eight species in North American tributaries of the Gulf of Mexico may be interpreted as the result of a recent Pleistocene origin for these species. Speciation dates in Micropterus range from 1.01 ± 0.32 to 11.17 ± 1.02 million years ago. Only one speciation event is dated to the Pleistocene, and rates of diversification are not significantly variable in Micropterus. The premise that the Pleistocene was an exceptional period of speciation in Micropterus is not supported. Instead, a Gulf Coast allopatric speciation model is proposed, and predicts periods of dynamic speciation driven by sea level fluctuations in the Late Miocene and Pliocene. The Pleistocene, however, was a period of significant intraspecific mitochondrial lineage diversification. The application of the Gulf Coast allopatric speciation model to the remaining aquatic fauna of the Gulf of Mexico coast in North America will rely on robust phylogenetic hypotheses and accurate age estimations of speciation events.     

The roles of wing color pattern and geography in the evolution of Neotropical Preponini butterflies

Diversification rates and evolutionary trajectories are known to be influenced by phenotypic traits and the geographic history of the landscapes that organisms inhabit. One of the most conspicuous traits in butterflies is their wing color pattern, which has been shown to be important in speciation. The evolution of many taxa in the Neotropics has also been influenced by major geological events. Using a dated, species‐level molecular phylogenetic hypothesis for Preponini, a colorful Neotropical butterfly tribe, we evaluated whether diversification rates were constant or varied through time, and how they were influenced by color pattern evolution and biogeographical events. We found that Preponini originated approximately 28 million years ago and that diversification has increased through time consistent with major periods of Andean uplift. Even though some clades show evolutionarily rapid transitions in coloration, contrary to our expectations, these shifts were not correlated with shifts in diversification. Involvement in mimicry with other butterfly groups might explain the rapid changes in dorsal color patterns in this tribe, but such changes have not increased species diversification in this group. However, we found evidence for an influence of major Miocene and Pliocene geological events on the tribe''s evolution. Preponini apparently originated within South America, and range evolution has since been dynamic, congruent with Andean geologic activity, closure of the Panama Isthmus, and Miocene climate variability.     

Phylogeny of Saxifraga section Saxifraga subsection Arachnoideae (Saxifragaceae) and the origin of low elevation shade-dwelling species

Saxifraga section Saxifraga subsection Arachnoideae is a lineage of 12 species distributed mainly in the European Alps. It is unusual in terms of ecological diversification by containing both high elevation species from exposed alpine habitats and low elevation species from shady habitats such as overhanging rocks and cave entrances. Our aims are to explore which of these habitat types is ancestral, and to identify the possible drivers of this remarkable ecological diversification. Using a Hybseq DNA-sequencing approach and a complete species sample we reconstructed and dated the phylogeny of subsection Arachnoideae. Using Landolt indicator values, this phylogenetic tree was used for the reconstruction of the evolution of temperature, light and soil pH requirements in this lineage. Diversification of subsection Arachnoideae started in the late Pliocene and continued through the Pleistocene. Both diversification among and within clades was largely allopatric, and species from shady habitats with low light requirements are distributed in well-known refugia. We hypothesize that low light requirements evolved when species persisting in cold-stage refugia were forced into marginal habitats by more competitive warm-stage vegetation. While we do not claim that such competition resulted in speciation, it very likely resulted in adaptive evolution.     

First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

What can hybrid zones tell us about speciation? The case of Heliconius erato and H. himera (Lepidoptera: Nymphalidae)

To understand speciation we need to study the genetics and ecology of intermediate cases where interspecific hybridization still occurs. Two closely related species of Heliconius butterflies meet this criterion: Heliconius himera is endemic to dry forest and thorn scrub in southern Ecuador and northern Peru, while its sister species, H. erato , is ubiquitous in wet forest throughout south and central America. In three known zones of contact, the two species remain distinct, while hybrids are found at low frequency. Collections in southern Ecuador show that the contact zone is about 5 km wide, half the width of the narrowest clines between colour pattern races of H. erato. The narrowness of this dine argues that very strong selection (s ≅ 1) is maintaining the parapatric distributions of these two species. The zone is closely related with a habitat transition from wet to dry forest, which suggests that the narrow zone of parapatry is maintained primarily by ecological adaptation. Selection on colour pattern loci, assortative mating and hybrid inviability may also be important. The genetics of hybrids between the two species shows that the major gene control of pattern elements is similar to that found in previous studies of H. erato races, and some of the loci are homologous. This suggests that similar genetic processes are involved in the morphological divergence of species and races. Evidence from related Heliconius supports a hypothesis that ecological adaptation is the driving force for speciation in the group.     

The combined role of dispersal and niche evolution in the diversification of Neotropical lizards

Ecological requirements and environmental conditions can influence diversification across temporal and spatial scales. Understanding the role of ecological niche evolution under phylogenetic contexts provides insights on speciation mechanisms and possible responses to future climatic change. Large‐scale phyloclimatic studies on the megadiverse Neotropics, where biomes with contrasting vegetation types occur in narrow contact, are rare. We integrate ecological and biogeographic data with phylogenetic comparative methods, to investigate the relative roles of biogeographic events and niche divergence and conservatism on the diversification of the lizard genus Kentropyx Spix, 1825 (Squamata: Teiidae), distributed in South American rainforests and savannas. Using five molecular markers, we estimated a dated species tree, which recovered three clades coincident with previously proposed species groups diverging during the mid‐Miocene. Biogeography reconstruction indicates a role of successive dispersal events from an ancestral range in the Brazilian Shield and western Amazonia. Ancestral reconstruction of climatic tolerances and niche overlap metrics indicates a trend of conservatism during the diversification of groups from the Amazon Basin and Guiana Shield, and a strong signal of niche divergence in the Brazilian Shield savannas. Our results suggest that climatic‐driven divergence at dynamic forest‐savanna borders might have resulted in adaptation to new environmental niches, promoting habitat shifts and shaping speciation patterns of Neotropical lizards. Dispersal and ecological divergence could have a more important role in Neotropical diversification than previously thought.     

Tests of pleistocene speciation in montane grasshoppers (genus Melanoplus) from the sky islands of western North America

There has a been a resurgence of debate on whether the Pleistocene glaciations inhibited speciation. This study tests a model of Pleistocene speciation, estimating the phylogenetic relationships and divergence times of 10 species of montane grasshoppers, genus Melanoplus, using 1300 bp of the mitochondrial gene cytochrome oxidase I (COI). Based on average pairwise distances (corrected for multiple substitutions using Kimura's two-parameter model), all species appear to have originated within the Pleistocene. Sequence divergences between species are less than 4%, corresponding to divergence times less than 1.7 million years ago. Branching patterns among the species suggest that speciation was associated with more than one glacial-interglacial cycle. A likelihood-ratio test rejected a model of simultaneous species origins, the predicted branching pattern if species arose from the fragmentation of a widespread ancestor. These grasshoppers live in an area that was previously glaciated and, as inhabitants of the northern Rocky Mountain sky islands, underwent latitudinal and probably altitudinal shifts in distribution in response to climatic fluctuations. Given the repeated distributional shifts and range overlap of the taxa, there most likely has been ample opportunity for population mixing. However, despite periodic glacial cycles, with more than 10 major glaciations over the past million years and climatic fluctuations over as short a time scale as 10(3) to 10(4) years, the dynamic history of the Pleistocene did not preclude speciation. Although relationships among some taxa remain unresolved, these grasshopper species, even with their recent origins, exhibit genetic coherence and monophyletic or paraphyletic gene trees. The frequency of glacial cycles suggests that the speciation process must have been extremely rapid. These species of grasshoppers are morphologically very similar, differing primarily in the shape of the male genitalia. These characters are posited to be under sexual selection, may play an important role in reproductive isolation, and are known to diverge rapidly. This suggests the rapidity of evolution of reproductive isolation may determine whether species divergences occurred during the Pleistocene glaciations.     

啄木鸟科内谱系生物地理学新探

啄木鸟群体共包含254种,但早期对啄木鸟的分类学研究局限于形态学、行为学和鸣唱法等方向,分子研究结果又彼此存在差异,因此啄木鸟分类系统目前还有争议。分子数据的大量累积为啄木鸟系统发育的研究提供了条件。利用从GenBank获取的ND2基因序列,对啄木鸟科179种啄木鸟采用NJ法、MP法以及ML法分别进行了分子进化树构建,同时用NJ法重建了Picumnus属分子进化树。此外又利用cytb基因序列采用NJ法对Dryobates属和Picoides属进行分子进化树重建。我们的结果显示,姬啄木亚科Picumnus属新大陆物种彼此的进化关系仍需进一步研究,而该属唯一一种旧大陆种P.innominatus,种内遗传差异在种间水平依然显著,这说明其种内遗传差异一直被低估了。在啄木鸟亚科内,Hemicircus属和Nesoctites属为姐妹关系,而Geocolaptes olivaceus插入了Campethera属支系中,并非与其构成姐妹群体。此外,一些姐妹群在新旧大陆上的分布,推测这种分布可能与第四纪冰川有关。啄木鸟从旧大陆向新大陆迁移,在第四纪冰川作用下不断发生分布范围的扩张与收缩,最终形成现存的分布格局。     

Mutualistic mimicry enhances species diversification through spatial segregation and extension of the ecological niche space

Species richness varies among clades, yet the drivers of diversification creating this variation remain poorly understood. While abiotic factors likely drive some of the variation in species richness, ecological interactions may also contribute. Here, we examine one class of potential contributors to species richness variation that is particularly poorly understood: mutualistic interactions. We aim to elucidate large‐scale patterns of diversification mediated by mutualistic interactions using a spatially explicit population‐based model. We focus on mutualistic Müllerian mimicry between conspicuous toxic prey species, where convergence in color patterns emerges from predators' learning process. To investigate the effects of Müllerian mimicry on species diversification, we assume that some speciation events stem from shifts in ecological niches, and can also be associated with shift in mimetic color pattern. Through the emergence of spatial mosaics of mimetic color patterns, Müllerian mimicry constrains the geographical distribution of species and allows different species occupying similar ecological niches to exist simultaneously in different regions. Müllerian mimicry and the resulting spatial segregation of mimetic color patterns thus generate more balanced phylogenetic trees and increase overall species diversity. Our study sheds light on complex effects of Müllerian mimicry on ecological, spatial, and phylogenetic diversification.     

Ice sheets promote speciation in boreal birds

The premise that Pleistocene ice ages played an important role in generating present-day species diversity has been challenged by genetic data indicating that most of the youngest terrestrial species on Earth coalesced long before major glacial advances. However, study has been biased towards faunas distributed at low latitudes that were not directly fragmented by advancing ice sheets. Using mitochondrial sequence divergence and a molecular clock, we compared the coalescence times of pairs of avian species belonging to superspecies complexes from the high-latitude boreal forest with those of sub-boreal and tropical avifaunas of the New World. Remarkably, all coalescence events in boreal superspecies date to the Pleistocene, providing direct evidence that speciation was commonly initiated during recent glacial periods. A pattern of endemism in boreal superspecies plausibly links the timing of divergence to the fragmentation of the boreal forest by ice sheets during the Mid- and Late Pleistocene. In contrast to the boreal superspecies, only 56% of sub-boreal and 46% of tropical superspecies members coalesced during the Pleistocene, suggesting that avifaunas directly fragmented by ice sheets experienced rapid rates of diversification, whereas those distributed farther south were affected to a lesser extent. One explanation for the absence of pre-Pleistocene superspecies in boreal avifaunas is that strong selection pressures operated in boreal refugia, causing superspecies members to achieve ecological differentiation at an accelerated rate.     

Wing pattern evolution and the origins of mimicry among North American admiral butterflies (Nymphalidae: Limenitis)

The evolution of wing pattern diversity in butterflies has emerged as a model system for understanding the origins and maintenance of adaptive phenotypic novelty. Admiral butterflies (genus Limenitis) are an attractive system for studying wing pattern diversity because mimicry is common among the North American species and hybrid zones occur wherever mimetic and non-mimetic wing pattern races meet. However, the utility of this system has been limited because the evolutionary relationships among these butterflies remain unclear. Here I present a robust species-level phylogeny of Limenitis based on 1911 bp of two mitochondrial genes (COI and COII) and 904 bp of EF1-alpha for all five of the Nearctic species/wing pattern races, the majority of the Palearctic species, and three outgroup genera; Athyma, Moduza (Limenitidini), and Neptis (Limenitidinae: Neptini). Maximum-likelihood and Bayesian analyses indicate that the North American species are a well-supported, monophyletic lineage that is most closely related to the widespread, Palearctic, Poplar admiral (L. populi). Within North America, the Viceroy (L. archippus) is the basal lineage while the relationships among the remaining species are not well resolved. A combined maximum-likelihood analysis, however, indicates that the two western North America species (L. lorquini and L. weidemeyerii) are sister taxa and closely related to the wing pattern subspecies of the polytypic Limenitis arthemis species complex. These results are consistent with (1) an ancestral host-shift to Salicaceae by the common ancestor of the Poplar admiral and the Nearctic admiral lineage, (2) a single colonization of the Nearctic, and (3) a subsequent radiation of the North American forms leading to at least three independent origins of mimicry.     

The polymorphism, mimicry, and host plant relationships of Hypolimnas butterflies

The 30 species of Hypolimnas butterflies recognized are divided amongst six pattern morphism classes: monomorphic, weakly dimorphic, simple dual, unimodal polymorphic, partial 9-limited polymorphic, and multiple ♀ dual polymorphic. Eighty per cent of Hypolimnas species belong to one or other of the first three categories. About 60% of the species show evidence of mimicry, some 40% being mimetic in both sexes. In the Oriental region the great majority of mimetic species within the genus are modelled on Euptoea , whilst those of the Afrotropical region are mostly modelled on Amauris. The host plants recorded for the genus strongly suggest that the Urticaceae are primitive for the group, but eleven other families of flowering plants are recorded as hosts. The apparent non-utilization of Urticaceae by H. misippus (L.), the most widespread species of the group, is discussed in terms of its possible ecological and phylogenetic relationships with H. bolina (L.). A check-list of butterflies feeding on Ipomoea (a host plant for both H. bolina and H. misippus) , and on other Convolvulaceae, is included.     

Rapid diversification associated with ecological specialization in Neotropical Adelpha butterflies

Rapid diversification is often associated with morphological or ecological adaptations that allow organisms to radiate into novel niches. Neotropical Adelpha butterflies, which comprise over 200 species and subspecies, are characterized by extraordinary breadth in host plant use and wing colour patterns compared to their closest relatives. To examine the relationship between phenotypic and species diversification, we reconstructed the phylogenetic history of Adelpha and its temperate sister genus Limenitis using genomewide restriction‐site‐associated DNA (RAD) sequencing. Despite a declining fraction of shared markers with increasing evolutionary distance, the RAD‐Seq data consistently generated well‐supported trees using a variety of phylogenetic methods. These well‐resolved phylogenies allow the identification of an ecologically important relationship with a toxic host plant family, as well as the confirmation of widespread, convergent wing pattern mimicry throughout the genus. Taken together, our results support the hypothesis that evolutionary innovations in both larvae and adults have permitted the colonization of novel host plants and fuelled adaptive diversification within this large butterfly radiation.     

Molecular phylogenetics and biogeography of the Neocellia Series of Anopheles mosquitoes in the Oriental Region

Molecular studies of population divergence and speciation across the Oriental Region are sparse, despite the region’s high biodiversity and extensive Pliocene and Pleistocene environmental change. A molecular phylogenetic study of the Neocellia Series of Anopheles mosquitoes was undertaken to identify patterns of diversification across the Oriental Region and to infer the role of Pleistocene and Pliocene climatic change. A robust phylogeny was constructed using CO2 and ND5 mitochondrial genes and ITS2 and D3 nuclear ribosomal markers. Bayesian analysis of mitochondrial genes was used to date divergence events. The repeated contraction and expansion of forest habitat resulting from Pleistocene climatic fluctuations appears to have had a substantial impact on intraspecific diversification, but has not driven speciation within this group. Primarily early to mid Pliocene speciation was detected within the Annularis Group, whereas speciation within the Maculatus and Jamesii Groups occurred during the mid and late Pliocene. Both allopatric divergence driven by late Pliocene environmental changes and ecological adaptation, involving altitudinal replacement and seasonality, are likely to have influenced speciation in the Maculatus Group.