Skip to main content

Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae)



One of the most perplexing questions in evolutionary biology is why some lineages diversify into many species, and others do not. In many cases, ecological opportunity has played an important role, leading to diversification along trophic or habitat-based axes. The Goodeidae (Teleostomi: Cyprinodontiformes) are a family of freshwater fishes with two subfamilies: Goodeinae (42 species, viviparous, heterogeneous habitats, Mesa Central of Mexico) and Empetrichthyinae (4 species, oviparous, homogeneous habitats, Great Basin of the United States). These discrepant sets of characteristics and their sister-group relationship make the goodeids amenable to a comparative study of diversification. We gathered lateral body images from more than 1600 specimens of all extant species in the family. Geometric morphometric, and phylogenetic comparative analyses were used to address whether higher species diversity correlates with higher rates of morphological shape evolution and whether there are differences in functional/habitat modules between the two subfamilies.


This study recovered a higher rate of overall body shape evolution in the Goodeinae that is nearly double in magnitude compared to the Empetrichthyinae. A modularity test indicated that the Goodeinae displayed elevated rates of morphological evolution in comparison to the Empetrichthyinae when only trunk (locomotor) regions were compared between subfamilies. No significant differences in evolutionary shape rates were recovered when the trophic (head) regions were compared between subfamilies.


These results support the hypothesis that Mexican goodeids radiated via an ecological opportunity scenario into a wide-array of novel habitats in the island-like Mesa Central as evidenced by their high rate of shape evolution, relative to the Empetrichthyinae. This study quantitatively unraveled the drivers of evolution and eliminated trophic specialization as a driving force within the Goodeidae.


A combination of phylogenetic and morphometric data, and phylogenetic comparative analyses were used to examine body shape rate evolution within the Goodeidae. Results support the hypothesis that species in the subfamily Goodeinae on the central Mexican plateau had a higher rate of body shape evolution relative to its sister subfamily Empetrichthyinae in the Great Basin suggesting that the Goodeinae diversified via an ecological opportunity scenario along habitat, rather than trophic axes.


Disproportionate species richness among clades is one of the most interesting patterns in evolutionary biology, with some clades being exceedingly diverse while others are relatively depauperate [1,2,3]. The drivers of speciation and diversification that result in clade imbalance have been the subject of an ongoing discussion within the literature for decades [4,5,6,7,8,9]. Although species richness might be expected to vary purely by stochastic processes, many other ideas have been put forth to explain discrepant patterns. First, clade age is believed to be an important factor, whereas older clades are expected to have higher species richness due to the greater length of time for diversification and speciation to occur [10, 11]. Second, it has been shown that differential diversification rates are a common explanation for clade disparity [3, 12, 13]. Differential speciation and extinction rates are often correlated with phenomena such as key innovations, adaptive radiations, which can increase diversification rates, and contrasting geologic and climatic histories, which can increase extinction rates [14,15,16,17]. In fact, the majority of the most heavily studied examples of adaptive radiation are related to trophic and/or habitat specializations [18,19,20,21,22].

It is generally expected that species-rich clades should harbor higher levels of phenotypic diversity in comparison to less speciose clades [23, 24]. Morphological disparity and speciation may be linked, and one hypothesis suggests that clades with higher rates of phenotypic evolution may be able to reach into novel ecological trait space, leading to an increase in diversification [24,25,26]. The idea that some organisms are morphologically more versatile than others leading to replacement of the latter is a key concept regarding ecological opportunity and adaptive zones [23, 27, 28].

The freshwater fish family Goodeidae (Order: Cyprinodontiformes) [29] is an ideal group to address the process of diversification. The family is found within the Great Basin of the southwestern United States and the central Mexican highlands (Fig. 1). This disjunct distributional pattern is unique to Goodeidae, and found in no other primary freshwater fish groups. Increasing desiccation of the Sonoran Desert during the Tertiary is hypothesized to have divided the ancestral goodeid lineage into two disjunct subfamilies [30, 31], the Goodeinae [29] and the Empetrichthyinae [32].

Fig. 1
figure 1

Distribution of the two subfamilies of the Goodeidae based on vouchered museum records (, May 2017)

The subfamily Goodeinae is endemic to Mexico, with approximately 18 genera and 42 extant species, with the highest diversity occurring in the geographic area known as the Mesa Central, a relatively depauperate, isolated highland plateau [33,34,35]. This region has been subjected to substantial volcanic and tectonic activity since the beginning of the Miocene, leading to intricate hydrological systems [36, 37], which have likely contributed to higher isolation and speciation in this area [35]. All species in the subfamily are viviparous, and embryos of all species possess a vascular rectal structure, known as a trophotaeniae, for nutrient absorption [38,39,40]. Across the Goodeinae, the species possess trophic ecologies ranging from strictly carnivorous (Alloophorus robustus and Allodontichthys tamazulae) to herbivorous (Goodea atripinnis), however, the majority of species possess similar diets, are omnivorous, and occupy different niches within the trophic spectrum between these two dietary extremes [35, 41, 42]. Species of Goodeinae inhabit lakes, creeks, marshes, canals, and large rivers [35], with some species being habitat specialists (i.e. springs only) and others being more generalist in terms of their habitat preferences.

The sister group to the Goodeinae, the subfamily Empetrichthyinae is much less diverse than the Goodeinae. Two genera, Crenichthys and Empetrichthys, and three species (and multiple subspecies) currently occupy the Great Basin of the United States. Three additional taxa have gone extinct within the last century [42,43,44]. All species of Empetrichthyinae are oviparous, are opportunistic omnivores and utilize similar niches in springs and pools of the Great Basin of the United States [44, 45].

The disparities between the two subfamilies within the Goodeidae offer an excellent opportunity to test for contributions to lineage diversification. In addition to the distinct morphological and life-history differences, one important advantage of the chosen taxa is that they are sister lineages, which therefore allows for the removal of the clade-age effect from the comparison. The objectives of this paper are two-fold. First, to determine whether there are differences in rates of multivariate morphological evolution between the two clades, and whether those differences reflect the hypothesis that the Goodeinae, due to their higher level of diversity and range of morphological variation, reach into novel areas of morphological trait space. Second, to test for differences in rates of evolution between phenotypic modules for the two groups. As stated earlier, trophic specialization and occupancy of new niches are most often associated with elevated speciation rates in other groups, particularly in adaptive radiation scenarios, which has been hypothesized for the Goodeinae on the Mesa Central [31]. In this comparison, it is expected that the head regions of the Goodeinae will vary at higher evolutionary rates if trophic specialization is the selective force (trophic diversification hypothesis). Alternatively, if the trunk region of the Goodeinae display elevated evolutionary rates relative to the Empetrichthyinae, this would indicate that the observed patterns are more likely the result of the occupancy of novel habitats and the associated hydrological constraints during the radiation (habitat diversification hypothesis).


