Skip to main content

Advertisement

The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico

Abstract

Background

Mesoamerica is a remarkable region with a high geological and ecological complexity. Within northern Mesoamerica, the biotic province of the Sierra Madre del Sur (SMS) in southwestern Mexico harbors exceptionally high avian endemism and diversity. Herein, we searched for spatially and temporally concordant phylogeographic patterns, in four bird genera from three distinct avian orders co-distributed across Mesoamerica and investigated their causes through hypothesis testing regarding historical processes. Selected species include endemic and differentiated populations across the montane forests of Mesoamerica, and particularly within the SMS.

Results

We gathered mitochondrial DNA sequences for at least one locus from 177 individuals across all species. We assessed genetic structure, demographic history, and defined a framework for the coalescent simulations used in biogeographic hypothesis testing temporal and spatial co-variance. Our analyses suggested shared phylogeographic breaks in areas corresponding to the SMS populations, and between the main montane systems in Mesoamerica, with the Central Valley of Oaxaca and the Nicaragua Depression being the most frequently shared breaks among analyzed taxa. Nevertheless, dating analyses and divergence patterns observed were consistent with the hypothesis of broad vicariance across Mesoamerica derived from mechanisms operating at distinct times across taxa in the SMS.

Conclusions

Our study provides a framework for understanding the evolutionary origins and historical factors enhancing speciation in well-defined regions within Mesoamerica, indicating that the evolutionary history of extant biota inhabiting montane forests is complex and often idiosyncratic.

Background

Comparison of phylogeographic patterns across co-distributed species provides valuable insights on the probable congruence of processes that may have driven intraspecific diversification in particular regions [1,2,3,4]. Several phylogeographic studies of individual bird species have shown that genetic structure is frequently associated to discontinuous ranges and geographical barriers that restrict admixture between populations, leading to genetic differentiation (see [5,6,7,8]) which may have proceeded independently or jointly with the divergence in morphological, behavioral or ecological traits (see [9,10,11,12,13]). However, a relatively lesser effort has been dedicated to provide robust biogeographic hypotheses seeking to explain shared distributional patterns, even when it may be expected that species sharing geographic areas should show congruent spatio-temporal patterns of differentiation [1, 14,15,16,17,18,19].

Situated at the northernmost Neotropics, Mesoamerica possess one of the highest levels of endemism and species diversity, but also one of the most globally threatened biota and high rates of deforestation [20, 21]. The complex geological history of this region, as well as the cyclic changes in vegetation and climate, pose it as a challenging area for biogeographic and evolutionary studies; in addition, constant orogenic processes have promoted a highly broken topography characterized by highland isolated patches of humid montane forest between 600 and 3000 m, which includes both humid pine-oak forest and cloud forest [22,23,24]. This mosaic-like landscape has been associated to centers of diversification along elevational gradients, in which both a high species richness and endemism have evolved for the last 2 Myr, likely as a result of both environmental and geological complexity, as well as Pleistocene climatic fluctuations [6, 15, 18, 22, 25,26,27,28,29], which is supported by relatively recent intraspecific differentiation processes in several groups of organisms (see [27, 30]), therefore explaining the existence of numerous endemic species in different taxonomic groups, including birds [31,32,33,34,35].

The Sierra Madre del Sur (SMS) biogeographic province is included within the Mesoamerican highlands, in the Mexican Transition Zone (MTZ, [36]). This isolated mountain range is mostly surrounded by dry lowlands of the Isthmus of Tehuantepec (IT) to the east, the Pacific slope to the south and west, the Balsas Depression and the central Trans-Mexican Volcanic Belt to the north (Fig. 1 [38,39,40,41];). The SMS with an overall extension of ca. 56,729 km2, is located along the Mexican Pacific slope (Jalisco to Oaxaca, 16–18°N, 95–102°W) [42]. The region hosts 622 avian species, 29 (4.6%) of which are classified as semi-endemic (species breeding elsewhere but for which the entire species’ population spends the wintering season within the political limits of Mexico [43]), 15 (2.4%) as quasi-endemic (species whose geographic range extends no more than 35,000 km2 to a neighboring country due to habitat continuity, [43]), and 54 (8.6%) are endemic to Mexico; in addition, 134 (21.5%) species are considered in some category of risk [34]. Furthermore, as in other Mesoamerican bird taxa [44,45,46,47,48], the SMS holds endemic and well-differentiated subspecies, suggesting the importance of isolation for the evolution of intraspecific variation [17, 49].

Fig. 1
figure1

