Ecological divergence and conservatism: spatiotemporal patterns of niche evolution in a genus of livebearing fishes (Poeciliidae: Xiphophorus)
© Culumber and Tobler. 2016
Received: 2 July 2015
Accepted: 13 January 2016
Published: 19 February 2016
Ecological factors often have a strong impact on spatiotemporal patterns of biodiversity. The integration of spatial ecology and phylogenetics allows for rigorous tests of whether speciation is associated with niche conservatism (constraints on ecological divergence) or niche divergence. We address this question in a genus of livebearing fishes for which the role of sexual selection in speciation has long been studied, but in which the potential role of ecological divergence during speciation has not been tested.
By combining reconstruction of ancestral climate tolerances and disparity indices, we show that the earliest evolutionary split in Xiphophorus was associated with significant divergence for temperature variables. Niche evolution and present day niches were most closely associated with each species’ geographic distribution relative to a biogeographic barrier, the Trans-Mexican Volcanic Belt. Tests for similarity of the environmental backgrounds of closely related species suggested that the relative importance of niche conservatism and divergence during speciation varied among the primary clades of Xiphophorus. Closely related species in the two swordtail clades exhibited higher levels of niche overlap than expected given environmental background similarity indicative of niche conservatism. In contrast, almost all species of platyfish had significantly divergent niches compared to environmental backgrounds, which is indicative of niche divergence.
The results suggest that the relative importance of niche conservatism and divergence differed among the clades of Xiphophorus and that traits associated with niche evolution may be more evolutionarily labile in the platyfishes. Our results ultimately suggest that the taxonomic scale of tests for conservatism and divergence could greatly influence inferences of their relative importance in the speciation process.
KeywordsPhyloclimatic modeling Geographic distribution Disparity Thermal biology
Ecological factors play a critical role in evolutionary processes by providing sources of selection that drive microevolutionary change and by imposing constraints that limit organismal performance [1–3]. Accordingly, they can shape phenotypic evolution and functional diversification and set the distributional limits for populations and species [4–6]. Ecological factors can also directly drive speciation. Ecological speciation models suggest that divergent natural selection drives the emergence of reproductive isolation as a byproduct of adaptation [7, 8]. Speciation in this context is accompanied by divergence in niche occupation along at least one axis of multidimensional niche space [9, 10]. Even when non-ecological factors drive the initial emergence of reproductive isolation, multifarious selection can stabilize and accelerate speciation processes  and lead to the secondary accumulation of ecological differences among lineages [12–15]. Alternatively, constraints imposed by ecological factors can also play a role in speciation, if evolutionary lineages become geographically isolated during periods of environmental change but retain ancestral niche affinities (niche conservatism ). In this case, vicariant speciation occurs when ancestral distributional ranges shift and break up (e.g., along elevational gradients during periods of warming or orogeny), and incipient species fail to adapt to novel environmental conditions that would facilitate the maintenance of gene flow [16, 17]. The relative contributions of niche divergence and niche conservatism during speciation processes should consequently affect the functional diversification of lineages and shape the degree of niche similarity between members of a phylogenetic group. However, the relative importance of these factors during diversification remains largely unknown. The integration of ecological modeling and biogeographic and phylogenetic analyses allows for comparative analyses of niche evolution across broad spatial scales, which can provide insight into the relationship between ecological differentiation and diversification across large taxonomic groups [18–20].
The fish family Poeciliidae has risen as a prominent model for investigating phenotypic evolution and speciation in response to sexual selection [21, 22]. Accordingly, we have a fairly sophisticated understanding of how female preferences have shaped the evolution of sexually selected traits [23–27]. There is also an emerging understanding of how natural selection mediated by physiochemical stressors and biotic interactions can drive adaptation and speciation in select groups (e.g., [28–31]). Nonetheless, the importance of ecological divergence at broader taxonomic scales and the relative contributions of ecological and sexual selection in speciation remain poorly understood . Here we focus on one of the most diverse poeciliid genera to conduct the first broad comparative study of ecological niche divergence in livebearing fishes.
Species of the genus Xiphophorus are well known as a model for investigating melanoma  and sexual selection (e.g., [34–36]). The evolutionary relationships among species have been the topic of considerable investigation and debate due to disagreement among phylogenies constructed from different data types (morphological vs. molecular: [26, 35, 37–39]), but the genus is generally comprised of three main, monophyletic clades: southern swordtails, northern swordtails, and platyfishes (Fig. 1a-c, Fig. 2a, Additional file 1). Discordance between mitochondrial and nuclear DNA phylogenies has prompted interest in hybrid speciation [40–43] and provided evidence of extensive historical gene flow among species . The role of prezygotic reproductive isolation in speciation of Xiphophorus has been widely discussed and tested largely in the context of female preferences for and against sexually selected traits and chemical cues [27, 40, 45–50], but the extent of ecological divergence and its potential role in promoting speciation has not been investigated. Here, we used ecological niche modeling in a phylogenetic context to test for the roles of niche conservatism and divergence during diversification and to examine patterns of niche evolution across the major clades within Xiphophorus.
Specifically, we tested the prediction that if niche evolution (and potentially patterns of speciation) is associated with niche conservatism, then sister and closely related species within each of the major clades of Xiphophorus should have more similar niches than expected under a null model. Alternatively, if speciation is associated with niche divergence, then closely related species should have less similar niches than expected under a null model. We first used a hypothesis-testing approach, in which a null distribution of expected niche overlap values was generated by randomly sampling points within the distributions of two focal species, in order to compare to the observed niche overlap between those species as calculated from their niche models. Our results indicated contrasting patterns of niche conservatism and divergence in the swordtails and platyfishes, respectively. To further characterize patterns of niche evolution across the genus, we employed phylogenetically controlled analyses of multivariate (principal component analysis) and univariate niche space (ancestral tolerances and accumulation of disparity). These additional analyses allowed us to illustrate patterns of niche similarities and differences across the genus and to examine cases of divergence and convergence upon similar niche characteristics for single environmental variables. Based on the results of our hypothesis testing, we predicted that species would tend to cluster with other members of their clade in multivariate space and for tolerances of single environmental variables. However, following evidence of reduced niche similarity among platyfishes, we also predicted that platyfishes would exhibit greater variation along niche dimensions than either swordtail clade.
Niche overlap & phylogeny
Niche identity and background
Across the 64 pairwise species comparisons made within the three major clades of Xiphophorus, only two had equivalent niches (Additional file 2). Niche identity tests showed that X. cortezi had an equivalent niche to that of X. pygmaeus and X. multilineatus, though the latter two did not have equivalent niches. Background tests revealed that niche conservatism was intermediate between species in each of the two major swordtail clades, with 17/28 (61 %) and 9/21 (43 %) of comparisons exhibiting significantly greater niche overlap than the null expectation (in at least one direction) in the northern and southern swordtails, respectively. The more restrictive comparison of backgrounds only between sister species showed a strong pattern of conservatism, as all comparisons exhibited greater niche overlap than the null expectation. In contrast to the swordtails, the platyfishes tended to have much lower overlap than expected based on the available background environment. Only 2/15 (13 %) tests indicated some degree of conservatism. Most importantly, no sister species of platyfishes exhibited niche overlap greater than expected based on the available background conditions.
Ancestral tolerance & phylogenetic PCA
Weighted means of the predicted niche occupancy for all environmental variables used in niche modeling and for all species of Xiphophorus included in our analyses. Worldclim variables related to temperature are reported in degrees Celsius and precipitation variables in centimeters. The Hydro1k variables retained for analyses were compound topographic index (CTI), aspect ratio (Aspect), and flow accumulation (FA). The species are grouped by their primary evolutionary clade and by their position in the phylogenetic tree (Fig. 2a). Species present in Fig. 2a but not in the table were not niche modeled due to too few occurrence points. Letters behind the platyfish species names indicate their position relative to the Trans-Mexican Volcanic Belt (TMVB) as north (N) or south (S). All northern swordtails are north and all southern swordtails south of the TMVB, respectively
As with the ATPs, the phylogenetic PCA indicated a general separation between species north and south of the TMVB, where the distributions of species south of the TMVB were influenced most by precipitation and minimum temperature in the coldest month, and those of species north of the TMVB were influenced by the remaining temperature variables (Fig. 2b-c). The phylogenetic MANOVA confirmed significant differences in the niches of species on either side of the TMVB (F 4,16: 38.4, P < 0.001). An additional pattern supported by the pPCA was the difference of the platyfishes X. maculatus and X. milleri from the rest of the platyfish clade, indicating a strong divergence in platyfish niches. Variation in maximum temperature and precipitation seasonality was related primarily to niche differences within clades, and high precipitation seasonality was particularly characteristic for the distribution of the sister species, X. mixei and X. monticolus.
Disparity plots for precipitation variables generally showed a pattern of positive MDI indicating greater evolution within subclades (e.g., bio15; Fig. 5). Precipitation seasonality (bio15) and precipitation in the warmest quarter (bio18) had the strongest signature of evolution within subclades (Additional file 4). Specifically, four nodes ranging from 0.67-0.74 relative time (~1.5 MYA), and spread across all three major clades of Xiphophorus, exhibited significant divergence for precipitation in the warmest quarter (P < 0.05). All precipitation variables showed large increases in disparity around 2.5 MYA, with three of four precipitation variables reaching their highest levels of disparity at that time point.
Evolutionary diversification of swordtails and platyfishes has been markedly shaped by sexual selection [27, 40, 50] and widespread introgressive hybridization ; yet, the importance and extent of ecological divergence in the genus Xiphophorus has remained largely unstudied. Our analysis of ecological niches across this group indicated contrasting patterns of niche evolution among the three major lineages. Hypothesis testing revealed a strong tendency of niche conservatism as indicated by significant similarity in niches compared to background environments among closely related species of the northern and southern swordtail clades. In contrast, species in the platyfish clade have undergone a high degree of niche divergence. These general patterns were well supported in all analyses and were discernible in tests of multivariate niche space (e.g., phylogenetic PCA) as well as in single environmental variables (e.g., ancestral tolerances). Furthermore, patterns of niche variation in Xiphophorus were strongly bisected by the Trans-Mexican Volcanic Belt, illustrating an important association between biogeography and patterns of niche evolution at broader taxonomic scales. This pattern was particularly evident in platyfishes on either side of the TMVB, which showed greater niche similarities to and shared ancestral tolerances with swordtails exhibiting a similar geographic distribution, rather than with phylogenetically related platyfishes. The consistency of the results from multiple analyses of niche evolution in Xiphophorus indicate that both divergence and constraint have accompanied diversification, which highlights avenues for future research on the roles of ecological adaptation and constraint in speciation. Our results also support previous studies in insects, amphibians, birds, and mammals in underscoring the importance of taxonomic scale for making inferences about divergence and conservatism in niche evolution [19, 51]. For example, a focus on either swordtail group alone would have suggested high niche conservatism during the diversification of Xiphophorus fishes, while a focus on platyfishes would have provided strong evidence towards niche divergence. In reality, the relative importance of niche conservatism and divergence during diversification can vary even among closely related species.
Has niche divergence accompanied speciation?
The extent and importance of ecological divergence during speciation remains an open question in ecology and evolution . Studies of climatic niche evolution have reported evidence of both divergence [20, 53] and conservatism [19, 51, 54, 55]. We uncovered a significant correlation between phylogenetic distance and niche overlap. This pattern is consistent with the accumulation of ecological differences over evolutionary time. However, this does not in itself allow for inferences to be made about the roles of niche conservatism and divergence during diversification. Examining the outliers from this correlation indicated that nodes exhibiting greater niche overlap than expected based on their age (i.e. potential conservatism) only occurred within swordtail clades. The nodes exhibiting less overlap than expected (i.e., potential divergence) mostly belonged to the platyfish clade (4 of 5 platyfish nodes had reduced overlap). Potential niche conservatism in swordtails and niche divergence in platyfishes was also supported by other analyses. Background similarity tests indicated that closely related species had higher levels of niche overlap than expected by chance, with 63 % of swordtail sister species exhibiting significant niche conservatism. Conservatism remained high at broader taxonomic scales within these two clades as 61 % (northern swordtails) and 43 % (southern swordtails) of pairwise comparisons indicated significantly greater overlap compared to a null model, and only 2 (7 %) and 3 (14 %) tests, respectively, indicated significantly divergent niches. However, niche identity tests revealed that only two pairwise comparisons within northern swordtails had equivalent niches. Low niche equivalency indicates that while swordtails are characterized by niche conservatism (niches are significantly more similar than their environmental backgrounds), speciation has been accompanied by shifts in climatic niches. In contrast to the swordtails, platyfishes exhibited no equivalent niches and almost no conservatism. Among comparisons restricted to sister species of platyfishes, none of the tests indicated niche conservatism. Across all pairwise comparisons of platyfishes, only 2/15 (13 %) comparisons had niche overlap values greater than expected, but 9/15 (60 %) exhibited significant niche divergence with less niche overlap compared to the null model. Our results for the two swordtail clades tend to coincide with previous work showing that 80-90 % of sister species exhibit niche conservatism in analyses of Mexican butterflies, birds, and mammals . In that context, it becomes evident that diversification in platyfishes is characterized by much greater niche diversification, perhaps suggesting that adaptation to divergent environmental conditions is an important driver of platyfish speciation. Our results are similar to findings on other North American fishes, which generally suggest that some clades are characterized by niche conservatism while others are not . Similarly, asymmetries in the degree of niche divergence and conservatism have also been observed in other non-fish systems, albeit at broader taxonomic scales than those in this study. For example, estimates of niche conservatism based on geographic overlap of terrestrial vertebrates and conservatism in habitat affinities in marine invertebrates have each exhibited variation among very divergent taxa [57, 58]. Our data show that such asymmetries among taxa can be apparent at much smaller geographic and taxonomic scales.
The cause for the contrasting evidence for niche divergence and conservatism in different clades of Xiphophorus fishes remains unclear. Given that the platyfishes broadly overlap the geographic distributions of swordtails and that both swordtail clades are characterized by high conservatism suggests that intrinsic rather extrinsic factors may play a role in niche evolution for these respective groups. In general, there may be three non-mutually exclusive hypotheses that could explain the observed pattern: (1) The documented differences among clades may be an artifact of incomplete taxon sampling. The only taxa not included in our analyses are those that are highly endemic with too few, independent occurrence points to conduct niche modeling. It is unlikely that we have underestimated divergence in the two swordtail clades, because our analyses were only missing one species each from the northern (X. continens) and the southern (X. kallmani) swordtail clades. Xiphophorus continens is closely related to X. pygmaeus [39, 44], and like all other northern swordtails, it occurs in the Rio Panuco drainage . Xiphophorus kallmani is endemic to the Laguna de Catemaco, which is within the geographic distribution of X. hellerii . Our analyses were missing three species of platyfishes. Xiphophorus gordoni and X. meyeri are closely related to X. couchianus comprising a single, monophyletic clade. All three species occur in tributaries of the Rio Grande basin and have restricted ranges in an arid region of northern Mexico that is geographically removed from the remainder of the genus . Accordingly, the three species may inhabit similar environments and exhibit similar niches to one another, but as a group are likely divergent from other platyfishes. Finally, X. andersi is a more basal platyfish [39, 44] that occurs south of the TMVB, but inhabits a high-elevation, headwater habitat that is distinct from its closest relatives X. maculatus and X. xiphidium . (2) Differences in the degree of niche conservatism between the swordtails and platyfishes could be a consequence of the environmental variables included in this study. Niche modeling focused on parameters that characterize climatic niches based on distributional patterns. Even in the absence of niche divergence in climatic niche space among northern and southern swordtail species, ecological divergence may have occurred along other biotic or abiotic niche dimensions . This is a particularly intriguing problem for the northern swordtails. All northern swordtail species occur in the Rio Panuco drainage, and some species occur sympatrically within the same habitats . This provides an opportunity for future studies to test for potential ecological differentiation among species along other niche axes. (3) The relative importance of niche divergence and conservatism during diversification may be related to the strength of sexual selection. Swordtails commonly have exaggerated male traits, including elaborate caudal extensions  and enlarged dorsal fins . In addition, there is a wealth of evidence for the presence of female mate choice for a variety of male traits in swordtails [25, 34–36, 62, 63], but see . Nonetheless, it is difficult to say whether platyfishes are characterized by weaker sexual selection. Even though platyfishes generally do not have the sword ornament, females exhibit preferences for swords , enlarged dorsal fins [64, 66], and platyfishes tend to have a much greater variety of genetic pattern and color traits compared to swordtails [67, 68]. Whether differences in niche divergence and conservatism within Xiphophorus is due to a constraining effect of sexual selection on niche evolution in swordtails remains to be investigated, but other studies have found evidence of natural and sexual selection reinforcing or constraining one another during diversification [69–71].
Biogeography and climatic niches
A consistent trend observed in multiple analyses was the broad separation in the niches of species on either side of the TMVB. For example, phyloclim analyses revealed that ancestral tolerances of platyfishes overlapped those of both swordtail clades. Specifically, the two southern platyfish, X. maculatus and X. milleri, tended to have tolerances similar to the southern swordtails (for example, Fig. 4a). Geographic isolation during speciation likely facilitates the accumulation of ecological differences. Biogeography and associated vicariance is therefore important for understanding evolutionary patterns of species niches [72, 73]. Studies in plants have suggested that biogeographic breaks in species distributions promote diversification and increases in trait disparity [20, 74], and our results generally agree with such patterns. The TMVB is a major biogeographic break for freshwater fishes [75, 76] and was associated with strong divergence between the niches of the ancestral southern swordtails and the northern swordtails as well as divergence within the platyfishes. In addition to ancestral tolerances, this pattern was also recovered in the phylogenetic PCA and MANOVA analyzing variation in multivariate niche space (Fig. 2c).
Previous work on the role of biogeography in niche evolution in Mexican butterflies, birds, and mammals found considerable niche conservatism between sister species on either side of the Isthmus of Tehuantepec, an arid region separating more moist, montane regions [19, 77]. Though a direct comparison to those studies is difficult here, particularly because sister species of Xiphophorus generally do not inhabit both sides of the TMVB, our data revealed significant differences in the niches of species living on opposite sides of the TMVB even after controlling for phylogenetic distance. If niche conservatism were a strong factor for niche evolution of sister species, we would have expected X. maculatus to have stronger niche similarity with X. xiphidium, and X. couchianus to have similarity with X. variatus, their closest phylogenetic relatives. Instead, X. maculatus and X. couchianus exhibited significant similarity only with their closest geographic congeners. Consequently, major biogeographic barriers like TMVB may be critical in shaping the evolution of ecological niches.
Lineages in unique environments
Despite high degrees of niche conservatism in swordtails, niche identity comparisons revealed that closely related species almost never have equivalent niches. Analyses of the accumulation of disparity and ancestral tolerances provided insight into specific instances of when and for what variables closely related species have diverged, and several examples of divergence and convergence were also apparent in multivariate space (phylogenetic PCA): (1) Xiphophorus monticolus and X. mixei have diverged from their sister species, X. clemenciae, for all precipitation variables despite close geographic proximity. (2) The northernmost platyfishes, X. couchianus and X. xiphidium, exhibited divergent tolerances for several temperature and precipitation variables compared to the rest of the platyfishes and were clustered in a unique area of multivariate space. (3) The niches of the northern swordtail clade containing X. pygmaeus and its closest relatives have diverged towards higher maximum temperatures (bio5; Additional file 3). (4) Most interestingly, members of each of the main clades (X. alvarezi in the southern swordtails, X. malinche in the northern swordtails, and X. evelynae in the platyfishes) have converged on similar high-elevation niches characterized by cooler maximum temperatures (bio5; Additional file 3). The convergence between X. evelynae and X. malinche was particularly apparent in the pPCA, and the strong divide in niche space between species north and south of the TMVB is likely the reason that X. alvarezi did not cluster with these two species, even though it was clearly divergent from its closest relatives (X. mayae and X. signum).
Species that occupy distinct portions of niche space will provide unique models to understand the relative roles of adaptation and plasticity that underlie differences in observed niches based on distributional data. This remains a major task for the field, because, adaptive divergence cannot be inferred from correlative modeling approaches alone. One well-documented case in Xiphophorus suggests that species are adapted to local climatic conditions. Replicated hybrid zones between X. birchmanni and X. malinche occur along an elevational and thermal gradient in multiple streams of the Rio Panuco drainage . Physiological tolerances to thermal extremes in wild-caught and common-garden fish together with gene expression analyses demonstrated the role of physiological adaptation to local thermal environments as a key factor in the maintenance of the geographic structure of hybrid zones .
It remains unclear whether adaptation to novel niche space promotes vicariance or whether vicariance followed by adaptation accelerates or facilitates the maintenance of reproductive isolation. A study on North American Monkeyflowers suggests that speciation is often initiated by small populations occupying novel niche dimensions , typically through colonization of the new environment followed by reproductive isolation (budding speciation). Budding speciation could therefore be one mechanism by which ecological divergence occurs. The uplift of mountain ranges, which slowly isolates upstream populations from downstream congeners, may be a special case of budding speciation. Such isolation by orogeny has been implicated for some species of Xiphophorus  and may explain replicated instances of divergence in novel niche dimensions.
The relationship between niche evolution and speciation remains widely studied, because conservatism and divergence might equally promote the emergence of new species . In our study, the degree of niche conservatism observed across the genus Xiphophorus varied markedly among clades. An equally interesting and important question stems from this pattern: what gives rise to such asymmetry in niche evolution within a single genus? Additional studies are needed to better understand the relative frequency of niche divergence and conservatism during speciation and the extent to which heterogeneous patterns of niche evolution occur in other taxa as we observed in Xiphophorus. Future studies will need to combine analyses of species distributions and ecological niches with functional assays of relevant physiological and morphological traits that ultimately shape organismal fitness. This is a substantial task considering that a persistent problem for fishes and other aquatic organisms is to model the functional ecological and physiological niche based on landscape level environmental data that has more predictable effects on organismal performance. Whereas temperature variation can give rise to explicit testable hypotheses related to temperature adaptation, the effects of many variables that are currently available and commonly used on ecological niche modeling for aquatic organisms, including various precipitation metrics, are often difficult to interpret, and most parts of the world do not have the necessary infrastructure nor historical data records to consider the use of more relevant environmental variables, such as dissolved oxygen concentrations, in models that span large spatial and taxonomic scales.
Distributional patterns and occurrence records
The most basal divergence in Xiphophorus gave rise to the southern swordtails and the common ancestor of platyfishes and northern swordtails . The southern swordtails have the most southern range and occur from southern Mexico to northern Guatemala and Honduras. There is considerable variation in the geographic ranges of individual species. Within southern swordtails, species of the X. clemenciae clade (X. clemenciae, X. mixei, and X. monticolus) are only found in upland reaches of the Rio Coatzacoalcos in Mexico. Four additional species (X. alvarezi, X. mayae, X. signum, and X. kallmani) are scattered across several drainages from southern Mexico through northern Honduras, but like members of the X. clemenciae clade, each has a relatively restricted range. In contrast to all other southern swordtails, X. hellerii is widely distributed from the Trans-Mexican Volcanic Belt (TMVB) to Honduras. The northern swordtails are found to the north of the TMVB and are almost entirely restricted to the Rio Panuco drainage. Each species has a relatively restricted distribution in different sub-drainages, or they segregate along elevational gradients within sub-drainages . The TMVB therefore broadly defines the geographic separation between the northern and southern swordtails. Interestingly, the platyfishes span both sides of the TMVB. The platyfish, X. maculatus, is basal to other species of platyfishes and is widely distributed to the south of the TMVB, suggesting that southern Mexico was a likely origin for the platyfishes. Two additional platyfishes inhabit the region south of the TMVB; X. milleri occurs in the Laguna de Catemaco and its tributaries (sympatric with X. kallmani also endemic to Catemaco), and X. andersi is known from a single locality in the state of Veracruz. More recently derived platyfishes (X. couchianus, X. gordoni, and X. meyeri) have colonized spring habitats of the Rio Grande drainage in the northern Mexican states of Coahuila and Nuevo Leon, which are geographically isolated from other species of Xiphophorus. The remaining platyfishes north of the TMVB vary in their distributions, ranging from highly restricted headwaters of a single basin (X. evelynae) to numerous streams of a single drainage (X. xiphidium), and being widely distributed from the Rio Nautla to the Rio Soto La Marina (X. variatus). Herein, we follow the vernacular of Kallman and Kazianis  in referring to northern and southern platyfishes based on their geographic distributions relative to the TMVB.
We obtained species occurrence data for all Xiphophorus species from the Fishnet2 Portal (http://fishnet2.net), the literature [59, 60, 81–86], and field collections in Mexico and Honduras conducted by the authors. All occurrence records were vetted against the literature to ensure that occurrences for each species fell within their documented, native distributions [59, 60], and non-native records were removed from the dataset. For each species, we also removed occurrence records that were <1 km apart to match the spatial resolution of environmental variables (see below).
In order to obtain a time-calibrated topology for phlyoclimatic analyses, we used the maximum clade credibility topology of Jones et al. , which was derived from genome-wide restriction-site associated DNA sequencing. This phylogeny includes all 26 described swordtail and platyfish species. Prior to time calibration, taxa with too few occurrence records (<5) for ecological niche modeling were manually pruned using Mesquite v2.91 . Divergence times were estimated using the ‘chronos’ command in the package ape  in R. In this method divergence times are estimated with a penalized likelihood approach, and a correlated substitution rate was employed . Generating a time-calibrated phylogeny for Xiphophorus is difficult due to a lack of fossil calibration points. Discordance between the mitochondrial and nuclear phylogenies in Xiphophorus, likely due to historical gene flow and introgressive hybridization , also precludes use of a mitochondrial phylogeny and molecular clock dating. Although the relative divergence times within the genus Xiphophorus is most important for the analyses and interpretations herein, the estimated absolute divergence times are in general agreement with several lines of evidence.
The timing of geological events in Mexico is relatively well known, including the closure of the TMVB, a key geographic barrier to dispersal. The final uplift of the TMVB towards Punta El Morro on the Gulf Coast of Mexico occurred approximately 5 million years ago (MYA). Prior studies of Mexican fishes have used the timing of this geological event to calibrate divergence times in cichlids and other poeciliids [75, 76, 90], due to the fact that it represents a major biogeographic break in the distributions of many fish taxa, including Xiphophorus . We therefore constrained the age of the initial divergence in Xiphophorus to a minimum of 4.5 MYA and a maximum of 5 MYA, assuming correlated substitution rates among clades. This approach placed the divergence of X. milleri, endemic to the Laguna de Catemaco, between 1.5–2 MYA (Additional file 1), which is consistent with estimated geological dates for the lake . The estimated age of X. milleri was also consistent with the estimated age of another livebearer endemic to the lake, Poeciliopsis catemaco, as calculated using a molecular clock approach with mitochondrial DNA .
Ecological niche modeling
We assembled a set of coverages for 24 environmental variables from the Worldclim (http://www.worldclim.org) and Hydro1k (https://lta.cr.usgs.gov/HYDRO1K) project databases. Data layers were downloaded at 30 arc-seconds (~1 km2) resolution. The Worldclim dataset is composed of nineteen bioclimatic variables related to temperature and precipitation, and the Hydro1k dataset is composed of five variables related to hydrography [93, 94]. Variables in these datasets have been successfully used in previous studies analyzing distributions of freshwater fish species with ecological niche modeling, including several poeciliids [79, 95–97]. All data layers were clipped to the area of study encompassing the native distribution of the species included in this study (see Fig. 1a-c).
Niche evolution analyses
Niche models were used to calculate niche overlap among all pairwise species based on Schoener’s D and Warren’s I  in R using the package phyloclim (http://www.christophheibl.de/Rpackages), which was used for all subsequent analyses unless otherwise stated. Both D and I range from 0 (no overlap) to 1 (complete overlap). Because low niche overlap may be a correlated with phylogenetic distance as a consequence of neutral evolution rather than niche differentiation through divergent selection, we tested for phylogenetic signal in patterns of niche overlap using the age.range.correlation (ARC) command and used SPSS v17.0 (IBM Corp.) to generate 95 % confidence intervals around the correlation. This procedure applies a linear regression of node age to niche overlap to test whether there is a relationship between time since divergence and niche overlap at all nodes of the phylogeny. No correlation would indicate a random association between node age and niche overlap, whereas either a significant negative or positive correlation indicate phylogenetic signal in niche evolution. While the correlation itself does not indicate niche conservatism or divergence, examining outliers can help identify nodes that exhibit greater or lesser niche overlap than expected for their age.
In order to explicitly evaluate the influence of niche conservatism and divergence during Xiphophorus diversification, we tested hypotheses about niche overlap using the niche identity and background tests . These tests quantify niche similarity between species in the context of a null model, enabling inferences of niche conservatism and divergence. Because tests for niche conservatism are most appropriate for investigating potential shifts in ecological niche dimensions among closely related species, we restricted our tests to pairwise comparisons within each of the three primary clades. Niche identity tests evaluate whether the niches of a given pair of species are identical. However, a rejection of niche identity between two species may be driven by differences in available environmental conditions within the distributions of the two species. Background similarity tests provide an additional comparison by testing whether observed niche overlap between two species is greater than the overlap expected based on the environmental conditions available to each species (i.e., niche conservatism). Alternatively, observed overlap may be less than that expected base on the environmental background of the species (i.e., niche divergence). Background similarity analyses generate random “pseudolocality” points within the observed distributions of two species and then test whether the observed niche overlap is more or less than would be expected given the available environmental backgrounds of the two species. As in previous studies [104–106], we used 100 replicates and a confidence interval of 0.95 to assess significance.
In order to trace the evolutionary history of species niches through time, we produced predicted niche occupancy (PNO) profiles . For a given species, this approach uses probability distributions output by MAXENT to quantify the occupancy of that species across all niche space (with the probabilities summing to 1). This process was repeated to generate PNOs for each species and all environmental variables. The PNO profiles were then used to reconstruct the ancestral tolerances for each environmental variable. These ancestral tolerances can be plotted similar to a phylogenetic tree and allow inference of divergent and convergent evolution along a particular niche axis. The ancestral tolerances were constructed using a maximum likelihood approach under a model of Brownian motion evolution and using 100 random samples from the PNO profiles of all taxa as implemented in phyloclim. In contrast to prior methods that reconstructed ancestral tolerances with respect to only extreme (maximum or minimum) or mean values of environmental variables, this approach reconstructs the entire distribution of ancestral environmental tolerances and can better account for intraspecific niche variation during ancestral state reconstruction . Finally, we calculated the weighted means of the PNOs for each species and environmental variable. Because ancestral tolerance analyses produce a separate plot for every environmental variable, we used a phylogenetic principal components analysis (pPCA: ) on the PNO weighted means to summarize occupancy of multivariate niche space across the genus using the package phytools . The weighted means of predicted niche occupancy (Table 1), ancestral tolerance plots (Fig. 4 and Additional file 3), and a preliminary PCA all indicated that species niches exhibited relatively little variation in hydrographic variables, rendering them uninformative for explaining differences among species. These variables were therefore excluded from the final PCA. To examine patterns of niche evolution associated with the TMVB biogeographic divide, we conducted a phylogenetic multivariate analysis of variance (MANOVA) using the first four phylogenetic principal components obtained from the phlogenetic PCA. The analysis was run in R with the package geiger  by comparing observed effects of position relative to the TMVB (north or south) to a null distribution of the F-statistic generated from 1000 simulations under a Brownian motion model.
As a final approach to examine niche evolution, we constructed relative disparity plots . This approach quantifies the extent to which disparity in the niche is distributed within or among subclades and illustrates temporal patterns in the accumulation of niche disparity among clades. The x-axis of a disparity plot indicates relative evolutionary time ranging from 0 (root of the phylogeny) to 1 (present day). The y-axis describes the observed disparity and the disparity expected under an unconstrained, Brownian model of evolution. Disparity was calculated as the average squared pairwise distance between taxa with respect to the mean predicted niche occupancy for each interior node and standardized by the disparity observed across the entire clade. We then calculated a disparity index based on environmental variables following the definition of the Morphological Disparity Index (MDI: [110, 111]), which allows comparison of the observed disparity values to values expected under a Brownian motion model of evolution, using 1000 simulations and a confidence interval of 0.95. Negative MDI values ─ where disparity is less than the expected value under the Brownian model on the disparity plots ─ indicate that disparity is distributed among subclades, consistent with conservatism within more inclusive clades and divergence among more inclusive clades. Positive MDI values indicate that disparity is within subclades, consistent with divergence among subclades. In order to test whether any nodes exhibited MDI that deviated significantly from its simulated null distribution, we calculated the rank of the observed disparity value for the node within the simulated values as described in Swenson . Disparity and MDI analyses were conducted in R using the package geiger .
We would like to thank R. Morehouse and R. Raghavan for assistance during niche modeling, as well as A. Hope and D. Shepard for helpful comments that improved the manuscript. This work was supported by an American Philosophical Society Franklin Research Grant to ZWC and National Science Foundation grant (IOS-1121832) to MT. Publication of this article was funded in part by the Kansas State University Open Access Publishing Fund.
Open AccessThis 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.
- Niewiarowski PH, Roosenburg W. Reciprocal transplant reveals sources of variation in growth rates of the lizard Sceloporus undulatus. Ecology. 1993;74(7):1992–2002.View ArticleGoogle Scholar
- Angert AL, Schemske DW. The evolution of species’ distributions: reciprocal transplants across the elevation ranges of Mimulus cardinalis and M. lewisii. Evolution. 2005;59(8):1671–84.PubMedView ArticleGoogle Scholar
- Barrett RDH, Hoekstra HE. Molecular spandrels: tests of adaptation at the genetic level. Nat Rev Genet. 2011;12(11):767–80.PubMedView ArticleGoogle Scholar
- Sexton JP, McIntyre PJ, Angert AL, Rice KJ. Evolution and ecology of species range limits. Ann Rev Ecol Evol Syst. 2009;40:415–36. doi:https://doi.org/10.1146/annurev.ecolsys.110308.120317.View ArticleGoogle Scholar
- Arbour JH, López-Fernández H. Adaptive landscape and functional diversity of Neotropical cichlids: implications for the ecology and evolution of Cichlinae (Cichlidae; Cichliformes). J Evol Biol. 2014;27(11):2431–42.PubMedView ArticleGoogle Scholar
- Bossu CM, Near TJ. Ecological constraint and the evolution of sexual dichromatism in darters. Evolution. 2015;69:1219–31.PubMedView ArticleGoogle Scholar
- Schluter D. Ecology and the origin of species. Trends Ecol Evol. 2001;16(7):372–80.PubMedView ArticleGoogle Scholar
- Rundle HD, Nosil P. Ecological speciation. Ecol Lett. 2005;8(3):336–52. doi:https://doi.org/10.1111/j.1461-0248.2004.00715.x.View ArticleGoogle Scholar
- Cox SC, Prys-Jones RP, Habel JC, Amakobe BA, Day JJ. Niche divergence promotes rapid diversification of East African sky island white-eyes (Aves: Zosteropidae). Mol Ecol. 2014;23(16):4103–18.PubMed CentralPubMedView ArticleGoogle Scholar
- Gómez-Rodríguez C, Baselga A, Wiens JJ. Is diversification rate related to climatic niche width? Global Ecol Biogeogr. 2014;24(4):383–95.View ArticleGoogle Scholar
- Chevin L-M, Decorzent G, Lenormand T. Niche dimensionality and the genetics of ecological speciation. Evolution. 2014;68(5):1244–56.PubMedView ArticleGoogle Scholar
- Owens IPF, Bennett PM, Harvey PH. Species richness among birds: body size, life history, sexual selection or ecology? Proc R Soc B. 1999;266(1422):933–9.PubMed CentralView ArticleGoogle Scholar
- Arnqvist G, Edvardsson M, Friberg U, Nilsson T. Sexual conflict promotes speciation in insects. Proc Natl Acad Sci. 2000;97(19):10460–4.PubMed CentralPubMedView ArticleGoogle Scholar
- Mank JE. Mating preferences, sexual selection and patterns of cladogenesis in ray-finned fishes. J Evol Biol. 2007;20(2):597–602.PubMedView ArticleGoogle Scholar
- Burke JM, Arnold ML. Genetics and the fitness of hybrids. Annu Rev Genet. 2001;35:31–52.PubMedView ArticleGoogle Scholar
- Wiens JJ. What is speciation and how should we study it? Am Nat. 2004;163(6):914–23.PubMedView ArticleGoogle Scholar
- Shepard DB, Burbrink FT. Lineage diversification and historical demography of a sky island salamander, Plethodon ouachitae, from the Interior Highlands. Mol Ecol. 2008;17(24):5315–35.PubMedView ArticleGoogle Scholar
- Pearman PB, Guisan A, Broennimann O, Randin CF. Niche dynamics in space and time. Trends Ecol Evol. 2008;23(3):149–58.PubMedView ArticleGoogle Scholar
- Warren DL, Glor RE, Turelli M. Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution. Evolution. 2008;62(11):2868–83. doi:https://doi.org/10.1111/j.1558-5646.2008.00482.x.PubMedView ArticleGoogle Scholar
- Evans MEK, Smith SA, Flynn RS, Donoghue MJ. Climate, niche evolution, and diversification of the “Bird-Cage” evening primroses (Oenothera, Sections Anogra and Kleinia). Am Nat. 2009;173(2):225–40.PubMedView ArticleGoogle Scholar
- Evans JP, Pilastro A, Schlupp I. Ecology and evolution of poeciliid fishes. Chicago: University of Chicago Press; 2011.View ArticleGoogle Scholar
- Pollux BJA, Meredith RW, Springer MS, Garland T, Reznick DN. The evolution of the placenta drives a shift in sexual selection in livebearing fish. Nature. 2014;513(7517):233–6.PubMedView ArticleGoogle Scholar
- Kodric-Brown A. Female preference and sexual selection for male coloration in the guppy (Poecilia reticulata). Behav Ecol Sociobiol. 1985;17(3):199–205.View ArticleGoogle Scholar
- Ptacek MB, Childress MJ, Petersen JA, Tomasso AO. Phylogenetic evidence for the gain and loss of a sexually selected trait in sailfin mollies. ISRN Zool. 2011;2011.Google Scholar
- Robinson DM, Tudor MS, Morris MR. Female preference and the evolution of an exaggerated male ornament: the shape of the preference function matters. Anim Behav. 2011;81(5):1015–21.View ArticleGoogle Scholar
- Meyer A. The evolution of sexually selected traits in male swordtail fishes (Xiphophorus: Poeciliidae). Heredity. 1997;79(3):329–37.View ArticleGoogle Scholar
- Wong BBM, Rosenthal GG. Female disdain for swords in a swordtail fish. Am Nat. 2006;167(1):136–40.PubMedView ArticleGoogle Scholar
- Tobler M, Palacios M, Chapman LJ, Mitrofanov I, Bierbach D, Plath M, et al. Evolution in extreme environments: replicated phenotypic differentiation in livebearing fish inhabiting sulfidic springs. Evolution. 2011;65:2213–28.PubMedView ArticleGoogle Scholar
- Bassar RD, Ferriere R, Lopez-Sepulcre A, Marshall MC, Travis J, Pringle CM, et al. Direct and indirect ecosystem effects of evolutionary adaptation in the Trinidadian guppy (Poecilia reticulata). Am Nat. 2012;180(2):167–85.PubMedView ArticleGoogle Scholar
- Ingley SJ, Billman EJ, Belk MC, Johnson JB. Morphological divergence driven by predation environment within and between species of Brachyrhaphis fishes. PLoS One. 2014;9(2), e90274.PubMed CentralPubMedView ArticleGoogle Scholar
- Martin RA, Riesch R, Heinen-Kay JL, Langerhans RB. Evolution of male coloration during a post-pleistocene radiation of Bahamas mosquitofish (Gambusia hubbsi). Evolution. 2014;68(2):397–411.PubMedView ArticleGoogle Scholar
- Endler JA. Integrative commentary on ecology and evolution of poeciliid fishes. In: Evans JP, Pilastro A, Schlupp I, editors. Ecology and evolution of poeciliid fishes. Chicago: University of Chicago Press; 2011. p. 301–10.Google Scholar
- Schartl M. Evolution of Xmrk: an oncogene, but also a speciation gene? Bioessays. 2008;30(9):822–32.PubMedView ArticleGoogle Scholar
- Ryan MJ, Hews DK, Wagner Jr WE. Sexual selection on alleles that determine body size in the swordtail Xiphophorus nigrensis. Behav Ecol Sociobiol. 1990;26(4):231–7.View ArticleGoogle Scholar
- Meyer A, Morrissey JM, Schartl M. Recurrent origin of a sexually selected trait in Xiphophorus fishes inferred from a molecular phylogeny. Nature. 1994;368(6471):539–42.PubMedView ArticleGoogle Scholar
- Kingston JJ, Rosenthal GG, Ryan MJ. The role of sexual selection in maintaining a colour polymorphism in the pygmy swordtail, Xiphophorus pygmaeus. Anim Behav. 2003;65:735–43.View ArticleGoogle Scholar
- Rauchenberger M, Kallman KD, Moritzot DC. Monophyly and geography of the Rio Pánuco basin swordtails (genus Xiphophorus) with descriptions of four new species. Am Mus Novitates. 1990;2975.Google Scholar
- Marcus JM, McCune AR. Ontogeny and phylogeny in the northern swordtail clade of Xiphophorus. Syst Biol. 1999;48(3):491–522.View ArticleGoogle Scholar
- Kang JH, Schartl M, Walter RB, Meyer A. Comprehensive phylogenetic analysis of all species of swordtails and platies (Pisces: Genus Xiphophorus) uncovers a hybrid origin of a swordtail fish, Xiphophorus monticolus, and demonstrates that the sexually selected sword originated in the ancestral lineage of the genus, but was lost again secondarily. BMC Evol Biol. 2013;13(1):25.PubMed CentralPubMedView ArticleGoogle Scholar
- Meyer A, Salzburger W, Schartl M. Hybrid origin of a swordtail species (Teleostei : Xiphophorus clemenciae) driven by sexual selection. Mol Ecol. 2006;15(3):721–30.PubMedView ArticleGoogle Scholar
- Jones JC, Perez-Sato J-A, Meyer A. A phylogeographic investigation of the hybrid origin of a species of swordtail fish from Mexico. Mol Ecol. 2012;21(11):2692–712.PubMedView ArticleGoogle Scholar
- Schumer M, Cui R, Boussau B, Walter R, Rosenthal G, Andolfatto P. An evaluation of the hybrid speciation hypothesis for Xiphophorus clemenciae based on whole genome sequences. Evolution. 2012;67(4):1155–68.PubMed CentralPubMedView ArticleGoogle Scholar
- Jones JC, Fan S, Franchini P, Schartl M, Meyer A. The evolutionary history of Xiphophorus fish and their sexually selected sword: a genome-wide approach using restriction site-associated DNA sequencing. Mol Ecol. 2013;22(11):2986–3001.PubMedView ArticleGoogle Scholar
- Cui R, Schumer M, Kruesi K, Walter R, Andolfatto P, Rosenthal G. Phylogenomics reveals extensive reticulate evolution in Xiphophorus fishes. Evolution. 2013;67(8):2166–79.PubMedView ArticleGoogle Scholar
- Crapon De Caprona M, Ryan MJ. Conspecific mate recognition in swordtails, Xiphophorus nigrensis and X. pygmaeus (Poeciliidae): olfactory and visual cues. Anim Behav. 1990;39(2):290–6.View ArticleGoogle Scholar
- Basolo A. Female preference for male sword length in the green swordtail, Xiphophorus helleri (Pisces, Poeciliidae). Anim Behav. 1990;40:339–49.View ArticleGoogle Scholar
- McLennan DA, Ryan MJ. Responses to conspecific and heterospecific olfactory cues in the swordtail Xiphophorus cortezi. Anim Behav. 1997;54(5):1077–88.PubMedView ArticleGoogle Scholar
- Hankison SJ, Morris MR. Sexual selection and species recognition in the pygmy swordtail, Xiphophorus pygmaeus: conflicting preferences. Behav Ecol Sociobiol. 2002;51(2):140–5.View ArticleGoogle Scholar
- Rosenthal GG, Wagner WE, Ryan MJ. Secondary reduction of preference for the sword ornament in the pygmy swordtail Xiphophorus nigrensis (Pisces: Poeciliidae). Anim Behav. 2002;63(1):37–45.View ArticleGoogle Scholar
- Morris MR, Moretz JA, Farley K, Nicoletto P. The role of sexual selection in the loss of sexually selected traits in the swordtail fish Xiphophorus continens. Anim Behav. 2005;69:1415–24.View ArticleGoogle Scholar
- Kozak KH, Wiens JJ. Does niche conservatism promote speciation? A case study in North American salamanders. Evolution. 2006;60(12):2604–21.PubMedView ArticleGoogle Scholar
- Nosil P, Harmon LJ, Seehausen O. Ecological explanations for (incomplete) speciation. Trends Ecol Evol. 2009;24(3):145–56.PubMedView ArticleGoogle Scholar
- Knouft JH, Losos JB, Glor RE, Kolbe JJ. Phylogenetic analysis of the evolution of the niche in lizards of the Anolis sagrei group. Ecology. 2006;87(7):S29–38.PubMedView ArticleGoogle Scholar
- Crisp MD, Arroyo MTK, Cook LG, Gandolfo MA, Jordan GJ, McGlone MS, et al. Phylogenetic biome conservatism on a global scale. Nature. 2009;458(7239):754–6.PubMedView ArticleGoogle Scholar
- Kozak KH, Wiens JJ. Niche conservatism drives elevational diversity patterns in Appalachian salamanders. Am Nat. 2010;176(1):40–54.PubMedView ArticleGoogle Scholar
- McNyset KM. Ecological niche conservatism in North American freshwater fishes. Biol J Linn Soc. 2009;96(2):282–95. doi:https://doi.org/10.1111/j.1095-8312.2008.01121.x.View ArticleGoogle Scholar
- Smith BT, Bryson RW, Houston DD, Klicka J. An asymmetry in niche conservatism contributes to the latitudinal species diversity gradient in New World vertebrates. Ecol Lett. 2012;15(11):1318–25.PubMedView ArticleGoogle Scholar
- Hopkins MJ, Simpson C, Kiessling W. Differential niche dynamics among major marine invertebrate clades. Ecol Lett. 2014;17(3):314–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller RR, Minckley WL, Norris SM, Gach MH. Freshwater fishes of Mexico. Chicago: University of Chicago Press; 2005.Google Scholar
- Kallman KD, Kazianis S. The genus Xiphophorus in Mexico and Central America. Zebrafish. 2006;3:271–85.PubMedView ArticleGoogle Scholar
- Winemiller KO, Fitzgerald DB, Bower LM, Pianka ER. Functional traits, convergent evolution, and periodic tables of niches. Ecol Lett. 2015;18(8):737–51.PubMedView ArticlePubMed CentralGoogle Scholar
- Rosenthal GG, Evans CS. Female preference for swords in Xiphophorus helleri reflects a bias for large apparent size. Proc Natl Acad Sci. 1998;95(8):4431–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Morris MR, Rios-Cardenas O, Brewer J. Variation in mating preference within a wild population influences the mating success of alternative mating strategies. Anim Behav. 2010;79(3):673–8.View ArticleGoogle Scholar
- MacLaren RD, Gagnon J, He R. Female bias for enlarged male body and dorsal fins in Xiphophorus variatus. Behav Processes. 2011;87(2):197–202.PubMedView ArticleGoogle Scholar
- Basolo AL. Female preference predates the evolution of the sword in swordtail fish. Science. 1990;250(4982):808–10.PubMedView ArticleGoogle Scholar
- Culumber ZW, Rosenthal GG. Mating preferences do not maintain the tailspot polymorphism in the platyfish, Xiphophorus variatus. Behav Ecol. 2013;24(6):1286–91.View ArticleGoogle Scholar
- Kallman KD. The platyfish, Xiphophorus maculatus. Handbook of genetics. Springer; 1975. p. 81-132.Google Scholar
- Basolo AL. Genetic linkage and color polymorphism in the southern platyfish (Xiphophorus maculatus): a model system for studies of color pattern evolution. Zebrafish. 2006;3(1):65–83.PubMedView ArticleGoogle Scholar
- Maan ME, Seehausen O. Ecology, sexual selection and speciation. Ecol Lett. 2013;14(6):591–602.View ArticleGoogle Scholar
- Safran RJ, Scordato ESC, Symes LB, Rodríguez RL, Mendelson TC. Contributions of natural and sexual selection to the evolution of premating reproductive isolation: a research agenda. Trends Ecol Evol. 2013;28(11):643–50.PubMedView ArticleGoogle Scholar
- Chenoweth SF, Appleton NC, Allen SL, Rundle HD. Genomic evidence that sexual selection impedes adaptation to a novel environment. Curr Biol. 2015;25(14):1860–6.PubMedView ArticleGoogle Scholar
- Pearson RG, Raxworthy CJ. The evolution of local endemism in Madagascar: Watershed versus climatic gradient hypothesis evaluated by null biogeographic models. Evolution. 2009;63(4):959–67. doi:https://doi.org/10.1111/j.1558-5646.2008.00596.x.PubMedView ArticleGoogle Scholar
- Van Rensburg BJ, Levin N, Kark S. Spatial congruence between ecotones and range-restricted species: implications for conservation biogeography at the sub-continental scale. Divers Distrib. 2009;15(3):379–89.View ArticleGoogle Scholar
- Moore BR, Donoghue MJ. Correlates of diversification in the plant clade Dipsacales: geographic movement and evolutionary innovations. Am Nat. 2007;170(S2):S28–55.PubMedView ArticleGoogle Scholar
- Agorreta A, Domínguez-Domínguez O, Reina RG, Miranda R, Bermingham E, Doadrio I. Phylogenetic relationships and biogeography of Pseudoxiphophorus (Teleostei: Poeciliidae) based on mitochondrial and nuclear genes. Mol Phylogen Evol. 2013;66(1):80–90.View ArticleGoogle Scholar
- Hulsey CD, Garci’ FJ, Johnson YS, Hendrickson DA, Near TJ. Temporal diversification of Mesoamerican cichlid fishes across a major biogeographic boundary. Mol Phylogen Evol. 2004;31(2):754–64.View ArticleGoogle Scholar
- Peterson AT, Soberón J, Sánchez-Cordero V. Conservatism of ecological niches in evolutionary time. Science. 1999;285(5431):1265–7.PubMedView ArticleGoogle Scholar
- Culumber ZW, Fisher HS, Tobler M, Mateos M, Barber PH, Sorenson MD, et al. Replicated hybrid zones of Xiphophorus swordtails along an elevational gradient. Mol Ecol. 2011;20(2):342–56.PubMedView ArticleGoogle Scholar
- Culumber ZW, Shepard DB, Coleman SW, Rosenthal GG, Tobler M. Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus). J Evol Biol. 2012;25(9):1800–14.PubMedView ArticleGoogle Scholar
- Grossenbacher DL, Veloz SD, Sexton JP. Niche and range size patterns suggest that speciation begins in small, ecologically diverged populations in North American monkeyflowers (Mimulus spp.). Evolution. 2014;68(5):1270–80.PubMedView ArticleGoogle Scholar
- Meyer MK, Schartl M. Xiphophorus mayae, a new species of swordtail from Guatemala (Teleostei: Poeciliidae). Ichthyol Explor Freshwat. 2002;13(1):59–66.Google Scholar
- Kallman KD, Walter RB, Morizot DC, Kazianis S. Two new species of Xiphophorus (Poeciliidae) from the isthmus of tehuantepec, oaxaca, Mexico with a discussion of the distribution of the X. clemenciae clade. Am Mus Novitates. 2004;3441:1–34.View ArticleGoogle Scholar
- Gutierrez-Rodriguez C, Shearer AE, Morris MR, de Queiroz K. Phylogeography and monophyly of the swordtail fish species Xiphophorus birchmanni (Cyprinodontiformes, Poeciliidae). Zool Scr. 2008;37(2):129–39. doi:https://doi.org/10.1111/j.1463-6409.2007.00311.x.View ArticleGoogle Scholar
- Fernandez AA, Morris MR. Mate choice for more melanin as a mechanism to maintain a functional oncogene. Proc Natl Acad Sci. 2008;105(36):13503–7. doi:https://doi.org/10.1073/pnas.0803851105.PubMed CentralPubMedView ArticleGoogle Scholar
- Fernandez AA, Bowser PR. Selection for a dominant oncogene and large male size as a risk factor for melanoma in the Xiphophorus animal model. Mol Ecol. 2010;19(15):3114–23.PubMed CentralPubMedView ArticleGoogle Scholar
- Summers K, Roney KE, Da Silva J, Capraro G, Cuthbertson BJ, Kazianis S, et al. Divergent patterns of selection on the DAB and DXB MHC class II loci in Xiphophorus fishes. Genetica. 2009;135(3):379–90.PubMedView ArticleGoogle Scholar
- Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 2.75. 2011. URL http://mesquiteproject.org. 2011.
- Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20(2):289–90.PubMedView ArticleGoogle Scholar
- Paradis E. Molecular dating of phylogenies by likelihood methods: a comparison of models and a new information criterion. Mol Phylogen Evol. 2013;67(2):436–44.View ArticleGoogle Scholar
- Obregón-Barboza H, Contreras-Balderas S, de Lourdes Lozano-Vilano M. The fishes of northern and central Veracruz, Mexico. Hydrobiologia. 1994;286(2):79–95.View ArticleGoogle Scholar
- West RC. Surface configuration and associated geology of Middle America. In: Wauchop RaW RC, editor. Handbook of Middle American Indians. Austin: University of Texas Press; 1964.Google Scholar
- Mateos M, Sanjur OI, Vrijenhoek RC. Historical biogeography of the livebearing fish genus Poeciliopsis (Poeciliidae: Cyprinodontiformes). Evolution. 2002;56(5):972–84.PubMedView ArticleGoogle Scholar
- Kozak KH, Graham CH, Wiens JJ. Integrating GIS-based environmental data into evolutionary biology. Trends Ecol Evol. 2008;23(3):141–8.PubMedView ArticleGoogle Scholar
- Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25(15):1965–78.View ArticleGoogle Scholar
- Dominguez-Dominguez O, Martinez-Meyer E, Zambrano L, De Leon GPP. Using ecological-niche modeling as a conservation tool for freshwater species: live-bearing fishes in central Mexico. Conserv Biol. 2006;20(6):1730–9.PubMedView ArticleGoogle Scholar
- Chen PF, Wiley EO, McNyset KM. Ecological niche modeling as a predictive tool: silver and bighead carps in North America. Biol Invasions. 2007;9(1):43–51.View ArticleGoogle Scholar
- Costa GC, Schlupp I. Biogeography of the Amazon molly: ecological niche and range limits of an asexual hybrid species. Global Ecol Biogeogr. 2010;19(4):442–51.Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190(3-4):231–59.View ArticleGoogle Scholar
- Pearson RG, Raxworthy CJ, Nakamura M, Peterson AT. Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar. J Biogeogr. 2007;34(1):102–17.View ArticleGoogle Scholar
- Elith J, Phillips SJ, Hastie T, Dudek M, Chee YE, Yates CJ. A statistical explanation of MaxEnt for ecologists. Divers Distrib. 2011;17(1):43–57.View ArticleGoogle Scholar
- Warren DL, Seifert SN. Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria. Ecol Appl. 2011;21(2):335–42.PubMedView ArticleGoogle Scholar
- Swets JA. Measuring the accuracy of diagnostic systems. Science. 1988;240(4857):1285–93.PubMedView ArticleGoogle Scholar
- Liu CR, Berry PM, Dawson TP, Pearson RG. Selecting thresholds of occurrence in the prediction of species distributions. Ecography. 2005;28(3):385–93.View ArticleGoogle Scholar
- Couvreur TLP, Porter-Morgan H, Wieringa JJ, Chatrou LW. Little ecological divergence associated with speciation in two African rain forest tree genera. BMC Evol Biol. 2011;11(1):296.PubMed CentralPubMedView ArticleGoogle Scholar
- Hidalgo-Galiana A, Sánchez-Fernández D, Bilton DT, Cieslak A, Ribera I. Thermal niche evolution and geographical range expansion in a species complex of western Mediterranean diving beetles. BMC Evol Biol. 2014;14(1):187.PubMed CentralPubMedView ArticleGoogle Scholar
- Lira-Noriega A, Peterson AT. Range-wide ecological niche comparisons of parasite, hosts and dispersers in a vector-borne plant parasite system. J Biogeogr. 2014;41(9):1664–73.View ArticleGoogle Scholar
- Revell LJ. Size-correction and principal components for interspecific comparative studies. Evolution. 2009;63(12):3258–68.PubMedView ArticleGoogle Scholar
- Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.View ArticleGoogle Scholar
- Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24(1):129–31.PubMedView ArticleGoogle Scholar
- Harmon LJ, Schulte JA, Larson A, Losos JB. Tempo and mode of evolutionary radiation in iguanian lizards. Science. 2003;301(5635):961–4.PubMedView ArticleGoogle Scholar
- Kozak KH, Larson A, Bonett RM, Harmon LJ. Phylogenetic analysis of ecomorphological divergence, community structure, and diversification rates in dusky salamanders (Plethodontidae: Desmognathus). Evolution. 2005;59(9):2000–16.PubMedView ArticleGoogle Scholar
- Swenson NG. Comparative methods and phylogenetic signal. Functional and Phylogenetic Ecology in R. Springer; 2014. p. 147-71.Google Scholar