Genetic architecture of a key reproductive isolation trait differs between sympatric and non-sympatric sister species of Lake Victoria cichlids

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.

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.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.

Genomic architecture of adaptive radiation and hybridization in Alpine whitefish

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

In this genomic study on Alpine whitefish radiations, the authors reveal details on the genetic architecture underlying the repeated eco-morphological diversification and the role of hybridization in the evolution of endemic whitefish species.

Genetic architecture of adaptive radiation across two trophic levels

Evolution of trophic diversity is a hallmark of adaptive radiation. Yet, transitions between carnivory and herbivory are rare in young adaptive radiations. Haplochromine cichlid fish of the African Great Lakes are exceptional in this regard. Lake Victoria was colonized by an insectivorous generalist and in less than 20 000 years, several clades of specialized herbivores evolved. Carnivorous versus herbivorous lifestyles in cichlids require many different adaptations in functional morphology, physiology and behaviour. Ecological transitions in either direction thus require many traits to change in a concerted fashion, which could be facilitated if genomic regions underlying these traits were physically linked or pleiotropic. However, linkage/pleiotropy could also constrain evolvability. To investigate components of the genetic architecture of a suite of traits that distinguish invertivores from algae scrapers, we performed quantitative trait locus (QTL) mapping using a second-generation hybrid cross. While we found indications of linkage/pleiotropy within trait complexes, QTLs for distinct traits were distributed across several unlinked genomic regions. Thus, a mixture of independently segregating variation and some pleiotropy may underpin the rapid trophic transitions. We argue that the emergence and maintenance of associations between the different genomic regions underpinning co-adapted traits that evolved and persist against some gene flow required reproductive isolation.

The genetic basis of a novel reproductive strategy in Sulawesi ricefishes: How modularity and a low number of loci shape pelvic brooding

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

Drainage-structuring of ancestral variation and a common functional pathway shape limited genomic convergence in natural high- and low-predation guppies

Studies of convergence in wild populations have been instrumental in understanding adaptation by providing strong evidence for natural selection. At the genetic level, we are beginning to appreciate that the re-use of the same genes in adaptation occurs through different mechanisms and can be constrained by underlying trait architectures and demographic characteristics of natural populations. Here, we explore these processes in naturally adapted high- (HP) and low-predation (LP) populations of the Trinidadian guppy, Poecilia reticulata. As a model for phenotypic change this system provided some of the earliest evidence of rapid and repeatable evolution in vertebrates; the genetic basis of which has yet to be studied at the whole-genome level. We collected whole-genome sequencing data from ten populations (176 individuals) representing five independent HP-LP river pairs across the three main drainages in Northern Trinidad. We evaluate population structure, uncovering several LP bottlenecks and variable between-river introgression that can lead to constraints on the sharing of adaptive variation between populations. Consequently, we found limited selection on common genes or loci across all drainages. Using a pathway type analysis, however, we find evidence of repeated selection on different genes involved in cadherin signaling. Finally, we found a large repeatedly selected haplotype on chromosome 20 in three rivers from the same drainage. Taken together, despite limited sharing of adaptive variation among rivers, we found evidence of convergent evolution associated with HP-LP environments in pathways across divergent drainages and at a previously unreported candidate haplotype within a drainage.

Identification of a novel sex determining chromosome in cichlid fishes that acts as XY or ZW in different lineages

Sex determination systems are highly conserved among most vertebrates with genetic sex determination, but can be variable and evolve rapidly in some. Here, we study sex determination in a clade with exceptionally high sex chromosome turnover rates. We identify the sex determining chromosomes in three interspecific crosses of haplochromine cichlid fishes from Lakes Victoria and Malawi. We find evidence for different sex determiners in each cross. In the Malawi cross and one Victoria cross the same chromosome is sex-linked but while females are the heterogametic sex in the Malawi species, males are the heterogametic sex in the Victoria species. This chromosome has not previously been reported to be sex determining in cichlids, increasing the number of different chromosomes shown to be sex determining in cichlids to 12. All Lake Victoria species of our crosses are less than 15,000 years divergent, and we identified different sex determiners among them. Our study provides further evidence for the diversity and evolutionary flexibility of sex determination in cichlids, factors which might contribute to their rapid adaptive radiations.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10750-021-04560-7.

Movement of transposable elements contributes to cichlid diversity

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Genetic architecture of a key reproductive isolation trait differs between sympatric and non-sympatric sister species of Lake Victoria cichlids

One hallmark of the East African cichlid radiations is the rapid evolution of reproductive isolation that is robust to full sympatry of many closely related species. Theory predicts that species persistence and speciation in sympatry with gene flow are facilitated if loci of large effect or physical linkage (or pleiotropy) underlie traits involved in reproductive isolation. Here, we investigate the genetic architecture of a key trait involved in behavioural isolation, male nuptial coloration, by crossing two sister species pairs of Lake Victoria cichlids of the genus Pundamilia and mapping nuptial coloration in the F2 hybrids. One is a young sympatric species pair, representative of an axis of colour motif differentiation, red-dorsum versus blue, that is highly recurrent in closely related sympatric species. The other is a species pair representative of colour motifs, red-chest versus blue, that are common in allopatric but uncommon in sympatric closely related species. We find significant quantitative trait loci (QTLs) with moderate to large effects (some overlapping) for red and yellow in the sympatric red-dorsum × blue cross, whereas we find no significant QTLs in the non-sympatric red-chest × blue cross. These findings are consistent with theory predicting that large effect loci or linkage/pleiotropy underlying mating trait differentiation could facilitate speciation and species persistence with gene flow in sympatry.