a Map of Mesoamerica showing the distribution of montane regions for which genetic data were analyzed in the present study, and summary of the main geographic groups in which they were divided; based on [37] and Natural Earth data (http://naturalearthdata.com). Highlighted in orange is the SMS region and main geographical breaks are shown with red lines. Evaluated barriers in this study are numbered as follows: 1) Balsas Depression, 2) Río Verde drainage, 3) Oaxaca’s Central Valley, 4) Isthmus of Tehuantepec, 5) Nicaragua Depression, and 6) Panama Channel Fault Zone (PCFZ). b Sierra Madre del Sur mountain range (SMS, ca. 56,729 km2) highlighted in orange, showing the broken topography and main mountain peaks. The red dot depicts Oaxaca’s Central Valley, whereas the red line depicts the Río Verde drainage

Isolated montane forest patches are generally thought to act as promoters for divergence between populations of several species, which have served as models for the understanding of speciation processes in island-like settings [50,51,52,53,54]. These fragmented habitats influence the population dynamics and demography, given that resident populations typically show reduced effective population sizes [18, 50, 55,56,57,58,59,60], they are prone to a reduced adaptive potential and increased extinction risk due to the stochastic effects of genetic drift [61,62,63,64,65,66,67]. Therefore, the study of evolutionary history of co-distributed species in these habitats is highly important for conservation purposes.

Herein, we compared phylogeographical structure patterns among four co-distributed and not closely-related Neotropical bird taxa with well-differentiated populations restricted to the SMS. We tested the following biogeographic hypotheses in an ABC framework: (a) a southern origin (southern Central America) followed by a northward expansion; (b) a Mexican origin and range expansion towards southern Central America; and (c) a northern Central American origin with north and southward dispersal. We also aimed to determine whether phylogeographic structure and divergence times in each taxon are congruent among taxa, therefore revealing patterns of simultaneous diversification. Specifically, we used mitochondrial datasets in a comparative phylogeography framework to test: (1) The presence of geographically-structured genetic variation; (2) the existence of spatially-congruent genetic variation across analyzed species; (3) if biogeographic history and similar diversification scenarios support vicariance, dispersal or population admixture. In addition, we also examined temporal patterns of diversification through hABC simulations, testing for either synchronous or multiple pulses for diversification in bird assemblages across the SMS region.

Results

Phylogenetic and population structure analyses

All of our phylogenetic and genetic assignment analyses supported previous findings on the genetic structure for all of the taxa (Fig. 2 and Fig. 3). Genetic structure is clearly delimited by landscape breaks corresponding to the main montane regions in Mexico and Central America. Nevertheless, we found that the detected genetic structure in these taxa is not coincident with the accepted taxonomy at the intraspecific level.

Fig. 2
figure2

Geographical distribution of sampled haplotypes for (a) Aulacorhynchus and (b) Chlorospingus. c and d depict median-joining network of concatenated mtDNA loci for Aulacorhynchus and Chlorospingus, numbers on lines depict mutational steps between haplotypes, gray dots represent median vectors inferred for the data. e and f depict maximum clade credibility trees using BEAST with branch support (BI/bootstrap). Values between brackets indicate the 95% highest posterior densities (HDP) of the estimated times of divergence events (in Myr). Nodes that have no represented time frame of diversification are those whose lower and upper bounds of the HDP interval had posterior probabilities inferior to 0.5

Fig. 3
figure3

Geographical distribution of sampled haplotypes for (a) Cardellina and (b) Eupherusa species. c and d depict median-joining network of concatenated mtDNA loci for Cardellina and Eupherusa respectively, numbers on lines depict mutational steps between haplotypes, gray dots represent median vectors inferred for the data. e and f depict maximum clade credibility trees using BEAST with branch support (BI/bootstrap). Values between brackets indicate the 95% highest posterior densities (HDP) of the estimated times of divergence events (in Myr). Nodes that have no represented time frame of diversification are those whose lower and upper bounds of the HDP interval had posterior probabilities inferior to 0.5

We found four distinctive lineages in Aulacorhynchus. Two of these, SMS (restricted to the SMS highlands) and EMx-NCA (widespread in Mesoamerica, including SMO, TM, Chiapas, and northern Central America) are differentiated by 27 mutational steps (Fig. 2c). The other two are distributed in southern Central America from Costa Rica to Panama. Phylogenetic analyses support this structure with ML bootstrap values and BI posterior probabilities above 80 and 0.9 respectively (Fig. 2e). We further observed two subclades in SMS samples, each representing Oaxaca and Guerrero populations.

For Chlorospingus, we observed six distinct mtDNA lineages, four of which included Mexican samples (Fig. 2d). Individuals representing the SMS (Sierra de Miahuatlán, Oaxaca and Guerrero) were recovered in a separated group from samples from northern Oaxaca, which clustered with samples from the eastern Sierra Madre Oriental highlands (SMO); both groups are separated by 26 mutational steps. A third group corresponded to the Chimalapas (Oax.) and northern Chiapas region, while the fourth group corresponded exclusively to the TM, separated from the SMS by 32 mutations. This genetic structure is supported by BI, as well as by genetic assignment analyses with overall high posterior probabilities. Some nodes however, showed low bootstrap and posterior probability values (Fig. 2f). As observed in Aulacorhynchus, samples from the states of Oaxaca and Guerrero were recovered as distinct subclades within the SMS.

In Cardellina, we found three lineages representing C. rubra separated by 11 to 24 mutational steps (Fig. 3c). Cardellina versicolor, the sister species, is differentiated by more than 30 mutational steps from the C. rubra group. Within C. rubra, all recovered lineages coincide with previous findings, with one group found in the SMOc, another in the TMVB, and a third one in the SMS. We found 8 samples from a locality in central Guerrero (Carrizal de Bravo) from the SMS to be grouped in the TMVB clade, which also includes samples from localities of Jalisco and Michoacán. This paraphyletic relationship was recovered in all analyses, receiving strong support from both bootstrap (> 80) and Bayesian posterior probabilities (1; Fig. 3e). Nevertheless, genetic assignment analysis yielded a distinct result for these Guerrero samples, particularly for three of them (jk04–349, jk04–351, and jk04–353), which clustered with the SMS group; the remaining five clustered with the TMVB group.

Finally, in Eupherusa-Thalurania, we found five lineages corresponding to the currently recognized species. T. ridgwayi group is separated by 39 mutational steps from all of the Eupherusa samples; species from the SMS E. cyanoprhys and E. poliocerca haplotypes are separated by 22 mutational steps; and finally, species E. nigriventris and E. eximia are separated by 50 mutations (Fig. 3d). Topologies from BI, ML (Fig. 3f) recovered high posterior probabilities and bootstrap support values (> 0.8, > 75 respectively). All analyses recovered a monophyletic Eupherusa genus, where the SMS is represented by E. poliocerca and its sister species E. cyanophrys, and a clade from eastern Mexico to Central America, composed of E. eximia sister to E. nigriventris. Finally, for each of these species, we did not find evidence supporting further intraspecific genetic structure, even when their distribution, as in E. poliocerca (Oaxaca-Guerrero) and E. eximia (E. e. nelsoni, E. e. eximia, and E. e. egregia), has been described to be disjunct [28, 68, 69].

Hierarchical AMOVAs (Table 1) for the analyzed species showed that variation was explained by differences among groups. In this case, the fixation indices among groups (FCT) were not significant, whereas fixation indices among populations within groups (FSC) were. Genetic diversity varied among datasets. Overall, we obtained high haplotype diversity (Hd > 0.5) and low nucleotide diversity (Table 2). The only low value of Hd was observed in E. cyanophrys. Fu’s FS and Tajima’s D values were not significant, excepting TMVB and SMS in Cardellina; however, we observed a tendency towards negative values (Table 2).

Table 1 Hierarchical AMOVAs for each matrix analyzed showing as sources of variation FCT, FSC and FST respectively
Table 2 Summary statistics of each analyzed species populations, based on K groups observed in phylogenetic and Geneland analyses, using concatenated mtDNA. N = number of sequences, Hd = haplotype diversity, π = nucleotide diversity. Fu’s FS and Tajima’s D values with an asterisk (*) represent significant values for their respective confidence intervals (< 0.02 and <  0.05). Columns with “n. a.” represent null genetic variation owing to the presence of only one haplotype. Region abbreviations are as follows: TMVB = Trans Mexican Volcanic Belt, SMS = Sierra Madre del Sur, SMOc = Sierra Madre Occidental, SMO = Sierra Madre Oriental, NCA = northern Central America, SCA = southern Central America, EMx-NCA = eastern Mexico-northern Central America, CR-NP = Costa Rica-northern Panama, SP = southern Panama

Overall, pairwise FST values showed significant high genetic differentiation between analyzed groups. However, probably due to small sample sizes, some comparisons showed no significant values, as in the SP cluster of Aulacorhynchus (Table 3), as well as in comparisons between E. nigriventris vs T. ridgwayi, and E. nigriventris vs E. cyanophrys clusters (Table 4). The lowest FST value obtained occurred between Cardellina clusters from the TMVB vs SMS (Table 4).

Table 3 Pairwise FST fixation index used to measure genetic structuring in Aulacorhynchus and Chlorospingus. Abbreviations correspond to defined geographical groups: EMx-NCA (eastern Mexico-northern Central America), CR-NP (Costa Rica-northern Panama), SP (southern Panama), SMS (Sierra Madre del Sur), SMO (Sierra Madre Oriental), NCA (northern Central America), SCA (southern Central America). An asterisk (*) represents significant values (< 0.05)
Table 4 Pairwise FST fixation index used to measure genetic structuring in Cardellina and Eupherusa. Abbreviations correspond to defined geographical groups: SMOc (Sierra Madre Occidental), TMVB (Trans Mexican Volcanic Belt), SMS (Sierra Madre del Sur). An asterisk (*) represents significant values (< 0.05)

Divergence time estimation and coalescent-based estimation of population history

Even when divergence dates varied, analyses showed that divergence between populations occurred during the Pleistocene, mainly within the last 2 Myr. The best-supported scenario for four population groups in Aulacorhynchus was scenario 1 (Additional file 5: Figure S5), with a posterior probability (PP) of 0.57 and a 95% confidence interval (CI) of 0.5584–0.5816 (Fig. 4a), which suggested that SMS, EMx-NCA, and CR-NP populations arose from an ancestral SP population, following a northwards divergence pattern. Confidence in scenario choice tested through type I and II errors, suggested the selection of a highly accurate scenario (Additional file 3). Posterior parameter estimates (Table 5) indicated that initial divergence (t4) from southern Panama northwards across the lowlands and narrowest region in Central America (Panama Channel Fault Zone) started during the Pliocene at 4.3 Myr (95% CI: 4–4.8 Myr), followed by successive northward expansion and further vicariance of populations, and the final divergence between eastern (EMx-NCA) and southwestern (SMS) Mexican populations (t2) about 1.3 Myr (95% CI: 420 kyr – 2.5 Myr) during the Pleistocene. When the lineages of Guerrero and Oaxaca were considered as distinct populations (Additional file 5: Figure S6), the best scenario was scenario 2 (PP of 0.677, 95% CI: 0.668–0.687; Fig. 4b). Type I and II errors indicated moderate confidence in the selected scenario (Additional file 3). Under this scenario, diversification patterns were similar to the latter, but the most recent divergence event occurred between Oaxaca and Guerrero populations (Additional file 4 and Additional file 5). Mean posterior parameter estimates (Table 5) indicated that the first divergence (t4) event is dated around 3.3 Myr (95% CI: 2.5–6.1 Myr), while the most recent divergence time (t1) occurred about 204 kyr (95% CI: 12.8–730 kyr).

Fig. 4
figure4

Best-supported biogeographic scenarios of diversification in analyzed species. a scenario 1 of diversification for Aulacorhynchus. b Scenario 2 of diversification for Aulacorhynchus considering as distinct populations lineages from Oaxaca and Guerrero. c Scenario 1 of diversification of Chlorospingus species. d best-supported scenario in Chlorospingus considering lineages of Guerrero and Oaxaca as distinct populations. e Scenario 2 of diversification in Eupherusa species including T. ridgwayi. f Scenario 1 of diversification in Eupherusa excluding T. rigdwayi. g Scenario 3 of diversification in Cardellina species

Table 5 Mean time posterior parameter estimates for the best supported scenarios of analyzed species. Estimates are based on the 1% of simulated datasets closest to the observed values. Time (t) is set in years

The best scenario for Chlorospingus when six groups are considered (Additional file 5: Figure S7), was scenario 1 (PP = 0.519, CI: 0.4965–0.5427; Fig. 4c). Confidence in scenario choice was high given the values for type I and II errors (Additional file 3). This scenario suggested an initial divergence (t6) of Mexican and Central American lineages during the Pliocene at 5.8 Myr (95% CI: 2.5–9.5 Myr (Table 5), diversification of Mexican lineages proceeded westward, with the latter divergence event (t2) between NCA and SCA at 1.3 Myr (95% CI: 138 kyr – 2.7 Myr; Additional file 4 and Additional file 5). When lineages from Oaxaca and Guerrero were considered as distinct populations, scenario 1 (Additional file 5: Figure S8) received the better support (PP = 0.6, 95% CI: 0.5909–0.625; Fig. 4d), with high confidence (Additional file 3). Scenario 1 suggested diversification events similar to the latter, mean parameter estimates (Table 5) suggested an initial divergence time (t6) at 2 Myr (95% CI: 1.1–5.3 Myr) between Mexican and Central American lineages, a divergence event (t2) estimated to occur after the Nicaragua Depression formation between Central American populations at 1 Myr (95% CI: 288 kyr – 2.1 Myr), and a final divergence event (t1) in the SMS lineage originating Guerrero and Oaxaca populations during a interglacial period in the pre-Illinoian stage at 656 kyr (95% CI: 50 kyr – 1.6 Myr; Additional file 4 and Additional file 5).

Best-supported scenario in Cardellina was scenario 3 (PP = 0.4414, CI: 0.43–0.45; Fig. 4g), with high confidence (Additional file 3). Scenario 3 (Additional file 5: Figure S9) predicts an initial northward dispersal event from northern Central America (C. versicolor) population through the Isthmus of Tehuantepec at t4 (~ 1.3 Myr); followed by vicariance events driven by the effect of Mexican highlands at t3 (~ 426 kyr) between SMOc and an hypothetical ancestral population that at t2 (~ 260 kyr) is fragmented into the TMVB and SMS populations (Table 5). This scenario also supports a recent change in the TMVB effective population size during the Wisconsin Glacial stage at 40 kyr (95% CI: 3500 yr – 122 kyr; Additional file 4 and Additional file 5).

The best scenario for Eupherusa including T. ridgwayi (Additional file 5: Figure S10) was scenario 2 (PP = 0.712, 95% CI: 0.7–0.72; Fig. 3e). Type I and II errors indicate high confidence in scenario choice (Additional file 3). Scenario 2 suggested divergence of T. ridgwayi and a hypothetical ancestral population during the late Pliocene at (t4) 2.7 Myr (95% CI: 2.1–3.4 Myr), in time (t3) and before major climatic oscillations in the mid-Pleistocene, the ancestral Eupherusa lineage splits into two hypothetical ancestral lineages, one that further splits at t2 (~ 1.3 Myr) giving rise to the populations in eastern Mexico and Central America (E. eximiaE.nigriventris), and the second which splits into populations in the SMS (E. poliocerca-E. cyanophrys) at 716 kyr (95% CI: 123 kyr – 1.5 Myr; Table 5). When T. ridgwayi is excluded (Additional file 5: Figure S11), the best-supported was scenario 1 (PP = 0.868, 95% CI: 0.859–0.877; Fig. 4f), which is similar to the previous scenario. Confidence in scenario choice was high given low values of type I and II errors (Additional file 3). Mean parameter estimates of scenario 1 (Table 5) indicated an initial divergence at t3 (~ 1.7 Myr) of an ancestral Eupherusa lineage into two ancestral populations. One diverged into eastern Mexican and Central American populations at t2 (~ 1.2 Myr), whereas the other splits into southwestern Mexican populations at 688 kyr (95% CI: 116 kyr – 1.5 Myr; Additional file 4 and Additional file 5).

Test of simultaneous divergence

The hABC analyses including four population pairs spanning along two identified putative barriers (Balsas Depression and Oaxaca’s Valleys), showed a relatively strong posterior probability support for four independent divergence events (PPΨ = 4 = 0.417) and the highest mode model-averaged posterior estimate of the hyperparameter Ω = 0.1415 (95% HDP interval = 0.0035–0.2998), in comparison with posterior probability of two (PPΨ = 2 = 0.282) or three (PPΨ = 3 = 0.275) divergence events (Fig. 5a).

Fig. 5
figure5

Joint posterior probability for E(τ) (average divergence time) and Ω (variance in divergence times/average divergence time) for (a) Four population pairs spanning across southwestern Mexico and eastern-central Mexico through Central America. b Two population pairs within SMS. c Three population pairs within SMS

Subsequently, we tested subsets of the data set for simultaneous divergence. First we tested two taxon pairs within SMS (Guerrero and Oaxaca populations separated by the Río Verde drainage) of Aulacorhynchus and Chlorospingus, obtaining high support for a single divergence event (PPΨ = 1 = 0.6749), along with a mode value of Ω consistent with synchronous divergence (Ω = 0.000475) and an Ω = 0, which is included within the 95% HDP intervals (Fig. 5b), whereas when including the Eupherusa taxon pair (cyanophrys-poliocerca), we found low support for a single divergence event, suggesting two divergence events across these taxa (PPΨ = 2 = 0.5436) and a Ω mode of 0.05 consistent with asynchronous divergence (Fig. 5).

When evaluating southwestern SMS against eastern populations spanning Oaxaca’s Central Valley barrier, there was a relatively strong support for a single divergence event between these populations (PPΨ = 1 = 0.4536) and Ω = 0.0022, also indicating support for synchronous diversification between these taxa occurring at ~ 1.06 Myr (Fig. 6d). Finally, the evaluation of population pairs distributed along Central America revealed a very high posterior probability of one divergence event (PPΨ = 1 = 0.7228), and a Ω hyperparameter value = 0 (95% HDP interval = 0.0–0.1167) across the Nicaraguan Depression lowlands between Costa Rica and Nicaragua occurring ~ 0.924 Myr (Fig. 6e).

Fig. 6
figure6

Joint posterior probability for E(τ) (average divergence time) and Ω (variance in divergence times/average divergence time) for (d) Three population pairs spanning across southwestern Mexico and eastern Mexico through Central America. e Three population pairs in Central America

Discussion

The aim of comparative phylogeography is to detect concordant patterns among co-distributed species, based on the idea that they may share common histories [1, 14, 16]. Thus, finding similar genetic patterns among different species may suggest similar historical processes influencing major biological components in a given region. Nevertheless, as shown in the results presented herein, complete concordance is commonly non-existent, which is probably a result of intrinsic characteristics of species (e.g., feeding habits) triggering different responses to particular events, therefore emphasizing the role of idiosyncratic events in the recent evolutionary history of genetic loci and bird taxa in the region [14, 58].

Furthermore, altitudinal gradients may shape broad biogeographical patterns and potentially population genetic variation but its influence is complex and variable among species [70,71,72,73]. Models show that cold and humid conditions were extended along the Mexican highlands during the Pleistocene; nevertheless the dynamics under which each taxon responded to the climatic cycles depends on several factors (e. g., time of establishment in an area, environmental preferences, dispersal capabilities [74]). Our study also indicates that at least two unrelated bird genera (Aulacorhynchus & Chlorospingus) with distinct natural histories but which inhabit roughly the same regions and wide elevational range along the humid montane forests of Mesoamerica, present highly similar phylogeographic structure, including a significant split between the highlands in Guerrero and Oaxaca within the SMS, which may be expected due to discontinuous distribution ranges, as well as to their presumably reduced dispersal abilities. The pattern observed in Cardellina and Eupherusa do not fully agree with the expectation that the degree of isolation between montane ranges is responsible for shaping the genetic differentiation between populations, as has been reported in previous studies of montane birds (see [45, 75]). The analyses in these two taxa evidenced exclusive haplotypes from Oaxaca, but genetic differentiation across the Río Verde drainage, as seen in Aulacorhynchus and Chlorospingus does not occur. Regarding Cardellina, despite the fact that the TMVB and the SMS montane ranges are separated by a significant geographic barrier such as the Balsas Depression (see [7]), incomplete lineage sorting between both regions is reported. These differences may also be due to the fact that both taxa inhabit different elevational ranges along the region: Eupherusa ranging to the lower limits of the cloud-forest and Cardellina ranging in the upper reaches of the humid montane forest, suggesting that isolation and connectivity cycles in the vegetation types might have a differential impact on the genetic pool of birds depending on the elevational range, which might explain the observed patterns [74].

Phylogenetic and population structure analyses

Our results, as several others from the Neotropics (see [11, 13, 18, 76, 77]), found shared phylogeographic breaks at which genetic differentiation occurs among widespread Mesoamerican highland bird species. We observed major breaks, as seen in Aulacorhynchus and Chlorospingus, between Chiapas and Mexican-Central American lineages. This has been also supported by other bird species [75], plants [18], rodents [20, 35], insects [51] and a snake [78]. Moreover, phylogeography shows isolation between populations inhabiting the mountains of southwestern Mexico (SMS) and northern Oaxaca (see Chlorospingus) evidencing that populations are clearly sorted within Oaxaca across its Central Valley, and more importantly morphological variation of these populations has also been documented as in E. eximia nelsoni [79]. Isolation between Central American highland populations (NCA-SCA) has been likely promoted by the Nicaraguan Depression (see [13, 77, 80]). The most relevant break in our study, occurred within the SMS region, between the states of Oaxaca and Guerrero highlands, which are separated by the Río Verde drainage, this break is supported by studies in plants [33], lizards [81, 82], frogs [31] and even in lowland species, such as iguanas [83].

The patterns observed in the haplotype networks mostly reveal a consistent match between haplogroups and subspecies, as well as species boundaries (see Eupherusa). The overall structures in the haplotype networks resemble an expected pattern for ancient divergence among the sky-islands where reciprocally monophyletic groups occur on each sky-island with no mixing (see Aulacorhynchus and Chlorospingus), and is consistent with estimated divergence and the gradual cooling occurring in the Gelasian stage of late Pliocene-early Pleistocene [84]; whereas in Eupherusa it is difficult to discern specific events for each divergence estimate due to wide confidence intervals. Cardellina is an exception which resembles a star-shaped network of a post glacial sky-island divergence [53] which is consistent with an interglacial period during the Pre-Illinoian stage [84, 85] and the estimated divergence in BEAST and DIYABC analyses.

AMOVA results highlighted the existence of genetic structure for all studied species. Most molecular variation was found among groups; nevertheless, the fixation indices (FSC) among populations within groups were significant, as in Aulacorhynchus, Chlorospingus and Cardellina. This pattern suggests the existence of further structured populations in at least one of the defined clusters, as supported by haplotype networks (Fig. 2b and Fig. 3d).

We obtained overall high values of haplotype diversity (Hd) and low nucleotide diversity (π) for the majority of analyzed groups within each taxon. These patterns have been attributed to population growth following demographic bottlenecks enhancing retention of novel mutations [86, 87], which is congruent in sky-island populations supporting lower nucleotide variation given reduced gene flow to other populations [52, 88, 89]. Population expansions may have been promoted by climatic oscillations during the Pleistocene; nevertheless, we found significant support for population expansion only in Cardellina populations from the TMVB, and SMS. For the rest of the taxa, Tajima’s D and Fu’s FS values suggest that the observed nucleotide polymorphism is selectively neutral even when they tend towards negative figures, thus favoring demographic stability during the Pleistocene. Pairwise comparisons of FST values for each taxon are high and mostly significant, which is expected given their disjunct geographic distribution. Overall, all populations concerned are allopatric and several lineages possess phenotypic and phylogenetic identity, thus the FST fixation indices support the idea that these populations have followed their own evolutionary trajectory for a long time (Tables 3, 4). Despite the FST pairwise comparison between SMS and TMVB populations was low; according to [90] these populations have great genetic differentiation. Moreover, some populations, despite allopatry, may be geographically close enough to show evidence of intermittent connections in response to environmental changes [91].

Divergence time estimation and coalescent-based estimation of population history

Our results highlight the importance of Pliocene-Pleistocene events in promoting the intra-specific genetic structure in the analyzed species, as reported in other studies. Divergence between SCA and northern populations is concordant with the formation of the Nicaragua Depression (ca. 4 Myr [92];), furthermore, divergence between NCA and Mexican lineages post-date the time frame of the formation of the volcanic arcs in Central America (Miocene-Pliocene or before [92];), the IT dry plains (ca. 6–3 Myr [38];), orogeny of SMS (ca. 35–20 Myr), and the formation of the TMVB as well as its period of volcanism in central Mexico (ca. 16–7 Myr [39,40,41, 93, 94];). Our results suggest that these geographic breaks may have acted indirectly as semi-permeable barriers to dispersal, leading to an accumulation of processes involving isolation and expansion [95], as the Cardellina analyses revealed.

DIYABC scenarios allowed us to infer evolutionary processes affecting populations. Although estimated divergence dates varied among populations of the analyzed taxa, major splits occurred during the Pleistocene. Scenarios from Aulacorhynchus, Chlorospingus, and Cardellina, along with their respective pairwise FST comparisons, support the hypothesis of sequential northward population range expansions from Central America with ancestral populations being split by vicariant events, such as fracture and changes in altitudinal ranges of montane forests linked to climatic oscillations, therefore promoting geographical isolation. Although our analyses support a southern origin, the pattern observed in Eupherusa + Thalurania ridgwayi suggest alternative scenarios as also probable. Based on previous evidence supporting T. ridgwayi as a member of the Eupherusa genus [96, 97], we included samples of T. ridgwayi in the diversification scenario comparisons, and we were unable to discern if diversification according to the best-supported scenario occurred from two potential regions that may have served as refugia (one in northwestern Mexico, another in southern Mexico-Central America), or whether SMS (E. poliocerca-E. cyanoprhys) and Central American lineages of Eupherusa come from a migration event from northwestern Mexico. Further studies should focus on the reconstruction of ancestral areas of distribution to clarify this issue (see [50]).

Overall, dates from the DIYABC analyses suggest more recent diversification times than our BEAST analyses, thus highlighting the effect that distinct assumed mutation rates have an effect of temporal frame estimation [98]; nevertheless, these values overlap their respective 95% HDP confidence intervals. Even in vicariant events, differences in time estimates can also be due to variance in the coalescent process and related to the demography of each species [18, 50], however we acknowledge the inconveniences of single locus based analyses and that distinct calibration for divergence time estimation may yield potential errors, thus future surveys incorporating more data should be able to obtain stronger and more accurate estimates of timing and synchrony of divergence.

Test of simultaneous divergence

Even when species ranges overlap across their distributional areas, some divergence time estimation frames are not congruent with each other, therefore, asynchrony in diversification within Mesoamerica could be depicted as a result of Pleistocene climatic oscillations which promoted rapid pulses of diversification and speciation rates across multiple lineages [13, 18, 19]. Our analyses including four population pairs detected non-simultaneous divergence events within the Pleistocene, pointing to the existence of several temporal frames where diversification and fixation of genetic differences within the same montane ranges occurred among distinct lineages. Our results also highlight a high probability of simultaneous divergence occurrence between two taxa (Aulacorhynchus and Chlorospingus) in the SMS region across the Río Verde drainage lowland barrier, with a mean divergence date of 0.9339 Myr, which 95% HDP confidence intervals falls within the estimated divergence dates between the same populations. When Eupherusa population pair (cyanophrys-poliocerca) was added, estimations supported two diversification events, thus implying that even when the three taxa share their overall distributional ranges in the SMS, they do not fully share a community-wide historical pattern across the landscape.

Even when we do not fully out rule any orogenic effect over the analyzed species, major changes to population structure should be promoted by the onset of climatic oscillations in the Pleistocene, which have been invoked to explain patterns of genetic differentiation in the increasing rates of speciation and diversification in Neotropical biota, given the evidence of a downward altitudinal migration of the forest line which resulted in a montane forest-like dominated landscape in the lowlands with further retreat when temperatures raised [99,100,101]. Although for the C. rubra lineage the observed strong signature of expansion in TMVB and SMS populations due to Pleistocene conditions and forest migrations resulted in the admixture of two isolated lineages; our results reveal that demographic dynamics for other species persisted in isolated regions throughout the Pleistocene allowing the evolution of high endemism, which could be an overall result of high effective population sizes.

Genetic endemism as a result of isolation of intraspecific lineages, as seen within the SMS should be considered in conservation measures given that they preserve historical components and maintain the species capability of adaptation [3, 66].

Conclusions

Our survey provides a framework for understanding the complex geographic history of Mesoamerica and its role in promoting lineage diversification in birds. From our analyses we recover multiple independently evolving lineages restricted to montane areas across the region, our results also highlight that congruence may be difficult to occur given idiosyncratic life histories of species; nevertheless, co-distributed taxa will share more common history than a random scenario could predict. Montane regions within Mesoamerica are strong drivers of lineage diversification, and as such our results point into major areas of great genetic diversity and endemism, such as isolated highlands of the SMS in the Mexican states of Oaxaca and Guerrero. Given ancient development of major montane systems, climatic oscillations and orogeny derived conditions were key drivers of diversification between these lineages as splits were calculated to occur mainly within the Pleistocene.

Methods

Sampling and sequence data

We selected four resident and co-distributed avian taxa from the montane forests of Mesoamerica, all of which represent different orders and families (Apodiformes: Trochilidae; Piciformes: Ramphastidae; and Passeriformes: Parulidae and Passerellidae), as well as life histories. Selected taxa include Toucanets (genus Aulacorhynchus), the Red and Red-faced Warbler (genus Cardellina), Bush-tanagers (genus Chlorospingus), and hummingbirds (genus Eupherusa-Thalurania). These taxa include at least one differentiated population restricted to the SMS for which we have previous genetic and morphometric information [5, 28, 68, 69, 96, 97, 102,103,104,105,106,107,108,109]. These studies have also estimated similar Pliocene-Pleistocene divergence times for all of the analyzed taxa.

Phylogenetic and population structure analyses

For each taxon, we retrieved mitochondrial DNA (mtDNA) sequences deposited in GenBank (https://www.ncbi.nlm.nih.gov/genbank/), and conducted analyses using alignments of concatenated mtDNA loci (Additional file 1). We estimated nucleotide substitution models and partition schemes for each species in PARTITIONFINDER [110], using the Bayesian Information Criterion (BIC) for model selection. Given the10 bp frameshift overlap in the ATPase 8 and 6 genes and the tRNA-Lys gene, we analyzed these as a single genetic region [111, 112]. Resultant partition schemes and model substitution parameters were used for conducting a phylogenetic reconstruction using the Bayesian inference approach (BI) for each species, as implemented in MRBAYES 3.2 [113]. We ran two independent searches using four Markov-Chains MonteCarlo for 107 generations sampling every 1000 generations. Convergence across runs was evaluated using two methods: a) the examination of the standard deviation of split frequencies (with acceptance values < 0.01); and b) by verification of parameter estimates in TRACER v1.6 [114], based on acceptable effective sample sizes (ESS values > 200). After checking for convergence, the first 25% of the generated trees were discarded as burn-in and the remaining 75% were kept to calculate posterior probabilities. In addition, we also estimated phylogenetic trees using maximum likelihood (ML) criteria, as implemented in RAXMLGUI 1.5b1 [115,116,117], using the GTRCAT model, and estimated nodal support via 1000 bootstrap iterations using the selected partition. We selected closely-related taxa as outgroups, as suggested by previous published studies: Aulacorhynchus albivitta [108], Chlorospingus flavopectus phaeocephalus [103], Cardellina versicolor [5], and Thalurania ridgwayi [96, 97].

We used NETWORK 4.6.1.1 (Fluxus Engineering, www.fluxus-engineering.com) to visualize the relationships among haplotypes by constructing networks using the median-joining algorithm [118, 119], assigning equal weights to all variable sites and an epsilon parameter with default values (ε = 0).

To explicitly test for phylogeographic structure in each taxon avoiding a priori criteria to delineate populations, we used a Bayesian model-based clustering algorithm, as implemented in GENELAND 4.0.3 [120, 121], which assigns samples to clusters (K) according to both geographical adjacency and genetic similarity through simulations with the Reversible Jump Markov Chain MonteCarlo (RJMCMC) algorithm. We varied the maximum number of expected clusters (Kmax) from 2 to 10, following possible sub structuring of populations within the main montane ranges in Mesoamerica. We performed 10 independent runs of 10,000 iterations each, a thinning value of 1000 and a 10% burn-in. Best results were selected according to highest posterior probabilities.

Genetic structure was assessed through hierarchical analysis of molecular variance (AMOVA) using pairwise differences, based on the number of clusters (K) obtained by GENELAND for each taxon. Following these results, samples of Aulacorhynchus were divided into four population groups (K = 4): 1) Sierra Madre del Sur (SMS) samples from the Mexican states Guerrero and Oaxaca, 2) eastern Mexico and north Central America (EMx-NCA), 3) Costa Rica and northern Panama (CR-NP), and 4) southern Panama (SP); for Chlorospingus, we obtained K = 6: 1) SMS (Guerrero-Oaxaca), 2) northern Chiapas (including the Chimalapas region in Oaxaca), 3) northeastern Oaxaca and Sierra Madre Oriental (SMO), 4) Los Tuxtlas Massif (TM), 5) northern Central America (NCA), and 6) southern Central America (SCA). Samples of Cardellina were divided into K = 4: 1) Trans Mexican Volcanic Belt (TMVB), 2) SMS, 3) Sierra Madre Occidental (SMOc), and 4) northern Central America (NCA-C. versicolor). Finally, the Eupherusa-Thalurania samples, where GENELAND clustering followed current species delimitation: 1) E. cyanophrys (Sierra de Miahuatlán, Oaxaca), 2) E. poliocerca (Guerrero-Oaxaca), 3) E. eximia (N Oaxaca-N Central America), 4) E. nigriventris (Costa Rica-Panama), and 5) Thalurania rigdwayi (Jalisco).

