The genomic timeline of cichlid fish diversification across continents

The macroevolutionary dynamics of pharyngognathy in fishes fail to support the key innovation hypothesis

Key innovations, traits that provide species access to novel niches, are thought to be a major generator of biodiversity. One commonly cited example of key innovation is pharyngognathy, a set of modifications to the pharyngeal jaws found in some highly species-rich fish clades such as cichlids and wrasses. Here, using comparative phylogenomics and phylogenetic comparative methods, we investigate the genomic basis of pharyngognathy and the impact of this innovation on diversification. Whole genomes resolve the relationships of fish clades with this innovation and their close relatives, but high levels of topological discordance suggest the innovation may have evolved fewer times than previously thought. Closer examination of the topology of noncoding elements accelerated in clades with the pharyngognathy innovation reveals hidden patterns of shared ancestry across putatively independent transitions to pharyngognathy. When our updated phylogenomic relationships are used alongside large-scale phylogenetic and ecological datasets, we find no evidence pharyngognathy consistently modifies the macroevolutionary landscape of trophic ecology nor does it increase diversification. Our results highlight the necessity of incorporating genomic information in studies of key innovation.

Genome assembly at chromosome scale with telomere ends for Pearlspot, Etroplus suratensis

The pearlspot, Etroplus suratensis is a climate resilient cichlid fish that exhibits unusual adaptation to salinity. The fish is able to complete full life cycle in diverse salinity habitats ranging from fresh water to marine environments. High-quality primary and phased genome assemblies were generated for pearlspot fish using PacBio HiFi and Arima HiC sequencing technologies, for the first time. The primary assembly is highly contiguous with contig N50 length of 36 Mb. The final assembly is of 1.247 Gb with N50 length of 51.57 Mb and 98% of the genome length anchored to 24 chromosomes. The genome was assessed to be 99.9% complete based on BUSCO evaluation and was predicted to contain 52.96% repeat elements. We have predicted 27,192 protein encoding genes, of which 21,580 were functionally annotated. The genome offers an invaluable resource to understand adaptation of pearlspot fish to diverse salinity habitats.

Replicated Functional Evolution in Cichlid Adaptive Radiations

AbstractAdaptive radiations highlight the mechanisms by which species and traits diversify and the extent to which these patterns are predictable. We used 1,110 high-speed videos of suction feeding to study functional and morphological diversification in 300 cichlid species from three African Great Lake radiations of varying ages (Victoria, Malawi, and Tanganyika) and an older, spatially dispersed continental radiation in the Neotropics. Among African radiations, standing diversity was reflective of time. Morphological and functional variance in Lake Victoria, the youngest radiation, was a subset of that within Lake Malawi, which itself was nested within the older Tanganyikan radiation. However, functional diversity in Neotropical cichlids was often lower than that in Lake Tanganyika, despite being much older. These two radiations broadly overlapped, but each diversified into novel trait spaces not found in the youngest lake radiations. Evolutionary rates across radiations were inversely related to age, suggesting extremely rapid trait evolution at early stages, particularly in lake radiations. Despite this support for early bursts, other patterns of trait diversity were inconsistent with expectations of adaptive radiations. This work suggests that cichlid functional evolution has played out in strikingly similar fashion in different radiations, with contingencies eventually resulting in lineage-specific novelties.

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

The genetics of niche-specific behavioral tendencies in an adaptive radiation of cichlid fishes

Behavior is critical for animal survival and reproduction, and possibly for diversification and evolutionary radiation. However, the genetics behind adaptive variation in behavior are poorly understood. In this work, we examined a fundamental and widespread behavioral trait, exploratory behavior, in one of the largest adaptive radiations on Earth, the cichlid fishes of Lake Tanganyika. By integrating quantitative behavioral data from 57 cichlid species (702 wild-caught individuals) with high-resolution ecomorphological and genomic information, we show that exploratory behavior is linked to macrohabitat niche adaptations in Tanganyikan cichlids. Furthermore, we uncovered a correlation between the genotypes at a single-nucleotide polymorphism upstream of the AMPA glutamate-receptor regulatory gene cacng5b and variation in exploratory tendency. We validated this association using behavioral predictions with a neural network approach and CRISPR-Cas9 genome editing.

