Exon-based Phylogenomics and the Relationships of African Cichlid Fishes: Tackling the Challenges of Reconstructing Phylogenies with Repeated Rapid Radiations.

African cichlids (subfamily: Pseudocrenilabrinae) are among the most diverse vertebrates, and their propensity for repeated rapid radiation has made them a celebrated model system in evolutionary research. Nonetheless, despite numerous studies, phylogenetic uncertainty persists, and riverine lineages remain comparatively underrepresented in higher-level phylogenetic studies. Heterogeneous gene histories resulting from incomplete lineage sorting (ILS) and hybridization are likely sources of uncertainty, especially during episodes of rapid speciation. We investigate the relationships of Pseudocrenilabrinae and its close relatives while accounting for multiple sources of genetic discordance using species tree and hybrid network analyses with hundreds of single-copy exons. We improve sequence recovery for distant relatives, thereby extending the taxonomic reach of our probes, with a hybrid reference guided/de novo assembly approach. Our analyses provide robust hypotheses for most higher-level relationships and reveal widespread gene heterogeneity, including in riverine taxa. ILS and past hybridization are identified as the sources of genetic discordance in different lineages. Sampling of various Blenniiformes (formerly Ovalentaria) adds strong phylogenomic support for convict blennies (Pholidichthyidae) as sister to Cichlidae and points to other potentially useful protein-coding markers across the order. A reliable phylogeny with representatives from diverse environments will support ongoing taxonomic and comparative evolutionary research in the cichlid model system. [African cichlids; Blenniiformes; Gene tree heterogeneity; Hybrid assembly; Phylogenetic network; Pseudocrenilabrinae; Species tree.].

Evolution, inactivation and loss of short wavelength-sensitive opsin genes during the diversification of Neotropical cichlids.

Natural variation in the number, expression and function of sensory genes in an organism's genome is often tightly linked to different ecological and evolutionary forces. Opsin genes, which code for the first step in visual transduction, are ideal models for testing how ecological factors such as light environment may influence visual system adaptation. Neotropical cichlid fishes are a highly ecologically diverse group that evolved in a variety of aquatic habitats, including black (stained), white (opaque) and clear waters. We used cross-species exon capture to sequence Neotropical cichlid short wavelength-sensitive (SWS) opsins, which mediate ultraviolet (UV) to blue visual sensitivity. Neotropical cichlid SWS1 opsin (UV-sensitive) underwent a relaxation of selective constraint during the early phases of cichlid diversification in South America, leading to pseudogenization and loss. Conversely, SWS2a (blue-sensitive) experienced a burst of episodic positive selection at the base of the South American cichlid radiation. This burst coincides with SWS1 relaxation and loss, and is consistent with findings in ecomorphological studies characterizing a period of extensive ecological divergence in Neotropical cichlids. We use ancestral sequence reconstruction and protein modelling to investigate mutations along this ancestral branch that probably modified SWS2a function. Together, our results suggest that variable light environments played a prominent early role in shaping SWS opsin diversity during the Neotropical cichlid radiation. Our results also illustrate that long-term evolution under light-limited conditions in South America may have reduced visual system plasticity; specifically, early losses of UV sensitivity may have constrained the evolutionary trajectory of Neotropical cichlid vision.

Exon-based phylogenomics strengthens the phylogeny of Neotropical cichlids and identifies remaining conflicting clades (Cichliformes: Cichlidae: Cichlinae).

The phenotypic, geographic, and species diversity of cichlid fishes have made them a group of great interest for studying evolutionary processes. Here we present a targeted-exon next-generation sequencing approach for investigating the evolutionary relationships of cichlid fishes (Cichlidae), with focus on the Neotropical subfamily Cichlinae using a set of 923 primarily single-copy exons designed through mining of the Nile tilapia (Oreochromis niloticus) genome. Sequence capture and assembly were robust, leading to a complete dataset of 415 exons for 139 species (147 terminals) that consisted of 128 Neotropical species, six African taxa, and five Indo-Malagasy cichlids. Gene and species trees were calculated using alternative partitioning schemes and reconstruction methods. In general, all methods yielded similar topologies to previously hypothesized relationships within the Cichlinae and clarified several relationships that were previously poorly supported or in conflict. Additional work will be needed to fully resolve all aspects of Cichlinae phylogeny. Overall, this approach yielded a well-resolved phylogeny of Neotropical cichlids that will be of utility for future assessments of the evolutionary and ecological processes within this diverse group of fishes. Furthermore, the general methodology employed here of exon targeting and capture should be applicable to any group of organisms with the availability of a reference genome.

Accelerated Evolution and Functional Divergence of the Dim Light Visual Pigment Accompanies Cichlid Colonization of Central America.

