Speciational history of North American Haemorhous finches (Aves: Fringillidae) inferred from multilocus data

We investigated species relationships and timing of speciation in North American Haemorhous finches by using a mitochondrial phylogeographic approach combined with a multilocus species tree reconstruction. Haemorhous purpureus and H. cassinii were strongly supported as sister taxa, and H. mexicanus was sister to H. purpureus + H. cassinii. Our divergence times indicated that diversification within Haemorhous occurred progressively from the Late Miocene into the Pleistocene. Our inferred pattern of speciation demonstrates the complexity of the origins of North American birds, and provides additional evidence that a single cause for speciation in closely related North American birds, such as Late Pleistocene glacial–interglacial cycles, is unlikely.     

Lineage divergence and speciation in the Web‐toed Salamanders (Plethodontidae: Hydromantes) of the Sierra Nevada,California

Peripatric speciation and the importance of founder effects have long been controversial, and multilocus sequence data and coalescent methods now allow hypotheses of peripatric speciation to be tested in a rigorous manner. Using a multilocus phylogeographical data set for two species of salamanders (genus Hydromantes) from the Sierra Nevada of California, hypotheses of recent divergence by peripatric speciation and older, allopatric divergence were tested. Phylogeographical analysis revealed two divergent lineages within Hydromantes platycephalus, which were estimated to have diverged in the Pliocene. By contrast, a low‐elevation species, Hydromantes brunus, diverged from within the northern lineage of H. platycephalus much more recently (mid‐Pleistocene), during a time of major climatic change in the Sierra Nevada. Multilocus species tree estimation and coalescent estimates of divergence time, migration rate, and growth rate reject a scenario of ancient speciation of H. brunus with subsequent gene flow and introgression from H. platycephalus, instead supporting a more recent divergence with population expansion. Although the small, peripheral distribution of H. brunus suggests the possibility of peripatric speciation, the estimated founding population size of the species was too large to have allowed founder effects to be important in its divergence. These results provide evidence for both recent speciation, most likely tied to the climatic changes of the Pleistocene, and older lineage divergence, possibly due to geological events, and add to evidence that Pleistocene glacial cycles were an important driver of diversification in the Sierra Nevada.     

Evolutionary history of the scorpionfly Dicerapanorpa magna (Mecoptera,Panorpidae)

Climate changes can have fundamental impacts on the distributional patterns of montane species, and range shifts frequently lead to allopatric divergence followed by the establishment of secondary contact zones. Many European and North American organisms have retreated to southern refugia during glacial periods and colonized northward during postglacial periods, but little is known about the evolutionary response of cold‐adapted insects to Pleistocene climate changes in eastern Asia. The scorpionfly Dicerapanorpa magna (Chou), with cold temperate habitat preference and weak dispersal ability, provides a good model system to explore how climate changes have influenced the distribution and divergence of cold‐adapted insects in eastern Asia. This study reconstructed the demographic dynamics and evolutionary history of D. magna with phylogeographic approaches, and predicted the species’ suitable areas under the Last Glacial Maximum (LGM) and current scenarios with the ecological niche modelling analysis. The mitochondrial cytochrome c oxidase subunit I resolved three phylogenetic lineages in D. magna dating back to Pleistocene, corresponding well with the geographically isolated Qinling, Bashan and Minshan Mountains. The ecological niche modelling recovered the suitable habitats for D. magna were the Qinling and Bashan Mountains under LGM and current conditions. The three lineages of D. magna might be in a process of incipient speciation, and likely derived their current distribution from separate glacial origins, followed by vicariance and divergence.     

CALIBRATING DIVERGENCE TIMES ON SPECIES TREES VERSUS GENE TREES: IMPLICATIONS FOR SPECIATION HISTORY OF APHELOCOMA JAYS

Estimates of the timing of divergence are central to testing the underlying causes of speciation. Relaxed molecular clocks and fossil calibration have improved these estimates; however, these advances are implemented in the context of gene trees, which can overestimate divergence times. Here we couple recent innovations for dating speciation events with the analytical power of species trees, where multilocus data are considered in a coalescent context. Divergence times are estimated in the bird genus Aphelocoma to test whether speciation in these jays coincided with mountain uplift or glacial cycles. Gene trees and species trees show general agreement that diversification began in the Miocene amid mountain uplift. However, dates from the multilocus species tree are more recent, occurring predominately in the Pleistocene, consistent with theory that divergence times can be significantly overestimated with gene‐tree based approaches that do not correct for genetic divergence that predates speciation. In addition to coalescent stochasticity, Haldane's rule could account for some differences in timing estimates between mitochondrial DNA and nuclear genes. By incorporating a fossil calibration applied to the species tree, in addition to the process of gene lineage coalescence, the present approach provides a more biologically realistic framework for dating speciation events, and hence for testing the links between diversification and specific biogeographic and geologic events.     

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.     