Phylogenetic analysis

The results of the phylogenetic analysis are similar to those of [33] (Fig. 2), but this present study included additional empetrichthyine taxa not included by those authors. The phylogenetic analysis resulted in a sister group relationship between Goodeinae and Empetrichthyinae. A monophyletic clade was inferred that consisted of all the subspecies of Crenichthys baileyi and Crenichthys nevadae, and this clade was sister to Empetrichthys latos, the sole extant species in this genus. Profundulidae is supported as the sister family to Goodeidae. The split of Empetrichthyinae and Goodeinae is estimated at 18.02 Mya (14.3–22.17 Mya). The tribes within Goodeinae diverged around 14.05 Mya (11.62–16.79). The separation of the genera Empetrichthys and Crenichthys is estimated at 6.88 Mya (4.35–9.48 Mya) (Table 1).

Fig. 2
figure 2

Time-calibrated phylogeny generated in BEAST using cytochrome b data from Doadrio and Dominguez [33], and additional taxa of Empetrichthyinae sequenced in this study. Fossil calibrations are shown at each node and referenced in Table 3

Table 1 Major divergence time estimates with the Goodeidae from the BEAST dating analysis

Goodeidae variation in body shape

The tests for overall body shape differences between the two subfamilies revealed significant differences compared to a Brownian motion model of evolution (Table 2). Specifically, the estimated rate of body shape evolution for the Goodeinae was almost twice as high as that for the Empetrichthyinae (σ 2E = 9.84 × 10− 4 vs. σ 2G = 1.90 × 10− 3).

Table 2 Evolutionary shape rate results and Modular Evolutionary Rate Results, with associated p values for significance testing. Test statistics include σ2R for the evolutionary shape rate ratio between the two subfamilies, Rmult is the module shape rate ratio between different modules, and the subscripts "E" and "G" represent Empetrichthyinae and Goodeinae, respectively

Modularity within the Goodeidae

When the morphometric data were further divided into functional/locomotor traits (head and trunk landmarks), only the trunk region showed a significant difference between the subfamilies, (σ2RHead = 1.564 p = 0.150, σ2RTrunk = 2.510 p = 0.021), indicating that these phenotype modules are evolving at different evolutionary rates within the Goodeidae. The head region did not show any significant differences in evolutionary rate between the two subfamilies, suggesting that the Empetrichthyinae and Goodeinae have not diversified along trophic axes.

The results from a comparison of only evolutionary trunk shape rate change, corresponding to the habitat diversification hypothesis, indicate that the taxa with Goodeinae displayed over two and a half times more trunk shape evolution than the Empetrichthyinae (σ 2E = 5.56 × 10− 3 σ 2G = 1.39 × 10− 2).

Body shape diversity within Goodeidae

The goodeids display a wide diversity of overall body shapes. When the extant taxa were plotted in phylomorphospace, the subfamily Goodeinae reaches into novel areas of morphospace for both overall (Fig. 3a) and trunk body shape plots (Fig. 3b). However, there is overlap among spring inhabiting species from both the Empetrichthyinae (Empetrichthys latos, Crenichthys baileyi, and C. nevadae) and Goodeinae (Allodontichthys hubbsi, A. polylepis, and A. tamazulae), which seem to independently have converged on a similar body shape in allopatry. Both PCs (1 and 2) showed the most shape change in the caudal peduncle and in the placement of dorsal and anal fin. PC1 (varianceOverall = 33%, varianceTrunk = 47%) was associated with a widening or compressing of the trunk area, with negative values associated with a widening of the dorsal lateral axis, and positive values are associated with a compression of the area. Positive PC2 (varianceOverall = 18%, varianceTrunk = 21%) values are associated with an elongation of the caudal peduncle and negative PC2 values are associated with a shortening of the caudal peduncle.

Fig. 3
figure 3

Phylomorphospace plots for all goodeid species from the time-calibrated phylogeny. Colors correspond to the seven clades recovered in Fig. 2. In particular, purple colors represent the Empetrichthyinae, whereas other colors represent the other groups of the Goodeinae. Phylomorphospace plots based on; a) tail shape, and b) overall body shape


Disparity of species diversity within and between clades has been heavily studied [1, 2, 10, 46]. As stated earlier, most studies of diversification have focused on trophic specializations [47,48,49] or habitat-based influences [11, 50,51,52,53] as the main drivers of diversification. Disentangling these phenomena (trophic- or habitat-based) has been challenging, but we now have the analytical tools to differentiate between these influences [54, 55].

Ecological opportunity or the evolution of a key innovation are two of the more common triggers of a diversification event in adaptive radiation scenarios [23, 55]. Historically, much of what is known about adaptive radiations was based on a few select and well-studied examples in isolated habitats such as oceanic islands, post-glacial lakes, and other unique and isolated environments [23, 56, 57]. However, the number of adaptive radiation examples has expanded dramatically in recent years, such that adaptive radiations no longer appear to be uncommon phenomena [58, 59].

Until recently, the correlation of morphological diversity and increased disparity, although a hallmark of adaptive radiations [27, 60, 61], has been difficult to test. Net differences in species diversification rates [62,63,64], examination of univariate phenotypic traits [25, 65, 66], and quantification of trophic based morphological differences [67,68,69] between clades have been utilized most often to test for influences on species diversification. These approaches ignore whole body shape or morphological associations as a whole (modules) and cannot provide the clearest picture of diversification. Linking innovations to increased speciation rates represent one of the strongest approaches for studying diversifications and adaptive radiations.