We also assessed genetic diversity of each population group within each species through the estimation of haplotype diversity (Hd), and nucleotide diversity (π). Genetic divergence between groups was measured using pairwise FST fixation index [122, 123], interpreting the results following guides in [90]. Significance of FST tests was assessed using 1000 permutations. To test for evidence of recent demographic changes in the selected species, we estimated historical demographic dynamics through the calculation of Fu’s FS statistic [124] and Tajima’s D statistic [125] neutrality tests. Significance of these tests (p <  0.02 in the case of FS statistic) was calculated through 1000 simulations. All analyses were conducted in ARLEQUIN 3.1 [126].

Divergence time estimation

We used a Bayesian MCMC-based approach to calculate divergence times among haplogroups within each species independently using BEAST v1.8.4 [127]. For each selected taxon, we first tested whether our dataset fits either to a strict clock model or to a relaxed clock model. We performed selection tests through the stepping-stone method (SS [128];), as implemented in MRBAYES 3.2 [113]. Given our partitioning models, the mean marginal likelihood of a strict clock performed better than a relaxed clock model for all taxa (Additional file 2). Thus, chains were run using a UPGMA starting tree, under a strict clock with substitution models (according to results from PARTITIONFINDER) for 108 generations sampling every 1000 steps, using a Yule speciation process [129, 130] with no topological constraints, and discarded the first 25% as burn-in. We used two approximations to convert branch lengths into time: (A) a body mass-corrected molecular clock rate of 0.0042 (min = 0.0011, max = 0.0158) subs/site/Myr [131], and (B) an uncorrelated lognormal relaxed clock substitution rate fixed at 0.01 average subs/site and SD = 0.003 [132,133,134] for each locus calculated in BEAST v1.8.4 [127]. Adjustment of the body mass-corrected and calculated clocks was evaluated using Bayes factors (BF), calculated from the marginal likelihoods from path sampling (PS) and stepping-stone (SS) methods in BEAST v1.8.4 [127, 135, 136]. Each marginal likelihood was estimated through 100 path steps with a Beta distribution (0.3, 1.0). We considered a 3 Log ml difference as strong evidence [137] against the null hypothesis (Additional file 2). We checked for stationarity using TRACER v1.6 [114]. Node ages are presented as mean heights and 95% credibility interval values with a posterior probability limit of 0.5 and resumed in a maximum clade credibility tree (MCCT). All of these were generated using TREEANNOTATOR v1.8.4 [127], and trees were visualized in FIGTREE v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/).

