Understanding explosive diversification through cichlid fish genomics

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Emotional contagion and prosocial behaviour in fish: An evolutionary and mechanistic approach

In this review, we consider the definitions and experimental approaches to emotional contagion and prosocial behaviour in mammals and explore their evolutionary conceptualisation for studying their occurrence in the evolutionarily divergent vertebrate group of ray-finned fish. We present evidence for a diverse set of fish phenotypes that meet definitional criteria for prosocial behaviour and emotional contagion and discuss conserved mechanisms that may account for some preserved social capacities in fish. Finally, we provide some considerations on how to address the question of interdependency between emotional contagion and prosocial response, highlighting the importance of recognition processes, decision-making systems, and ecological context for providing evolutionary explanations.

Genetic and phenotypic diversification in a widespread fish, the Sailfin Molly (Poecilia latipinna)

Widespread species often experience significant environmental clines over the area they naturally occupy. We investigated a widespread livebearing fish, the Sailfin molly (Poecilia latipinna) combining genetic, life-history, and environmental data, asking how structured populations are. Sailfin mollies can be found in coastal freshwater and brackish habitats from roughly Tampico, Veracruz in Mexico to Wilmington, North Carolina, in the USA. In addition, they are found inland on the Florida peninsula. Using microsatellite DNA, we genotyped 168 individuals from 18 populations covering most of the natural range of the Sailfin molly. We further determined standard life-history parameters for both males and females for these populations. Finally, we measured biotic and abiotic parameters in the field. We found six distinct genetic clusters based on microsatellite data, with very strong indication of isolation by distance. However, we also found significant numbers of migrants between adjacent populations. Despite genetic structuring we did not find evidence of cryptic speciation. The genetic clusters and the migration patterns do not match paleodrainages. Life histories vary between populations but not in a way that is easy to interpret. We suggest a role of humans in migration in the sailfin molly, for example in the form of a ship channel that connects southern Texas with Louisiana which might be a conduit for fish migration.

The succession of ecological divergence and reproductive isolation in adaptive radiations

Adaptive radiation is a major source of biodiversity but the way in which known components of ecological opportunity, ecological differentiation, and reproductive isolation underpin such biodiversity patterns remains elusive. Much is known about the evolution of ecological differentiation and reproductive isolation during single speciation events, but exactly how those processes scale up to complete adaptive radiations is less understood. Do we expect complete reproductive barriers between newly formed species before the ecological differentiation continues, or does proper species formation occur much later, long after the ecological diversification? Our goal is to improve our mechanistic understanding of adaptive radiations by analyzing an individual-based model that includes a suite of mechanisms that are known to contribute to biodiversity. The model includes variable biogeographic settings, ecological opportunities, and types of mate choice, which makes several different scenarios of an adaptive radiation possible. We find that evolving clades tend to exploit ecological opportunities early whereas reproductive barriers evolve later, demonstrating a decoupling of ecological differentiation and species formation. In many cases, we also find a long-term trend where assortative mating associated with ecological traits is replaced by sexual selection of neutral display traits as the primary mechanism for reproductive isolation. Our results propose that reticulate phylogenies are likely common and stem from initially low reproductive barriers, rather than the previously suggested idea of repeated hybridization events between well-separated species.

Genetic and developmental divergence in the neural crest programme between cichlid fish species

Neural crest (NC) is a vertebrate-specific embryonic progenitor cell population at the basis of important vertebrate features such as the craniofacial skeleton and pigmentation patterns. Despite the wide-ranging variation of NC-derived traits across vertebrates, the contribution of NC to species diversification remains underexplored. Here, leveraging the adaptive diversity of African Great Lakes' cichlid species, we combined comparative transcriptomics and population genomics to investigate the evolution of the NC genetic programme in the context of their morphological divergence. Our analysis revealed substantial differences in transcriptional landscapes across somitogenesis, an embryonic period coinciding with NC development and migration. This included dozens of genes with described functions in the vertebrate NC gene regulatory network, several of which showed signatures of positive selection. Among candidates showing between-species expression divergence, we focused on teleost-specific paralogs of the NC-specifier sox10 (sox10a and sox10b) as prime candidates to influence NC development. These genes, expressed in NC cells, displayed remarkable spatio-temporal variation in cichlids, suggesting their contribution to inter-specific morphological differences. Finally, through CRISPR/Cas9 mutagenesis, we demonstrated the functional divergence between cichlid sox10 paralogs, with the acquisition of a novel skeletogenic function by sox10a. When compared to the teleost models zebrafish and medaka, our findings reveal that sox10 duplication, although retained in most teleost lineages, had variable functional fates across their phylogeny. Altogether, our study suggests that NC-related processes - particularly those controlled by sox10s - might be involved in generating morphological diversification between species and lays the groundwork for further investigations into mechanisms underpinning vertebrate NC diversification.

Spatial population genetic structure of Caquetaia kraussii (Steindachner, 1878) evidenced by species-specific microsatellite loci in the middle and low basin of the Cauca River, Colombia