Diet and habitat as determinants of intestine length in fishes

Fish biologists have long assumed a link between intestinal length and diet, and relative gut length or Zihler's index are often used to classify species into trophic groups. This has been done for specific fish taxa or specific ecosystems, but not for a global fish dataset. Here, we assess these relationships across a dataset of 468 fish species (254 marine, 191 freshwater, and 23 that occupy both habitats) in relation to body mass and fish length. Herbivores had significantly relatively stouter bodies and longer intestines than omni- and faunivores. Among faunivores, corallivores had longer intestines than invertivores, with piscivores having the shortest. There were no detectable differences between herbivore groups, possibly due to insufficient understanding of herbivorous fish diets. We propose that reasons for long intestines in fish include (i) difficult-to-digest items that require a symbiotic microbiome, and (ii) the dilution of easily digestible compounds with indigestible material (e.g., sand, wood, exoskeleton). Intestinal indices differed significantly between dietary groups, but there was substantial group overlap. Counter-intuitively, in the largest dataset, marine species had significantly shorter intestines than freshwater fish. These results put fish together with mammals as vertebrate taxa with clear convergence in intestine length in association with trophic level, in contrast to reptiles and birds, even if the peculiar feeding ecology of herbivorous fish is probably more varied than that of mammalian herbivores.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11160-024-09853-3.

Toward transparent taxonomy: an interactive web-tool for evaluating competing taxonomic arrangements

Informative and consistent taxonomy above the species level is essential to communication about evolution, biodiversity and conservation, and yet the practice of taxonomy is considered opaque and subjective by non-taxonomist scientists and the public alike. While various proposals have tried to make the basis for the ranking and inclusiveness of taxa more transparent and objective, widespread adoption of these ideas has lagged. Here, we present TaxonomR, an interactive online decision-support tool to evaluate alternative taxonomic classifications. This tool implements an approach that quantifies the criteria commonly used in taxonomic treatments and allows the user to interactively manipulate weightings for different criteria to compare scores for taxonomic groupings under those weights. We use the butterfly taxon Argynnis to demonstrate how different weightings applied to common taxonomic criteria result in fundamentally different genus-level classifications that are predominantly used in different continents and geographic regions. These differences are objectively compared and quantified using TaxonomR to evaluate the kinds of criteria that have been emphasized in earlier classifications, and the nature of the support for current alternative taxonomic arrangements. The main role of TaxonomR is to make taxonomic decisions transparent via an explicit prioritization scheme. TaxonomR is not a prescriptive application. Rather, it aims to be a tool for facilitating our understanding of alternative taxonomic classifications that can, in turn, potentially support global harmony in biodiversity assessments through evidence-based discussion and community-wide resolution of historically entrenched taxonomic tensions.

Host adaptive radiation is associated with rapid virus diversification and cross-species transmission in African cichlid fishes

Adaptive radiations are generated through a complex interplay of biotic and abiotic factors. Although adaptive radiations have been widely studied in the context of animal and plant evolution, little is known about how they impact the evolution of the viruses that infect these hosts, which in turn may provide insights into the drivers of cross-species transmission and hence disease emergence. We examined how the rapid adaptive radiation of the cichlid fishes of African Lake Tanganyika over the last 10 million years has shaped the diversity and evolution of the viruses they carry. Through metatranscriptomic analysis of 2,242 RNA sequencing libraries, we identified 121 vertebrate-associated viruses among various tissue types that fell into 13 RNA and 4 DNA virus groups. Host-switching was commonplace, particularly within the Astroviridae, Metahepadnavirus, Nackednavirus, Picornaviridae, and Hepacivirus groups, occurring more frequently than in other fish communities. A time-calibrated phylogeny revealed that hepacivirus diversification was not constant throughout the cichlid radiation but accelerated 2-3 million years ago, coinciding with a period of rapid cichlid diversification and niche packing in Lake Tanganyika, thereby providing more closely related hosts for viral infection. These data depict a dynamic virus ecosystem within the cichlids of Lake Tanganyika, characterized by rapid virus diversification and frequent host jumping, and likely reflecting their close phylogenetic relationships that lower the barriers to cross-species virus transmission.