A multilocus perspective on the speciation history of a North American aridland toad (Anaxyrus punctatus)

Interpretations of phylogeographic patterns can change when analyses shift from single gene-tree to multilocus coalescent analyses. Using multilocus coalescent approaches, a species tree and divergence times can be estimated from a set of gene trees while accounting for gene-tree stochasticity. We utilized the conceptual strengths of a multilocus coalescent approach coupled with complete range-wide sampling to examine the speciation history of a broadly distributed, North American warm-desert toad, Anaxyrus punctatus. Phylogenetic analyses provided strong support for three major lineages within A. punctatus. Each lineage broadly corresponded to one of three desert regions. Early speciation in A. punctatus appeared linked to late Miocene-Pliocene development of the Baja California peninsula. This event was likely followed by a Pleistocene divergence associated with the separation of the Chihuahuan and Sonoran Deserts. Our multilocus coalescent-based reconstruction provides an informative contrast to previous single gene-tree estimates of the evolutionary history of A. punctatus.     

Shifting distributions and speciation: species divergence during rapid climate change

Questions about how shifting distributions contribute to species diversification remain virtually without answer, even though rapid climate change during the Pleistocene clearly impacted genetic variation within many species. One factor that has prevented this question from being adequately addressed is the lack of precision associated with estimates of species divergence made from a single genetic locus and without incorporating processes that are biologically important as populations diverge. Analysis of DNA sequences from multiple variable loci in a coalescent framework that (i) corrects for gene divergence pre-dating speciation, and (ii) derives divergence-time estimates without making a priori assumptions about the processes underlying patterns of incomplete lineage sorting between species (i.e. allows for the possibility of gene flow during speciation), is critical to overcoming the inherent logistical and analytical difficulties of inferring the timing and mode of speciation during the dynamic Pleistocene. Estimates of species divergence that ignore these processes, use single locus data, or do both can dramatically overestimate species divergence. For example, using a coalescent approach with data from six loci, the divergence between two species of montane Melanoplus grasshoppers is estimated at between 200,000 and 300,000 years before present, far more recently than divergence estimates made using single-locus data or without the incorporation of population-level processes. Melanoplus grasshoppers radiated in the sky islands of the Rocky Mountains, and the analysis of divergence between these species suggests that the isolation of populations in multiple glacial refugia was an important factor in promoting speciation. Furthermore, the low estimates of gene flow between the species indicate that reproductive isolation must have evolved rapidly for the incipient species boundaries to be maintained through the subsequent glacial periods and shifts in species distributions.     

Metabolic Rate and Climatic Fluctuations Shape Continental Wide Pattern of Genetic Divergence and Biodiversity in Fishes

Taxonomically exhaustive and continent wide patterns of genetic divergence within and between species have rarely been described and the underlying evolutionary causes shaping biodiversity distribution remain contentious. Here, we show that geographic patterns of intraspecific and interspecific genetic divergence among nearly all of the North American freshwater fish species (>750 species) support a dual role involving both the late Pliocene-Pleistocene climatic fluctuations and metabolic rate in determining latitudinal gradients of genetic divergence and very likely influencing speciation rates. Results indicate that the recurrent glacial cycles caused global reduction in intraspecific diversity, interspecific genetic divergence, and species richness at higher latitudes. At the opposite, longer geographic isolation, higher metabolic rate increasing substitution rate and possibly the rapid accumulation of genetic incompatibilities, led to an increasing biodiversity towards lower latitudes. This indicates that both intrinsic and extrinsic factors similarly affect micro and macro evolutionary processes shaping global patterns of biodiversity distribution. These results also indicate that factors favouring allopatric speciation are the main drivers underlying the diversification of North American freshwater fishes.     

Glacial cycles as an allopatric speciation pump in north‐eastern American freshwater fishes

