Phylogenomic analysis of Lake Malawi cichlid fishes: Further evidence that the three-stage model of diversification does not fit

African cichlid fishes: morphological data and taxonomic insights from a genus-level survey of supraneurals, pterygiophores, and vertebral counts (Ovalentaria, Blenniiformes, Cichlidae, Pseudocrenilabrinae)

Background

The iconic freshwater cichlid fishes (Cichlidae) comprise about 1750 validly named species and hundreds more that are known, but not yet described and named. Cichlids are an important source of protein for millions of people on several continents, are model organisms in studies of evolution, speciation, ecology, development, behaviour and physiology and are popular as aquarium fishes. Yet, comparative studies of cichlid internal anatomy are rare. Even their osteology has not been taxonomically surveyed. The cichlid postcranial skeleton has been especially neglected.

New information

Here, I provide the first survey in cichlids of the considerable variation in numbers of vertebrae, supraneurals and dorsal- and anal-fin supports (pterygiophores), as well as the patterns with which the pterygiophores insert between the neural or haemal spines. The study includes some 1700 specimens of nearly 400 cichlid species. Focusing on the largest subfamily, the African cichlids or Pseudocrenilabrinae, the survey furnishes data from species in all but one of its 166 genera. Limited data from species in the other cichlid subfamilies (Etroplinae, Ptychochrominae and Cichlinae) and from the related leaffishes, Polycentridae, are also presented. Key examples of pterygiophore insertion patterns from throughout the range of variation are illustrated and discussed. Detailed analytical tables and all raw data are provided in supplementary files.A bizarre specialisation in Cyprichromis is noted, evidently for the first time. Uniquely in this Lake Tanganyikan genus, five to seven anal pterygiophores are abdominal in position, located anterior to the anal fin and inserting toward or between successive pairs of pleural ribs.Taxonomic changes: The most speciose tribe of African cichlids, currently known as Haplochromini, is correctly called Pseudocrenilabrini. Based chiefly on the molecular phylogenetic findings of other workers, I propose four pseudocrenilabrine subtribes, one occurring in rivers and three endemic to Lake Malawi. I also re-assign the Lake Tanganyikan tribe Tropheini as another subtribe of Pseudocrenilabrini, in line with numerous molecular studies placing tropheines firmly within this tribe. The remaining genera of Pseudocrenilabrini remain incertae sedis in this tribe pending clarification of their phylogenetic relationships.The character complex here surveyed is a promising source of taxonomically and phylogenetically informative characteristics distinguishing or uniting cichlid taxa at multiple hierarchical levels, from species through subfamily. This reference set of novel character data can also provide information for palaeontological studies of African cichlids. These attributes are skeletal features potentially available for study in well preserved fossils and may help determine their correct taxonomic placement.

Preying on cyprinid snout warts (pearl organs) as a novel and peculiar habit in the Lake Malawi cichlid Docimodus evelynae

Cichlid fishes in the African Great Lakes have undergone explosive speciation, acquiring markedly varying ecologies and diets. There are multiple lineages of scale-eating cichlids, and their natural history and evolutionary ecology is only partially understood. We examined the feeding habit of Docimodus evelynae, a known scale eater, in Lake Malawi. The stomach contents of young individuals mainly consisted of unknown 1 mm hard, white warts (> 30%). To clarify the origin of these warts, we conducted an X-ray fluorometer analysis, and found they were rich in sulphur but low in silicon and calcium, suggesting they were epidermal tissues. Histological and morphological analyses revealed they were multicellular and cup-shaped. These characteristics matched only those of the pearl organs of the coexisting cyprinid Labeo cylindricus. DNA was extracted from the warts found in the stomach of five D. evelynae individuals, followed by PCR using primers targeting the partial COI gene of L. cylindricus. The resulting sequences exhibited 98% similarity to those of L. cylindricus. Pearl organs, never reported as a primary food for fish, could offer a substantial nutritional source based on calorific calculations. Understanding how this peculiar diet is foraged is essential for full comprehension of the food-web structure in this lake.

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.

Differences among reciprocal hybrids of Labeotropheus

