Phylogenomics of a putatively convergent novelty: did hypertrophied lips evolve once or repeatedly in Lake Malawi cichlid fishes?

BACKGROUND: Phylogenies provide critical information about convergence during adaptive radiation. To test whether there have been multiple origins of a distinctive trophic phenotype in one of the most rapidly radiating groups known, we used ultra-conserved elements (UCEs) to examine the evolutionary affinities of Lake Malawi cichlids lineages exhibiting greatly hypertrophied lips. RESULTS: The hypertrophied lip cichlids Cheilochromis euchilus, Eclectochromis ornatus, Placidochromis "Mbenji fatlip", and Placidochromis milomo are all nested within the non-mbuna clade of Malawi cichlids based on both concatenated sequence and single nucleotide polymorphism (SNP) inferred phylogenies. Lichnochromis acuticeps that exhibits slightly hypertrophied lips also appears to have evolutionary affinities to this group. However, Chilotilapia rhoadesii that lacks hypertrophied lips was recovered as nested within the species Cheilochromis euchilus. Species tree reconstructions and analyses of introgression provided largely ambiguous patterns of Malawi cichlid evolution. CONCLUSIONS: Contrary to mitochondrial DNA phylogenies, bifurcating trees based on our 1024 UCE loci supported close affinities of Lake Malawi lineages with hypertrophied lips. However, incomplete lineage sorting in Malawi tends to render these inferences more tenuous. Phylogenomic analyses will continue to provide powerful inferences about whether phenotypic novelties arose once or multiple times during adaptive radiation.

Replicated divergence in cichlid radiations mirrors a major vertebrate innovation.

Decoupling of the upper jaw bones--jaw kinesis--is a distinctive feature of the ray-finned fishes, but it is not clear how the innovation is related to the extraordinary diversity of feeding behaviours and feeding ecology in this group. We address this issue in a lineage of ray-finned fishes that is well known for its ecological and functional diversity--African rift lake cichlids. We sequenced ultraconserved elements to generate a phylogenomic tree of the Lake Tanganyika and Lake Malawi cichlid radiations. We filmed a diverse array of over 50 cichlid species capturing live prey and quantified the extent of jaw kinesis in the premaxillary and maxillary bones. Our combination of phylogenomic and kinematic data reveals a strong association between biting modes of feeding and reduced jaw kinesis, suggesting that the contrasting demands of biting and suction feeding have strongly influenced cranial evolution in both cichlid radiations.

Next generation phylogeography of cave and surface Astyanax mexicanus.

The loss of traits is a commonly observed evolutionary pattern in cave organisms, but due to extensive morphological convergence, inferring relationships between cave and surface populations can be difficult. For instance, Astyanax mexicanus (the blind Mexican cavefish) is thought to have repeatedly lost its eyes following colonization of cave environments, but the number of evolutionarily independent invasions of this species into caves remains unclear. Because of these repeated losses, it has become a model organism for studying the genetic basis of phenotypic trait loss. Here we reconstruct a high-resolution phylogeography for A. mexicanus inferred from both mitochondrial DNA and several thousand single nucleotide polymorphisms. We provide novel insight into the origin of cave populations from the Sabinos and Río Subterráneo caves and present evidence that the Sabinos cave population is part of a unique cave lineage unrelated to other A. mexicanus cave populations. Our results indicate A. mexicanus cave populations have at least four independent origins.

Reconciling gene trees of eastern North American minnows.

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.

Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids.

Phenotypically diverse Lake Malawi cichlids exhibit similar genomes. The extensive sharing of genetic polymorphism among forms has both intrigued and frustrated biologists trying to understand the nature of diversity in this and other rapidly evolving systems. Shared polymorphism might result from hybridization and/or the retention of ancestrally polymorphic alleles. To examine these alternatives, we used new genomic tools to characterize genetic differentiation in widespread, geographically structured populations of Labeotropheus fuelleborni and Metriaclima zebra. These phenotypically distinct species share mitochondrial DNA (mtDNA) haplotypes and show greater mtDNA differentiation among localities than between species. However, Bayesian analysis of nuclear single nucleotide polymorphism (SNP) data revealed two distinct genetic clusters corresponding perfectly to morphologically diagnosed L. fuelleborni and M. zebra. This result is a function of the resolving power of the multi-locus dataset, not a conflict between nuclear and mitochondrial partitions. Locus-by-locus analysis showed that mtDNA differentiation between species (F(CT)) was nearly identical to the median single-locus SNP F(CT). Finally, we asked whether there is evidence for gene flow at sites of co-occurrence. We used simulations to generate a null distribution for the level of differentiation between co-occurring populations of L. fuelleborni and M. zebra expected if there was no hybridization. The null hypothesis was rejected for the SNP data; populations that co-occur at rock reef sites were slightly more similar than expected by chance, suggesting recent gene flow. The coupling of numerous independent markers with extensive geographic sampling and simulations utilized here provides a framework for assessing the prevalence of gene flow in recently diverged species.