Coalescent-based estimation of population history

We tested for different hypotheses of population divergence and admixture using an approximate Bayesian computation (ABC) approach in DIYABC 2.1.0 [138]. We conducted initial simulations using 15 competing evolutionary scenarios per species (Additional file 3). Evolutionary scenarios for each species were built considering results of the phylogenetic and GENELAND analyses, biogeographic diversification scenarios proposed in previous studies (see [68, 106, 109]), and geographical breaks invoked as responsible of diversification events. We followed recommendations in [50, 139, 140] for the assessment of appropriate priors, as well as to select the highly-informative summary statistics from the simulations resembling datasets similar to the empirical ones. When phylogenies suggested the existence of two genetic groups within the SMS (p. e. distinctive Guerrero and Oaxaca lineages, as in Aulacorhynchus and Chlorospingus), we contrasted the same scenarios considering those clades as different populations in the analyses. All competing scenarios were eliminated in successive rounds, where the preferred hypotheses were the ones with highest posterior probabilities.

From the selection process described above, we tested three final models for Aulacorhynchus, Chlorospingus, and Cardellina; and two for Eupherusa-Thalurania (Additional file 5). For each matrix analyzed, we simulated 1 × 106 datasets and obtained summary statistics per scenario in each simulation. Based on rates reported previously for birds [141,142,143,144,145], we used an HKY mutation model with a uniform prior distribution, a mean mutation rate with a gamma distribution set to 2.0 × 10− 8 substitutions/site/year (min = 1.6 × 10− 8, max = 2.9 × 10− 8 substitutions/site/year) The posterior probabilities of competing scenarios were computed using a logistic regression on the 1% of simulated datasets closest to the observed data. The selected scenario was that with highest probability value and a non-overlapping 95% confidence interval. For the best supported scenario, we performed a model-checking procedure using a Principal Components Analysis (PCA) on test statistics to visualize the fit between simulated and observed datasets. Confidence on the chosen scenario was assessed though the calculation of type I and type II error rate, from 500 simulated pseudo-observed datasets (PODs) generated with the data of the best-supported scenario [146, 147]. Point estimates for demographic and temporal parameters were obtained by local linear regression on the 1% of simulations closest to the observed dataset for the best-supported scenario [138]. Divergence time obtained from DIYABC output was transformed assuming a conservative two-year generation time, which has been previously used in similar species groups [66, 148].

Test of simultaneous divergence