Evidence for selfing in a vertebrate from whole-genome sequencing

A growing number of recent genomic studies report asexual parthenogenetic reproduction in a wide range of taxa, including vertebrate species from the reptile, bird, and fish lineages. Yet, self-fertilization (selfing) has been recorded only in a single vertebrate, the mangrove killifish Kryptolebias marmoratus In cichlid fishes, sex determination is notably diverse and can be influenced by the environment, and sequential hermaphroditism has been reported for some species. Here, we present evidence for a case of facultative selfing in the cichlid fish Benitochromis nigrodorsalis, which is otherwise known as biparentally reproducing ovophilic mouthbrooder from Western Africa. Our laboratory observations revealed that a wild-caught individual produced repeatedly viable offspring in absence of a mating partner. By analyzing genome-wide single-nucleotide polymorphism (SNP) data, we compare that individual and two of its offspring to shed light on its reproductive mode. First, our results confirm uniparental reproduction. Second, overall heterozygosity is reduced in the offspring compared with outbred individuals. Retained maternal heterozygosity in the offspring is ∼51%, which is close to the theoretically expected value of a heterozygosity reduction of 50% by selfing. Heterozygosity patterns along individual chromosomes do not point to alternative parthenogenetic reproductive mechanisms like automixis by terminal or central fusion. Facultative selfing may represent an adaptive strategy ensuring reproduction when mating partners are absent and, hence, contribute to the cichlids' enormous evolutionary success.

The cichlid pharyngeal jaw novelty enhances evolutionary integration in the feeding apparatus

The modified pharyngeal jaw system of cichlid fishes is widely viewed as a key innovation that substantially facilitated the evolutionary exuberance of this iconic evolutionary radiation. We conduct comparative phylogenetic analyses of integration, disparity, and rate of evolution among feeding-related, skeletal structures in Neotropical cichlids and North American centrarchids, which lack the specialized pharyngeal jaw. Contrasting evolutionary patterns in these two continental radiations, we test a classic decoupling hypothesis. Specifically, we ask whether the modified pharyngeal jaw in cichlids resulted in enhanced evolutionary independence of the oral and pharyngeal jaws, leading to increased diversity of trophic structures. Contrary to this prediction, we find significantly stronger evolutionary integration between the oral and pharyngeal jaws in cichlids compared to centrarchids, although the two groups do not differ in patterns of integration within each jaw system. Further, though we find no significant differences in disparity, centrarchids show faster rates of morphological evolution. Our results suggest that the modified pharyngeal jaw resulted in less evolutionary independence and slower rates of evolution within the feeding system. Thus, we raise the possibility that the cichlid novelty enhances feeding performance, but does not prompt increased morphological diversification within the feeding apparatus, as has long been thought.

Mitogenomic Characterization of Cameroonian Endemic Coptodon camerunensis (Cichliformes: Cichlidae) and Matrilineal Phylogeny of Old-World Cichlids