Function of a key morphological innovation: fusion of the cichlid pharyngeal jaw.

The pharyngeal jaw of cichlids may represent a key innovation that facilitated their unparalleled trophic divergence. In cichlids, 'fusion' of the lower pharyngeal jaw (LPJ) results from suturing between the two lower ceratobranchials. To examine, what novel abilities a more extensively fused pharyngeal jaw may confer, the function of LPJ suturing was examined in Heroine cichlids. Greater LPJ suturing, pharyngeal jaw splitting under compression and the forces used to crush molluscs in the wild suggest increased LPJ fusion in the trophically polymorphic Herichthys minckleyi operates to strengthen the pharyngeal jaw. Among Heroine cichlid species, the presence of an external LPJ suture and feeding specialization on molluscs was evolutionarily quite variable, but greater LPJ fusion estimated from the amount of external suturing was highly correlated with molluscivory. Throughout cichlid diversification, increased pharyngeal jaw fusion via suturing has likely helped to reinforce the LPJ during pharyngeal processing thereby facilitating the ability of cichlids to exploit durable prey.

Mitochondrial genome primers for Lake Malawi cichlids.

Resolving the evolutionary history of rapidly diversifying lineages like the Lake Malawi Cichlid Flock demands powerful phylogenetic tools. Although this clade of over 500 species of fish likely diversified in less than two million years, the availability of extensive sequence data sets, such as complete mitochondrial genomes, could help resolve evolutionary patterns in this group. Using a large number of newly developed primers, we generated whole mitochondrial genome sequences for 14 Lake Malawi cichlids. We compared sequence divergence across protein-coding regions of the mitochondrial genome and also compared divergence in the mitochondrial loci to divergence at two nuclear protein-coding loci, Mitfb and Dlx2. Despite the widespread sharing of haplotypes of identical sequences at individual loci, the combined use of all protein-coding mitochondrial loci provided a bifurcating phylogenetic hypothesis for the exemplars of major lineages within the Lake Malawi cichlid radiation. The primers presented here could have substantial utility for evolutionary analyses of mitochondrial evolution and hybridization within this diverse clade.

Temporal diversification of Mesoamerican cichlid fishes across a major biogeographic boundary.

The Mexican Neovolcanic Plateau sharply divides the vertebrate fauna of Mesoamerica where the climate of both the neotropics and temperate North America gradually blend. Only a few vertebrate groups such as the Heroine cichlids, distributed from South America to the Rio Grande in North America, are found both north and south of the Neovolcanic Plateau. To better understand the geography and temporal diversification of cichlids at this geologic boundary, we used mitochondrial DNA sequences of the cytochrome b (cyt b) gene to reconstruct the relationships of 52 of the approximately 80 species of Heroine cichlids in Mesoamerica. Our analysis suggests several cichlids in South America should be considered as part of the Mesoamerican Heroine clade because they and the cichlids north of the Isthmus of Panama are clearly supported as monophyletic with respect to all other Neotropical cichlids. We also recovered a group containing species in Paratheraps+Paraneetroplus+Vieja as the sister clade to Herichthys. Herichthys is the only cichlid clade north of the Mexican Plateau and it is monophyletic. Non-parametric rate smoothing of cichlid cyt b sequence resulted in an estimated divergence time of approximately 6 million years for Herichthys. This temporal diversification is concordant with divergence times estimated for anurans in the genus Bufo, a group that exhibits a similar geographic distribution. Our results indicate the 5-million-year-old extension of the Mexican Neovolcanic Plateau to the Gulf Coast of Mexico has strongly influenced the current transition between the vertebrate faunas of the Neotropics and Nearctic.