The disparity in species richness between the subfamilies Goodeinae and Empetrichthyinae offered a rare chance to test for contributors to lineage diversification. By utilizing geometric morphometrics, phylogenetic hypotheses, and phylogenetic comparative methods [51, 70], this study addressed two principle questions. First, does higher species diversity within the Goodeinae correlate with higher rates of morphological shape evolution in comparison to the Empetrichthyinae. The results from this study show a higher rate of overall body shape evolution in the Goodeinae that is almost double in magnitude than in the Empetrichthyinae. These results support the hypothesis that Mexican goodeids radiated in the Mesa Central, in a short time frame, as evidenced by their high rate of shape evolution, relative to the Empetrichthyinae. Additionally, the phylomorphospace plots show that many members of the Goodeinae reach into novel areas of morphological trait space and are often not clustered phylogenetically, whereas, the Empetrichthyinae did not expand into novel trait space and did not diversify to the same extent as the Goodeinae. Alternatively, the evolution of viviparity also could have contributed to the adaptation of the Goodeinae in these novel environments, followed by morphological diversification. The viviparous life-history of the group could have freed them being linked to particular substrate or structure for flow regime for spawning and egg attachment, thereby allowing them to diversify in this region.

The Mesa Central, the center of Goodeinae diversification, is an otherwise faunally depauperate, isolated highland plateau that has experienced hydrological compartmentalization events due to volcanism and orogeny during the Pliocene and Pleistocene [71]. One other group of fishes has diversified in the Mesa Central (Chirostoma: Atherinopsidae) via an adaptive radiation [71,72,73]. The Mesa Central, and the availability of an unoccupied array of heterogeneous habitats resulted in species of Goodeinae being were able to inhabit varying hydrological regimes from slow to fast flows throughout the Rio Lerma basin [34]. Although this has been speculated as contributing to Goodeinae diversification [35, 74], this idea previously has not been tested, as other studies have focused on the patterns of diversification within the Goodeinae [37] or for select groups within the subfamily [75].

The second objective of this study specifically addressed whether there are differences in trophic/locomotor modules between the two subfamilies, thereby testing whether the differences in disparity between goodeid lineages is the result of trophic or habitat specialization. Phenotypic trait modules can give insights into the developmental, functional, and locomotory influences of morphological disparity. Contrasting what is generally expected, the modularity rate tests within the Goodeidae (Goodeinae+ Empetrichthyinae) recovered a much higher rate of trunk region evolution than the head region. In general, goodeid species in both subfamilies are generalists in terms of feeding niches, and the majority of the species seem to be opportunistic omnivores [34]. This could explain why head shape showed no significant difference in evolutionary shape rates between the two clades of goodeid fishes. In general, diversification along trophic axes is a common phenomenon, especially for fishes. For example, east African Rift Lake cichlids have long served as a model in our understanding of diversification and adaptive radiation [76]. The large range in both species and phenotypic diversity within cichlids has been hypothesized to be a result of utilization of different dietary components [77,78,79,80]. In several cichlid groups, morphological structures have been linked to trophic diversity, such as the pharyngeal jaw and certain muscular elements of the jaw [81, 82], and are believed to allow for increased diversification because of trophic partitioning [83]. Studies have investigated diversification shifts within the group [77, 84, 85], or have examined the correlation between morphology and ecomorphological variables [68] to explain species disparity. Similarly, pupfish radiations also have been linked to differences in head morphology by assigning functional variables a priori and comparing rates of diversification [26].

When evolutionary shape rates of the trunk region were compared between the subfamilies, the Goodeinae showed a relative rate increase of over two-fold, in comparison to the same trunk region of the Empetrichthyinae using the approach outlined in [58]. Goodeinae taxa have diversified into various aquatic systems (canals, streams, rivers, lakes, ditches, and outflows) [34], whereas the Empetrichthyinae are confined to homogeneous habitats including small springs and pools [86]. The higher rate of trunk region evolution in Goodeinae may in fact show a constraint release of the phenotypic liability of the trunk region due to the Goodeinae experiencing different hydrodynamic pressure than the Empetrichthyinae.

Similar to this study, other studies have recovered a strong link between hydrodynamic flow and body shape [87, 88, 89,90,91,92]. Therefore, for fish species occupying a wide range of different hydrodynamic regimes (i.e. lakes, rivers, streams, springs), such as the Goodeinae, the variation in overall body shape, and more specifically, the trunk region of the body is expected to increase relative to close relatives occupying more homogeneous habitats. The modularity data supports the idea that such variation may be the result of the functional constraints of life in new hydrologic regimes (ecological release, sensu 71), as a result of ecological opportunity on the Mesa Central, thereby allowing for an increase in ecological niche space. All species within the Empetrichthyinae occupy spring and pool habitats and show little variation in trunk morphology and a low rate of trunk shape evolution.


The results from this study support the hypothesis that the Goodeinae radiated via an ecological opportunity scenario (sensu 74), along habitat axes. The Goodeinae occupy the faunally depauperate, island-like Mesa Central. This is relevant because others have suggested that adaptive radiations primarily occur on islands or island–like habitats [75], although there are examples of radiations occurring under non-island like conditions [93]. The higher rate of trunk shape evolution for the Goodeinae supports the hypothesis that the Goodeinae ecological release following colonization of a diverse array of novel, unoccupied habitats during the Plio-Pleistocene facilitated directional selection on few key survival traits. It is likely that the diverse array of aquatic habitats and niches were either already occupied by competitors or not available for the Empetrichthyinae during the formation of the Great Basin, therefore both the species diversification and morphological diversification rates of the Empetrichthyinae were lower in comparison to the Goodeinae. Finally, the dated phylogeny suggests that there was rapid speciation and morphological diversification within the Goodeinae, lending additional support to the ecological opportunity hypothesis. Unlike most other studies of adaptive radiation, however, this study was able to disentangle the drivers of evolution and quantitatively eliminate trophic specialization as a driving force within the Goodeidae.


Phylogenetic analysis

