The taxonomic diversity of the cichlid fish fauna of ancient Lake Tanganyika, East Africa

Diversity, evolution, and emergence of fish viruses

The production of aquatic animals has more than doubled over the last 50 years and is anticipated to continually increase. While fish are recognized as a valuable and sustainable source of nutrition, particularly in the context of human population growth and climate change, the rapid expansion of aquaculture coincides with the emergence of highly pathogenic viruses that often spread globally through aquacultural practices. Here, we provide an overview of the fish virome and its relevance for disease emergence, with a focus on the insights gained through metagenomic sequencing, noting potential areas for future study. In particular, we describe the diversity and evolution of fish viruses, for which the majority have no known disease associations, and demonstrate how viruses emerge in fish populations, most notably at an expanding domestic-wild interface. We also show how wild fish are a powerful and tractable model system to study virus ecology and evolution more broadly and can be used to identify the major factors that shape vertebrate viromes. Central to this is a process of virus-host co-divergence that proceeds over many millions of years, combined with ongoing cross-species virus transmission.

Reference Genome Choice and Filtering Thresholds Jointly Influence Phylogenomic Analyses

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate the extent to which the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count (MAC) >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with MAC >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short-read genomic data for phylogenetic inference.

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.

Lack of host specialization despite selective host use in brood parasitic cuckoo catfish

Host-parasite dynamics involve coevolutionary arms races, which may lead to host specialization and ensuing diversification. Our general understanding of the evolution of host specialization in brood parasites is compromised by a restricted focus on bird and insect lineages. The cuckoo catfish (Synodontis multipunctatus) is an obligate parasite of parental care of mouthbrooding cichlids in Lake Tanganyika. Given the ecological and taxonomic diversity of mouthbrooding cichlids in the lake, we hypothesized the existence of sympatric host-specific lineages in the cuckoo catfish. In a sample of 779 broods from 20 cichlid species, we found four species parasitized by cuckoo catfish (with prevalence of parasitism of 2%-18%). All parasitized cichlids were from the tribe Tropheini, maternal mouthbrooders that spawn over a substrate (rather than in open water). Phylogenetic analysis based on genomic (ddRAD sequencing) and mitochondrial (Dloop) data from cuckoo catfish embryos showed an absence of host-specific lineages. This was corroborated by analyses of genetic structure and co-ancestry matrix. Within host species, parasitism was not associated with any individual characteristic we recorded (parent size, water depth), but was costly as parasitized parents carried smaller clutches of their own offspring. We conclude that the cuckoo catfish is an intermediate generalist and discuss costs, benefits and constraints of host specialization in this species and brood parasites in general.

Visual opsin gene expression evolution in the adaptive radiation of cichlid fishes of Lake Tanganyika

Tuning the visual sensory system to the ambient light is essential for survival in many animal species. This is often achieved through duplication, functional diversification, and/or differential expression of visual opsin genes. Here, we examined 753 new retinal transcriptomes from 112 species of cichlid fishes from Lake Tanganyika to unravel adaptive changes in gene expression at the macro-evolutionary and ecosystem level of one of the largest vertebrate adaptive radiations. We found that, across the radiation, all seven cone opsins-but not the rhodopsin-rank among the most differentially expressed genes in the retina, together with other vision-, circadian rhythm-, and hemoglobin-related genes. We propose two visual palettes characteristic of very shallow- and deep-water living species, respectively, and show that visual system adaptations along two major ecological axes, macro-habitat and diet, occur primarily via gene expression variation in a subset of cone opsin genes.

No evidence for sex chromosomes in natural populations of the cichlid fish Astatotilapia burtoni

Sex determination (SD) is not conserved among teleost fishes and can even differ between populations of the same species. Across the outstandingly species-rich fish family Cichlidae, more and more SD systems are being discovered. Still, the picture of SD evolution in this group is far from being complete. Lake Tanganyika and its affluent rivers are home to Astatotilapia burtoni, which belongs to the extremely successful East African cichlid lineage Haplochromini. Previously, in different families of an A. burtoni laboratory strain, an XYW system and an XY system have been described. The latter was also found in a second laboratory strain. In a laboratory-reared family descending from a population of the species' southern distribution, a second XY system was discovered. Yet, an analysis of sex chromosomes for the whole species distribution is missing. Here, we examined the genomes of 11 natural populations of A. burtoni, encompassing a wide range of its distribution, for sex-linked regions. We did not detect signs of differentiated sex chromosomes and also not the previously described sex chromosomal systems present in laboratory lines, suggesting different SD systems in the same species under natural and (long-term) artificial conditions. We suggest that SD in A. burtoni is more labile than previously assumed and consists of a combination of non-genetic, polygenic, or poorly differentiated sex chromosomes.