We tested whether genetic differentiation occurred simultaneously between the main geographical barriers in the analyzed clusters by performing a hierarchical approximate Bayesian computation (hABC) analysis as implemented in MSBAYES [4]. Populations inhabiting the montane regions throughout the study area showed similar patterns of genetic structuring even when they are not closely related, we may therefore expect inter-specific simultaneous divergence among regions if the processes driving diversification are common. We performed hABC analyses in 15 population pairs that spanned the same putative barriers to gene flow: between populations within the SMS in southwestern Mexico; between populations of the SMS and eastern populations across a lowland barrier in Oaxaca’s Central Valleys, and between populations spanning the Nicaraguan Depression in Central America (Fig. 1).

We used jModelTest [149, 150] to estimate a transition-transversion rate for each population pair, as required by MSBAYES. The analyses involved the estimation of a vector of summary statistics from the sequences utilized; afterwards 1 × 106 data sets were simulated under the specified multi-taxon model using prior distribution for the demographic parameters. The prior of the maximum possible number of divergence events (Ψ) was set to be equal to the number of lineage pairs tested. For the last stage of the MSBAYES analyses, we used the acceptance/rejection algorithm to approximate the posterior distribution for the hyper parameters that characterize the degree of variability of demographic and temporal parameters given the empirical data (e. g., τ, Ψ, Ω). We used Ω estimates to evaluate the support for each hypothesis: synchronous diversification is expected when Ω ≤ 0.01; asynchronous diversification if Ω > 0.01 [4, 18, 151, 152]. To estimate mean divergence times, we converted model-averaged E(τ) estimates (provided in MSBAYES as coalescent units of 4NAVE generations) to absolute time (Tdiv), using the equation Tdiv = E(τ) × (θAVE/μ) × g, where g is the generation time in years, and θAVE is the mean of the upper θ prior [151], and μ is the mutation rate. We used a mean mutation rate of 0.0042 substitutions/site/Myr prior [131], and an average generation time of 2 years [66, 148].

Availability of data and materials

All of the data generated or analyzed during this study are included in this published article [and its supplementary information files].

Abbreviations

CR-NP:

Costa Rica and northern Panama

EMx-NCA:

Eastern Mexico and north Central America

N Chiapas:

Northern Chiapas

NCA:

Northern Central America

SCA:

Southern Central America

SMO:

Northeastern Oaxaca and Sierra Madre Oriental

SMOc:

Sierra Madre Occidental

SMS:

Sierra Madre del Sur

SP:

Southern Panama

TM:

Tuxtlas Massif

TMVB:

Trans-Mexican Volcanic Belt