The adaptative responses and divergent evolution shown in the environments habited by the Cichlidae family allow to understand different biological properties, including fish genetic diversity and structure studies. In a zone that has been historically submitted to different anthropogenic pressures, this study assessed the genetic diversity and population structure of cichlid Caquetaia kraussii, a sedentary species with parental care that has a significant ecological role for its contribution to redistribution and maintenance of sedimentologic processes in its distribution area. This study developed de novo 16 highly polymorphic species-specific microsatellite loci that allowed the estimation of the genetic diversity and differentiation in 319 individuals from natural populations in the area influenced by the Ituango hydroelectric project in the Colombian Cauca River. Caquetaia kraussii exhibits high genetic diversity levels (Ho: 0.562-0.885; He: 0.583-0.884) in relation to the average neotropical cichlids and a three group-spatial structure: two natural groups upstream and downstream the Nechí River mouth, and one group of individuals with high relatedness degree, possibly independently formed by founder effect in the dam zone. The three genetic groups show recent bottlenecks, but only the two natural groups have effective population size that suggest their long-term permanence. The information generated is relevant not only for management programs and species conservation purposes, but also for broadening the available knowledge on the factors influencing neotropical cichlids population genetics.

Ontogeny and social context regulate the circadian activity patterns of Lake Malawi cichlids

Activity patterns tend to be highly stereotyped and critical for executing many different behaviors including foraging, social interactions, and predator avoidance. Differences in the circadian timing of locomotor activity and rest periods can facilitate habitat partitioning and the exploitation of novel niches. As a consequence, closely related species often display highly divergent activity patterns, suggesting that shifts from diurnal to nocturnal behavior, or vice versa, are critical for survival. In Africa's Lake Malawi alone, there are over 500 species of cichlids, which inhabit diverse environments and exhibit extensive phenotypic variation. We have previously identified a substantial range in activity patterns across adult Lake Malawi cichlid species, from strongly diurnal to strongly nocturnal. In many species, including fishes, ecological pressures differ dramatically across life-history stages, raising the possibility that activity patterns may change over ontogeny. To determine if rest-activity patterns change across life stages, we compared the locomotor patterns of six Lake Malawi cichlid species. While total rest and activity did not change between early juvenile and adult stages, rest-activity patterns did, with juveniles displaying distinct activity rhythms that are more robust than adults. One distinct difference between juveniles and adults is the emergence of complex social behavior. To determine whether social context is required for activity rhythms, we next measured locomotor behavior in group-housed adult fish. We found that when normal social interactions were allowed, locomotor activity patterns were restored, supporting the notion that social interactions promote circadian regulation of activity in adult fish. These findings reveal a previously unidentified link between developmental stage and social interactions in the circadian timing of cichlid activity.

Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria

Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.

Developmental plasticity and variability in the formation of egg-spots, a pigmentation ornament in the cichlid Astatotilapia calliptera

Vertebrate pigmentation patterns are highly diverse, yet we have a limited understanding of how evolutionary changes to genetic, cellular, and developmental mechanisms generate variation. To address this, we examine the formation of a sexually-selected male ornament exhibiting inter- and intraspecific variation, the egg-spot pattern, consisting of circular yellow-orange markings on the male anal fins of haplochromine cichlid fishes. We focus on Astatotilapia calliptera, the ancestor-type species of the Malawi cichlid adaptive radiation of over 850 species. We identify a key role for iridophores in initializing egg-spot aggregations composed of iridophore-xanthophore associations. Despite adult sexual dimorphism, aggregations initially form in both males and females, with development only diverging between the sexes at later stages. Unexpectedly, we found that the timing of egg-spot initialization is plastic. The earlier individuals are socially isolated, the earlier the aggregations form, with iridophores being the cell type that responds to changes to the social environment. Furthermore, we observe apparent competitive interactions between adjacent egg-spot aggregations, which strongly suggests that egg-spot patterning results mostly from cell-autonomous cellular interactions. Together, these results demonstrate that A. calliptera egg-spot development is an exciting model for investigating pigment pattern formation at the cellular level in a system with developmental plasticity, sexual dimorphism, and intraspecific variation. As A. calliptera represents the ancestral bauplan for egg-spots, these findings provide a baseline for informed comparisons across the incredibly diverse Malawi cichlid radiation.

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.

Mitogenomes do not substantially improve phylogenetic resolution in a young non-model adaptive radiation of freshwater gastropods

BACKGROUND

Species flocks in ancient lakes, and particularly those arising from adaptive radiation, make up the bulk of overall taxonomic and morphological diversity in these insular ecosystems. For these mostly young species assemblages, classical mitochondrial barcoding markers have so far been key to disentangle interspecific relationships. However, with the rise and further development of next-generation sequencing (NGS) methods and mapping tools, genome-wide data have become an increasingly important source of information even for non-model groups.

RESULTS

Here, we provide, for the first time, a comprehensive mitogenome dataset of freshwater gastropods endemic to Sulawesi and thus of an ancient lake invertebrate species flock in general. We applied low-coverage whole-genome sequencing for a total of 78 individuals including 27 out of the 28 Tylomelania morphospecies from the Malili lake system as well as selected representatives from Lake Poso and adjacent catchments. Our aim was to assess whether mitogenomes considerably contribute to the phylogenetic resolution within this young species flock. Interestingly, we identified a high number of variable and parsimony-informative sites across the other 'non-traditional' mitochondrial loci. However, although the overall support was very high, the topology obtained was largely congruent with previously published single-locus phylogenies. Several clades remained unresolved and a large number of species was recovered polyphyletic, indicative of both rapid diversification and mitochondrial introgression.