Allopatric speciation may be the principal mechanism generating new species. Yet, it remains difficult to judge the generality of this process because few studies have provided evidence that geographic isolation has triggered the development of reproductive isolation over multiple species of a regional fauna. Here, we first combine results from new empirical data sets (7 taxa) and published literature (9 taxa) to show that the eastern Great Lakes drainage represents a multispecies suture zone for glacial lineages of freshwater fishes with variable levels of genetic divergence. Second, we performed amplified fragment length polymorphism analyses among four pairs of lineages. Results indicate that lineages with relatively deep levels of mtDNA 5′ COI (barcode) sequence divergence (>2%) developed strong reproductive barriers, while lineages with lower levels of divergence show weaker reproductive isolation when found in sympatry. This suggests that a threshold of 2% sequence divergence at mtDNA could be used as a first step to flag cryptic species in North American freshwater fishes. By describing different levels of divergence and reproductive isolation in different co‐occurring fishes, we offer strong evidence that allopatric speciation has contributed significantly to the diversification of north‐eastern American freshwater fishes and confirm that Pleistocene glacial cycles can be viewed as a ‘speciation pump’ that played a predominant role in generating biodiversity.     

Nitric oxide (NO): a key player in the senescence of Medicago truncatula root nodules

The ectomycorrhizal fungus Tricholoma populinum is host-specific with Populus species. T. populinum has wind-dispersed progagules and may be capable of long-distance dispersal. In this study, we tested the hypothesis of a panmictic population between Scandinavia and North America. DNA sequences from five nuclear loci were used to assess phylogeographic structure and nucleotide divergence between continents. Tricholoma populinum was composed of Scandinavian and North American lineages with complete absence of shared haplotypes and only one shared nucleotide mutation. Divergence of these lineages was estimated at approx. 1.7-1.0 million yr ago (Ma), which occurred after the estimated divergence of host species Populus tremula and Populus balsamifera/Populus trichocarpa at 5 Ma. Phylogeographic structure was not observed within Scandinavian or North American lineages of T. populinum. Intercontinental divergence appears to have resulted from either allopatric isolation; a recent, rare long-distance dispersal founding event followed by genetic drift; or the response in an obligate mycorrhizal fungus with a narrow host range to contractions and expansion of host distribution during glacial and interglacial episodes within continents. Understanding present genetic variation in populations is important for predicting how obligate symbiotic fungi will adapt to present and future changing climatic conditions.     

Molecular evidence for Pleistocene glacial cycles driving diversification of a North American desert spider, Agelenopsis aperta

The influence of historical climatic vs. geological changes on species diversification patterns was investigated in a widely distributed North American desert spider, Agelenopsis aperta (Araneae: Agelenidae), with particular reference to Pleistocene glacial cycles and earlier patterns of mountain building. Levels of sequence divergence obtained from the mitochondrial gene, cytochrome oxidase I, dated to the Pleistocene, eliminating Rocky Mountain orogeny as a cause of diversification, as orogeny ended 4 million years ago. The results of phylogenetic and network analyses showed the presence of three geographically defined clades, which were consistent with the presence of at least three glacial refugia: (i) east of the Rocky Mountains; (ii) between the Rocky Mountains and Sierra Nevadas; and (iii) west of the Sierra Nevadas. In addition, populations within the Rocky Mountains exhibited significantly lower genetic diversity than populations east of the Rocky Mountains and the haplotypes found within the Rockies were a subset of eastern haplotypes. These patterns suggest that a post-Pleistocene range expansion occurred out of an eastern glacial refugium into the Rocky Mountains. Examination of phylogeographical studies of other North American desert taxa indicated that mountain building explained diversification patterns more effectively for some taxa but Pleistocene climate change was more important for others, including A. aperta.     

Comparative Nucleotide Diversity Across North American and European Populus Species