Current evidence suggests that hybridization played a crucial role in the early evolution and diversification of the species flocks of cichlid fishes in the African Great Lakes. Nonetheless, evidence for hybridization in the extant cichlid fauna is scant, suggesting that hybridization is rare in the modern era, perhaps enforced by natural or sexual selection acting against F1 hybrids. Additionally, most experimental studies of hybridization perform a hybrid cross in one direction, ignoring the reciprocal hybrid. In this study, we perform reciprocal crosses between sympatric congeners from Lake Malaŵi, Labeotropheus fuelleborni and L. trewavasae, in order to compare the body shape and coloration of males of both of these hybrids, as well as to examine how these hybrids fare during both inter- and intrasexual interactions. We found that L. trewavasae-sired hybrid males are intermediate to the parental species both morphologically and chromatically, while the reciprocal L. fuelleborni-sired hybrids are likely transgressive hybrids. Males of these transgressive hybrids also fare poorly during our mate choice experiments. While female L. trewavasae reject them as possible mates, male L. trewavasae do not make a distinction between them and conspecific males. Selection against transgressive F1 hybrids as observed in our crossing experiments may help explain why contemporary hybridization in Lake Malaŵi cichlids appears to be rare.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10750-022-05092-4.

Phylogenomic analyses show repeated evolution of hypertrophied lips among Lake Malawi cichlid fishes

Cichlid fishes have repeatedly evolved an astounding diversity of trophic morphologies. For example, hypertrophied lips have evolved multiple times in both African and Neotropical cichlids and could have even evolved convergently within single species assemblages such as African Lake Malawi cichlids. However, the extremely high diversification rate in Lake Malawi cichlids and extensive potential for hybridization has cast doubt on whether even genome-level phylogenetic reconstructions could delineate if these types of adaptations have evolved once or multiple times. To examine the evolution of this iconic trait using protein-coding and noncoding single nucleotide polymorphisms (SNPs), we analyzed the genomes of 86 Lake Malawi cichlid species, including 33 de novo resequenced genomes. Surprisingly, genome-wide protein-coding SNPs exhibited enough phylogenetic informativeness to reconstruct interspecific and intraspecific relationships of hypertrophied lip cichlids, although noncoding SNPs provided better support. However, thinning of noncoding SNPs indicated most discrepancies come from the relatively smaller number of protein-coding sites and not from fundamental differences in their phylogenetic informativeness. Both coding and noncoding reconstructions showed that several “sand-dwelling” hypertrophied lip species, sampled intraspecifically, form a clade interspersed with a few other nonhypertrophied lip lineages. We also recovered Abactochromis labrosus within the rock-dwelling “mbuna” lineage, starkly contrasting with the affinities of other hypertrophied lip taxa found in the largely sand-dwelling “nonmbuna” component of this radiation. Comparative analyses coupled with tests for introgression indicate there is no widespread introgression between the hypertrophied lip lineages and taken together suggest this trophic phenotype has likely evolved at least twice independently within-lake Malawi.

The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of East African cichlid fishes

Color patterns are often linked to the behavioral and morphological characteristics of an animal, contributing to the effectiveness of such patterns as antipredatory strategies. Species-rich adaptive radiations, such as the freshwater fish family Cichlidae, provide an exciting opportunity to study trait correlations at a macroevolutionary scale. Cichlids are also well known for their diversity and repeated evolution of color patterns and body morphology. To study the evolutionary dynamics between color patterns and body morphology, we used an extensive dataset of 461 species. A phylogenetic supertree of these species shows that stripe patterns evolved ~70 times independently and were lost again ~30 times. Moreover, stripe patterns show strong signs of correlated evolution with body elongation, suggesting that the stripes' effectiveness as antipredatory strategy might differ depending on the body shape. Using pedigree-based analyses, we show that stripes and body elongation segregate independently, indicating that the two traits are not genetically linked. Their correlation in nature is therefore likely maintained by correlational selection. Lastly, by performing a mate preference assay using a striped CRISPR-Cas9 mutant of a nonstriped species, we show that females do not differentiate between striped CRISPR mutant males and nonstriped wild-type males, suggesting that these patterns might be less important for species recognition and mate choice. In summary, our study suggests that the massive rates of repeated evolution of stripe patterns are shaped by correlational selection with body elongation, but not by sexual selection.

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 Cautionary Note on the Use of Genotype Callers in Phylogenomics