CONCLUSIONS

This once again illustrates that, despite the higher number of characters available, mitogenomes behave like a single locus and thus can only make a limited contribution to resolving species boundaries, particularly when introgression events are involved.

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Dynamic co-evolution of transposable elements and the piRNA pathway in African cichlid fishes

East African cichlid fishes have diversified in an explosive fashion, but the (epi)genetic basis of the phenotypic diversity of these fishes remains largely unknown. Although transposable elements (TEs) have been associated with phenotypic variation in cichlids, little is known about their transcriptional activity and epigenetic silencing. Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed an expansion of piwil1 genes in Lake Malawi cichlids, likely driven by PiggyBac TEs. The expanded piwil1 copies have signatures of positive selection and retain amino acid residues essential for catalytic activity. Furthermore, the gonads of African cichlids express a Piwi-interacting RNA (piRNA) pathway that target TEs. We define the genomic sites of piRNA production in African cichlids and find divergence in closely related species, in line with fast evolution of piRNA-producing loci. Our findings suggest dynamic co-evolution of TEs and host silencing pathways in the African cichlid radiations. We propose that this co-evolution has contributed to cichlid genomic diversity.

Influence of ambient water coloration on habitat and conspecific choice in the female Lake Malawi cichlid, Metriaclima zebra

Female cichlid fish living in African great lakes are known to have sensory systems that are adapted to ambient light environments. These sensory system adaptations are hypothesized to have influenced the evolution of the diverse male nuptial coloration. In rock-dwelling Lake Malawi mbuna cichlids, however, the extent to which ambient light environments influence female sensory systems and potentially associated male nuptial coloration remains unknown. Yet, the ubiquitous blue flank coloration and UV reflection of male mbuna cichlids suggest the potential impacts of the blue-shifted ambient light environment on these cichlid's visual perception and male nuptial coloration in the shallow water depth in Lake Malawi. In the present study, we explored whether and how the sensory bias of females influences intersexual communication in the mbuna cichlid, Metriaclima zebra. A series of choice experiments in various light environments showed that M. zebra females 1) have a preference for the blue-shifted light environment, 2) prefer to interact with males in blue-shifted light environments, 3) do not show a preference between dominant and subordinate males in full-spectrum, long-wavelength filtered, and short-wavelength filtered light environments, and 4) show a "reversed" preference for subordinate males in the UV-filtered light environment. These results suggest that the visual perception of M. zebra females may be biased to the ambient light spectra in their natural habitat by local adaptation and that this sensory bias may influence the evolution of blue and UV reflective patterns in male nuptial coloration.

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.

Epigenetics and the evolution of form: Experimental manipulation of a chromatin modification causes species-specific changes to the craniofacial skeleton

A central question in biology is the molecular origins of phenotypic diversity. While genetic changes are key to the genotype-phenotype relationship, alterations to chromatin structure and the physical packaging of histone proteins may also be important drivers of vertebrate divergence. We investigate the impact of such an epigenetic mechanism, histone acetylation, within a textbook example of an adaptive radiation. Cichlids of Lake Malawi have adapted diverse craniofacial structures, and here we investigate how histone acetylation influences morphological variation in these fishes. Specifically, we assessed the effect of inhibiting histone deacetylation using the drug trichostatin A (TSA) on developing facial structures. We examined this during three critical developmental windows in two cichlid species with alternate adult morphologies. Exposure to TSA during neural crest cell (NCC) migration and as postmigratory NCCs proliferate in the pharyngeal arches resulted in significant changes in lateral and ventral shape in Maylandia, but not in Tropheops. This included an overall shortening of the head, widening of the lower jaw, and steeper craniofacial profile, all of which are paedomorphic morphologies. In contrast, treatment with TSA during early chondrogenesis did not result in significant morphological changes in either species. Together, these data suggest a sensitivity to epigenetic alterations that are both time- and species-dependent. We find that morphologies are due to nonautonomous or potentially indirect effects on NCC development, including in part a global developmental delay. Our research bolsters the understanding that proper histone acetylation is essential for early craniofacial development and identifies a species-specific robustness to developmental change. Overall, this study demonstrates how epigenetic regulation may play an important role in both generating and buffering morphological variation.

Developmental and environmental plasticity in opsin gene expression in Lake Victoria cichlid fish

In many organisms, sensory abilities develop and evolve according to the changing demands of navigating, foraging, and communication across different environments and life stages. Teleost fish inhabit heterogeneous light environments and exhibit a large diversity in visual system properties among species. Cichlids are a classic example of this diversity; visual system variation is generated by different tuning mechanisms that involve both genetic factors and phenotypic plasticity. Here, we document the developmental progression of visual pigment gene expression in Lake Victoria cichlids and test if these patterns are influenced by variation in light conditions. We reared two sister species of Pundamilia to adulthood in two distinct visual conditions that resemble the light environments that they naturally inhabit in Lake Victoria. We also included interspecific first-generation hybrids. We focused on the four opsins that are expressed in Pundamilia adults (using real-time quantitative polymerase chain reaction (RT-qPCR)) (SWS2B, SWS2A, RH2A, and LWS) at 17 time points. We find that opsin expression profiles progress from shorter-wavelength sensitive opsins to longer-wavelength sensitive opsins with increasing age, in both species and their hybrids. The developmental trajectories of opsin expression also responded plastically to the visual conditions. Developmental and environmental plasticity in opsin expression may provide an important stepping stone in the evolution of cichlid visual system diversity.