Phylogenomics of trophically diverse cichlids disentangles processes driving adaptive radiation and repeated trophic transitions

Cichlid fishes of the tribe Tropheini are a striking case of adaptive radiation, exemplifying multiple trophic transitions between herbivory and carnivory occurring in sympatry with other established cichlid lineages. Tropheini evolved highly specialized eco-morphologies to exploit similar trophic niches in different ways repeatedly and rapidly. To better understand the evolutionary history and trophic adaptations of this lineage, we generated a dataset of 532 targeted loci from 21 out of the 22 described Tropheini species. We resolved the Tropheini into seven monophyletic genera and discovered one to be polyphyletic. The polyphyletic genus, Petrochromis, represents three convergent origins of the algae grazing trophic specialization. This repeated evolution of grazing may have been facilitated by adaptive introgression as we found evidence for gene flow among algae grazing genera. We also found evidence of gene flow among algae browsing genera, but gene flow was restricted between herbivorous and carnivorous genera. Furthermore, we observed no evidence supporting a hybrid origin of this radiation. Our molecular evolutionary analyses suggest that opsin genes likely evolved in response to selection pressures associated with trophic ecology in the Tropheini. We found surprisingly little evidence of positive selection in coding regions of jaw-shaping genes in this trophically diverse lineage. This suggests low degrees of freedom for further change in these genes, and possibly a larger role for regulatory variation in driving jaw adaptations. Our study emphasizes Tropheini cichlids as an important model for studying the evolution of trophic specialization and its role in speciation.

Somewhere I belong: phylogeny and morphological evolution in a species-rich lineage of ectoparasitic flatworms infecting cichlid fishes

A substantial portion of biodiversity has evolved through adaptive radiation. However, the effects of explosive speciation on species interactions remain poorly understood. Metazoan parasites infecting radiating host lineages could improve our knowledge because of their intimate host relationships. Yet limited molecular, phenotypic and ecological data discourage multivariate analyses of evolutionary patterns and encourage the use of discrete characters. Here, we assemble new molecular, morphological and host range data widely inferred from a species-rich lineage of parasites (Cichlidogyrus, Platyhelminthes: Monogenea) infecting cichlid fishes to address data scarcity. We infer a multimarker (28S/18S rDNA, ITS1, COI mtDNA) phylogeny of 58 of 137 species and characterize major lineages through synapomorphies inferred from mapping morphological characters. We predict the phylogenetic position of species without DNA data through shared character states, a morphological phylogenetic analysis, and a classification analysis with support vector machines. Based on these predictions and a cluster analysis, we assess the systematic informativeness of continuous characters, search for continuous equivalents for discrete characters, and suggest new characters for morphological traits not analysed to date. We also model the attachment/reproductive organ and host range evolution using the data for 136 of 137 described species and multivariate phylogenetic comparative methods (PCMs). We show that discrete characters not only can mask phylogenetic signals, but also are key for characterizing species groups. Regarding the attachment organ morphology, a divergent evolutionary regime for at least one lineage was detected and a limited morphological variation indicates host and environmental parameters affecting its evolution. However, moderate success in predicting phylogenetic positions, and a low systematic informativeness and high multicollinearity of morphological characters call for a revaluation of characters included in species characterizations.

Tracing evolutionary decoupling of oral and pharyngeal jaws in cichlid fishes

Evolutionary innovations can facilitate diversification if the novel trait enables a lineage to exploit new niches or by expanding character space. The elaborate pharyngeal jaw apparatus of cichlid fishes is often referred to as an evolutionary "key innovation" that has promoted the spectacular adaptive radiations in these fishes. This goes back to the idea that the structural and functional independence of the oral and pharyngeal jaws for food capturing and food processing, respectively, permitted each jaw type to follow independent evolutionary trajectories. This "evolutionary decoupling" is thought to have facilitated novel trait combinations and, hence, ecological specialization, ultimately allowing more species to coexist in sympatry. Here, we test the hypotheses of evolutionary decoupling of the oral and pharyngeal jaws in the massive adaptive radiation of cichlid fishes in African Lake Tanganyika. Based on phylogenetic comparative analyses of oral jaw morphology and lower pharyngeal jaw shape across most of the ∼240 cichlid species occurring in that lake, we show that the two jaws evolved coupled along the main axes of morphological variation, yet most other components of these trait complexes evolved largely independently over the course of the radiation. Further, we find limited correlations between the two jaws in both overall divergence and evolutionary rates. Moreover, we show that the two jaws were evolutionary decoupled at a late stage of the radiation, suggesting that decoupling contributed to micro-niche partitioning and the associated rapidly increasing trophic diversity during this phase.