Most mitochondrial DNA sequences of cytochrome b (cytb) were obtained from Genbank (~ 1140 bp, for the Goodeinae [33] and outgroup taxa (Additional file 1: Table S1). Additional sequences were generated in this study and tissue samples for these species were obtained from colleagues, natural history collections, or natural resource agencies. Sequence data were pruned to include one representative individual for each species. For the subfamily Empetrichthyinae, tissue samples were obtained, and DNA was extracted from fin clips using the DNeasy Blood and Tissue kit (Qiagen) following the manufacturer’s instructions. The cytb gene was amplified via polymerase chain reaction (PCR) with the primers L14724 and H1595 from [94].

Each 25 μl cytb PCR reaction consisted of the following: 0.75 μl MgCl2, 2.5 μl 10× buffer, 0.5 μl NTP, 0.5 μl of each 10 μM primer, 0.25 μl Taq, 1 μl of DNA template, and 19–20 μl of water. The thermal cycling parameters for cytb were as follows: 94 °C for 2 min, 27 cycles of 94 °C for 45 s, 54 °C for 30s, 72 °C for 1 min, and a final extension of 74 °C for 10 min. DNA sequencing was performed by the Beckman Coulter Genomics Facility (Danvers, MA). Alignment and editing of sequences was performed using Geneious 9.1.4 [95] Sequence data were submitted to GenBank (Additional file 1: Table S1).

Beast v1.8 [96] was used to generate a time-calibrated phylogeny, using the General Time Reversible (GTR) model, a relaxed lognormal clock, and a tree prior of Yule Process [97]. Divergence time estimations were conducted using eight fossil calibrations, lognormal priors, and minimum time estimates based on the oldest known fossils of select ingroup and outgroup taxa (Table 3). Two separate runs consisting of 40 million generations each were run sampling every 1000 generations. Tracer V1.5 [98] was used to determine effective sampling size (ESS) of all parameters and to determine that the analysis reached stationary. All ESS parameters were ˃200, indicating that the parameters robustly sampled the posterior distribution. A maximum clade credibility tree using a 20% burn-in was created with the program Tree Annotator V1.8.3 [96] and was used as the backbone for phenotypic diversification analyses.

Table 3 Fossil calibrations and parameters used in BEAST to construct a dated phylogeny

Geometric Morphometrics

Images of museum specimens (N = 1669) of all species (40 species and 4 subspecies) were taken from the left lateral side using a Nikon D5000 SLR camera and a 50 mm macro lens. Museum abbreviations follow [99]. Bent, warped, or juvenile specimens were not utilized to limit ontogenetic effects and non-biological variation. Eighteen homologous, two dimensional landmarks were used to summarize body shape [100] using the program TPSdig version 2.12 [101]. The landmarks (Fig. 4a) included; 1) anterior tip of the snout, 2) opening of mouth, 3) posterior edge of jaw, 4) posterior edge of the neurocranium, 5) upper edge of eye, 6) posterior edge of eye, 7) ventral edge of eye, 8) anterior edge of eye, 9) anterior edge of the dorsal fin, 10) posterior edge of the dorsal fin, 11) dorsal edge of the caudal fin, 12) ventral edge of the caudal fin, 13) anterior insertion of the anal fin, 14) posterior insertion of anal fin, 15) intersection of gill opening and ventral margin of body, 16) posterior most edge of neurocranium, 17) upper insertion of the pectoral fin, and 18) lower insertion of the pectoral fin.

Fig. 4
figure 4

a Distribution of landmarks used for geometric morphometric analyses and evolutionary shape rate tests, b modularity rate test comparing head and trunk morphologies of the Goodeinae [G] versus the Empetrichythinae [E], and c) modularity rate test by body region testing the trophic diversification hypothesis versus habitat diversification hypothesis by subfamily

Rate of body shape evolution analyses

We used the multivariate rates methodology developed by [102] to examine rates of body shape evolution using all of the landmarks. This approach utilizes geometric morphometric data and a phylogeny to test for differences in evolutionary shape rates (σ2). This method uses multivariate distances and a Brownian motion model of evolution that assumes that variance increases with time. The distance-based approach has also been shown to be statistically robust when increasing trait dimensions, which is unavoidable in multivariate data [103].

A Procrustes superimposition analysis was carried out in the program MorphoJ [104], which corrects for scaling, rotation, translation, and size biases. All other analyses were performed in R [105] using the packages ‘geomorph’, ‘phytools’, and ‘geiger’. Mean shape data from the Procrustes data were calculated for each species for the shape analyses. The time-calibrated phylogeny was pruned to include only the taxa for which morphological data were available using the ‘drop.tip’ function in the package ‘phytools’. The multivariate technique for testing evolutionary rates (σ2) using the shape data collected and the time-calibrated phylogeny was implemented [102]. The evolutionary shape rate for each subfamily was recovered (σ2G, σ2E), along with a ratio for Goodeidae and Empetrichthyinae (σ2R), which represents relative differences in evolutionary shape rate between the two subfamilies. Phylogenetic simulation of 9999 iterations was used to test the null hypothesis that one single rate is found between the two subfamilies.

We used the procedures in [106] and [102] to test for multivariate differences between sets of phenotypic traits (modules). Modules are sets of traits that show highly mingled connections to each other, but are loosely connected to other sets of traits [107, 108]. Modules have long been thought to vary in evolutionary rates; for example, the steady swimming hypothesis predicts caudal peduncle ratio differences among individuals in fast vs. slow moving aquatic systems [87, 88].

We conducted a modularity test to specifically test for differences in trophic/locomotor modules (head vs. trunk regions), corresponding to the trophic and habitat diversification hypotheses (Fig. 4b). The head region consisted of landmarks 1–8, 15–18, whereas the trunk region included landmarks 9–14 (Fig. 4c). This test was accomplished using a two-step approach incorporating the Procrustes aligned shape data considered by module, along with a time-calibrated phylogeny under a Brownian model of evolution. First, we compared evolutionary shape rates between the two phenotypic modules (head vs. trunk) (Fig. 4b) to determine whether the modules are evolving at different rates, and therefore represent independent anatomical units. Second, we then compared evolutionary shape rates within each module by subfamily to determine if there were rate differences at the subfamily level within each module. The procedures follow [106] and were carried out in the package ‘geomorph’. Landmarks were divided by morphological regions, and then tested using relative rate ratios.