A protocol for whole-brain Ca2+ imaging in Astyanax mexicanus, a model of comparative evolution

In this protocol, we describe a comparative approach to study the evolution of brain function in the Mexican tetra, Astyanax mexicanus. We developed surface fish and two independent populations of cavefish with pan-neuronal expression of the Ca2+ sensor GCaMP6s. We describe a methodology to prepare samples and image activity across the optic tectum and olfactory bulb.

Differential Expression of miRNAs, lncRNAs, and circRNAs between Ovaries and Testes in Common Carp (Cyprinus carpio)

Female common carp grow faster than male individuals, implying that rearing females could be more profitable in aquaculture. Non-coding RNAs (ncRNAs) serve as versatile regulators with multiple functions in diverse biological processes. However, the roles of ncRNAs in the sex differentiation of common carp are less studied. In this study, we investigated the expression profiles of ncRNAs, including miRNAs, lncRNAs, and circRNAs, in the gonads to comprehend the roles of ncRNAs in sex differentiation in common carp. A substantial number of differentially expressed (DE) ncRNAs in ovaries and testes were identified. Some miRNAs, notably miR-205, miR-214, and miR-460-5p, might modulate hormone synthesis and thus maintain sex. A novel miRNA, novel_158, was predicted to suppress the expression of foxl3. DE lncRNAs were associated with oocyte meiosis, GnRH signaling pathways, and steroid biosynthesis, while DE circRNA target genes were enriched in the WNT signaling pathway and MAPK signaling pathway. We also analyzed ncRNA-mRNA interactions to shed light on the crosstalk between competing endogenous RNAs (ceRNAs), which is the critical mechanism by which lncRNAs and circRNAs function. Some lncRNAs and circRNAs may be able to competitively bind novel_313, a new miRNA, and thus regulate hsd17β3. Our research will provide a valuable resource for understanding the genetic basis of gonadal differentiation and development in common carp.

Phenotypic rate and state are decoupled in response to river-to-lake transitions in cichlid fishes

Geographic access to isolated ecosystems is an important catalyst of adaptive radiation. Cichlid fishes repeatedly colonized rift, crater, and volcanic lakes from surrounding rivers. We test the "lake effect" on the phenotypic rate and state across 253 cichlid species. The rate of evolution was consistently higher (~10-fold) in lakes, and consistent across different dimensions of the phenotype. Rate shifts tended to occur coincident with or immediately following river-to-lake transitions, generally resulting in 2- to 5-fold faster rates than in the founding riverine lineage. By contrast, river- and lake-dwelling cichlids exhibit considerable overlap in phenotypes, generally with less disparity in lakes, but often different evolutionary optima. Taken together, these results suggest that lake radiations rapidly expand into niches largely already represented by ancestral riverine lineages, albeit in different frequencies. Lakes may provide ecological opportunity via ecological release (e.g., from predators/competitors) but need not be coupled with access to novel ecological niches.

Genome editing in East African cichlids and tilapias: state-of-the-art and future directions

African cichlid fishes of the Cichlidae family are a group of teleosts important for aquaculture and research. A thriving research community is particularly interested in the cichlid radiations of the East African Great Lakes. One key goal is to pinpoint genetic variation underlying phenotypic diversification, but the lack of genetic tools has precluded thorough dissection of the genetic basis of relevant traits in cichlids. Genome editing technologies are well established in teleost models like zebrafish and medaka. However, this is not the case for emerging model organisms, such as East African cichlids, where these technologies remain inaccessible to most laboratories, due in part to limited exchange of knowledge and expertise. The Cichlid Science 2022 meeting (Cambridge, UK) hosted for the first time a Genome Editing Workshop, where the community discussed recent advances in genome editing, with an emphasis on CRISPR/Cas9 technologies. Based on the workshop findings and discussions, in this review we define the state-of-the-art of cichlid genome editing, share resources and protocols, and propose new possible avenues to further expand the cichlid genome editing toolkit.

Teleost genomic repeat landscapes in light of diversification rates and ecology

Repetitive DNA make up a considerable fraction of most eukaryotic genomes. In fish, transposable element (TE) activity has coincided with rapid species diversification. Here, we annotated the repetitive content in 100 genome assemblies, covering the major branches of the diverse lineage of teleost fish. We investigated if TE content correlates with family level net diversification rates and found support for a weak negative correlation. Further, we demonstrated that TE proportion correlates with genome size, but not to the proportion of short tandem repeats (STRs), which implies independent evolutionary paths. Marine and freshwater fish had large differences in STR content, with the most extreme propagation detected in the genomes of codfish species and Atlantic herring. Such a high density of STRs is likely to increase the mutational load, which we propose could be counterbalanced by high fecundity as seen in codfishes and herring.

Genes, Morphology, Performance, and Fitness: Quantifying Organismal Performance to Understand Adaptive Evolution