Next-generation-sequencing genotype callers are commonly used in studies to call variants from newly sequenced species. However, due to the current availability of genomic resources, it is still common practice to use only one reference genome for a given genus, or even one reference for an entire clade of a higher taxon. The problem with traditional genotype callers, such as the one from GATK, is that they are optimized for variant calling at the population level. However, when these callers are used at the phylogenetic level, the consequences for downstream analyses can be substantial. Here, we performed simulations to compare the performance between the genotype callers of GATK and ATLAS, and present their differences at various phylogenetic scales. We show that the genotype caller of GATK substantially underestimates the number of variants at the phylogenetic level, but not at the population level. We also found that the accuracy of heterozygote calls declines with increasing distance to the reference genome. We quantified this decline and found that it is very sharp in GATK, while ATLAS maintains high accuracy even at moderately divergent species from the reference. We further suggest that efforts should be taken towards acquiring more reference genomes per species, before pursuing high-scale phylogenomic studies. [ATLAS; efficiency of SNP calling; GATK; heterozygote calling; next-generation sequencing; reference genome; variant calling.].

Neural innervation as a potential trigger of morphological color change and sexual dimorphism in cichlid fish

Many species change their coloration during ontogeny or even as adults. Color change hereby often serves as sexual or status signal. The cellular and subcellular changes that drive color change and how they are orchestrated have been barely understood, but a deeper knowledge of the underlying processes is important to our understanding of how such plastic changes develop and evolve. Here we studied the color change of the Malawi golden cichlid (Melanchromis auratus). Females and subordinate males of this species are yellow and white with two prominent black stripes (yellow morph; female and non-breeding male coloration), while dominant males change their color and completely invert this pattern with the yellow and white regions becoming black, and the black stripes becoming white to iridescent blue (dark morph; male breeding coloration). A comparison of the two morphs reveals that substantial changes across multiple levels of biological organization underlie this polyphenism. These include changes in pigment cell (chromatophore) number, intracellular dispersal of pigments, and tilting of reflective platelets (iridosomes) within iridophores. At the transcriptional level, we find differences in pigmentation gene expression between these two color morphs but, surprisingly, 80% of the genes overexpressed in the dark morph relate to neuronal processes including synapse formation. Nerve fiber staining confirms that scales of the dark morph are indeed innervated by 1.3 to 2 times more axonal fibers. Our results might suggest an instructive role of nervous innervation orchestrating the complex cellular and ultrastructural changes that drive the morphological color change of this cichlid species.

Gene markers for exon capture and phylogenomics in ray-finned fishes

Gene capture coupled with the next-generation sequencing has become one of the preferred methods of subsampling genomes for phylogenomic studies. Many exon markers have been developed in plants, sharks, frogs, reptiles, fishes, and others, but no universal exon markers have been tested in ray-finned fishes. Here, we identified a suite of "single-copy" protein-coding sequence (CDS) markers through comparing eight fish genomes, and tested them empirically in 83 species (33 families and nine orders or higher clades: Acipenseriformes, Lepisosteiformes, Elopomorpha, Osteoglossomorpha, Clupeiformes, Cypriniformes, Gobiaria, Carangaria, and Eupercaria; sensu Betancur et al. 2013). Sorting the markers according to their completeness and phylogenetic decisiveness in taxa tested resulted in a selection of 4,434 markers, which were proven to be useful in reconstructing phylogenies of the ray-finned fishes at different taxonomic levels. We also proposed a strategy of refining baits (probes) design a posteriori based on empirical data. The markers that we have developed may greatly enrich the batteries of exon markers for phylogenomic study in ray-finned fishes.

Pleiotropic jaw morphology links the evolution of mechanical modularity and functional feeding convergence in Lake Malawi cichlids

Complexity in how mechanistic variation translates into ecological novelty could be critical to organismal diversification. For instance, when multiple distinct morphologies can generate the same mechanical or functional phenotype, this could mitigate trade-offs and/or provide alternative ways to meet the same ecological challenge. To investigate how this type of complexity shapes diversity in a classic adaptive radiation, we tested several evolutionary consequences of the anterior jaw four-bar linkage for Lake Malawi cichlid trophic diversification. Using a novel phylogenetic framework, we demonstrated that different mechanical outputs of the same four jaw elements are evolutionarily associated with both jaw protrusion distance and jaw protrusion angle. However, these two functional aspects of jaw protrusion have evolved independently. Additionally, although four-bar morphology showed little evidence for attraction to optima, there was substantial evidence of adaptive peaks for emergent four-bar linkage mechanics and jaw protrusion abilities among Malawi feeding guilds. Finally, we highlighted a clear case of two cichlid species that have -independently evolved to graze algae in less than 2 Myr and have converged on similar jaw protrusion abilities as well as four-bar linkage mechanics, but have evolved these similarities via non-convergent four-bar morphologies.