Nucleotide polymorphisms in two North American balsam poplars (Populus trichocarpa Torr. & Gray and P. balsamifera L.; section Tacamahaca), and one Eurasian aspen (P. tremula L.; section Populus) were compared using nine loci involved in defense, stress response, photoperiodism, freezing tolerance, and housekeeping. Nucleotide diversity varied among species and was highest for P. tremula (θ(w) = 0.005, π(T) = 0.007) as compared to P. balsamifera (θ(w) = 0.004, π(T) = 0.005) or P. trichocarpa (θ(w) = 0.002, π(T) = 0.003). Across species, the defense and the stress response loci accounted for the majority of the observed level of nucleotide diversity. In general, the studied loci did not deviate from neutral expectation either at the individual locus (non-significant normalized Fay and Wu's H) or at the multi-locus level (non-significant HKA test). Using molecular clock analysis, section Tacamahaca probably shared a common ancestor with section Populus approximately 4.5 million year ago. Divergence between the two closely related balsam poplars was about 0.8 million years ago, a pattern consistent with an isolation-with-migration (IM) model. As expected, P. tremula showed a five-fold higher substitution rate (2 × 10(-8) substitution/site/year) compared to the North American species (0.4 × 10(-8) substitution/site/year), probably reflecting its complex demographic history. Linkage disequilibrium (LD) varied among species with a more rapid decay in the North American species (<400 bp) in comparison to P. tremula (?400 bp). The similarities in nucleotide diversity pattern and LD decay of the two balsam poplar species likely reflects the recent time of their divergence.     

Phylogenetic assemblage structure of North American trees is more strongly shaped by glacial–interglacial climate variability in gymnosperms than in angiosperms

How fast does biodiversity respond to climate change? The relationship of past and current climate with phylogenetic assemblage structure helps us to understand this question. Studies of angiosperm tree diversity in North America have already suggested effects of current water–energy balance and tropical niche conservatism. However, the role of glacial–interglacial climate variability remains to be determined, and little is known about any of these relationships for gymnosperms. Moreover, phylogenetic endemism, the concentration of unique lineages in restricted ranges, may also be related to glacial–interglacial climate variability and needs more attention. We used a refined phylogeny of both angiosperms and gymnosperms to map phylogenetic diversity, clustering and endemism of North American trees in 100‐km grid cells, and climate change velocity since Last Glacial Maximum together with postglacial accessibility to recolonization to quantify glacial–interglacial climate variability. We found: (1) Current climate is the dominant factor explaining the overall patterns, with more clustered angiosperm assemblages toward lower temperature, consistent with tropical niche conservatism. (2) Long‐term climate stability is associated with higher angiosperm endemism, while higher postglacial accessibility is linked to to more phylogenetic clustering and endemism in gymnosperms. (3) Factors linked to glacial–interglacial climate change have stronger effects on gymnosperms than on angiosperms. These results suggest that paleoclimate legacies supplement current climate in shaping phylogenetic patterns in North American trees, and especially so for gymnosperms.     

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.     

EFFECTS OF CLIMATIC AND GEOLOGICAL PROCESSES DURING THE PLEISTOCENE ON THE EVOLUTIONARY HISTORY OF THE NORTHERN CAVEFISH,AMBLYOPSIS SPELAEA (TELEOSTEI: AMBLYOPSIDAE)

Climatic and geological processes associated with glaciation cycles during the Pleistocene have been implicated in influencing patterns of genetic variation and promoting speciation of temperate flora and fauna. However, determining the factors promoting divergence and speciation is often difficult in many groups because of our limited understanding of potential vicariant barriers and connectivity between populations. Pleistocene glacial cycles are thought to have significantly influenced the distribution and diversity of subterranean invertebrates; however, impacts on subterranean aquatic vertebrates are less clear. We employed several hypothesis‐driven approaches to assess the impacts of Pleistocene climatic and geological changes on the Northern Cavefish, Amblyopsis spelaea, whose current distribution occurs near the southern extent of glacial advances in North America. Our results show that the modern Ohio River has been a significant barrier to dispersal and is correlated with patterns of genetic divergence. We infer that populations were isolated in two refugia located north and south of the Ohio River during the most recent two glacial cycles with evidence of demographic expansion in the northern isolate. Finally, we conclude that climatic and geological processes have resulted in the formation of cryptic forms and advocate recognition of two distinct phylogenetic lineages currently recognized as A. spelaea.     

New mitochondrial DNA data affirm the importance of Pleistocene speciation in North American birds

The timing of origin of modern North American bird species in relation to Pleistocene glaciations has long been the topic of significant discussion and disagreement. Recently, Klicka and Zink (1997) and Avise and Walker (1998) enlivened this debate by using calibrated molecular distance values to estimate timing of speciations. Here we use new molecular studies to test their conclusions. Molecular distance values for 39 pairs of proven sister species, 27 of which are based on new data, alter the currently perceived pattern that avian species splits occurred mainly in the Pliocene and early-mid-Pleistocene. Mitochondrial DNA divergence values for this set of taxa showed a skewed distribution pointing toward relatively young speciation times, in contrast to the pattern presented by Klicka and Zink (1997) for 35 sister plus non-sister species pairs. Our pattern was not significantly different from that of Avise and Walker (1998) for "intraspecific phylogroups," some of which are species. We conclude that the entire Pleistocene, including the last two glacial cycles (<250,000 years ago), was important in speciations of modern North American birds. A substantial number of speciations were both initiated and completed in the last 250,000 years. Simultaneously, many taxa began to diverge in the Pleistocene but their speciations are not yet complete (per Avise and Walker 1998). The suggestion that durations of speciations average two million years is probably a substantial overestimate.     

Quaternary refugia are associated with higher speciation rates in mammalian faunas of the Western Palaearctic

Quaternary climate changes have contributed to shape the biogeographic distribution of extant species. The combination of climatic niche conservatism and glacial‐interglacial cycles forced many species to retract their range limits for surviving under the advance of Pleistocene ice‐sheets. Refugia offered geographical opportunities for species to retreat, persist and, later on, begin recolonization processes under favourable environmental conditions. Here we explore the hypothesis that refugia have been not only crucial for the survival of multiple species but also acted as speciation centres for Western Palaearctic mammals. We define ‘recurrent massive refugia’ as those geographical regions that have historically accumulated the highest levels of co‐occurring species for several Quaternary cycles. Our assemblage‐level analyses identify the existence of refugia within the Iberian Peninsula and the Atlantic French margin that were recurrently selected by most mammals. The topographic heterogeneity, climatic stability and microhabitat availability of these refugial areas may have offered suitable habitat conditions for multiple species during different climatic events over time. Using a Bayesian analysis of macroevolutionary mixtures we detected that the higher level of divergence and accumulative evolutionary changes in mammals of the Western Palaearctic are found in refugia. The continuous retractions and expansions of species’ ranges during the Pleistocene promoted temporal changes in the composition and richness of communities in this biogeographic region. The reorganization of ecological composition driven by cyclical climatic events may have favoured the emergence of biotic interactions and ecological responses conducive to novel selective pressures. Our findings suggest, first, that multiple climatic changes in the form of glacial‐interglacial transitions during the Quaternary have left a detectable imprint on the observed geographical patterns of species richness in mammalian faunas of the Western Palaearctic, and second, highlight the importance of refugia for the preservation of species (‘museums’) and as centers of speciation and endemism (‘cradles’) as well.     

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia,multidirectional expansion,and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range‐wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north‐to‐south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial–interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.     

Division within the North American boreal forest: Ecological niche divergence between the Bicknell's Thrush (Catharus bicknelli) and Gray‐cheeked Thrush (C. minimus)

Sister species that diverged in allopatry in similar environments are expected to exhibit niche conservatism. Using ecological niche modeling and a multivariate analysis of climate and habitat data, I test the hypothesis that the Bicknell's Thrush (Catharus bicknelli) and Gray‐cheeked Thrush (C. mimimus), sister species that breed in the North American boreal forest, show niche conservatism. Three tree species that are important components of breeding territories of both thrush species were combined with climatic variables to create niche models consisting of abiotic and biotic components. Abiotic‐only, abiotic+biotic, and biotic‐only models were evaluated using the area under the curve (AUC) criterion. Abiotic+biotic models had higher AUC scores and did not over‐project thrush distributions compared to abiotic‐only or biotic‐only models. From the abiotic+biotic models, I tested for niche conservatism or divergence by accounting for the differences in the availability of niche components by calculating (1) niche overlap from ecological niche models and (2) mean niche differences of environmental values at occurrence points. Niche background similarity tests revealed significant niche divergence in 10 of 12 comparisons, and multivariate tests revealed niche divergence along 2 of 3 niche axes. The Bicknell's Thrush breeds in warmer and wetter regions with a high abundance of balsam fir (Abies balsamea), whereas Gray‐cheeked Thrush often co‐occurs with black spruce (Picea mariana). Niche divergence, rather than conservatism, was the predominant pattern for these species, suggesting that ecological divergence has played a role in the speciation of the Bicknell's Thrush and Gray‐cheeked Thrush. Furthermore, because niche models were improved by the incorporation of biotic variables, this study validates the inclusion of relevant biotic factors in ecological niche modeling to increase model accuracy.