Cichlids encompass one of the most diverse groups of fishes in South and Central America, and show extensive variation in life history, morphology, and colouration. While studies of visual system evolution in cichlids have focussed largely on the African rift lake species flocks, Neotropical cichlids offer a unique opportunity to investigate visual system evolution at broader temporal and geographic scales. South American cichlid colonization of Central America has likely promoted accelerated rates of morphological evolution in Central American lineages as they encountered reduced competition, renewed ecological opportunity, and novel aquatic habitats. To investigate whether such transitions have influenced molecular evolution of vision in Central American cichlids, we sequenced the dim-light rhodopsin gene in 101 Neotropical cichlid species, spanning the diversity of the clade. We find strong evidence for increased rates of evolution in Central American cichlid rhodopsin relative to South American lineages, and identify several sites under positive selection in rhodopsin that likely contribute to adaptation to different photic environments. We expressed a Neotropical cichlid rhodopsin protein invitro for the first time, and found that while its spectral tuning properties were characteristic of typical vertebrate rhodopsin pigments, the rate of decay of its active signalling form was much slower, consistent with dim light adaptation in other vertebrate rhodopsins. Using site-directed mutagenesis combined with spectroscopic assays, we found that a key amino acid substitution present in some Central American cichlids accelerates the rate of decay of active rhodopsin, which may mediate adaptation to clear water habitats.

A targeted next-generation sequencing toolkit for exon-based cichlid phylogenomics.

Cichlid fishes (family Cichlidae) are models for evolutionary and ecological research. Massively parallel sequencing approaches have been successfully applied to study relatively recent diversification in groups of African and Neotropical cichlids, but such technologies have yet to be used for addressing larger-scale phylogenetic questions of cichlid evolution. Here, we describe a process for identifying putative single-copy exons from five African cichlid genomes and sequence the targeted exons for a range of divergent (>tens of millions of years) taxa with probes designed from a single reference species (Oreochromis niloticus, Nile tilapia). Targeted sequencing of 923 exons across 10 cichlid species that represent the family's major lineages and geographic distribution resulted in a complete taxon matrix of 564 exons (649,549 bp), representing 559 genes. Maximum likelihood and Bayesian analyses in both species tree and concatenation frameworks yielded the same fully resolved and highly supported topology, which matched the expected backbone phylogeny of the major cichlid lineages. This work adds to the body of evidence that it is possible to use a relatively divergent reference genome for exon target design and successful capture across a broad phylogenetic range of species. Furthermore, our results show that the use of a third-party laboratory coupled with accessible bioinformatics tools makes such phylogenomics projects feasible for research groups that lack direct access to genomic facilities. We expect that these resources will be used in further cichlid evolution studies and hope the protocols and identified targets will also be useful for phylogenetic studies of a wider range of organisms.

Colonization and/or mitochondrial selective sweeps across the North Atlantic intertidal assemblage revealed by multi-taxa approximate Bayesian computation.

Intertidal and subtidal communities of the western and eastern coasts of the North Atlantic Ocean were greatly affected by Pleistocene glaciations, with some taxa persisting on both coasts, and others recolonizing after being extirpated on one coast during the Last Glacial Maximum. In the original spirit of comparative phylogeography, we conducted a comparative analysis using mtDNA sequence data and a hierarchical approximate Bayesian computation (ABC) approach for testing these two scenarios across 12 intertidal and subtidal coastal invertebrates spanning the North Atlantic to determine the temporal dynamics of species membership of these two ephemeral communities. Conditioning on a low gene-flow model, our results suggested that a colonization or mitochondrial selective sweep history was predominant across all taxa, with only the bivalve mollusc Mytilus edulis showing a history of trans-Atlantic persistence. Conditioning on a high gene-flow model weakened the support for this assemblage-level demographic history. The predominance of a colonization-type history also highlights concerns about analyses based on single-locus data where genetic hitchhiking may be incorrectly inferred as colonization. In conclusion, driving factors in shifting species range distributions and membership of ephemeral coastal communities could be species-specific environmental tolerances, species interactions, and/or stochastic demographic extinction. Through a re-examination of a long-standing question of North Atlantic phylogeography, we highlight the flexibility and statistical honesty of using a model-based ABC approach.

Molecular resolution of the systematics of a problematic group of fishes (Teleostei: Osmeridae) and evidence for morphological homoplasy.

Relationships among the species of Northern Hemisphere smelts (family Osmeridae) have long been debated in the fish systematics literature. Eight independent studies based on morphological characters failed to reach any consensus on osmerid interrelationships. We reconstruct the osmerid phylogeny based on DNA sequence data from three mitochondrial (cytb, 16S, 12S) and three nuclear (ITS2, S71, RAG1) gene regions from multiple individuals of the 14 species in 6 genera, using the Japanese ayu (Plecoglossus altivelis) as the outgroup. Analyses with different combinations of nuclear and mitochondrial datasets yielded a generally well-resolved phylogeny of the genera that conflicts with previous hypotheses of osmerid interrelationships, and Shimodaira-Hasegawa tests suggest our topology with the current molecular dataset is significantly better than earlier reconstructions. In addition, mapping 114 morphological characters used in previous studies onto our phylogeny shows widespread homoplasy, which is likely the source of the systematic disagreement produced in earlier works.