The mitogenomic evolution of old-world cichlids is still largely incomplete in Western Africa. In this present study, the complete mitogenome of the Cameroon endemic cichlid, Coptodon camerunensis, was determined by next-generation sequencing. The mitogenome was 16,557 bp long and encoded with 37 genes (13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a control region). The C. camerunensis mitogenome is AT-biased (52.63%), as exhibited in its congener, Coptodon zillii (52.76% and 53.04%). The majority of PCGs start with an ATG initiation codon, except COI, which starts with a GTG codon and five PCGs and ends with the TAA termination codon and except seven PCGs with an incomplete termination codon. In C. camerunensis mitogenome, most tRNAs showed classical cloverleaf secondary structures, except tRNA-serine with a lack of DHU stem. Comparative analyses of the conserved blocks of two Coptodonini species control regions revealed that the CSB-II block was longer than other blocks and contained highly variable sites. Using 13 concatenated PCGs, the mitogenome-based Bayesian phylogeny easily distinguished all the examined old-world cichlids. Except for Oreochromini and Coptodinini tribe members, the majority of the taxa exhibited monophyletic clustering within their respective lineages. C. camerunensis clustered closely with Heterotilapia buttikoferi (tribe Heterotilapiini) and had paraphyletic clustering with its congener, C. zillii. The Oreochromini species also displayed paraphyletic grouping, and the genus Oreochromis showed a close relationship with Coptodinini and Heterotilapiini species. In addition, illustrating the known distribution patterns of old-world cichlids, the present study is congruent with the previous hypothesis and proclaims that prehistoric geological evolution plays a key role in the hydroclimate of the African continent during Mesozoic, which simultaneously disperses and/or colonizes cichlids in different ichthyological provinces and Rift Lake systems in Africa. The present study suggests that further mitogenomes of cichlid species are required, especially from western Africa, to understand their unique evolution and adaptation.

Expanded olfactory system in ray-finned fishes capable of terrestrial exploration

BACKGROUND

Smell abilities differ greatly among vertebrate species due to distinct sensory needs, with exceptional variability reported in the number of olfactory genes and the size of the odour-processing regions of the brain. However, key environmental factors shaping genomic and phenotypic changes linked to the olfactory system remain difficult to identify at macroevolutionary scales. Here, we investigate the association between diverse ecological traits and the number of olfactory chemoreceptors in approximately two hundred ray-finned fishes.

RESULTS

We found independent expansions producing large gene repertoires in several lineages of nocturnal amphibious fishes, generally able to perform active terrestrial exploration. We reinforced this finding with on-purpose genomic and transcriptomic analysis of Channallabes apus, a catfish species from a clade with chemosensory-based aerial orientation. Furthermore, we also detected an augmented information-processing capacity in the olfactory bulb of nocturnal amphibious fishes by estimating the number of cells contained in this brain region in twenty-four actinopterygian species.

CONCLUSIONS

Overall, we report a convergent genomic and phenotypic magnification of the olfactory system in nocturnal amphibious fishes. This finding suggests the possibility of an analogous evolutionary event in fish-like tetrapod ancestors during the first steps of the water-to-land transition, favouring terrestrial adaptation through enhanced aerial orientation.

Genomics of cold adaptations in the Antarctic notothenioid fish radiation

Numerous novel adaptations characterise the radiation of notothenioids, the dominant fish group in the freezing seas of the Southern Ocean. To improve understanding of the evolution of this iconic fish group, here we generate and analyse new genome assemblies for 24 species covering all major subgroups of the radiation, including five long-read assemblies. We present a new estimate for the onset of the radiation at 10.7 million years ago, based on a time-calibrated phylogeny derived from genome-wide sequence data. We identify a two-fold variation in genome size, driven by expansion of multiple transposable element families, and use the long-read data to reconstruct two evolutionarily important, highly repetitive gene family loci. First, we present the most complete reconstruction to date of the antifreeze glycoprotein gene family, whose emergence enabled survival in sub-zero temperatures, showing the expansion of the antifreeze gene locus from the ancestral to the derived state. Second, we trace the loss of haemoglobin genes in icefishes, the only vertebrates lacking functional haemoglobins, through complete reconstruction of the two haemoglobin gene clusters across notothenioid families. Both the haemoglobin and antifreeze genomic loci are characterised by multiple transposon expansions that may have driven the evolutionary history of these genes.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Cenozoic colonisation of the Indian Ocean region by the Australian freshwater originating glassperch family Ambassidae (Teleostei)