To understand the complexities of morphological evolution, we must understand the relationships between genes, morphology, performance, and fitness in complex traits. Genomicists have made tremendous progress in finding the genetic basis of many phenotypes, including a myriad of morphological characters. Similarly, field biologists have greatly advanced our understanding of the relationship between performance and fitness in natural populations. However, the connection from morphology to performance has primarily been studied at the interspecific level, meaning that in most cases we lack a mechanistic understanding of how evolutionarily relevant variation among individuals affects organismal performance. Therefore, functional morphologists need methods that will allow for the analysis of fine-grained intraspecific variation in order to close the path from genes to fitness. We suggest three methodological areas that we believe are well suited for this research program and provide examples of how each can be applied within fish model systems to build our understanding of microevolutionary processes. Specifically, we believe that structural equation modeling, biological robotics, and simultaneous multi-modal functional data acquisition will open up fruitful collaborations among biomechanists, evolutionary biologists, and field biologists. It is only through the combined efforts of all three fields that we will understand the connection between evolution (acting at the level of genes) and natural selection (acting on fitness).

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.

Support for a radiation of free-living flatworms in the African Great Lakes region and the description of five new Macrostomum species

BACKGROUND

The African Great Lakes have long been recognized as an excellent location to study speciation. Most famously, cichlid fishes have radiated in Lake Tanganyika and subsequently spread into Lake Malawi and Lake Victoria, where they again radiated. Other taxa have diversified in these lakes, such as catfish, ostracods, gastropods, and Monegenean gill parasites of cichlids. However, these radiations have received less attention, and the process leading to their speciation in this unique region remains to be further explored. Here I present evidence that suggests a radiation of Macrostomum flatworms has occurred in the African Great Lakes region, offering a good opportunity for such investigations.

RESULTS

Recent field work has revealed a monophyletic clade of 16 Macrostomum flatworms that have, to date, only been collected from Lake Tanganyika. Additionally, a species collected from Lake Malawi was found nested within this clade. Molecular phylogenetic analysis, largely based on transcriptome data, suggests that this clade underwent rapid speciation, possibly due to a large habitat diversity in the lake. I also observed significant differences in the sperm morphology of these flatworms compared to those of species found outside Lake Tanganyika and Lake Malawi. These included the elongation of an anterior structure, a reduction in the size of the lateral sperm bristles, and changes in relative proportions. I propose functional hypotheses for these changes in sperm design, and formally describe Macrostomum gracilistylum sp. nov from Lake Malawi and its sister species Macrostomum crassum sp. nov., Macrostomum pellitum sp. nov., Macrostomum longispermatum sp. nov., and Macrostomum schäreri sp. nov., from Lake Tanganyika.

CONCLUSIONS

The available evidence is consistent with the hypothesis that Macrostomum flatworms have radiated in Lake Tanganyika and subsequently spread to Lake Malawi. However, whether this represents a bona fide adaptive radiation still needs to be determined. Therefore, the African Great Lakes are promising targets for further research into flatworm diversity and speciation.

Ancient standing genetic variation facilitated the adaptive radiation of Lake Victoria cichlids

Cichlid fishes are textbook examples of explosive speciation and adaptive radiation, providing a great opportunity to understand how the genomic substrate yields extraordinary species diversity. Recently, we performed comparative genomic analyses of three Lake Victoria cichlids to reveal the genomic substrates underlying their rapid speciation and adaptation. We found that long divergent haplotypes derived from large-scale standing genetic variation, which originated before the adaptive radiation of Lake Victoria cichlids, may have contributed to their rapid diversification. In addition, the present study on genomic data from other East African cichlids suggested the reuse of alleles that may have originated in the ancestral lineages of Lake Tanganyika cichlids during cichlid evolution. Therefore, our results highlight that the primary factor that could drive repeated adaptive radiation across East African cichlids was allelic reuse from standing genetic variation to adapt to their own specific environment. In this report, we summarize the main results and discuss the evolutionary mechanisms of cichlids, based on our latest findings.