Dissecting a potential spandrel of adaptive radiation: Body depth and pectoral fin ecomorphology coevolve in Lake Malawi cichlid fishes

The evolution of body shape reflects both the ecological factors structuring organismal diversity as well as an organism's underlying anatomy. For instance, body depth in fishes is thought to determine their susceptibility to predators, attractiveness to mates, as well as swimming performance. However, the internal anatomy influencing diversification of body depth has not been extensively examined, and changes in body depth could arise as a by-product of functional changes in other anatomical structures. Using an improved phylogenetic hypothesis for a diverse set of Lake Malawi cichlid fishes, we tested the evolutionary association between body depth and the height of the pectoral girdle. To refine the functional importance of the observed substantial correlation, we also tested the coevolution of pectoral girdle height and pectoral fin area. The extensive coevolution of these traits suggests body depth in fishes like the Lake Malawi cichlids could diverge simply as a by-product of being tightly linked to ecomorphological divergence in other functional morphological structures like the pectoral fins.

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

BACKGROUND

Phylogenies provide critical information about convergence during adaptive radiation. To test whether there have been multiple origins of a distinctive trophic phenotype in one of the most rapidly radiating groups known, we used ultra-conserved elements (UCEs) to examine the evolutionary affinities of Lake Malawi cichlids lineages exhibiting greatly hypertrophied lips.

RESULTS

The hypertrophied lip cichlids Cheilochromis euchilus, Eclectochromis ornatus, Placidochromis "Mbenji fatlip", and Placidochromis milomo are all nested within the non-mbuna clade of Malawi cichlids based on both concatenated sequence and single nucleotide polymorphism (SNP) inferred phylogenies. Lichnochromis acuticeps that exhibits slightly hypertrophied lips also appears to have evolutionary affinities to this group. However, Chilotilapia rhoadesii that lacks hypertrophied lips was recovered as nested within the species Cheilochromis euchilus. Species tree reconstructions and analyses of introgression provided largely ambiguous patterns of Malawi cichlid evolution.

CONCLUSIONS

Contrary to mitochondrial DNA phylogenies, bifurcating trees based on our 1024 UCE loci supported close affinities of Lake Malawi lineages with hypertrophied lips. However, incomplete lineage sorting in Malawi tends to render these inferences more tenuous. Phylogenomic analyses will continue to provide powerful inferences about whether phenotypic novelties arose once or multiple times during adaptive radiation.

Anterior-posterior gene expression differences in three Lake Malawi cichlid fishes with variation in body stripe orientation

Morphological differentiation among closely related species provides opportunities to study mechanisms shaping natural phenotypic variation. Here, we address variation in the orientation of melanin-colored body stripes in three cichlid species of the tribe Haplochromini. Melanochromis auratus displays a common pattern of dark, straight horizontal body stripes, whereas in Aristochromis christyi and Buccochromis rhoadesii, oblique stripes extend from the anterior dorsal to the posterior mid-lateral trunk. We first validated a stably reference gene, and then, investigated the chromatophore distribution in the skin by assessing the expression levels of the iridophore and melanophore marker genes, ltk and slc24a5, respectively, as well as pmel, a melanophore pigmentation marker gene. We found anterior-posterior differences in the expression levels of the three genes in the oblique-striped species. The higher anterior expression of ltk, indicates increased iridophore density in the anterior region, i.e., uneven horizontal distribution of iridophores, which coincides with the anterior dorsalization of melanophore stripe in these species. The obliqueness of the horizontal body stripes might be a result of distinct migratory or patterning abilities of melanophores in anterior and posterior stripe regions which could be reflected by variation in the expression of genes involved in melanophore patterning. To address this, we investigated anterior-posterior expression levels of a primary set of candidate target genes with known functions in melanophore migration and stripe patterning in the adult zebrafish, and their related gene regulatory network. Among these genes, those with differences in anterior-posterior expression showed only species-specific differential expression, e.g., sdf1a, col14a1a, ifitm5, and agpat3, with the exception of fbxw4/hagoromo (differentially expressed in an oblique-and the straight-striped species). In summary, distinct anterior-posterior gradients in iridophore density found to be more similar characteristic between the two oblique-striped species. Furthermore, the species-specific differential expression of genes involved in stripe patterning might also implicate distinct molecular processes underlying the obliqueness of body stripe in two closely related cichlid species.