We examined the phylogeny and biogeography of the glassperch family Ambassidae (Teleostei), which is widely distributed in the freshwater, brackish and marine coastal habitats across the Indo-West Pacific region. We first built a comprehensive time-calibrated phylogeny of Ambassidae using five genes. We then used this tree to reconstruct the evolution of the salinity preference and ancestral areas. Our results indicate that the two largest genera of Ambassidae, Ambassis and Parambassis, are each not monophyletic. The most recent common ancestor of Ambassidae was freshwater adapted and lived in Australia about 56 million years ago. Three independent freshwater-to-marine transitions are inferred, but no marine-to-freshwater ones. To explain the distribution of ambassids, we hypothesise two long-distance marine dispersal events from Australia. A first event was towards Southeast Asia during the early Cenozoic, followed by a second one towards Africa during mid-Cenozoic. The phylogenetic signal associated with the salinity adaptation of these events was not detected, possibly because of the selective extinction of intermediate marine lineages. The Ambassidae shares two characteristics with other freshwater fish groups distributed in continental regions surrounding the Indian Ocean: They are too young to support the hypothesis that their distribution is the result of the fragmentation of Gondwana, but they did not retain the phylogenetic signal of their marine dispersal.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

East African cichlid fishes

Cichlid fishes are a very diverse and species-rich family of teleost fishes that inhabit lakes and rivers of India, Africa, and South and Central America. Research has largely focused on East African cichlids of the Rift Lakes Tanganyika, Malawi, and Victoria that constitute the biodiversity hotspots of cichlid fishes. Here, we give an overview of the study system, research questions, and methodologies. Research on cichlid fishes spans many disciplines including ecology, evolution, physiology, genetics, development, and behavioral biology. In this review, we focus on a range of organismal traits, including coloration phenotypes, trophic adaptations, appendages like fins and scales, sensory systems, sex, brains, and behaviors. Moreover, we discuss studies on cichlid phylogenies, plasticity, and general evolutionary patterns, ranging from convergence to speciation rates and the proximate and ultimate mechanisms underlying these processes. From a methodological viewpoint, the last decade has brought great advances in cichlid fish research, particularly through the advent of affordable deep sequencing and advances in genetic manipulations. The ability to integrate across traits and research disciplines, ranging from developmental biology to ecology and evolution, makes cichlid fishes a fascinating research system.

Species diversity and phylogeography of the Australoheros autrani group (Teleostei, Cichlidae) in eastern Brazil

Cichlid fishes are an important model system in evolutionary biology, primarily because of their exceptional diversity. However, while some cichlid assemblages, such as the ones of the African Great Lakes, have received considerable attention, others are not well studied, including many riverine species. Here, we focus on the Australoheros autrani species group and first report a new record of Australoheros in the upper Paranaíba River drainage, extending the known distribution range of this genus. Through Bayesian inference and maximum likelihood phylogenetic analyses of the mitochondrial cytochrome b gene of these specimens as well as available sequences, we assigned the newly discovered population to Australoheros barbosae. We corroborate the monophyly of the A. autrani species group and the presence of three species in the upper/middle Paraíba do Sul River basin as well as molecular diagnostic characters for each. Finally, we provide evidence for a recent expansion of A. barbosae.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10228-022-00888-9.

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 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 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.

The impact of paleoclimatic changes on body size evolution in marine fishes

General rules are useful tools for understanding how organisms evolve. Cope’s rule (tendency to increase in size over evolutionary time) and Bergmann’s rule (tendency to grow to larger sizes in cooler climates) both relate to body size, an important factor that affects the biology, ecology, and physiology of organisms. These rules are well studied in endotherms but remain poorly understood among ectotherms. Here, we show that paleoclimatic changes strongly shaped the trajectory of body size evolution in tetraodontiform fishes. Their body size evolution is explained by both Cope’s and Bergmann’s rules, highlighting the impact of paleoclimatic changes on aquatic organisms, which rely on their environment for temperature regulation and are likely more susceptible than terrestrial vertebrates to climatic changes.

