Discussions aimed at resolution of the Tree of Life are most often focused on the interrelationships of major organismal lineages. In this study, we focus on the resolution of some of the most apical ...branches in the Tree of Life through exploration of the phylogenetic relationships of darters, a species-rich clade of North American freshwater fishes. With a near-complete taxon sampling of close to 250 species, we aim to investigate strategies for efficient multilocus data sampling and the estimation of divergence times using relaxed-clock methods when a clade lacks a fossil record. Our phylogenetic data set comprises a single mitochondrial DNA (mtDNA) gene and two nuclear genes sampled from 245 of the 248 darter species. This dense sampling allows us to determine if a modest amount of nuclear DNA sequence data can resolve relationships among closely related animal species. Darters lack a fossil record to provide age calibration priors in relaxed-clock analyses. Therefore, we use a near-complete species-sampled phylogeny of the perciform clade Centrarchidae, which has a rich fossil record, to assess two distinct strategies of external calibration in relaxed-clock divergence time estimates of darters: using ages inferred from the fossil record and molecular evolutionary rate estimates. Comparison of Bayesian phylogenies inferred from mtDNA and nuclear genes reveals that heterospecific mtDNA is present in approximately 12.5% of all darter species. We identify three patterns of mtDNA introgression in darters: proximal mtDNA transfer, which involves the transfer of mtDNA among extant and sympatric darter species, indeterminate introgression, which involves the transfer of mtDNA from a lineage that cannot be confidently identified because the introgressed haplotypes are not clearly referable to mtDNA haplotypes in any recognized species, and deep introgression, which is characterized by species diversification within a recipient clade subsequent to the transfer of heterospecific mtDNA. The results of our analyses indicate that DNA sequences sampled from single-copy nuclear genes can provide appreciable phylogenetic resolution for closely related animal species. A well-resolved near-complete species-sampled phylogeny of darters was estimated with Bayesian methods using a concatenated mtDNA and nuclear gene data set with all identified heterospecific mtDNA haplotypes treated as missing data. The relaxed-clock analyses resulted in very similar posterior age estimates across the three sampled genes and methods of calibration and therefore offer a viable strategy for estimating divergence times for clades that lack a fossil record. In addition, an informative rank-free clade-based classification of darters that preserves the rich history of nomenclature in the group and provides formal taxonomie communication of darter clades was constructed using the mtDNA and nuclear gene phylogeny. On the whole, the appeal of mtDNA for phylogeny inference among closely related animal species is diminished by the observations of extensive mtDNA introgression and by finding appreciable phylogenetic signal in a modest sampling of nuclear genes in our phylogenetic analyses of darters.
Eastern North America is the location of the world's most species-rich temperate freshwater fish fauna. Hypotheses regarding the geographic and temporal scale of teleost diversification in this ...region have not been broadly investigated using absolute divergence time estimates among the constituent lineages. This study used time-calibrated molecular phylogenies estimated from mitochondrial and nuclear genes to investigate the temporal and geographic signatures of diversification within barcheek darters, a clade of allopatrically distributed species endemic to the Eastern Highlands. Results from divergence time estimates using an uncorrelated lognormal model suggest that the barcheek darters are an ancient group with a crown node estimate of 16.3 mya, 95% highest posterior density (HPD): 12.4, 20.5, and the clade is characterized by substantial intraspecific divergence times within several species. In particular, the Caney Fork endemic Etheostoma basilare comprises five strongly supported and deeply divergent clades with a most recent common ancestor estimated at 8.0 mya, 95% HPD: 5.6, 10.7. These results are concordant with the hypothesis that geologically stable areas of eastern North America have facilitated both the generation and preservation of lineages across a substantial breadth of evolutionary time, and that allopatric speciation in darters has occurred at much smaller spatial scales than previously realized.
Interspecific divergence along a benthic to pelagic habitat axis is ubiquitous in freshwater fishes inhabiting lentic environments. In this study, we examined the influence of this habitat axis on ...the macroevolution of a diverse, lotic radiation using mtDNA and nDNA phylogenies for eastern North America's most species-rich freshwater fish clade, the open posterior myodome (OPM) cyprinids. We used ancestral state reconstruction to identify the earliest benthic to pelagic transition in this group and generated fossil-calibrated estimates of when this shift occurred. This transition could have represented evolution into a novel adaptive zone, and therefore, we tested for a period of accelerated lineage accumulation after this historical habitat shift.
Ancestral state reconstructions inferred a similar and concordant region of our mtDNA and nDNA based gene trees as representing the shift from benthic to pelagic habitats in the OPM clade. Two independent tests conducted on each gene tree suggested an increased diversification rate after this inferred habitat transition. Furthermore, lineage through time analyses indicated rapid early cladogenesis in the clade arising after the benthic to pelagic shift.
A burst of diversification followed the earliest benthic to pelagic transition during the radiation of OPM cyprinids in eastern North America. As such, the benthic/pelagic habitat axis has likely influenced the generation of biodiversity across disparate freshwater ecosystems.
Cichlid fishes are classic examples of adaptive radiation because of their putative tendency to explosively diversify after invading novel environments. To examine whether ecological opportunity ...increased diversification (speciation minus extinction) early in a species-rich cichlid radiation, we determined if Heroine cichlids experienced a burst of diversification following their invasion of Central America.
We first reconstructed the Heroine phylogeny and determined the basal node to use as the root of Central American Heroine diversification. We then examined the influence of incomplete taxon sampling on this group's diversification patterns. First, we added missing species randomly to the phylogeny and assessed deviations from a constant rate of lineage accumulation. Using a range of species numbers, we failed to recover significant deviations from a pure-birth process and found little support for an early burst of diversification. Then, we examined patterns of lineage accumulation as nodes were increasingly truncated. We assumed that as we removed more recently diverged lineages that sampling would become more complete thereby increasing the power to detect deviations from a pure-birth model. However, truncation of nodes provided even less support for an early burst of diversification.
Contrary to expectations, our analyses suggest Heroine cichlids did not undergo a burst of diversification when they invaded from South America. Throughout their history in Central America, Heroine cichlids appear to have diversified at a constant rate.
The Greater Antilles could have provided a land-bridge between South and Central America for many organisms including cichlid fishes. Using both a coalescent and concatenated phylogenetic species ...tree of 10 genes, we reconstructed the phylogeny of 22 heroine cichlid lineages to test whether the cichlids of the Greater Antilles are the sister group to all heroine cichlids in Central America. Display omitted
► Ten gene species tree of Central American and Caribbean cichlids. ► Individual gene trees topologically very discordant. ► The concatenated and coalescent-based tree differed in support and topologically. ► Greater Antillean cichlids phylogenetically nested within Central American heroines.
Heroine cichlids are major components of the fish faunas in both Central America and the Caribbean. To examine the evolutionary patterns of how cichlids colonized both of these regions, we reconstructed the phylogenetic relationships among 23 cichlid lineages. We used three phylogenetically novel nuclear markers (Dystropin b, Myomesin1, and Wnt7b) in combination with sequence data from seven other gene regions (Nd2, Rag1, Enc1, Sreb2, Ptr, Plagl2, and Zic1) to elucidate the species tree of these cichlids. The species examined represent major heroine lineages in South America, Central America, and the Greater Antilles. The individual gene trees of these groups were topologically quite discordant. Therefore, we combined the genetic partitions and inferred the species tree using both concatenation and a coalescent-based Bayesian method. The two resulting phylogenetic topologies were largely concordant but differed in two fundamental ways. First, more nodes in the concatenated tree were supported with substantial or 100% Bayesian posterior support than in the coalescent-based tree. Second, there was a minor, but biogeographically critical, topological difference between the concatenated and coalescent-based trees. Nevertheless, both analyses recovered topologies consistent with the Greater Antillean heroines being phylogenetically nested within the largely Central American heroine radiation. This study suggests that reconstructions of cichlid phylogeny and historical biogeography should account for the vagaries of individual gene histories.
Display omitted
► Seven gene species tree of eastern North America minnows ► Individual gene trees topologically very discordant ► Concatenated and coalescent-based trees differed slightly in support ...and topology ► Phenacobius more closely related to the more pelagic shiners than to the benthic chubs
Most eastern North American cyprinid fishes belong to a clade known as the “open posterior myodome” (OPM) minnows, but phylogenetic relationships within this clade have been difficult to ascertain. Previous attempts to resolve relationships among the generally benthic “chubs” and the more pelagic “shiners”, that constitute the majority of OPM minnows, have led to highly discordant phylogenetic hypotheses. To further examine relationships among the OPM minnows, we utilized both a concatenated Bayesian approach and a coalescent-based species tree method to analyze data from six protein coding nuclear loci (Enc1, Ptr, Ryr3, Sh3px3, Tbr1, and Zic1), as well as the mitochondrial locus (Cytb). We focused our analyses on the chub-like genus Phenacobius, a group that has drifted topologically between other benthic chubs and the more pelagic shiners, and also included exemplar taxa from 11 other OPM lineages. Individual gene trees were highly discordant regarding relationships within Phenacobius and across the OPM clade. The concatenated Bayesian analysis and coalescent-based species tree reconstruction recovered slightly different phylogenetic topologies. Additionally, the posterior support values for clades using the coalescent-based approach were consistently lower than the concatenated analysis. However, Phenacobius was resolved as monophyletic and as the sister lineage to Erimystax regardless of the combined data approach taken. Furthermore, Phenacobius+Erimystax was recovered as more closely related to the shiners we examined than to other chubs. Relationships within Phenacobius varied depending on the combined phylogenetic method utilized. Our results highlight the importance of multi-locus, coalescent-based approaches for resolving the phylogeny of diverse clades like the eastern North American OPM minnows.
The ability of Perciform fishes to protrude their jaw has likely been critical to the trophic diversification of this group, which includes approximately 20% of all vertebrates. The length of the ...ascending process of the premaxilla is thought to influence the maximum extent that cichlids and other Perciforms protrude their oral jaw. Using a combination of morphometrics, kinematics, and new phylogenetic hypotheses for 20 Heroine cichlid species, we tested the evolutionary relationship between the length of the premaxillary ascending process and maximum jaw protrusion. In this clade, the length of the ascending process of the premaxilla ranged from 11.6-32.7% with respect to standard length whereas maximum jaw protrusion ranged from 3.5-23.4% with respect to standard length. The evolutionary relationships among the Heroine cichlids obtained from the genetic partitions cytochrome b, S7, and RAG1 showed limited concordance. However, correlations between the length of the ascending process and maximum jaw protrusion were highly significant when examined as independent contrasts using all three topologies. Evolutionary change in the length of the ascending process of the premaxilla is likely critical for determining the amount of jaw protrusion in Perciform groups such as cichlid fishes.