Finally, we used the phylomorphospace approach [109, 110] to characterize the evolutionary patterns of body shape diversity within the Goodeidae (overall, and within phenotypic modules). This multivariate phenotypic method projects the phylogeny and the first two principal components of a principal components analysis of the aligned morphometric data to depict the evolutionary history of morphospace occupancy for the goodeid fishes.


  1. 1.

    Slowinski JB, Guyer C. Testing the stochasticity of patterns of organismal diversity: an improved null model. Amer Nat. 1989;134:907–21.

    Article  Google Scholar 

  2. 2.

    Heard SB. Patterns in phylogenetic tree balance with variable and evolving speciation rates. Evolution. 1996;50:2141–8.

    PubMed  Article  Google Scholar 

  3. 3.

    Ricklefs RE. Estimating diversification rates from phylogenetic information. Trends Ecol Evolut. 2007;22:601–10.

    Article  Google Scholar 

  4. 4.

    Thorson G. Bottom communities. In: Hedgpeth JW, editor. Treatise on marine ecology and paleoecology Geological Society of America memoirs, boulder, CO; 1957. p. 461–534.

    Google Scholar 

  5. 5.

    Fischer AG. Latitudinal variations in organic diversity. Evolution. 1992;14:64–81.

    Article  Google Scholar 

  6. 6.

    Hutchinson GE. Homage to Santa Rosalia or why are there so many kinds of animals? Amer Nat. 1959;870:145–59.

    Article  Google Scholar 

  7. 7.

    Foote M. Contributions of individual taxa to overall morphological disparity. Paleobiology. 1993;19:403–19.

    Article  Google Scholar 

  8. 8.

    Losos JB, Miles DB. Testing the hypothesis that a clade has adaptively radiated: iguanid lizard clades as a case study. Amer Nat. 2002;160:147–57.

    Article  Google Scholar 

  9. 9.

    Burbrink FT, Chen X, Myers EA, Brandeley MC, Pyron AR. Evidence for determinism in species diversification and contingency in phenotypic evolution during adaptive radiation. Proc R Soc Lond. 2012;279:4817–26.

    Article  Google Scholar 

  10. 10.

    McPeek MA, Brown JM. Clade age and not diversification rate explains species richness among animal taxa. Amer Nat. 2007;169:97–106.

    Article  Google Scholar 

  11. 11.

    Bloom DD, Weir JT, Piller KR, Lovejoy NR. Do freshwater fishes diversify faster than marine fishes? A test using state-dependent diversification analyses and molecular phylogenetics of new world silversides Atherinopsidae. Evolution. 2013;67:2040–57.

    PubMed  Article  Google Scholar 

  12. 12.

    Cook JA, Lessa EP. Are rates of diversification in subterranean south American tuco-tucos genus Ctenomys, Rodentia: Octodontidae unusually high? Evolution. 1998;52:1521–7.

    PubMed  Google Scholar 

  13. 13.

    Magallon S, Sanderson MJ. Absolute diversification rates in angiosperm clades. Evolution. 2001;55:1762–80.

    PubMed  Article  CAS  Google Scholar 

  14. 14.

    Erwin DH. A preliminary classification of evolutionary radiations. Hist Biol. 1992;2:133–47.

    Article  Google Scholar 

  15. 15.

    Foote M. Discordance and concordance between morphological and taxonomic diversity. Paleobiology. 1993;19:185–204.

    Article  Google Scholar 

  16. 16.

    Heard SB, Hauser DL. Key evolutionary innovations and their ecological mechanisms. Hist Biol. 1995;10:151–73.

    Article  Google Scholar 

  17. 17.

    Hunter JP. Key innovations and the ecology of macroevolution. Trends Ecol Evol. 1998;13:31–6.

    PubMed  Article  CAS  Google Scholar 

  18. 18.

    Greenwood PJ. Morphology, endemism and speciation in African cichlid fishes. Verhandlungen der Deutschen Zoologischen Gesellschaft. 1973;66:115–24.

    Google Scholar 

  19. 19.

    Losos JB. Ecomorphology, performance capability, and scaling of west Indian Anolis lizards: an evolutionary analysis. Ecol Monogr. 1990;60:369–88.

    Article  Google Scholar 

  20. 20.

    Chakrabarty P. Testing conjectures about morphological diversity in cichlids of lakes Malawi and Tanganyika. Copeia. 2005;2:359–73.

    Article  Google Scholar 

  21. 21.

    López-Fernández H, Arbour JH, Winemiller KO, Honeycutt HL. Testing for ancient adaptive radiations in Neotropical cichlid fishes. Evolution. 2013;67:1321–37.

    PubMed  Google Scholar 

  22. 22.

    Martin CH, Feinstein LC. Novel trophic niches drive variable progress toward ecological speciation within an adaptive radiation of pupfishes. Mol Ecol. 2014;23:1846–62.

    PubMed  Article  Google Scholar 

  23. 23.

    Schluter D. The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.

    Google Scholar 

  24. 24.

    Slater GJ, Price SA, Francesco S, Alfaro ME. Diversity versus disparity and the radiation of modern cetaceans. Proc R Soc Lond. 2010;1697:3097–104.

    Article  Google Scholar 

  25. 25.

    Parent CE, Crespi BJ. Ecological opportunity in adaptive radiation of Galapagos endemic land snails. Amer Nat. 2009;174:898–905.

    Article  Google Scholar 

  26. 26.

    Martin CH, Wainwright PC. Trophic novelty is linked to exceptional rates of morphological diversification in two adaptive radiations of Cyprinodon pupfish. Evolution. 2011;65:2197–227.

    PubMed  Article  Google Scholar 

  27. 27.

    Simpson GG. Tempo and mode in evolution Columbia University press; 1944.

    Google Scholar 

  28. 28.

    Vermeij GJ. Biological versatility and earth history. Proc Natl Acad Sci U S A. 1973;70:1936–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  29. 29.

    Jordan DS. A classification of fishes, including families and genera as far as known. Stanford: Stanford University Press; 1923.

    Google Scholar 

  30. 30.

    Parenti LR. A phylogenetic and biogeographic analysis of cyprinodontiform fishes Teleostei, Atherinomorpha. Bull Am Mus Nat Hist. 1981;168:341–547.

    Google Scholar 

  31. 31.

    Webb SA, Graves JA, Macias-Garcia C, Magurran AE, Foighil DÓ, Ritchie MG. Molecular phylogeny of the livebearing Goodeidae:Cyprinodontiformes. Mol Phylogenet Evol. 2004;30:527–44.

    PubMed  Article  CAS  Google Scholar 

  32. 32.

    Jordan DS, Evermann BW. The fishes of north and middle America. Bulletin of the United States National Museum. 1930;47:2183–3136.

    Google Scholar 

  33. 33.

    Doadrio I, Domínguez O. Phylogenetic relationships within the fish family Goodeidae based on cytochrome b sequence data. Mol Phylogenet Evol. 2004;31:416–30.

    PubMed  Article  CAS  Google Scholar 

  34. 34.

    Miller RR, Minckley WL, Norris SM. Freshwater fishes of Mexico. Chicago: The University of Chicago Press; 2005.

    Google Scholar 

  35. 35.

    Domínguez-Domínguez O, Pedraza-Lara C, Gurrola-Sánchez N, Pérez-Rodríguez R, Israde-Alcántara I, Garduño-Monroy VH, Brooks DR. Historical biogeography of the Goodeinae: Cyprinodontiformes. In: UribeMC GHJ, editor. Viviparous fishes II. New life publications; 2010. p. 33–74.

    Google Scholar 

  36. 36.

    Nieto-Samaniego AF, Ferrari L, Alaniz-Alvarez SA, Labarthe-Hernández G, Rosas-Elguera J. Variation of Cenozoic extension and volcanism across the southern sierra Madre occidental volcanic province. Geol Soc Am Bull. 1999;111:347–63.

    Article  Google Scholar 

  37. 37.

    Dominguez-Dominguez O, Doadrio I, Perez-Ponce de Leon G. Historical biogeography of some river basins in Central Mexico evidenced by their goodeine freshwater fishes: a preliminary hypothesis using secondary brooks parsimony analysis. J Biogeogr. 2006;33:1437–47.

    Article  Google Scholar 

  38. 38.

    Hubbs CL, Turner CL. Studies of the fishes of the order Cyprinodontes, XVI: a revision of the Goodeidae. Misc Publ Museum Zool. 1939;42:1–100.

    Google Scholar 

  39. 39.

    Wourms JP, Cohen DM. Trophotaeniae, embryonic adaptations, in the viviparous ophidioid fish, Oligopus longhursti: a study of museum specimens. J Morphol. 1975;147:385–401.

    Article  PubMed  Google Scholar 

  40. 40.

    Lombardi J, Wourms JP. Embryonic growth and trophotaenial development in goodeid fishes (Teleostei: Atheriniformes). J Morphol. 1988;197:193–208.

    PubMed  Article  Google Scholar 

  41. 41.

    Soto-Galera E, Paulo-Maya J, López-López E, Serna-Hernández JA, Lyons J. Change in fish fauna as indication of aquatic ecosystem condition in Río Grande de Morelia-Lago de Cuitzeo basin. Environ Manag. 1999;1:133–40.

    Google Scholar 

  42. 42.

    Lyons J, Gutierrez-Hernandez A, Diaz-Pardo E, Soto-Galera M, Medina-Nava M, Pineda-Lopez R. Development of a preliminary index of biotic integrity IBI based on fish assemblages to assess ecosystem condition in lakes in Mexico. Hydrobiologia. 2000;418:57–72. Minckley WL, Deacon JE. Southwestern fishes and the enigma of “endangered species”. Science. 1968; 159:1424-32

    Article  Google Scholar 

  43. 43.

    Grant EC, Riddle BR. Are the endangered springfish Crenichthys Hubbs and Poolfish Empetrichthys Gilbert fundulines or goodeids? A mitochondrial DNA assessment. Copeia 1995:209–12.

  44. 44.

    Williams JE. Threatened fishes of the world: Empetrichthys latos miller, 1948 (Cyprinodontidae). Environ Biol Fish. 1996;45:272.

    Article  Google Scholar 

  45. 45.

    Williams CD, Williams JE. Summer food habits of fishes from two springs in east-Central Nevada. Southwest Nat. 1982;27:437–45.

    Article  Google Scholar 

  46. 46.

    Stadler T. Mammalian phylogeny reveals recent diversification rate shifts. ProNat Acad Sci USA. 2011;108:6187–92.

    Article  Google Scholar 

  47. 47.

    Collar DC, O'Meara BC, Wainwright PC, Near TJ. Piscivory limits diversification of feeding morphology in centrarchid fishes. Evolution. 2009;63:1557–73.

    PubMed  Article  Google Scholar 

  48. 48.

    Thomas GH, Meiri S, Phillimore AB. Body size diversification in Anolis: novel environment and island effects. Evolution. 2009;63:2017–30.

    PubMed  Article  CAS  Google Scholar 

  49. 49.

    Strauss R. Allometry and functional feeding morphology in haplochromine cichlids. In: Echelle AA, Kornfield I, editors. Evolution of fish species flocks. University of Maine at Orono press; 1984. p. 217–29.

    Google Scholar 

  50. 50.

    Rüber L, Adams DC. Evolutionary convergence of body shape and trophic morphology in cichlids from Lake Tanganyika. J Evol Biol. 2001;14:325–32.

    Article  Google Scholar 

  51. 51.

    Bouton N, De Visser J, Barel CDN. Correlating head shape and ecological variables in rock-dwelling haplochromines Teleostei: Cichlidae from Lake Victoria. Biol J Linnean Soc. 2002;76:39–48.

    Article  Google Scholar 

  52. 52.

    Sabaj MH. Standard symbolic codes for institutional resource collections in herpetology and ichthyology: an Online Reference. Version 6.5. 2016. Electronically accessible at, American Society of Ichthyologists and Herpetologists, Washington, DC.

  53. 53.

    Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Geometric Morphometrics for Biologists: A Primer. New York: Elsevier Academic Press; 2004.

    Google Scholar 

  54. 54.

    Barbour CD. A biogeographical history of Chirostoma (Pisces: Atherinidae): a species flock from the Mexican plateau. Copeia. 1973;1973:533–66.

    Article  Google Scholar 

  55. 55.

    Echelle AA, Echelle AF. Evolutionary genetics of a ‘species flock:‘Atherinid fishes on the mesa central of México. In: Echelle AA, Kornfield I, editors. Evolution of fish species flocks. University of Maine Press; 1984. p. 94–109.

    Google Scholar 

  56. 56.

    Bloom DD, Piller KR, Lyons J, Mercado-Silva N, Medina-Nava M. Systematics and biogeography of the silverside tribe Menidiini based on the mitochondrial ND2 gene. Copeia. 2009;2009:408–17.

    Article  Google Scholar 

  57. 57.

    Ritchie MG, Hamill RM, Graves JA, Magurran AE, Webb SA, Macías Garcia C. Sex and differentiation: population genetic divergence and sexual dimorphism in Mexican goodeid fish. J Evol Biol. 2007;20:2048–55.

    PubMed  Article  CAS  Google Scholar 

  58. 58.

    Webb SA. Molecular systematics of the genus Allodontichthys (Cyprinodontiformes: Goodeidae). Rev Fish Biol Fish. 2002;12:193–205.

    Article  Google Scholar 

  59. 59.

    Muschick M, Indermaura A, Salzburger W. Convergent evolution within an adaptive radiation of cichlid fishes. Curr Biol. 2012;22:2362–8.

    PubMed  Article  CAS  Google Scholar 

  60. 60.

    Day JJ, Cotton JA, Barraclough TG. Tempo and mode of diversification of Lake Tanganyika cichlid fishes. PLoS One. 2008;

  61. 61.

    Kocher T. Adaptive evolution and explosive speciation: the cichlid fish model. Nature Rev Genet. 2004;5:288–98.

    PubMed  Article  CAS  Google Scholar 

  62. 62.

    Clabaut C, Bunje PM, Salzburger W, Meyer A. Geometric morphometric analyses provide evidence for the adaptive character of the Tanganyikan cichlid fish radiations. Evolution. 2007;61:560–5.

    PubMed  Article  Google Scholar 

  63. 63.

    Liem KF. Evolutionary strategies and morphological innovations: cichlid pharyngeal jaws. Syst Zool. 1974;20:425–41.

    Google Scholar 

  64. 64.

    Liem KF. Aquatic versus terrestrial feeding modes: possible impacts on the trophic ecology of vertebrates. Amer Zool. 1990;30:209–22.

    Article  Google Scholar 

  65. 65.

    Hulsey CD, Hollingsworth PR, Fordyce JA. Temporal diversification of central American cichlids. BMC Evol Biol. 2010;10:1–12.

    Article  CAS  Google Scholar 

  66. 66.

    Robinson BW, Wilson DS. Character release and displacement in fishes: a neglected literature. Amer Nat. 1994;1:596–627.

    Article  Google Scholar 

  67. 67.

    Hulsey CD, Mims MC, Parnell NF, Streelman JT. Comparative rates of lower jaw diversification in cichlid adaptive radiations. J Evol Biol. 2010;23:1456–67.

    PubMed  Article  CAS  Google Scholar 

  68. 68.

    McMahan CD, Chakrabarty P, Sparks JS, Smith WML, Davis MP. Temporal patterns of diversification across global cichlid biodiversity (Acanthomorpha: Cichlidae). PLoS One. 2013;

  69. 69.

    Denton JSS, Adams DC. A new phylogenetic test for comparing multiple high-dimensional evolutionary rates suggests interplay of evolutionary rates and modularity in lanternfishes (Myctophiformes; Myctophidae). Evolution. 2015;69:2425–40.

    PubMed  Article  Google Scholar 

  70. 70.

    Minckley WL, Marsh PC. Inland fishes of the greater southwest: Chronicle of a Vanishing Biota. Tucson: University of Arizona Press; 2009.

    Google Scholar 

  71. 71.

    Langerhans RB. Predictability of phenotypic differentiation across flow regimes in fishes. Integr Comp Biol. 2008;48:750–68.

    PubMed  Article  Google Scholar 

  72. 72.

    Foster K, Bower LM, Piller KR. Getting in shape: habitat-based morphological divergence for two sympatric fishes. Biol J Linnean Soc. 2015;114:152–62.

    Article  Google Scholar 

  73. 73.

    Webb PW. Fast-start performance and body form in seven species of teleost fish. J Exp Biol. 1978;4:211–26.

    Google Scholar 

  74. 74.

    Gatz AJ Jr. Ecological morphology of freshwater stream fishes. Tulane Stud Zool Bot. 1979;21:91–124.

    Google Scholar 

  75. 75.

    Gatz AJ Jr. Community organization in fishes as indicated by morphological features. Ecology. 1979:711–8.

  76. 76.

    Triantafyllou MS, Triantafyllou GS, Yue DKP. Hydrodynamics of fishlike swimming? Annu Rev Fluid Mech. 2000;32:33–53.

    Article  Google Scholar 

  77. 77.

    Tanentzap AJ, Brandt AJ, Smissen RD, Heenan PB, Fukami T, Lee WG. When do plant radiations influence community assembly? The importance of historical contingency in the race for niche space. New Phytol. 2015;207:468–79.

    PubMed  Article  Google Scholar 

  78. 78.

    Schmidt TR, Gold JR. Complete sequence of the mitochondrial cytochrome b gene in the Cherryfin shiner, Lythrurus roseipinnis (Teleostei: Cyprinidae). Copeia. 1993;3:880–3.

    Article  Google Scholar 

  79. 79.

    Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.

    PubMed  PubMed Central  Article  Google Scholar 

  80. 80.

    Drummond AJ, et al. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  81. 81.

    Kingman JFC. The coalescent. Stoch Process Appl. 1982;13:235–48.

    Article  Google Scholar 

  82. 82.

    Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Tracer v1.7, 2018.

  83. 83.

    Rohlf, FJ. tpsDig, digitize landmarks and outlines, 2.05: Department of Ecology and Evolution, State University of New York at Stony Brook. 2006.

  84. 84.

    Adams DC, James Rohlf F, Slice DE. Geometric morphometrics: ten years of progress following the ‘revolution’. Ital J Zool. 2004;71:5–16.

    Article  Google Scholar 

  85. 85.

    McArdle BH, Anderson MJ. Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology. 2001;82:290–7.

    Article  Google Scholar 

  86. 86.

    Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol. 2011;11:353–7.

    Article  Google Scholar 

  87. 87.

    Salzburger, W, Mack, T, Verheyen, E, Meyer, A. Out of Tanganyika: Genesis, explosive speciation, key-innovations and phylogeography of the haplochromine cichlid fishes. BMC Evolutionary Biology; 2005:17,

  88. 88.

    Simpson GG. The major features of evolution. New York: Columbia University Press; 1953.

    Google Scholar 

  89. 89.

    R Development Core Team, R: A language and environment for statistical computing. R D C team, Ed R foundation for statistical computing. 2011.

    Google Scholar 

  90. 90.

    Klingenberg CP. Morphological integration and developmental modularity. Annu Rev Ecol Syst. 2008;39:115–32.

    Article  Google Scholar 

  91. 91.

    Clune J, Mouret J, Lipson H. The evolutionary origins of modularity. Proc R Soc Lond B Biol Sci. 2013;280:2812–63.

    Article  Google Scholar 

  92. 92.

    Sidlauskas B. Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008;62:3135–56.

    PubMed  Article  Google Scholar 

  93. 93.

    Klingenberg CP, Gidaszewski NA. Testing and quantifying phylogenetic signals and homoplasy in morphometric data. Syst Biol. 2010;59:245–61.

    PubMed  Article  CAS  Google Scholar 

  94. 94.

    Rabosky DL. Ecological limits and diversification rate: alternative paradigms to explain the variation in species richness among clades and regions. Ecol Lett. 2009;12:735–43.

    PubMed  Article  Google Scholar 

  95. 95.

    Partridge L. Field and laboratory observations on the foraging and feeding techniques of blue tits Parus caeruleus and coal tits P. ater in relation to their habitats. Anim Behav. 1976;24:534–44.

    Article  Google Scholar 

  96. 96.

    Moreno E, Carrascal LM. Leg morphology and feeding postures in four Parus species: an experimental ecomorphological approach. Ecology. 1993;4:2037–44.

    Article  Google Scholar 

  97. 97.

    Bouton N, Witte F, Van Alphen J, Schenk A, Seehausen O. Local adaptations in populations of rock-dwelling haplochromines (Pisces: Cichlidae) from southern Lake Victoria. Proc R Soc Lond B. 1999;1417:335–60.

    Google Scholar 

  98. 98.

    Carr GD, Robichaux RH, Witter MS, Kyhos DW. Adaptive radiation of the Hawaiian silversword alliance (Compositae–Madiinae): a comparison with Hawaiian picture-winged Drosophila. In: Giddings LV, Kaneshiro KY, Anderson WW, editors. Genetics, speciation and the founder principle. Oxford: Oxford Univ Press; 1989. p. 79–97.

    Google Scholar 

  99. 99.

    Armbruster WS. Evolution of plant pollination systems: hypotheses and tests with the Neotropical vine Dalechampia. Evolution. 1993;12:1480–505.

    Article  Google Scholar 

  100. 100.

    Harrison MK, Crespi BJ. A phylogenetic test of ecomorphological adaptation in Cancer crabs. Evolution. 1999;99:961–5.

    Article  Google Scholar 

  101. 101.

    Stromberg CA. Decoupled taxonomic radiation and ecological expansion of open-habitat grasses in the Cenozoic of North America. Proc Natl Acad Sci. 2005;102:11980–4.

    PubMed  Article  CAS  Google Scholar 

  102. 102.

    Adams DC. Quantifying and comparing phylogenetic evolutionary rates for shape and other high-dimensional phenotypic data. Syst Biol. 2014;63:166–78.

    PubMed  Article  Google Scholar 

  103. 103.

    Harmon LJ, Weir J, Brock C, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–31.

    PubMed  Article  CAS  Google Scholar 

  104. 104.

    Yoder JB, Clancey E, Roches DS, Eastman JM, Gentry L, Godsoe W, Hagey TJ, Jochimsen D, Oswald BP, Robertson J, et al. Ecological opportunity and the origin of adaptive radiations. J Evol Biol. 2010;23:1581–96.

    PubMed  Article  CAS  Google Scholar 

  105. 105.

    Grant PR, Grant BR. How and why species multiply the radiation of Darwin’s finches. Princeton: University Press; 2008.

    Google Scholar 

  106. 106.

    Losos JB. Lizards in an evolutionary tree: ecology and adaptive radiation of anoles. Berkeley: University of California Press; 2009.

    Google Scholar 

  107. 107.

    Stroud JT, Losos JB. Ecological opportunity and adaptive radiation. Annu Rev Ecol Evol Syst. 2016;47:507–32.

    Article  Google Scholar 

  108. 108.

    Sturmbauer C. Explosive speciation in cichlid fishes of the African Great Lakes: a dynamic model of adaptive radiation. J Fish Biol. 1998;53:18–36.

    Article  Google Scholar 

  109. 109.

    Schluter D. Ecology and the origin of species. Trends Ecol Evol. 2001;16:372–80.

    PubMed  Article  CAS  Google Scholar 

  110. 110.

    Rabosky DL. Likelihood methods for detecting temporal shifts in diversification rates. Evolution. 2006;60:1152–64.

    PubMed  Article  Google Scholar 

  111. 111.

    Marchio EA, Piller KR. Cryptic diversity in a widespread live-bearing fish (Poeciliidae: Belonesox). Biol J Linn Soc. 2013;109:848–60.

    Article  Google Scholar 

  112. 112.

    Santini F, Harmon LJ, Carnevale G, Alfaro ME. Did genome duplication drive the origin of teleosts? A comparative study of diversification in ray-finned fishes. BMC Evol Biol. 2009;9:194.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  113. 113.

    Bickley, W.B Jr. Paleoenvironmental reconstruction of late quaternary lacustrine sediments (Seibold site) in southeastern North Dakota. 1970. M.S. thesis, University of North Dakota, Grand Forks.

  114. 114.

    Cvancara AM, Clayton L, Bickley WB Jr, Jacob AF, Ashworth AC, Brophy JA, Shay T, Delorme LD, Lammers GE. Paleolimnology of late quaternary deposits: Seibold site, North Dakota. Science. 1971;171:172–4.

    PubMed  Article  CAS  Google Scholar 

  115. 115.

    Uyeno T, Miller RR. Relationships of Empetrichthys erdisi, a Pliocene cyprinodontid fish from California, with re-marks on the Fundulinae and Cyprinodontinae. Copeia. 1962;1962:520–32.

    Article  Google Scholar 

  116. 116.

    Álvarez J, Arriola-Longoria J. Primer goodeido fósil procedente del Plioceno jalisciense (Pisces: Teleostomi). Boletín de la Sociedad de Ciencias Naturales de Jalisco. 1972;6:6–15.

    Google Scholar 

