- Research article
- Open Access
Disentangling the drivers of diversification in an imperiled group of freshwater fishes (Cyprinodontiformes: Goodeidae)
© The Author(s). 2018
- Received: 4 December 2017
- Accepted: 20 June 2018
- Published: 18 July 2018
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.
- Ecological opportunity
- Body shape
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–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–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–17]. In fact, the majority of the most heavily studied examples of adaptive radiation are related to trophic and/or habitat specializations [18–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–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 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–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 . 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–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 , 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–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 . 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).
Major divergence time estimates with the Goodeidae from the BEAST dating analysis
Split of Empetrichythinae and Goodeinae
Split Empetrichythys from Crenichthys
Split of Ilyodontini
Split of Girardinichthyini and Goodiini
Split of Characodon
Split of Ilyodon and Xenotaenia from Allodontichythys
Split of Goodea and Ataeniobius
Split of Allotoca
Split of Chapalichthyini
Goodeidae variation in body shape
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
9.839 × 10−4
1.900 × 10−3
2.368 × 10−3
3.730 × 10−3
LM = 9–14
5.567 × 10−3
1.397 × 10−2
A = LM 1–8, 15–18
B = LM 9–14
2.74 × 10−3
1.07 × 10− 3
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
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–49] or habitat-based influences [11, 50–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–64], examination of univariate phenotypic traits [25, 65, 66], and quantification of trophic based morphological differences [67–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 . One other group of fishes has diversified in the Mesa Central (Chirostoma: Atherinopsidae) via an adaptive radiation [71–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 . 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  or for select groups within the subfamily .
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 . 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 . 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–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 . Studies have investigated diversification shifts within the group [77, 84, 85], or have examined the correlation between morphology and ecomorphological variables  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 .
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 . Goodeinae taxa have diversified into various aquatic systems (canals, streams, rivers, lakes, ditches, and outflows) , whereas the Empetrichthyinae are confined to homogeneous habitats including small springs and pools . 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–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 , although there are examples of radiations occurring under non-island like conditions . 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.
Most mitochondrial DNA sequences of cytochrome b (cytb) were obtained from Genbank (~ 1140 bp, http://www.ncbi.nlm.nih.gov/Genbank/) for the Goodeinae  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 .
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  Sequence data were submitted to GenBank (Additional file 1: Table S1).
Fossil calibrations and parameters used in BEAST to construct a dated phylogeny
A) †Polymixia sp.
McMahan et al. 2013 
B) †Gymnogeophagus eocenicus
McMahan et al. 2013 
C) †Plesioheros and †Tremembichthys
McMahan et al. 2013 
Marchio and Piller 2013 
E) Poeciliidae and Anablepidae
Santini et al. 2009 
Jordan 1923; Uyeno and Miller 1962 
Alvarez and Arriola-Longoria 1972 
Rate of body shape evolution analyses
We used the multivariate rates methodology developed by  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 .
A Procrustes superimposition analysis was carried out in the program MorphoJ , which corrects for scaling, rotation, translation, and size biases. All other analyses were performed in R  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 . 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  and  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  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.
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.
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.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
- Slowinski JB, Guyer C. Testing the stochasticity of patterns of organismal diversity: an improved null model. Amer Nat. 1989;134:907–21.View ArticleGoogle Scholar
- Heard SB. Patterns in phylogenetic tree balance with variable and evolving speciation rates. Evolution. 1996;50:2141–8.PubMedView ArticleGoogle Scholar
- Ricklefs RE. Estimating diversification rates from phylogenetic information. Trends Ecol Evolut. 2007;22:601–10.View ArticleGoogle Scholar
- 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
- Fischer AG. Latitudinal variations in organic diversity. Evolution. 1992;14:64–81.View ArticleGoogle Scholar
- Hutchinson GE. Homage to Santa Rosalia or why are there so many kinds of animals? Amer Nat. 1959;870:145–59.View ArticleGoogle Scholar
- Foote M. Contributions of individual taxa to overall morphological disparity. Paleobiology. 1993;19:403–19.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- McPeek MA, Brown JM. Clade age and not diversification rate explains species richness among animal taxa. Amer Nat. 2007;169:97–106.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- Magallon S, Sanderson MJ. Absolute diversification rates in angiosperm clades. Evolution. 2001;55:1762–80.PubMedView ArticleGoogle Scholar
- Erwin DH. A preliminary classification of evolutionary radiations. Hist Biol. 1992;2:133–47.View ArticleGoogle Scholar
- Foote M. Discordance and concordance between morphological and taxonomic diversity. Paleobiology. 1993;19:185–204.View ArticleGoogle Scholar
- Heard SB, Hauser DL. Key evolutionary innovations and their ecological mechanisms. Hist Biol. 1995;10:151–73.View ArticleGoogle Scholar
- Hunter JP. Key innovations and the ecology of macroevolution. Trends Ecol Evol. 1998;13:31–6.PubMedView ArticleGoogle Scholar
- Greenwood PJ. Morphology, endemism and speciation in African cichlid fishes. Verhandlungen der Deutschen Zoologischen Gesellschaft. 1973;66:115–24.Google Scholar
- Losos JB. Ecomorphology, performance capability, and scaling of west Indian Anolis lizards: an evolutionary analysis. Ecol Monogr. 1990;60:369–88.View ArticleGoogle Scholar
- Chakrabarty P. Testing conjectures about morphological diversity in cichlids of lakes Malawi and Tanganyika. Copeia. 2005;2:359–73.View ArticleGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Schluter D. The ecology of adaptive radiation. Oxford: Oxford University Press; 2000.Google Scholar
- 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.View ArticleGoogle Scholar
- Parent CE, Crespi BJ. Ecological opportunity in adaptive radiation of Galapagos endemic land snails. Amer Nat. 2009;174:898–905.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Simpson GG. Tempo and mode in evolution Columbia University press; 1944.Google Scholar
- Vermeij GJ. Biological versatility and earth history. Proc Natl Acad Sci U S A. 1973;70:1936–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Jordan DS. A classification of fishes, including families and genera as far as known. Stanford: Stanford University Press; 1923.Google Scholar
- Parenti LR. A phylogenetic and biogeographic analysis of cyprinodontiform fishes Teleostei, Atherinomorpha. Bull Am Mus Nat Hist. 1981;168:341–547.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- Jordan DS, Evermann BW. The fishes of north and middle America. Bulletin of the United States National Museum. 1930;47:2183–3136.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- Miller RR, Minckley WL, Norris SM. Freshwater fishes of Mexico. Chicago: The University of Chicago Press; 2005.Google Scholar
- 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
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- Lombardi J, Wourms JP. Embryonic growth and trophotaenial development in goodeid fishes (Teleostei: Atheriniformes). J Morphol. 1988;197:193–208.PubMedView ArticleGoogle Scholar
- 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
- 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-32View ArticleGoogle Scholar
- 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.Google Scholar
- Williams JE. Threatened fishes of the world: Empetrichthys latos miller, 1948 (Cyprinodontidae). Environ Biol Fish. 1996;45:272.View ArticleGoogle Scholar
- Williams CD, Williams JE. Summer food habits of fishes from two springs in east-Central Nevada. Southwest Nat. 1982;27:437–45.View ArticleGoogle Scholar
- Stadler T. Mammalian phylogeny reveals recent diversification rate shifts. ProNat Acad Sci USA. 2011;108:6187–92.View ArticleGoogle Scholar
- Collar DC, O'Meara BC, Wainwright PC, Near TJ. Piscivory limits diversification of feeding morphology in centrarchid fishes. Evolution. 2009;63:1557–73.PubMedView ArticleGoogle Scholar
- Thomas GH, Meiri S, Phillimore AB. Body size diversification in Anolis: novel environment and island effects. Evolution. 2009;63:2017–30.PubMedView ArticleGoogle Scholar
- 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
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Sabaj MH. Standard symbolic codes for institutional resource collections in herpetology and ichthyology: an Online Reference. Version 6.5. 2016. Electronically accessible at http://www.asih.org/, American Society of Ichthyologists and Herpetologists, Washington, DC.
- Zelditch ML, Swiderski DL, Sheets HD, Fink WL. Geometric Morphometrics for Biologists: A Primer. New York: Elsevier Academic Press; 2004.Google Scholar
- Barbour CD. A biogeographical history of Chirostoma (Pisces: Atherinidae): a species flock from the Mexican plateau. Copeia. 1973;1973:533–66.View ArticleGoogle Scholar
- 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
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Webb SA. Molecular systematics of the genus Allodontichthys (Cyprinodontiformes: Goodeidae). Rev Fish Biol Fish. 2002;12:193–205.View ArticleGoogle Scholar
- Muschick M, Indermaura A, Salzburger W. Convergent evolution within an adaptive radiation of cichlid fishes. Curr Biol. 2012;22:2362–8.PubMedView ArticleGoogle Scholar
- Day JJ, Cotton JA, Barraclough TG. Tempo and mode of diversification of Lake Tanganyika cichlid fishes. PLoS One. 2008; doi.org/10.1371/journal.pone.0001730
- Kocher T. Adaptive evolution and explosive speciation: the cichlid fish model. Nature Rev Genet. 2004;5:288–98.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Liem KF. Evolutionary strategies and morphological innovations: cichlid pharyngeal jaws. Syst Zool. 1974;20:425–41.Google Scholar
- Liem KF. Aquatic versus terrestrial feeding modes: possible impacts on the trophic ecology of vertebrates. Amer Zool. 1990;30:209–22.View ArticleGoogle Scholar
- Hulsey CD, Hollingsworth PR, Fordyce JA. Temporal diversification of central American cichlids. BMC Evol Biol. 2010;10:1–12.View ArticleGoogle Scholar
- Robinson BW, Wilson DS. Character release and displacement in fishes: a neglected literature. Amer Nat. 1994;1:596–627.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- McMahan CD, Chakrabarty P, Sparks JS, Smith WML, Davis MP. Temporal patterns of diversification across global cichlid biodiversity (Acanthomorpha: Cichlidae). PLoS One. 2013; doi.org/10.1371/journal.pone.007116
- 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.PubMedView ArticleGoogle Scholar
- Minckley WL, Marsh PC. Inland fishes of the greater southwest: Chronicle of a Vanishing Biota. Tucson: University of Arizona Press; 2009.Google Scholar
- Langerhans RB. Predictability of phenotypic differentiation across flow regimes in fishes. Integr Comp Biol. 2008;48:750–68.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Webb PW. Fast-start performance and body form in seven species of teleost fish. J Exp Biol. 1978;4:211–26.Google Scholar
- Gatz AJ Jr. Ecological morphology of freshwater stream fishes. Tulane Stud Zool Bot. 1979;21:91–124.Google Scholar
- Gatz AJ Jr. Community organization in fishes as indicated by morphological features. Ecology. 1979:711–8.Google Scholar
- Triantafyllou MS, Triantafyllou GS, Yue DKP. Hydrodynamics of fishlike swimming? Annu Rev Fluid Mech. 2000;32:33–53.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Drummond AJ, et al. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Kingman JFC. The coalescent. Stoch Process Appl. 1982;13:235–48.View ArticleGoogle Scholar
- Rambaut A, Drummond AJ, Xie D, Baele G, Suchard MA. Tracer v1.7, 2018. http://tree.bio.ed.ac.uk/software/tracer/
- Rohlf, FJ. tpsDig, digitize landmarks and outlines, 2.05: Department of Ecology and Evolution, State University of New York at Stony Brook. 2006. http://life.bio.sunysb.edu/morph/index.html
- Adams DC, James Rohlf F, Slice DE. Geometric morphometrics: ten years of progress following the ‘revolution’. Ital J Zool. 2004;71:5–16.View ArticleGoogle Scholar
- McArdle BH, Anderson MJ. Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology. 2001;82:290–7.View ArticleGoogle Scholar
- Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol. 2011;11:353–7.View ArticleGoogle Scholar
- 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, https://doi.org/10.1186/1471-2148-5-17.
- Simpson GG. The major features of evolution. New York: Columbia University Press; 1953.Google Scholar
- 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
- Klingenberg CP. Morphological integration and developmental modularity. Annu Rev Ecol Syst. 2008;39:115–32.View ArticleGoogle Scholar
- Clune J, Mouret J, Lipson H. The evolutionary origins of modularity. Proc R Soc Lond B Biol Sci. 2013;280:2812–63.View ArticleGoogle Scholar
- Sidlauskas B. Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008;62:3135–56.PubMedView ArticleGoogle Scholar
- Klingenberg CP, Gidaszewski NA. Testing and quantifying phylogenetic signals and homoplasy in morphometric data. Syst Biol. 2010;59:245–61.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Moreno E, Carrascal LM. Leg morphology and feeding postures in four Parus species: an experimental ecomorphological approach. Ecology. 1993;4:2037–44.View ArticleGoogle Scholar
- 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
- 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
- Armbruster WS. Evolution of plant pollination systems: hypotheses and tests with the Neotropical vine Dalechampia. Evolution. 1993;12:1480–505.View ArticleGoogle Scholar
- Harrison MK, Crespi BJ. A phylogenetic test of ecomorphological adaptation in Cancer crabs. Evolution. 1999;99:961–5.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Adams DC. Quantifying and comparing phylogenetic evolutionary rates for shape and other high-dimensional phenotypic data. Syst Biol. 2014;63:166–78.PubMedView ArticleGoogle Scholar
- Harmon LJ, Weir J, Brock C, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–31.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Grant PR, Grant BR. How and why species multiply the radiation of Darwin’s finches. Princeton: University Press; 2008.Google Scholar
- Losos JB. Lizards in an evolutionary tree: ecology and adaptive radiation of anoles. Berkeley: University of California Press; 2009.Google Scholar
- Stroud JT, Losos JB. Ecological opportunity and adaptive radiation. Annu Rev Ecol Evol Syst. 2016;47:507–32.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Schluter D. Ecology and the origin of species. Trends Ecol Evol. 2001;16:372–80.PubMedView ArticleGoogle Scholar
- Rabosky DL. Likelihood methods for detecting temporal shifts in diversification rates. Evolution. 2006;60:1152–64.PubMedView ArticleGoogle Scholar
- Marchio EA, Piller KR. Cryptic diversity in a widespread live-bearing fish (Poeciliidae: Belonesox). Biol J Linn Soc. 2013;109:848–60.View ArticleGoogle Scholar
- 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. doi.org/10.1186/1471-2148-9-194 PubMedPubMed CentralView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Á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