References

  1. 1.

    Avise JC. Phylogeography: the history and formation of species. Cambridge: Harvard University Press; 2000.

  2. 2.

    Avise JC. Phylogeography: retrospect and prospect. J Biogeogr. 2008;36:3–15.

  3. 3.

    D’Amen M, Zimmermann NE, Pearman PB. Conservation of phylogeographic lineages under climate change. Glob Ecol Biogeogr. 2012;22:93–104.

  4. 4.

    Hickerson MJ, Stahl E, Takebayashi N. msBayes: Pipeline for testing comparative phylogeographic histories using hierarchical approximate Bayesian computation. BMC Bioinform. 2007;8:268. https://doi.org/10.1186/1471-2105-8-268.

  5. 5.

    Barrera-Guzmán AO, Milá B, Sánchez-González LA, Navarro-Sigüenza AG. Speciation in an avian complex endemic to the mountains of middle America (Ergaticus, Aves: Parulidae). Mol Phylogenet Evol. 2012;62(3):907–20. https://doi.org/10.1016/j.ympev.2011.11.020.

  6. 6.

    Hewitt GM. The genetic legacy of the quaternary ice ages. Nature. 2000;405:907–13.

  7. 7.

    Navarro-Sigüenza AG, Peterson AT, Nyari A, García-Deras GM, García-Moreno J. Phylogeography of the Buarremon brush-finch complex (Aves, Emberizidae) in Mesoamerica. Mol Phylogenet Evol. 2008;47:21–35.

  8. 8.

    Navarro-Sigüenza AG, García-Hernández MA, Peterson AT. A new species of brush-finch (Arremon; Emberizidae) from western Mexico. Wilson J Ornithol. 2013;125:443–684.

  9. 9.

    Cicero C, Koo MS. The role of niche divergence and phenotypic adaptation in promoting lineage diversification in the sage sparrow (Artemisiospiza belli, Aves: Emberizidae). Biol J Linn Soc. 2012;107:332–54.

  10. 10.

    González C, Ornelas JF, Gutiérrez-Rodríguez C. Selection and geographic isolation influence hummingbird speciation: genetic, acoustic and morphological divergence in the wedge-tailed sabrewing (Campylopterus curvipennis). BMC Evol Biol. 2011;11:38. https://doi.org/10.1186/1471-2148-11-38.

  11. 11.

    Milá B, Wayne RK, Fitze P, Smith TB. Divergence with gene flow and fine-scale phylogeographical structure in the wedge-billed woodcreeper, Glyphorynchus spirurus, a Neotropical rainforest bird. Mol Ecol. 2009;18:2979–95. https://doi.org/10.1111/j.1365-294X.2009.04251.x.

  12. 12.

    Pérez-Emán JL, Mumme RL, Jabloński PG. Phylogeography and adaptive plumaje evolution in central american subspecies of the slate-throated redstart (Myioborus miniatus). Ornithol Monogr. 2010;67:90–102.

  13. 13.

    Ortíz-Ramírez MF, Andersen MJ, Zaldívar-Riverón A, Ornelas JF, Navarro-Sigüenza AG. Geographic isolation drives divergence of uncorrelated genetic and song variation in the Ruddy-capped nightingale-thrush (Catharus frantzii; Aves: Turdidae). Mol Phylogenet Evol. 2016;94:74–86. https://doi.org/10.1016/j.ympev.2015.08.017.

  14. 14.

    Arbeláez-Cortés E. Filogeografía comparada: conceptos, métodos y patrones generales en aves Neotropicales. Acta Biol Colomb. 2012;17(1):19–38.

  15. 15.

    Arbeláez-Cortés E, Milá B, Navarro-Sigüenza AG. Multilocus analysis of intraspecific differentiation in three endemic bird species from the northern Neotropical dry forest. Mol Phylogenet Evol. 2014;70:236–377. https://doi.org/10.1016//j.mpev.2013.10.006.

  16. 16.

    Linder HP. On areas of endemism, with an example from the African Restionaceae. Syst Biol. 2001;50:892–912. https://doi.org/10.1080/106351501753462867.

  17. 17.

    Sánchez-González LA, Navarro-Sigüenza AG. History meets ecology: a geographical analysis of ecological restriction in the Neotropical humid montane forest avifaunas. Divers Distrib. 2009;15:1–11.

  18. 18.

    Ornelas JF, Sosa V, Soltis DE, Daza JM, González C, Soltis PS, Gutiérrez-Rodríguez C, de los Monteros AE, Castoe TA, Bell C, Ruiz-Sánchez E. Comparative phylogeographic analyses illustrate the complex evolutionary history of threatened cloud forests of northern Mesoamerica. PLoS One. 2013;8:e56283.

  19. 19.

    Barber BR, Klicka J. Two pulses of diversification across the isthmus of Tehuantepec in a montane Mexican bird fauna. Proc Biol Sci. 2010;277(1694):2675–81. https://doi.org/10.1098/rspb.2010.0343.

  20. 20.

    Harris D, Rogers DS, Sullivan J. Phylogeography of Peromyscus furvus (Rodentia; Muridae) based on cytochrome b sequence data. Mol Ecol. 2000;9:2129–35. https://doi.org/10.1046/j.1365-294x.2000.01135.x.

  21. 21.

    Muñoz-Villers LE, López-Blanco J. Land use/cover changes using Landsat TM/ETM images in a tropical and biodiverse montainous area of Central-Eastern Mexico. Int J Remote Sens. 2007;29:71–93.

  22. 22.

    Hernández-Baños BE, Peterson AT, Navarro-Sigüenza AG, Escalante-Pliego P. Bird faunas of the humid montane forests of Mesoamerica. Biogeographic patterns and priorities for conservation. Bird Conserv Int. 1995;5:251–77.

  23. 23.

    Luna-Vega I, Espinosa D, Contreras-Medina R. Biodiversidad de la Sierra Madre del Sur. Ciudad de México: Universidad Nacional Autónoma de México; 2016.

  24. 24.

    Morales R, Escalante T, Noguera-Urbano EA, Gámez N, Rodríguez-Tapia G. Conservation biogeography in the Mexican Mountain Component: bridging conservation and patterns of endemism. Therya. 2016;7(2):215–29. https://doi.org/10.12933/therya-16-329.

  25. 25.

    Antonelli A, Kissling WD, Flantua SG, Bermúdez MA, Mulch A, Muellner-Riehl AN, Kreft H, Linder HP, Badgley C, Fjeldså J, Fritz SA, Rahbek C, Herman F, Hooghiemstra H, Hoorn C. Geological and climatic influences on mountain biodiversity. Nat Geosci. 2018;11:718–25. https://doi.org/10.1038/s41561-018-0236-z.

  26. 26.

    Bertelli S, Szumik C, Goloboff PA, Giannini NP, Navarro-Sigüenza AG, Peterson AT, Cracraft J. Mexican land birds reveal complexity in fine-scale patterns of endemism. J Biogeogr. 2017;44(8):1836–46. https://doi.org/10.1111/jbi.12987.

  27. 27.

    Carnaval AC, Hickerson MJ, Haddad CFB, Rodrigues MT, Moritz C. Stability predicts genetic diversity in the Brazilian Atlantic forest hotspot. Science. 2009;323(5915):785–9. https://doi.org/10.1126/science.1166955.

  28. 28.

    García-Moreno J, Navarro-Sigüenza AG, Peterson AT, Sánchez-González LA. Genetic variation coincides with geographic structure in the common bush-tanager (Chlorospingus ophthalmicus) complex from Mexico. Mol Phylogenet Evol. 2004;33:186–96. https://doi.org/10.1016/j.ympev.2004.05.007.

  29. 29.

    Ornelas JF, Ruiz-Sánchez E, Sosa V. Phylogeography of Podocarpus matudae (Podocarpaceae): pre-quaternary relicts in northern Mesoamerican cloud forests. J Biogeogr. 2010;37:2384–96. https://doi.org/10.1111/j.1365-2699.2010.02372.x.

  30. 30.

    Halffter G, Llorente-Bousquets J, Morrone JJ. La perspectiva biogeográfica histórica. In: Capital natural de México vol. I: Conocimiento actual de la biodiversidad. México: CONABIO; 2008. p. 67–86.

  31. 31.

    Campbell JA, Blancas-Hernández JC, Smith EN. A new species of stream-breeding treefrog of the genus Charadrahyla (Hylidae) from the Sierra Madre del Sur of Guerrero, Mexico. Copeia. 2009;2:287–95.

  32. 32.

    Feria-Ortiz M, García-Vázquez UO. A new species of Plestiodon (Squamata: Scincidae) from Sierra Madre del Sur of Guerrero, México. Zootaxa. 2012;3339:57–68.

  33. 33.

    González-Villarreal LM. Two new species of oak (Fagaceae, Quercus sect. Lobatae) from the Sierra Madre del Sur, Mexico. Brittonia. 2003;55(1):49–60. https://doi.org/10.1663/0007-196X(2003)055[0049:TNSOOF]2.0.CO;2.

  34. 34.

    Navarro-Sigüenza AG, Blancas-Calva E, Almazán-Núñez RC, Hernández-Baños BE, García-Trejo EA, Peterson AT. Diversidad y endemismo de las aves de la Sierra Madre del Sur. In: Luna-Vega I, Espinosa D, Contreras-Medina R, editors. Biodiversidad de la Sierra Madre del Sur. Ciudad de México: Universidad Nacional Autónoma de México; 2016. p. 381–411.

  35. 35.

    Ordóñez-Garza N, Matson JO, Strauss RE, Bradley RD, Salazar-Bravo J. Patterns of phenotypic and genetic variation in three species of endemic Mesoamerican Peromyscus (Rodentia: Cricetidae). J Mammal. 2010;91(4):848–59. https://doi.org/10.1644/09-MAMM-A-167.1.

  36. 36.

    Morrone JJ. Halffter’s Mexican transition zone (1962-2014), cenocrons and evolutionary biogeography. J Zool Syst Evol Res. 2015;53:249–57. https://doi.org/10.1111/jzs.12098.

  37. 37.

    Inegi-Conabio-INE (Instituto Nacional de estadística y Geografía- Comisión Nacional para el Conocimiento y Uso de la Biodiversidad- Instituto Nacional de Ecología): Ecorregiones terrestres de México. Escala 1:1,000,000. Mexico City, 2008.

  38. 38.

    Barrier E, Velasquillo L, Chavez M, Gaulon R. Neotectonic evolution of the isthmus of Tehuantepec (southeastern Mexico). Tetonophysics. 1998;287(1–4):77–96. https://doi.org/10.1016/s0040-1951(98)80062-0.

  39. 39.

    Ferrari L, Orozco-Esquivel T, Manea V, Manea M. The dynamic history of the trans-Mexican Volcanic Belt and the Mexico subduction zone. Tectonophysics. 2012;522:122–49.

  40. 40.

    Morán-Zenteno DJ, Tolson G, Martínez-Serrano RG, Martiny B, Schaaf P, Silva-Romo G, Macías-Romo C, Alba-Aldave L, Hernández-Bernal MS, Solis-Pichardo GN. Tertiary arc-magnetism of the Sierra Madre del Sur, Mexico, and its transition of the volcanic activity of the trans-Mexican Volcanic Belt. J S Am Earth Sci. 1999;12(6):513–35.

  41. 41.

    Morán-Zenteno DJ, Martiny B, Tolson G, Solís-Pichardo G, Alba-Aldave L, del Sol H-BM, Macías-Romo C, Martínez-Serrano RG, Schaaf P, Silva-Romo G. Geocronología y características geoquímicas de las rocas magmáticas de la Sierra Madre del Sur. B Soc Geol Mex. 2000;53:27–58.

  42. 42.

    Morrone JJ. Biogeographic regionalization of the Sierra Madre del Sur province, Mexico. Rev Mex Biodivers. 2017;88:710–4. https://doi.org/10.1016/j.rmb.2017.07.012.

  43. 43.

    González-García F, Gómez de Silva H. Especies endémicas: Riqueza, patrones de distribución y retos para su conservación. In: Gómez de Silva H, Oliveras de Ita A, editors. Conservación de aves: Experiencias en México: CIPAMEX; 2003. p. 150–94.

  44. 44.

    Arbeláez-Cortés E, Roldán-Piña D, Navarro-Sigüenza AG. Multilocus phylogeography and morphology give insights into the recent evolution of a Mexican endemic songbird: Vireo hypochryseus. J Avian Biol. 2014;45(3):253–63. https://doi.org/10.1111/j.1600-048X.2013.00335.x.

  45. 45.

    Cortés-Rodríguez N, Hernández-Baños BE, Navarro-Sigüenza AG, Peterson AT, García-Moreno J. Phylogeography and population genetics of the amethyst-throated hummingbird (Lampornis amethystinus). Mol Phylogenet Evol. 2008;48:1–11. https://doi.org/10.1016/j.ympev.2008.02.005.

  46. 46.

    Klicka J, Spellman GM, Winker K, Chua V, Smith BT. A phylogeographic and population genetic analysis of a widespread, sedentary north American bird: the hairy woodpecker (Picoides villosus). Auk. 2011;128(2):346–62. https://doi.org/10.1525/auk.2011.10264.

  47. 47.

    Mason NA, Olvera-Vital A, Lovette IJ, Navarro-Sigüenza AG. Hidden endemism, deep polyphyly, and repeated dispersal across the isthmus of Tehuantepec: diversification of the white-collared seedeater complex (Thraupidae: Sporophila torqueola). Ecol Evol. 2018;8(3):1867–81. https://doi.org/10.1002/ece3.3799.

  48. 48.

    Núñez-Zapata J, Peterson AT, Navarro-Sigüenza AG. Pleistocene diversification and speciation of white-throated thrush (Turdus assimilis; Aves: Turdidae). J Ornithol. 2016;157(4):1073–85. https://doi.org/10.1007/s10336-016-1350-6.

  49. 49.

    Navarro-Sigüenza AG, Peterson AT. An alternative species taxonomy of the birds of Mexico. Biota Neotrop. 2004;4:1–32.

  50. 50.

    Chan LM, Brown JL, Yoder AD. Integrating statistical genetic and geospatial methods bring new power to phylogeography. Mol Phylogenet Evol. 2011;59:523–37.

  51. 51.

    Duennes MA, Lozier JD, Hines HM, Cameron SA. Geographical patterns of genetic divergence in the widespread Mesoamerican bumblebee Bombus ephippiatus (Hymenoptera: Apidae). Mol Phylogenet Evol. 2012;64(1):219–31. https://doi.org/10.1016/j.ympev.2012.03.018.

  52. 52.

    Masta S. Phylogeography of the jumping spider Habronattus pugillis (Araneae: Salticidae): recent vicariance of sky island populations? Evolution. 2000;54:1699–711.

  53. 53.

    McCormack JE, Bowen BS, Smith TB. Integrating paleoecology and genetics of bird populations in two sky island archipelagos. BMC Biol. 2008;6:28. https://doi.org/10.1186/1441-7007-6-28.

  54. 54.

    Tennessen JA, Zamudio KR. Genetic differentiation among mountain island populations of the striped plateau lizard, Sceloporus virgatus (Squamata: Phrynosomatidae). Copeia. 2008;2008(3):558–64. https://doi.org/10.1643/CG-06-038.

  55. 55.

    Endo Y, Nash M, Hoffmann AA, Slatyer R, Miller AD. Comparative phylogeography of alpine invertebrates indicates deep lineage diversification and historical refugia in the Australian Alps. J Biogeogr. 2015;42:89–102. https://doi.org/10.1111/jbi.12387.

  56. 56.

    Ornelas JF, González C, Hernández-Baños BE, García-Moreno J. Molecular and iridescent feather reflectance data reveal recent genetic diversification and phenotypic differentiation in a cloud forest hummingbird. Ecol Evol. 2016;6(4):1104–27. https://doi.org/10.1002/ece3.1950.

  57. 57.

    Oswald JA, Overcast I, Mauck WM, Andersen MJ, Smith BT. Isolation with asymmetric gene flow during the nonsynchronous divergence of dry forest birds. Mol Ecol. 2017;26(5):1386–400. https://doi.org/10.1111/mec.14013.

  58. 58.

    Qu Y, Lei F. Comparative phylogeography of two endemic birds of the Tibetan plateau, the white-rumped snow finch (Onychostruthus taczanowskii) and the Hume’s ground tit (Pseudopodoces humilis). Mol Phylogenet Evol. 2009;51:312–26. https://doi.org/10.1016/j.ympev.2009.01.013.

  59. 59.

    Qu Y, Lei F, Zhang R, Lu X. Comparative phylogeography of five avian species: implications for Pleistocene evolutionary history in the Qinghai-Tibetan plateau. Mol Ecol. 2010;19:338–51. https://doi.org/10.1111/j.1365-294X.2009.04445.x.

  60. 60.

    Soltis DE, Morris AB, McLachlan JS, Manos PS, Soltis PS. Comparative phylogeography of unglaciated eastern North America. Mol Ecol. 2006;15:4261–93. https://doi.org/10.1111/j.1365-294X.2006.03061.x.

  61. 61.

    Frankham R. Challenges and opportunities of genetic approaches to biological conservation. Biol Conserv. 2010;143:1919–27. https://doi.org/10.1016/j.biocon.2010.05.011.

  62. 62.

    Frankham R, Bradshaw CJA, Brook BW. Genetics in conservation management: revised recommendations for the 50/50 rules, red list criteria and population viability analyses. Biol Conserv. 2014;170:56–63. https://doi.org/10.1016/j.biocon.2013.12.036.

  63. 63.

    Gamble T, Colli GR, Rodrigues MT, Werneck FP, Simons AM. Phylogeny and cryptic diversity in geckos (Phyllopezus; Phyllodactylidae; Gekkota) from South America’s open biomes. Mol Phylogenet Evol. 2011;62:943–53. https://doi.org/10.1016/j.ympev.2011.11.033.

  64. 64.

    Moritz C, Faith DP. Comparative phylogeography and the identification of genetically divergent areas for conservation. Mol Ecol. 1998;7:419–29.

  65. 65.

    Pearse DE, Crandall KA. Beyond FST: analysis of population genetic data for conservation. Conserv Genet. 2004;5:585–602.

  66. 66.

    Qu Y, Song G, Gao B, Quan Q, Ericson PGP, Lei F. The influence of geological events on the endemism of east Asian birds studied through comparative phylogeography. J Biogeogr. 2015;42:179–92.

  67. 67.

    Yannic G, St-Laurent MH, Ortego J, Taillon J, Beauchemin A, Bernatchez L, Dussault C, Côté SD. Integrating ecological and genetic structure to define management units for caribou in eastern Canada. Conserv Genet. 2016;17:437–53. https://doi.org/10.1007/s10592-015-0795-0.

  68. 68.

    Hernández-Baños BE, Zamudio-Beltrán LE, García-Moreno J, Eguiarte LE. Phylogenetic relationships within the genus Eupherusa inferred from mtDNA sequences. Ornitol Neotrop. 2012;23:73–9.

  69. 69.

    Sánchez-González LA, Navarro-Sigüenza AG, Peterson AT, García-Moreno J. Taxonomy of Chlorospingus ophthalmicus in Mexico and northern Central America. Bull B O C. 2007;127(1):34–49.

  70. 70.

    Cicero C, Koo MS. The role of niche divergence and phenotypic adaptation in promoting lineage diversification in the sage sparrow (Artemisiospiza belli, Aves: Emberizidae). Biol J Linn Soc. 2012;107:332–54. https://doi.org/10.1111/j.1095-8312.2012.01942.x.

  71. 71.

    Sánchez-González LA, Navarro-Sigüenza AG, Krabbe NK, Fjeldså J, García-Moreno J. Diversification in the Andes: the Atlapetes brush-finches. Zool Scr. 2015;44(2):135–52. https://doi.org/10.1111/zsc.12088.

  72. 72.

    Múrria C, Bonada N, Arnedo MA, Prat N, Vogler AP. Higher β-and γ-diversity at species and genetic levels in headwaters than in mid-order streams in Hydropsyche (Trichoptera). Freshw Biol. 2013;58:2226–36.

  73. 73.

    Hahn T, Kettle CJ, Ghazoul J, Frei ER, Matter P, Pluess AR. Patterns of genetic variation across altitude in three plant species of semi-dry grasslands. PLoS One. 2012;7(8):e41608.

  74. 74.

    Mastretta-Yanes A, Moreno-Letelier A, Piñero D, Jorgensen TH, Emerson BC. Biodiversity in the Mexican highlands and the interaction of geology, geography and climate within the trans-Mexican Volcanic Belt. J Biogeogr. 2015;42:1586–600. https://doi.org/10.1111/jbi.12546.

  75. 75.

    Zamudio-Beltrán LE, Hernández-Baños BE. Genetic and morphometric divergence in the garnet-throated hummingbird Lamprolaima rhami (Aves: Trochilidae). Peer J. 2018;6:e5733. https://doi.org/10.7717/peerj.5733.

  76. 76.

    Cabanne G, Dhorta F, Sari E, Santos F, Miyaki C. Nuclear and mitochondrial phylogeography of the Atlantic forest endemic Xiphorhynchus fuscus (Aves: Dendrocolaptidae): biogeography and systematics implications. Mol Phylogenet Evol. 2008;49(3):760–73. https://doi.org/10.1016/j.ympev.2008.09.013.

  77. 77.

    Sánchez-González LA, Morrone JJ, Navarro-Sigüenza AG. Distributional patterns of the Neotropical montane forest avifaunas. Biol J Linn Soc. 2008;94(1):175–94. https://doi.org/10.1111/j.1095-8312.2008.00979.x.

  78. 78.

    Devitt TJ. Phylogeography of the Western Lyresnake (Trimorphodon biscutatus): testing aridland biogeographical hypotheses across the Neartic-Neotropical transition. Mol Ecol. 2006;15(14):4387–407. https://doi.org/10.1111/j.1365-294X.2006.03015.x.

  79. 79.

    Hernández-Baños BE. Filogenia del género Eupherusa Gould (Aves: Trochilidae). In: Doctoral dissertation. México: Facultad de Ciencias, Universidad Nacional Autónoma de México; 1998.

  80. 80.

    Vázquez-Miranda H, Navarro-Sigüenza AG, Omland KE. Phylogeography of the Rufous-Naped wren (Campylorhynchus rufinucha): speciation and hybridization in Mesoamerica. Auk. 2009;126(4):765–78. https://doi.org/10.1525/auk.2009.07048.

  81. 81.

    Flores-Villela O, Sánchez-H O. A new species of Abronia (Squamata:Anguidae) from the Sierra Madre del Sur of Guerrero, Mexico, with comments on Abronia deppii. Herpetologica. 2003;59(4):524–31. https://doi.org/10.1655/02-39.

  82. 82.

    Nieto-Montes de Oca A, Campbell JA, Flores-Villela O. A new species of Xenosaurus (Squamata: Xenosauridae) from the Sierra Madre del Sur of Oaxaca, Mexico. Herpetologica. 2001;57(1):32–47.

  83. 83.

    Zarza E, Reynoso VH, Emerson BC. Diversification in the northern neotropics: mitochondrial and nuclear DNA phylogeography of the iguana Ctenosaura pectinata and related species. Mol Ecol. 2008;17:3259–75. https://doi.org/10.1111/j.1365-294X.2008.03826.x.

  84. 84.

    Gibbard P, van Kolfschoten T. The Pleistocene and Holocene Epochs. In: Gradstein FH, Ogg JG, Smith AG, editors. A Geologic Time Scale, vol. 2005: Cambridge University Press; 2004. p. 441–52.

  85. 85.

    Lüthi D, Le Floch M, Bereiter B, Blunier T, Barnola JM, Siegenthaler U, Raynaud D, Jouzel J, Fischer H, Kawamura K, Stocker TF. High-resolution carbon dioxide concentration record 650,000-800,000 years before present. Nature. 2008;453(7193):379–82. https://doi.org/10.1038/nature06949.

  86. 86.

    Grant WS, Bowen BW. Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation. J Hered. 1998;89:415–26.

  87. 87.

    Watterson GA. Allele frequencies after a bottleneck. Theor Popul Biol. 1984;26:387–407.

  88. 88.

    Barber PH. Phylogeography of the canyon treefrog, Hyla arenicolor (cope) based on mitochondrial DNA sequence data. Mol Ecol. 1999;8(4):547–62. https://doi.org/10.1046/j.1365-294X.1999.00593.x.

  89. 89.

    Frankham R. Do island populations have less genetic variation than mainland populations? Heredity. 1997;78(3):311–27.

  90. 90.

    Hartl DL, Clark AG. Principles of population genetics (Vol. 116). Sunderland: Sinauer Associates; 1997.

  91. 91.

    Barker BS, Rodríguez-Robles JA, Aran VS, Montoya A, Waide RB, Cook JA. Sea level, topography, and island diversity: Phylogeography of the Puerto Rican red-eyed Coquí, Eleutherodactylus antillensis. Mol Ecol. 2012;21(24):6033–52. https://doi.org/10.1111/mec.12020.

  92. 92.

    Gutiérrez-García TA, Vázquez-Domínguez E. Consensus between genes and stones in the biogeographic and evolutionary history of Central America. Quat Res. 2013;79(3):311–24. https://doi.org/10.1016/j.yqres.2012.12.007.

  93. 93.

    Bryson RW, García-Vázquez UO, Riddle BR. Relative roles of Neogene vicariance and quaternary climate change on the historical diversification of bunchgrass lizards (Sceloporus scalaris group) in Mexico. Mol Phylogenet Evol. 2012;62:447–57. https://doi.org/10.1016/j.ympev.2011.10.014.

  94. 94.

    Verma SP, Carrasco-Núñez G. Reappraisal of the geology and geochemistry of Volcán Zamorano, Central Mexico: implications for discriminating the Sierra Madre occidental and the Mexican Volcanic Belt. Int Geol Rev. 2003;45:724–52.

  95. 95.

    Smith BT, McCormack JE, Cuervo AM, Hickerson MJ, Aleixo A, Cadena CD, Pérez-Emán J, Burney CW, Xie X, Harvey MG, Faircloth BC, Glenn TC, Derryberry EP, Prejean J, Fields S, Brumfield RT. The drivers of tropical speciation. Nature. 2014;515(7527):406–9. https://doi.org/10.1038/nature13687.

  96. 96.

    McGuire JA, Witt CC, Altshuler DL, Remsen JV. Phylogenetic systematics and biogeography of hummingbirds: Bayesian and maximum likelihood analyses of partitioned data and selection of an appropriate partitioning strategy. Syst Biol. 2007;56(5):837–56. https://doi.org/10.1080/10635150701656360.

  97. 97.

    McGuire JA, Witt CC, Remsen JV, Corl A, Rabosky DL, Altshuler DL, Dudley R. Molecular phylogenetics and the diversification of hummingbirds. Curr Biol. 2014;24(8):910–6. https://doi.org/10.1016/j.cub.2014.03.016.

  98. 98.

    Arbogast BS, Edwards SV, Wakeley J, Beerli P, Slowinski JB. Estimating divergence times from molecular data on phylogenetic and population genetic timescales. Annu RevEcol Evol Syst. 2002;33:707–40. https://doi.org/10.1146/annurev.ecolsys.33.010802.150500.

  99. 99.

    Caballero M, Lozano-García S, Vázquez-Selem L, Ortega B. Evidencias de cambio climático y ambiental en registros glaciales y en cuencas del centro de México durante el último máximo glacial. B Soc Geol Mex. 2010;62:359–77.

  100. 100.

    Lachniet MS, Vázquez-Selem L. Last glacial maximum equilibrium line altitudes in the circum-Caribbean (Mexico, Guatemala, Costa Rica, Colombia, and Venezuela). Quat Int. 2005;138:129–44. https://doi.org/10.1016/j.quaint.2005.02.010.

  101. 101.

    Ramírez-Barahona S, Eguiarte LE. The role of glacial cycles in promoting genetic diversity in the Neotropics: the case of cloud forests during the last glacial maximum. Ecol Evol. 2013;3(3):725–38. https://doi.org/10.1002/ece3.483.

  102. 102.

    Avendaño JE, Stiles FG, Cadena CD. A new subspecies of common bush-tanager (Chlorospingus flavopectus, Emberizidae) from the east slope of the Andes of Colombia. Ornitol Colomb. 2013;13:44–58.

  103. 103.

    Bonaccorso E, Navarro-Sigüenza AG, Sánchez-González LA, Peterson AT, García-Moreno J. Genetic differentiation of the Chlorospingus ophthalmicus complex in Mexico and Central America. J Avian Biol. 2008;39:311–21. https://doi.org/10.1111/j.2008.0908-8857.04233.x.

  104. 104.

    Bonaccorso E, Guayasamin JM, Peterson AT, Navarro-Sigüenza AG. Molecular phylogeny and systematics of Neotropical toucanets in the genus Aulacorhynchus. Zool Scr. 2011;40:336–49.

  105. 105.

    Lovette IJ, Pérez-Emán JL, Sullivan JP, Banks RC, Fiorentino I, Córdoba-Córdoba S, Echeverry-Galvis M, Barker FK, Burns KJ, Klicka J, Lanyon SM, Bermingham E. A comprehensive multilocus phylogeny for the wood-warblers and a revised classification of the Parulidae (Aves). Mol Phylogenet Evol. 2010;57:753–70. https://doi.org/10.1016/j.ympev.2010.07.018.

  106. 106.

    Maldonado-Sánchez D, Gutiérrez-Rodríguez C, Ornelas JF. Genetic divergence on the common bush-tanager Chlorospingus ophthalmicus (Aves: Emberizidae) throughout Mexican cloud forests: the role of geography, ecology and Pleistocene climatic fluctuations. Mol Phylogenet Evol. 2016;99:76–88. https://doi.org/10.1016/j.ympev.2016.03.014.

  107. 107.

    Navarro-Sigüenza AG, Peterson AT, López-Medrano E, Benítez-Díaz H. Species limits in Mesoamerican Aulacorhynchus toucanets. Wilson J Ornithol. 2001;113(4):363–72. https://doi.org/10.1676/0043-5643(2001)113[0363:SLIMAT]2.0.CO;2.

  108. 108.

    Puebla-Olivares F, Bonaccorso E, de los Monteros AE, Omland KE, Llorente-Busquets JE, Peterson, AT, Navarro-Sigüenza AG: Speciation in the emerald toucanet (Aulacorhynchus prasinus) complex. Auk 2008, 125: 39–50.

  109. 109.

    Schindel DE, Stoeckle MY, Milensky CM, Trizna M, Schmidt BK, Gebhard CA, Graves GR. Project description: DNA barcodes of bird species in the National Museum of Natural History, Smithsonian Institution, USA. ZooKeys. 2011;152:87–91. https://doi.org/10.3897/zookeys.152.2473.

  110. 110.

    Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29(6):1695–701.

  111. 111.

    Hunt JS, Bermingham E, Ricklefs RE. The molecular systematics and biogeography of Antillean thrashers, tremblers and mockingbirds (Aves: Mimidae). Auk. 2001;118(1):35–55.

  112. 112.

    Lovette IJ. Molecular phylogeny and plumage signal evolution in a trans Andean and circum Amazonian avian species complex. Mol Phylogenet Evol. 2005;32:512–23. https://doi.org/10.1016/j.ympev.2004.01.007.

  113. 113.

    Ronquist F, Huelsenbeck J, Teslenko M. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2011;61(3):539–42. https://doi.org/10.1093/sysbio/sys029.

  114. 114.

    Rambaut A, Suchard MA, Xie D, Drummond AJ: Tracer v1.6 (2014, available from http://beast.bio.ed.ac.uk/Tracer).

  115. 115.

    Silvestro D, Michalak I. raxmlGUI: a graphical front-end for RAxML. Org Divers Evol. 2012;12(4):335–7. https://doi.org/10.1007/s13127-011-0056-0.

  116. 116.

    Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90. https://doi.org/10.1093/bioinformatics/btl446.

  117. 117.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3. https://doi.org/10.1093/bioinformatics/btu033.

  118. 118.

    Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16(1):37–48.

  119. 119.

    Crandall KA, Templeton AR. Applications of intraspecific phylogenetics. In: Harvey PH, Leigh Brown AJ, Maynard Smith J, Nee S, editors. New uses for new phylogenies. Oxford: Oxford University Press; 1996. p. 81–99.

  120. 120.

    Guillot G, Estoup A, Mortier F, Cosson JF. A spatial statistical model for landscape genetics. Genetics. 2005;170:1261–80.

  121. 121.

    Guillot G, Mortier F, Estoup A. Geneland: a program for landscape genetics. Mol Ecol Notes. 2005;5:712–5.

  122. 122.

    Wright S. Evolution in Mendelian populations. Genetics. 1931;16(2):97–159.

  123. 123.

    Wright S. Evolution and genetics of populations. Vol. 4.Variability Within and Among Natural Populations. Chicago: Univ. of Chicago Press; 1978.

  124. 124.

    Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147(2):915–25.

  125. 125.

    Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123(3):585–95.

  126. 126.

    Excoffier L, Laval G, Schneider S. Arlequin ver 3.0: an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.

  127. 127.

    Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUTi and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.

  128. 128.

    Xie WG, Lewis PO, Fan Y, Kuo L, Chen MH. Improving marginal likelihood estimation for bayesian phylogenetic model selection. Syst Biol. 2011;60:150–60.

  129. 129.

    Gernhard T. The conditioned reconstructed process. J Theor Biol. 2008;253:769–78. https://doi.org/10.1016/j.jtbi.2008.04.005.

  130. 130.

    Yule GU. A Mathematical theory of evolution based on the conclusions of Dr. J C Willis. Philos Trans R Soc B. 1925;213:21–87. https://doi.org/10.1098/rstb.1925.0002.

  131. 131.

    Nabholz B, Lanfear R, Fuchs J. Body mass-corrected molecular rate for bird mitochondrial DNA. Mol Ecol. 2016;25:4438–49. https://doi.org/10.1111/mec.13780.

  132. 132.

    Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4(5):e88. https://doi.org/10.1371/journal.pbio.0040088.

  133. 133.

    Pons J, Ribera I, Bertranpetit J, Balke M. Nucleotide substitution rates for the full set of mitochondrial protein-coding genes in Coleoptera. Mol Phylogenet Evol. 2010;56:796–807. https://doi.org/10.1016/j.ympev.2010.02.007.

  134. 134.

    Tong KJ, Duchêne DA, Duchêne S, Geoghegan JL, Ho SYW. A comparison of methods for estimating substitution rates from ancient sequence data. BMC Evol Biol. 2018;18:70. https://doi.org/10.1186/s12862-018-1192-3.

  135. 135.

    Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko AV. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol. 2012;29(9):2157–67. https://doi.org/10.1093/molbev/mss084.

  136. 136.

    Baele G, Li WLS, Drummond AJ, Suchard MA, Lemey P. Accurate model selection of relaxed molecular clocks in Bayesian phylogenetics. Mol Biol Evol. 2013;30(2):239–43. https://doi.org/10.1093/molbev/mss243.

  137. 137.

    Raftery AE. Approximate Bayes factors and accounting for model uncertainty in generalized linear models. Biometrika. 1996;83:251–66.

  138. 138.

    Cornuet JM, Pudlo P, Veyssier J, Dehne-Garcia A, Gautier M, Leblois R, Martin JM, Estoup A. DIY ABC v2.0: a software to make approximate Bayesian computation inferences about population history using single nucleotide polymorphism, DNA sequence and microsatellite data. Bioinformatics. 2014;30(8):1187–9. https://doi.org/10.1093/bioinformatics/btt763.

  139. 139.

    Bertorelle G, Benazzo A, Mona S. ABC as a flexible framework to estimate demography over space and time: some cons, many pros. Mol Ecol. 2010;19:2609–25.

  140. 140.

    Knowles LL, Alvarado-Serrano DF. Exploring the population genetic consequences of the colonization process with spatio-temporally explicit models: insights from coupled ecological, demographic and genetic models in montane grasshoppers. Mol Ecol. 2010;19(17):3727–45. https://doi.org/10.1111/j.1365-294X.2010.04702.x.

  141. 141.

    Lerner HRL, Meyer M, James HF, Hofreiter M, Fleischer RC. Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers. Curr Biol. 2011;21:1838–44.

  142. 142.

    Lim HC, Sheldon FH. Multilocus analysis of the evolutionary dynamics of rainforest bird populations in Southeast Asia. Mol Ecol. 2011;20:3414–38.

  143. 143.

    Manthey JD, Klicka J, Spellman GM. Isolation-driven divergence: speciation in a widespread north American songbird (Aves: Certhiidae). Mol Ecol. 2011;20:4371–84.

  144. 144.

    Smith BT, Klicka J. The profound influence of the late Pliocene Panamanian uplift on the exchange, diversification, and distribution of New World birds. Ecography. 2010;33:333–42.

  145. 145.

    Weir JT, Schluter D. Calibrating the avian molecular clock. Mol Ecol. 2008;17:2321–8.

  146. 146.

    Cornuet JM, Ravigne V, Estoup A: Inference on population history and model checking using DNA sequence and microsatellite data with the software DIYABC (v1.0). BMC Bioinform 2010, 11: 401. doi: https://doi.org/10.1186/1471-2105-11-401.

  147. 147.

    Fontaine MC, Austerlitz F, Giraud T, Labbé F, Papura D, Richard-Cervera S, Delmotte F. Genetic signature of a range expansion and leap-frog event after the recent invasion of Europe by the grapevine downy mildew pathogen Plasmopara viticola. Mol Ecol. 2013;22(10):2771–86. https://doi.org/10.1111/mec.12293.

  148. 148.

    Licona-Vera Y, Ornelas JF, Wethington S, Bryan KB. Pleistocene range expansions promote divergence with gene flow between migratory and sedentary populations of Calothorax hummingbirds. Biol J Linn Soc. 2018;124:645–67.

  149. 149.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.

  150. 150.

    Guindon S, Gascuel O. A simple, fast and accurate method to estimate large phylogenies by maximum-likelihood. Syst Biol. 2003;52:696–704.

  151. 151.

    Bagley JC, Hickerson M, Johnson JB. Testing hypotheses of diversification in Panamanian frogs and freshwater fishes using hierarchical approximate Bayesian computation with model averaging. Diversity. 2018;10(4):120. https://doi.org/10.3390/d10040120.

  152. 152.

    Stone GN, Lohse K, Nicholls JA, Fuentes-Utrilla P, Sinclair F, Schoenrogge K, Csoka G, Melika G, Nieves-Aldrey JL, Pujade-Villar J, Tavakoli M, Askew RR, Hickerson MJ. Reconstructing community assembly in time and space reveals enemy scape in a western Paleartic insect community. Curr Biol. 2012;22(6):532–7. https://doi.org/10.1016/j.cub.2012.01.059.

Download references

Acknowledgments

This paper is a partial fulfillment for the obtention of a M. Sc. Degree of AR-M at the Posgrado en Ciencias Biológicas, UNAM. We thank all the scientific collectors that directly or indirectly contributed with the obtaining of samples used in this study. We also thank Alicia Mastretta-Yanes and Leopoldo Vázquez-Reyes and two anonymous reviewers for their helpful comments on the manuscript.

Funding

We thank CONACyT scholarship (851173), PAPIIT (IN215818), and CONACyT (152060), for funding that supported the study.

Author information

AR-M and AGN-S conceived and designed the study. AR-M conducted the analyses with input from LAS-G, CG, and AGN-S. AR-M wrote the manuscript with contributions from all of the authors, who read and approved the final manuscript.

Correspondence to Alberto Rocha-Méndez.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Genbank accession numbers. Samples obtained for the realization of analyses in this study.

Additional file 2: Marginal likelihoods for sequence evolution. Marginal likelihoods estimation for sequence fit to a strict or relaxed clock and comparison of evolutionary rate assumptions through Bayes factors (BFs).

Additional file 3: Hypothetical scenarios of diversification of 4 bird taxa in Mesoamerica. Scenarios tested with ABC for Aulacorhynchus, Chlorospingus, Cardellina, and Eupherusa-Thalurania with confidence scenario choice values (Type I and II errors).

Additional file 4: Posterior parameter estimates for the best supported scenarios. Population and divergence time estimates obtained with ABC for Aulacorhynchus, Chlorospingus, Cardellina, and Eupherusa-Thalurania. (DOCX 33 kb)

Additional file 5: Final tested hypothetical scenarios of diversification of 4 bird taxa in Mesoamerica. Competing demographic scenarios of Aulacorhynchus, Chlorospingus, Cardellina, and Eupherusa-Thalurania with posterior probabilities and model checking.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rocha-Méndez, A., Sánchez-González, L.A., González, C. et al. The geography of evolutionary divergence in the highly endemic avifauna from the Sierra Madre del Sur, Mexico. BMC Evol Biol 19, 237 (2019). https://doi.org/10.1186/s12862-019-1564-3

Download citation

Keywords

  • Approximate Bayesian computation
  • Coalescent
  • Comparative phylogeography
  • Mesoamerica
  • Mexico
  • Montane forest
  • mtDNA
  • Pleistocene
  • Sierra Madre del Sur