Multiple paths toward repeated phenotypic evolution in the spiny-leg adaptive radiation (Tetragnatha; Hawai'i)

The repeated evolution of phenotypes provides clear evidence for the role of natural selection in driving evolutionary change. However, the evolutionary origin of repeated phenotypes can be difficult to disentangle as it can arise from a combination of factors such as gene flow, shared ancestral polymorphisms or mutation. Here, we investigate the presence of these evolutionary processes in the Hawaiian spiny-leg Tetragnatha adaptive radiation, which includes four microhabitat-specialists or ecomorphs, with different body pigmentation and size (Green, Large Brown, Maroon, and Small Brown). We investigated the evolutionary history of this radiation using 76 newly generated low-coverage, whole-genome resequenced samples, along with phylogenetic and population genomic tools. Considering the Green ecomorph as the ancestral state, our results suggest that the Green ecomorph likely re-evolved once, the Large Brown and Maroon ecomorphs evolved twice and the Small Brown evolved three times. We found that the evolution of the Maroon and Small Brown ecomorphs likely involved ancestral hybridization events, while the Green and Large Brown ecomorphs likely evolved through novel mutations, despite a high rate of incomplete lineage sorting in the dataset. Our findings demonstrate that the repeated evolution of ecomorphs in the Hawaiian spiny-leg Tetragnatha is influenced by multiple evolutionary processes.

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.

Lifespans of Fur Color Morphs in Polymorphic Populations of the Mole Vole and the Hypothesis of Adaptive Polymorphism

Different lifespans were for the first time demonstrated for three (brown, bicolor, and black) fur color morphs in ten mole vole populations of the Volga, Ural, and Trans-Ural regions. With the longest lifespan of 5 years in the species, morphs that numerically dominate in a population can live 1-4 years longer than accompanying morphs. Spearman's correlation coefficient between the longest lifespan of the morphs and their proportion in the population was Rsp = 0.81 (p < 0.0001). A number of morphological and functional features were identified in the color morphs. The findings are of general biological significance, confirming the hypothesis of adaptive polymorphism. Evolutionary and ecological mechanisms whereby selective advantages develop in morphs (as probable ecomorphs) are possible to evaluate using the morphs as a natural model of the initial step of sympatric form development in different parts of the range.

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.

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.

Morphological and temporal variation in early embryogenesis contributes to species divergence in Malawi cichlid fishes

The cichlid fishes comprise the largest extant vertebrate family and are the quintessential example of rapid "explosive" adaptive radiations and phenotypic diversification. Despite low genetic divergence, East African cichlids harbor a spectacular intra- and interspecific morphological diversity, including the hyper-variable, neural crest (NC)-derived traits such as coloration and craniofacial skeleton. Although the genetic and developmental basis of these phenotypes has been investigated, understanding of when, and specifically how early, in ontogeny species-specific differences emerge, remains limited. Since adult traits often originate during embryonic development, the processes of embryogenesis could serve as a potential source of species-specific variation. Consequently, we designed a staging system by which we compare the features of embryogenesis between three Malawi cichlid species-Astatotilapia calliptera, Tropheops sp. 'mauve' and Rhamphochromis sp. "chilingali"-representing a wide spectrum of variation in pigmentation and craniofacial morphologies. Our results showed fundamental differences in multiple aspects of embryogenesis that could underlie interspecific divergence in adult adaptive traits. First, we identified variation in the somite number and signatures of temporal variation, or heterochrony, in the rates of somite formation. The heterochrony was also evident within and between species throughout ontogeny, up to the juvenile stages. Finally, the identified interspecific differences in the development of pigmentation and craniofacial cartilages, present at the earliest stages of their overt formation, provide compelling evidence that the species-specific trajectories begin divergence during early embryogenesis, potentially during somitogenesis and NC development. Altogether, our results expand our understanding of fundamental cichlid biology and provide new insights into the developmental origins of vertebrate morphological diversity.

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.

Mitogenomic Features and Evolution of the Nile River Dominant Tilapiine Species (Perciformes: Cichlidae)

To better understand the diversity and evolution of cichlids, we sequenced, assembled, and annotated the complete mitochondrial genomes of three Nile tilapiine species (Coptodon zillii, Oreochromis niloticus, and Sarotherodon galilaeus) dominating the Nile River waters. Our results showed that the general mitogenomic features were conserved among the Nile tilapiine species. The genome length ranged from 16,436 to 16,631 bp and a total of 37 genes were identified (two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and 1 control region). The ND6 was the only CDS that presented a negative AT skew and a positive GC skew. The most extended repeat sequences were in the D-loop followed by the pseudogenes (trnS^GCU). The ND5 showed relatively high substitution rates whereas ATP8 had the lowest substitution rate. The codon usage bias displayed a greater quantity of NNA and NNC at the third position and anti-bias against NNG. The phylogenetic relationship based on the complete mitogenomes and CDS was able to differentiate the three species as previously reported. This study provides new insight into the evolutionary connections between various subfamilies within cichlids while providing new molecular data that can be applied to discriminate between Nile tilapiine species and their populations.

Epigenetic divergence during early stages of speciation in an African crater lake cichlid fish

Epigenetic variation can alter transcription and promote phenotypic divergence between populations facing different environmental challenges. Here, we assess the epigenetic basis of diversification during the early stages of speciation. Specifically, we focus on the extent and functional relevance of DNA methylome divergence in the very young radiation of Astatotilapia calliptera in crater Lake Masoko, southern Tanzania. Our study focuses on two lake ecomorphs that diverged approximately 1,000 years ago and a population in the nearby river from which they separated approximately 10,000 years ago. The two lake ecomorphs show no fixed genetic differentiation, yet are characterized by different morphologies, depth preferences and diets. We report extensive genome-wide methylome divergence between the two lake ecomorphs, and between the lake and river populations, linked to key biological processes and associated with altered transcriptional activity of ecologically relevant genes. Such genes differing between lake ecomorphs include those involved in steroid metabolism, hemoglobin composition and erythropoiesis, consistent with their divergent habitat occupancy. Using a common-garden experiment, we found that global methylation profiles are often rapidly remodeled across generations but ecomorph-specific differences can be inherited. Collectively, our study suggests an epigenetic contribution to the early stages of vertebrate speciation.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Adaptive phenotypic plasticity is under stabilizing selection in Daphnia

The adaptive nature of phenotypic plasticity is widely documented. However, little is known about the evolutionary forces that shape genetic variation of plasticity within populations. Whether genetic variation in plasticity is driven by stabilizing or diversifying selection and whether the strength of such forces remains constant through time, remain open questions. Here, we address this issue by assessing the evolutionary forces that shape genetic variation in antipredator developmental plasticity of Daphnia pulex. Antipredator plasticity in D. pulex is characterized by the growth of a pedestal and spikes in the dorsal head region upon exposure to predator cue. We characterized genetic variation in plasticity using a method that describes the entire dorsal shape amongst >100 D. pulex strains recently derived from the wild. We observed the strongest reduction in genetic variation in dorsal areas where plastic responses were greatest, consistent with stabilizing selection. We compared mutational variation (Vm) to standing variation (Vg) and found that Vg/Vm is lowest in areas of greatest plasticity, again consistent with stabilizing selection. Our results suggest that stabilizing selection operates directly on phenotypic plasticity in Daphnia and provide a rare glimpse into the evolution of fitness-related traits in natural populations.

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.

Genome-wide species delimitation analyses of a silverside fish species complex in central Mexico indicate taxonomic over-splitting

BACKGROUND

Delimiting species across a speciation continuum is a complex task, as the process of species origin is not generally instantaneous. The use of genome-wide data provides unprecedented resolution to address convoluted species delimitation cases, often unraveling cryptic diversity. However, because genome-wide approaches based on the multispecies coalescent model are known to confound population structure with species boundaries, often resulting in taxonomic over-splitting, it has become increasingly evident that species delimitation research must consider multiple lines of evidence. In this study, we used phylogenomic, population genomic, and coalescent-based species delimitation approaches, and examined those in light of morphological and ecological information, to investigate species numbers and boundaries comprising the Chirostoma "humboltianum group" (family Atherinidae). The humboltianum group is a taxonomically controversial species complex where previous morphological and mitochondrial studies produced conflicting species delimitation outcomes. We generated ddRADseq data for 77 individuals representing the nine nominal species in the group, spanning their distribution range in the central Mexican plateau.

RESULTS

Our results conflict with the morphospecies and ecological delimitation hypotheses, identifying four independently evolving lineages organized in three geographically cohesive clades: (i) chapalae and sphyraena groups in Lake Chapala, (ii) estor group in Lakes Pátzcuaro and Zirahuén, and (iii) humboltianum sensu stricto group in Lake Zacapu and Lerma river system.

CONCLUSIONS

Overall, our study provides an atypical example where genome-wide analyses delineate fewer species than previously recognized on the basis of morphology. It also highlights the influence of the geological history of the Chapala-Lerma hydrological system in driving allopatric speciation in the humboltianum group.

Characterizing the genetic basis of trait evolution in the Mexican cavefish

Evolution in response to a change in ecology often coincides with various morphological, physiological, and behavioral traits. For most organisms little is known about the genetic and functional relationship between evolutionarily derived traits, representing a critical gap in our understanding of adaptation. The Mexican tetra, Astyanax mexicanus, consists of largely independent populations of fish that inhabit at least 30 caves in Northeast Mexico, and a surface fish population, that inhabit the rivers of Mexico and Southern Texas. The recent application of molecular genetic approaches combined with behavioral phenotyping have established A. mexicanus as a model for studying the evolution of complex traits. Cave populations of A. mexicanus are interfertile with surface populations and have evolved numerous traits including eye degeneration, insomnia, albinism, and enhanced mechanosensory function. The interfertility of different populations from the same species provides a unique opportunity to define the genetic relationship between evolved traits and assess the co-evolution of behavioral and morphological traits with one another. To define the relationships between morphological and behavioral traits, we developed a pipeline to test individual fish for multiple traits. This pipeline confirmed differences in locomotor activity, prey capture, and startle reflex between surface and cavefish populations. To measure the relationship between traits, individual F2 hybrid fish were characterized for locomotor behavior, prey-capture behavior, startle reflex, and morphological attributes. Analysis revealed an association between body length and slower escape reflex, suggesting a trade-off between increased size and predator avoidance in cavefish. Overall, there were few associations between individual behavioral traits, or behavioral and morphological traits, suggesting independent genetic changes underlie the evolution of the measured behavioral and morphological traits. Taken together, this approach provides a novel system to identify genetic underpinnings of naturally occurring variation in morphological and behavioral traits.

An atlas of fish genome evolution reveals delayed rediploidization following the teleost whole-genome duplication

Teleost fishes are ancient tetraploids descended from an ancestral whole-genome duplication that may have contributed to the impressive diversification of this clade. Whole-genome duplications can occur via self-doubling (autopolyploidy) or via hybridization between different species (allopolyploidy). The mode of tetraploidization conditions evolutionary processes by which duplicated genomes return to diploid meiotic pairing, and subsequent genetic divergence of duplicated genes (cytological and genetic rediploidization). How teleosts became tetraploid remains unresolved, leaving a fundamental gap in the interpretation of their functional evolution. As a result of the whole-genome duplication, identifying orthologous and paralogous genomic regions across teleosts is challenging, hindering genome-wide investigations into their polyploid history. Here, we combine tailored gene phylogeny methodology together with a state-of-the-art ancestral karyotype reconstruction to establish the first high-resolution comparative atlas of paleopolyploid regions across 74 teleost genomes. We then leverage this atlas to investigate how rediploidization occurred in teleosts at the genome-wide level. We uncover that some duplicated regions maintained tetraploidy for more than 60 million years, with three chromosome pairs diverging genetically only after the separation of major teleost families. This evidence suggests that the teleost ancestor was an autopolyploid. Further, we find evidence for biased gene retention along several duplicated chromosomes, contradicting current paradigms that asymmetrical evolution is specific to allopolyploids. Altogether, our results offer novel insights into genome evolutionary dynamics following ancient polyploidizations in vertebrates.

ceRNA crosstalk mediated by ncRNAs is a novel regulatory mechanism in fish sex determination and differentiation

Competing endogenous RNAs (ceRNAs) are vital regulators of gene networks in mammals. The involvement of noncoding RNAs (ncRNAs) as ceRNA in genotypic sex determination (GSD) and environmental sex determination (ESD) in fish is unknown. The Chinese tongue sole, which has both GSD and ESD mechanisms, was used to map the dynamic expression pattern of ncRNAs and mRNA in gonads during sex determination and differentiation. Transcript expression patterns shift during the sex differentiation phase, and ceRNA modulation occurs through crosstalk of differentially expressed long ncRNAs (lncRNAs), circular RNAs (circRNAs), microRNAs (miRNAs), and sex-related genes in fish. Of note was the significant up-regulation of a circRNA from the sex-determining gene dmrt1 (circular RNA dmrt1) and a lncRNA, called AMSDT (which stands for associated with male sex differentiation of tongue sole) in Chinese tongue sole testis. These two ncRNAs both share the same miRNA response elements with gsdf, which has an up-regulated expression when they bind to miRNA cse-miR-196 and concurrent down-regulated female sex-related genes to facilitate testis differentiation. This is the first demonstration in fish that ceRNA crosstalk mediated by ncRNAs modulates sexual development and unveils a novel regulatory mechanism for sex determination and differentiation.

Allopatric origin of sympatric whitefish morphs with insights on the genetic basis of their reproductive isolation

The European whitefish (Coregonus lavaretus) species complex is a classic example of recent adaptive radiation. Here, we examine a whitefish population introduced to northern Finnish Lake Tsahkal in the late 1960s, where three divergent morphs (viz. littoral, pelagic, and profundal feeders) were found 10 generations after. Using demographic modeling based on genomic data, we show that whitefish morphs evolved during a phase of strict isolation, refuting a rapid sympatric divergence scenario. The lake is now an artificial hybrid zone between morphs originated in allopatry. Despite their current syntopy, clear genetic differentiation remains between two of the three morphs. Using admixture mapping, we identify five SNPs associated with gonad weight variation, a proxy for sexual maturity and spawning time. We suggest that ecological adaptations in spawning time evolved in allopatry are currently maintaining partial reproductive isolation in the absence of other barriers to gene flow.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

Evolution of miRNA binding sites and regulatory networks in cichlids

The divergence of regulatory regions and gene regulatory network (GRN) rewiring is a key driver of cichlid phenotypic diversity. However, the contribution of miRNA binding site turnover has yet to be linked to GRN evolution across cichlids. Here, we extend our previous studies by analysing the selective constraints driving evolution of miRNA and transcription factor (TF) binding sites of target genes, to infer instances of cichlid GRN rewiring associated with regulatory binding site turnover. Comparative analyses identified increased species-specific networks that are functionally associated to traits of cichlid phenotypic diversity. The evolutionary rewiring is associated with differential models of miRNA and TF binding site turnover, driven by a high proportion of fast-evolving polymorphic sites in adaptive trait genes compared to subsets of random genes. Positive selection acting upon discrete mutations in these regulatory regions is likely to be an important mechanism in rewiring GRNs in rapidly radiating cichlids. Regulatory variants of functionally associated miRNA and TF binding sites of visual opsin genes differentially segregate according to phylogeny and ecology of Lake Malawi species, identifying both rewired e.g. clade-specific and conserved network motifs of adaptive trait associated GRNs. Our approach revealed several novel candidate regulators, regulatory regions and three-node motifs across cichlid genomes with previously reported associations to known adaptive evolutionary traits.

Evolutionary Divergence and Radula Diversification in Two Ecomorphs from an Adaptive Radiation of Freshwater Snails

(1) Background: Adaptive diversification of complex traits plays a pivotal role in the evolution of organismal diversity. In the freshwater snail genus Tylomelania, adaptive radiations were likely promoted by trophic specialization via diversification of their key foraging organ, the radula. (2) Methods: To investigate the molecular basis of radula diversification and its contribution to lineage divergence, we used tissue-specific transcriptomes of two sympatric Tylomelania sarasinorum ecomorphs. (3) Results: We show that ecomorphs are genetically divergent lineages with habitat-correlated abundances. Sequence divergence and the proportion of highly differentially expressed genes are significantly higher between radula transcriptomes compared to the mantle and foot. However, the same is not true when all differentially expressed genes or only non-synonymous SNPs are considered. Finally, putative homologs of some candidate genes for radula diversification (hh, arx, gbb) were also found to contribute to trophic specialization in cichlids and Darwin’s finches. (4) Conclusions: Our results are in line with diversifying selection on the radula driving Tylomelania ecomorph divergence and indicate that some molecular pathways may be especially prone to adaptive diversification, even across phylogenetically distant animal groups.

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.