Download references


We would like to thank the following natural history collections and their museum curators and staff for the loan of specimens, including Tulane University (Hank Bart and Justin Mann), the University of Michigan Museum of Zoology (Bill Fink and Doug Nelson), and the Field Museum (Caleb McMahan). We also would like to thank Brandon Senger, Kevin Guadalupe, and the Nevada Department of Wildlife for providing tissue samples of Crenichthys and Empetrichthys. We would like to thank Cooper Campbell for the goodeid illustrations.


Funding for this project was provided by the National Science Foundation (DEB 1354930) to KRP.

Availability of data and materials

Supplementary materials, data sets, and Genbank accession numbers are available in Dryad digital repository (, doi:

Author information




Both authors agreed with the structure of the analysis and agreed on the final manuscript, KRP provided funding, materials and intellectual structure and insight. KLF gathered the data and conducted the analyses. Both KLF and KRP conceptualized the study and wrote the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Kyle R. Piller.

Ethics declarations

Ethics approval and consent to participate

Samples used in this study either were obtained from natural history collections, natural resource agencies, or donated by colleagues. These samples were obtained from the field under the appropriate permits and followed the animal care protocols of these respective institutions.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional file

Additional file 1:

Table S1. Genbank accession numbers and locality information used to generated the time calibrated phylogeny. (XLSX 13 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( 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

Foster, K.L., Piller, K.R. Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae). BMC Evol Biol 18, 116 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Ecological opportunity
  • Diversification
  • Fishes
  • Mexico
  • Body shape
  • Phylogenetics