Body size is an important species trait, correlating with life span, fecundity, and other ecological factors. Over Earth’s geological history, climate shifts have occurred, potentially shaping body size evolution in many clades. General rules attempting to summarize body size evolution include Bergmann’s rule, which states that species reach larger sizes in cooler environments and smaller sizes in warmer environments, and Cope’s rule, which poses that lineages tend to increase in size over evolutionary time. Tetraodontiform fishes (including pufferfishes, boxfishes, and ocean sunfishes) provide an extraordinary clade to test these rules in ectotherms owing to their exemplary fossil record and the great disparity in body size observed among extant and fossil species. We examined Bergmann’s and Cope’s rules in this group by combining phylogenomic data (1,103 exon loci from 185 extant species) with 210 anatomical characters coded from both fossil and extant species. We aggregated data layers on paleoclimate and body size from the species examined, and inferred a set of time-calibrated phylogenies using tip-dating approaches for downstream comparative analyses of body size evolution by implementing models that incorporate paleoclimatic information. We found strong support for a temperature-driven model in which increasing body size over time is correlated with decreasing oceanic temperatures. On average, extant tetraodontiforms are two to three times larger than their fossil counterparts, which otherwise evolved during periods of warmer ocean temperatures. These results provide strong support for both Bergmann’s and Cope’s rules, trends that are less studied in marine fishes compared to terrestrial vertebrates and marine invertebrates.

The macroparasite fauna of cichlid fish from Nicaraguan lakes, a model system for understanding host-parasite diversification and speciation

The Nicaraguan lakes represent an ideal continent-island-like setting to study the colonization patterns of both fish and their parasites. The dominant fish fauna are cichlids, particularly the Midas cichlid species complex Amphilophus spp., a well-studied model for recent sympatric speciation. Here, we characterized the Midas cichlid macroparasite diversity in Nicaraguan lakes. We evaluated patterns of parasite diversity across host populations. Morphological and molecular analyses were conducted, revealing a macroparasite fauna composed by 37 taxa, including platyhelminths, nematodes, copepods, branchiurans, hirudineans and oribatids. Three invasive species are reported for the first time. The Midas cichlid was infected by 22 parasite taxa, 18 shared with other cichlids. Eight taxa conformed the core parasite fauna of the Midas cichlid. The large lakes had higher parasite diversity than the smaller and isolated crater lakes, although parasite infracommunity diversity was lower. Environmental factors along with the differential distribution of intermediate hosts, the potential resistance gained by their hosts after colonization of new lakes, competitive exclusion among parasites, or the introduction of exotic fish, may determine the observed pattern of parasite heterogeneous distribution. Our study provides a ground to explore the evolutionary history of both, hosts and parasites within the context of speciation and diversification processes.

Dactylogyridae 2022: a meta-analysis of phylogenetic studies and generic diagnoses of parasitic flatworms using published genetic and morphological data

Dactylogyridae is one of the most studied families of parasitic flatworms with more than 1000 species and 166 genera described to date including ecto- and endoparasites. Dactylogyrid monogeneans were suggested as model organisms for host-parasite macroevolutionary and biogeographical studies due to the scientific and economic importance of some of their host lineages. Consequently, an array of phylogenetic research into different dactylogyrid lineages has been produced over the past years but the last family-wide study was published 16 years ago. Here, we provide a meta-analysis of the phylogenetic relationships of Dactylogyridae including representatives of all genera with available molecular data (n=67). First, we investigate the systematic informativeness of morphological characters widely used to diagnose dactylogyrid genera through a parsimony analysis of the characters, character mapping, and phylogenetic comparative methods. Second, we provide an overview of the current state of the systematics of the family and its subfamilies, and summarise potentially poly- and paraphyletic genera. Third, we elaborate on the implications of taxonomic, citation, and confirmation bias in past studies. Fourth, we discuss host range, biogeographical, and freshwater-marine patterns. We found two well-supported macroclades which we assigned to the subfamilies Dactylogyrinae and Ancyrocephalinae. These subfamilies further include 16 well-supported clades with only a few synapomorphies that could be deduced from generic diagnoses in the literature. Furthermore, few morphological characters considered systematically informative at the genus level display a strong phylogenetic signal. However, the parsimony analysis suggests that these characters provide little information on the relationships between genera. We conclude that a strong taxonomic bias and low coverage of DNA sequences and regions limit knowledge on morphological and biogeographical evolutionary patterns that can be inferred from these results. We propose addressing potential citation and confirmation biases through a 'level playing field' multiple sequence alignment as provided by this study.

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