Speciation dynamics and extent of parallel evolution along a lake-stream environmental contrast in African cichlid fishes

Understanding the dynamics of speciation is a central topic in evolutionary biology. Here, we investigated how morphological and genomic differentiation accumulated along the speciation continuum in the African cichlid fish Astatotilapia burtoni. While morphological differentiation was continuously distributed across different lake-stream population pairs, we found that there were two categories with respect to genomic differentiation, suggesting a “gray zone” of speciation at ~0.1% net nucleotide divergence. Genomic differentiation was increased in the presence of divergent selection and drift compared to drift alone. The quantification of phenotypic and genetic parallelism in four cichlid species occurring along a lake-stream environmental contrast revealed parallel and antiparallel components in rapid adaptive divergence, and morphological convergence in species replicates inhabiting the same environments. Furthermore, we show that the extent of parallelism was higher when ancestral populations were more similar. Our study highlights the complementary roles of divergent selection and drift on speciation and parallel evolution.

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.

Mapping of quantitative trait loci underlying a magic trait in ongoing ecological speciation

Background

Telmatochromis temporalis is a cichlid fish endemic to Lake Tanganyika. The normal and dwarf morphs of this fish are a clear example of ongoing ecological speciation, and body size plays an important role in this speciation event as a magic trait. However, the genetic basis underlying this trait has not been studied.

Results

Based on double-digested restriction-site associated DNA (ddRAD) sequencing of a hybrid cross between the morphs that includes F0 male, F0 female, and 206 F2 individuals, we obtained a linkage map consisting of 708 ddRAD markers in 22 linkage groups, which corresponded to the previously reported Oreochromis niloticus chromosomes, and identified one significant and five suggestive quantitative trait loci (QTL) for body size. From the body-size distribution pattern, the significant and three of the five suggestive QTL are possibly associated with genes responsible for the difference in body size between the morphs.

Conclusions

The QTL analysis presented here suggests that multiple genes, rather than a single gene, control morph-specific body size. The present results provide further insights about the genes underlying the morph specific body size and evolution of the magic trait during ecological speciation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07908-4.

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

The morphology of fish scales has been investigated for > 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.

Trophic specialisation reflected by radular tooth material properties in an "ancient" Lake Tanganyikan gastropod species flock

Background

Lake Tanganyika belongs to the East African Great Lakes and is well known for harbouring a high proportion of endemic and morphologically distinct genera, in cichlids but also in paludomid gastropods. With about 50 species these snails form a flock of high interest because of its diversity, the question of its origin and the evolutionary processes that might have resulted in its elevated amount of taxa. While earlier debates centred on these paludomids to be a result of an intralacustrine adaptive radiation, there are strong indications for the existence of several lineages before the lake formation. To evaluate hypotheses on the evolution and radiation the detection of actual adaptations is however crucial. Since the Tanganyikan gastropods show distinct radular tooth morphologies hypotheses about potential trophic specializations are at hand.

Results

Here, based on a phylogenetic tree of the paludomid species from Lake Tanganyika and adjacent river systems, the mechanical properties of their teeth were evaluated by nanoindentation, a method measuring the hardness and elasticity of a structure, and related with the gastropods’ specific feeding substrate (soft, solid, mixed). Results identify mechanical adaptations in the tooth cusps to the substrate and, with reference to the tooth morphology, assign distinct functions (scratching or gathering) to tooth types. Analysing pure tooth morphology does not consistently reflect ecological specializations, but the mechanical properties allow the determination of eco-morphotypes.

Conclusion

In almost every lineage we discovered adaptations to different substrates, leading to the hypothesis that one main engine of the flock’s evolution is trophic specialization, establishing distinct ecological niches and allowing the coexistence of taxa.

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.

Brood parasitism of an open-water spawning cichlid by the cuckoo catfish

The cuckoo catfish Synodontis multipunctatus and S. grandiops are endemic to Lake Tanganyika and the only known nonavian vertebrates that exhibit obligate interspecific brood parasitism. Seven maternal mouth-brooding cichlid fish species are reported to be natural hosts of the parasitic catfish and share a common reproductive behaviour that the catfish exploits: cichlid females spawn eggs on the bottom, allowing the catfish female to place her eggs near the cichlid eggs, and the cichlid females collect the catfish eggs by mouth together with their own eggs. However, so far it has not been reported that the cuckoo catfish exploit different spawning behaviours. The genus Cyprichromis consists of five maternal mouth-brooding species endemic to Lake Tanganyika, most of which spawn and collect eggs in open water. This study reports that the cuckoo catfish also parasitizes the open-water spawning Cyprichromis coloratus, although it may not be a regular host.