A Chromosome-Level Genome Assembly of Mozambique Tilapia (Oreochromis mossambicus) Reveals the Structure of Sex Determining Regions

The Mozambique tilapia (Oreochromis mossambicus) is a fascinating taxon for evolutionary and ecological research. It is an important food fish and one of the most widely distributed tilapias. Because males grow faster than females, genetically male tilapia are preferred in aquaculture. However, studies of sex determination and sex control in O. mossambicus have been hindered by the limited characterization of the genome. To address this gap, we assembled a high-quality genome of O. mossambicus, using a combination of high coverage of Illumina and Nanopore reads, coupled with Hi-C and RNA-Seq data. Our genome assembly spans 1,007 Mb with a scaffold N50 of 11.38 Mb. We successfully anchored and oriented 98.6% of the genome on 22 linkage groups (LGs). Based on re-sequencing data for male and female fishes from three families, O. mossambicus segregates both an XY system on LG14 and a ZW system on LG3. The sex-patterned SNPs shared by two XY families narrowed the sex determining regions to ∼3 Mb on LG14. The shared sex-patterned SNPs included two deleterious missense mutations in ahnak and rhbdd1, indicating the possible roles of these two genes in sex determination. This annotated chromosome-level genome assembly and identification of sex determining regions represents a valuable resource to help understand the evolution of genetic sex determination in tilapias.

The cichlid oral and pharyngeal jaws are evolutionarily and genetically coupled

Evolutionary constraints may significantly bias phenotypic change, while "breaking" from such constraints can lead to expanded ecological opportunity. Ray-finned fishes have broken functional constraints by developing two jaws (oral-pharyngeal), decoupling prey capture (oral jaw) from processing (pharyngeal jaw). It is hypothesized that the oral and pharyngeal jaws represent independent evolutionary modules and this facilitated diversification in feeding architectures. Here we test this hypothesis in African cichlids. Contrary to our expectation, we find integration between jaws at multiple evolutionary levels. Next, we document integration at the genetic level, and identify a candidate gene, smad7, within a pleiotropic locus for oral and pharyngeal jaw shape that exhibits correlated expression between the two tissues. Collectively, our data show that African cichlid evolutionary success has occurred within the context of a coupled jaw system, an attribute that may be driving adaptive evolution in this iconic group by facilitating rapid shifts between foraging habitats, providing an advantage in a stochastic environment such as the East African Rift-Valley.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Dynamics of sex chromosome evolution in a rapid radiation of cichlid fishes

Sex is a fundamental trait determined by environmental and/or genetic factors, including sex chromosomes. Sex chromosomes are studied in species scattered across the tree of life, yet little is known about tempo and mode of sex chromosome evolution among closely related species. Here, we examine sex chromosome evolution in the adaptive radiation of cichlid fishes in Lake Tanganyika. Through the analysis of male and female genomes from 244 cichlid taxa (189 described species with 5 represented with two local variants/populations; 50 undescribed species) and of 396 multitissue transcriptomes from 66 taxa, we identify signatures of sex chromosomes in 79 taxa, involving 12 linkage groups. We find that Tanganyikan cichlids have the highest rates of sex chromosome turnover and heterogamety transitions known to date. We show that sex chromosome recruitment is not at random. Moreover convergently emerged sex chromosomes in cichlids support the “limited options” hypothesis of sex chromosome evolution.

Comparative scale morphology in the adaptive radiation of cichlid fishes (Perciformes: Cichlidae) from Lake Tanganyika

The morphology of fish scales has been investigated for &gt; 200 years, but research on evolutionary patterns of scale morphology is scarce. Here, we study scale morphology and its evolution in the adaptive radiation of cichlid fishes from Lake Tanganyika, which are known for their exceptional diversity in habitat use, feeding ecology and morphology. Based on a geometric morphometric approach on eight scales per specimen (covering different body regions), we quantify scale types and morphology across nearly all ~240 species of the cichlid adaptive radiation in Lake Tanganyika. We first show that scale type, shape and ctenii coverage vary along the body, which is probably attributable to adaptations to different functional demands on the respective scales. Our comparative analyses reveal that flank scale size is tightly linked to phylogeny, whereas scale shape and ctenii coverage can be explained only in part by phylogenetic history and/or our proxy for ecology (stable isotopes and body shape), suggesting an additional adaptive component. We also show that our measured scale characteristics can help to assign an individual scale to a taxonomic group or ecotype. Thus, our data may serve as a valuable resource for taxonomic studies and to interpret fossil finds.

Comparing Ultraconserved Elements and Exons for Phylogenomic Analyses of Middle American Cichlids: When Data Agree to Disagree

Choosing among types of genomic markers to be used in a phylogenomic study can have a major influence on the cost, design, and results of a study. Yet few attempts have been made to compare categories of next-generation sequence markers limiting our ability to compare the suitability of these different genomic fragment types. Here, we explore properties of different genomic markers to find if they vary in the accuracy of component phylogenetic trees and to clarify the causes of conflict obtained from different data sets or inference methods. As a test case, we explore the causes of discordance between phylogenetic hypotheses obtained using a novel data set of ultraconserved elements (UCEs) and a recently published exon data set of the cichlid tribe Heroini. Resolving relationships among heroine cichlids has historically been difficult, and the processes of colonization and diversification in Middle America and the Greater Antilles are not yet well understood. Despite differences in informativeness and levels of gene tree discordance between UCEs and exons, the resulting phylogenomic hypotheses generally agree on most relationships. The independent data sets disagreed in areas with low phylogenetic signal that were overwhelmed by incomplete lineage sorting and nonphylogenetic signals. For UCEs, high levels of incomplete lineage sorting were found to be the major cause of gene tree discordance, whereas, for exons, nonphylogenetic signal is most likely caused by a reduced number of highly informative loci. This paucity of informative loci in exons might be due to heterogeneous substitution rates that are problematic to model (i.e., computationally restrictive) resulting in systematic errors that UCEs (being less informative individually but more uniform) are less prone to. These results generally demonstrate the robustness of phylogenomic methods to accommodate genomic markers with different biological and phylogenetic properties. However, we identify common and unique pitfalls of different categories of genomic fragments when inferring enigmatic phylogenetic relationships.

The non-gradual nature of adaptive radiation

Adaptive radiation is a major source of biodiversity. Still, many aspects of this evolutionary process remain poorly understood. Our recent integrative examination of the cichlid adaptive radiation in African Lake Tanganyika provides new insights into the process of explosive diversification. The in-depth phylogenetic comparative analysis of nearly all species occurring in that lake permitted us to trace patterns of eco-morphological evolution throughout the phylogenetic history of the radiation and revealed that it occurred in a non-gradual manner, in the form of time-shifted bursts of accelerated evolution. The temporal pattern of divergence along different axes of morphological evolution provides empirical support for a scenario that potentially unifies two popular models of adaptive radiation, the "early burst model" and the "stages model".

Gene expression dynamics during rapid organismal diversification in African cichlid fishes

Changes in gene expression play a fundamental role in phenotypic evolution. Transcriptome evolutionary dynamics have so far mainly been compared among distantly related species and remain largely unexplored during rapid organismal diversification, in which gene regulatory changes have been suggested as particularly effective drivers of phenotypic divergence. Here we studied gene expression evolution in a model system of adaptive radiation, the cichlid fishes of African Lake Tanganyika. By comparing gene expression profiles of 6 different organs in 74 cichlid species representing all subclades of this radiation, we demonstrate that the rate of gene expression evolution varies among organs, transcriptome parts and the subclades of the radiation, indicating different strengths of selection. We found that the noncoding part of the transcriptome evolved more rapidly than the coding part, and that the gonadal transcriptomes evolved more rapidly than the somatic ones, with the exception of liver. We further show that the rate of gene expression change was not constant over the course of the radiation but accelerated at its later phase. Finally, we show that-at the per-gene level-the evolution of expression patterns is dominated by stabilizing selection.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.