An intronic transposon insertion associates with a trans-species color polymorphism in Midas cichlid fishes

Bilateral asymmetry in craniofacial structures and kinematics of feeding attacks in the scale-eating cichlid fish, Perissodus microlepis

Cichlid fishes are a textbook example for adaptive radiations, since they diversified into several hundred highly specialized species in each of three great East African lakes. Even scale-eating, an extremely specialized feeding mode, evolved independently multiple times in these radiations and in Lake Tanganyika alone, six endemic scale-eating species occupy this extremely specialized ecological niche. Perissodus microlepis went a step further, by evolving bilaterally asymmetrical heads with an intra-specific polymorphism where left- and right-headed morphs predominantly scrape scales from the opposite sides of their prey. While the bilateral asymmetry of scale-eating cichlids has been known, exactly which craniofacial features explain the laterality of the heads remained unclear. Here we aimed, by utilizing micro-computed tomography (μCT), to resolve this issue of how bilateral symmetry in the skeletal structure is broken in scale-eating Perissodus. Our 3D geometric morphometrics analysis clearly separated and identified the two groups of either left- or right-headed fish. In addition, we observed consistent asymmetric volume changes in the premaxilla, maxilla, and mandible of the craniofacial structures, where left-headed fish have larger jaw elements on the right side, and vice versa. The bimodality implies that the effect sizes of environmental factors might be minor while genetics might be responsible to a larger extent for the asymmetry observed in their head morphology. High-speed video analyses of attacks by asymmetrical morphotypes revealed that they utilize their asymmetrical mouth protrusion, as well as lateralized behavior, to re-orientate the gape towards the preferred side of their prey fish to more efficiently scrape scales.

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

BACKGROUND

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. We set out to bridge this gap and to understand the interactions between TEs and their cichlid hosts.

RESULTS

Here, we describe dynamic patterns of TE expression in African cichlid gonads and during early development. Orthology inference revealed strong conservation of TE silencing factors in cichlids, and 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 targets 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.

CONCLUSIONS

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.

Structural genomic variation and behavioral interactions underpin a balanced sexual mimicry polymorphism

How phenotypic diversity originates and persists within populations are classic puzzles in evolutionary biology. While balanced polymorphisms segregate within many species, it remains rare for both the genetic basis and the selective forces to be known, leading to an incomplete understanding of many classes of traits under balancing selection. Here, we uncover the genetic architecture of a balanced sexual mimicry polymorphism and identify behavioral mechanisms that may be involved in its maintenance in the swordtail fish Xiphophorus birchmanni. We find that ∼40% of X. birchmanni males develop a "false gravid spot," a melanic pigmentation pattern that mimics the "pregnancy spot" associated with sexual maturity in female live-bearing fish. Using genome-wide association mapping, we detect a single intergenic region associated with variation in the false gravid spot phenotype, which is upstream of kitlga, a melanophore patterning gene. By performing long-read sequencing within and across populations, we identify complex structural rearrangements between alternate alleles at this locus. The false gravid spot haplotype drives increased allele-specific expression of kitlga, which provides a mechanistic explanation for the increased melanophore abundance that causes the spot. By studying social interactions in the laboratory and in nature, we find that males with the false gravid spot experience less aggression; however, they also receive increased attention from other males and are disdained by females. These behavioral interactions may contribute to the maintenance of this phenotypic polymorphism in natural populations. We speculate that structural variants affecting gene regulation may be an underappreciated driver of balanced polymorphisms across diverse species.

Remarkably High Repeat Content in the Genomes of Sparrows: The Importance of Genome Assembly Completeness for Transposable Element Discovery

Transposable elements (TE) play critical roles in shaping genome evolution. Highly repetitive TE sequences are also a major source of assembly gaps making it difficult to fully understand the impact of these elements on host genomes. The increased capacity of long-read sequencing technologies to span highly repetitive regions promises to provide new insights into patterns of TE activity across diverse taxa. Here we report the generation of highly contiguous reference genomes using PacBio long-read and Omni-C technologies for three species of Passerellidae sparrow. We compared these assemblies to three chromosome-level sparrow assemblies and nine other sparrow assemblies generated using a variety of short- and long-read technologies. All long-read based assemblies were longer (range: 1.12 to 1.41 Gb) than short-read assemblies (0.91 to 1.08 Gb) and assembly length was strongly correlated with the amount of repeat content. Repeat content for Bell's sparrow (31.2% of genome) was the highest level ever reported within the order Passeriformes, which comprises over half of avian diversity. The highest levels of repeat content (79.2% to 93.7%) were found on the W chromosome relative to other regions of the genome. Finally, we show that proliferation of different TE classes varied even among species with similar levels of repeat content. These patterns support a dynamic model of TE expansion and contraction even in a clade where TEs were once thought to be fairly depauperate and static. Our work highlights how the resolution of difficult-to-assemble regions of the genome with new sequencing technologies promises to transform our understanding of avian genome evolution.

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.

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.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Transposable elements and xenobiotic resistance

Transposable elements or TEs are well known drivers of adaptive change in plants and animals but their role in insecticide resistance remains poorly documented. This review examines the potential role of transposons in resistance and identifies key areas where our understanding remains unclear. Despite well-known model systems such as upregulation of Drosophila Cyp6g1, many putative examples lack functional validation. The potential types of transposon-associated changes that could lead to resistance are reviewed, including changes in up-regulation, message stability, loss of function and alternative splicing. Where potential mechanisms appear absent from the resistance literature examples are drawn from other areas of biology. Finally, ways are suggested in which transgenic expression could be used to validate the biological significance of TE insertion. In the absence of such functional expression studies many examples of the association of TEs and resistance genes therefore remain as correlations.

SINE Insertion in the Pig Carbonic Anhydrase 5B (CA5B) Gene Is Associated with Changes in Gene Expression and Phenotypic Variation

Transposons are genetic elements that are present in mammalian genomes and occupy a large proportion of the pig genome, with retrotransposons being the most abundant. In a previous study, it was found that a SINE retrotransposon was inserted in the 1^st intron of the CA5B gene in pigs, and the present study aimed to investigate the SINE insertion polymorphism in this gene in different pig breeds. Polymerase chain reaction (PCR) was used to confirm the polymorphism in 11 pig breeds and wild boars), and it was found that there was moderate polymorphism information content in 9 of the breeds. Further investigation in cell experiments revealed that the 330 bp SINE insertion in the RIP-CA5B site promoted expression activity in the weak promoter region of this site. Additionally, an enhancer verification vector experiment showed that the 330 bp SINE sequence acted as an enhancer on the core promoter region upstream of the CA5B gene region. The expression of CA5B in adipose tissue (back fat and leaf fat) in individuals with the (SINE^+/+) genotype was significantly higher than those with (SINE^+/-) and (SINE^-/-) genotypes. The association analysis revealed that the (SINE^+/+) genotype was significantly associated with a higher back fat thickness than the (SINE^-/-) genotype. Moreover, it was observed that the insertion of SINE at the RIP-CA5B site carried ATTT repeats, and three types of (ATTT) repeats were identified among different individuals/breeds (i.e., (ATTT)₄, (ATTT)₆ and (ATTT)₉). Overall, the study provides insights into the genetic basis of adipose tissue development in pigs and highlights the role of a SINE insertion in the CA5B gene in this process.

The influence of transposable elements on animal colouration

Transposable elements (TEs) are mobile genetic sequences present within host genomes. TEs can contribute to the evolution of host traits, since transposition is mutagenic and TEs often contain host regulatory and protein coding sequences. We review cases where TEs influence animal colouration, reporting major patterns and outstanding questions. TE-induced colouration phenotypes typically arise via introduction of novel regulatory sequences and splice sites, affecting pigment cell development or pigment synthesis. We discuss if particular TE types may be more frequently involved in the evolution of colour variation in animals, given that examples involving long terminal repeat (LTR) elements appear to dominate. Currently, examples of TE-induced colouration phenotypes in animals mainly concern model and domesticated insect and mammal species. However, several influential recent examples, coupled with increases in genome sequencing, suggest cases reported from wild species will increase considerably.

Molecular parallelism in the evolution of a master sex-determining role for the anti-Mullerian hormone receptor 2 gene (amhr2) in Midas cichlids

The evolution of sex chromosomes and their differentiation from autosomes is a major event during genome evolution that happened many times in several lineages. The repeated evolution and lability of sex-determination mechanisms in fishes makes this a well-suited system to test for general patterns in evolution. According to current theory, differentiation is triggered by the suppression of recombination following the evolution of a new master sex-determining gene. However, the molecular mechanisms that establish recombination suppression are known from few examples, owing to the intrinsic difficulties of assembling sex-determining regions (SDRs). The development of forward-genetics and long-read sequencing have generated a wealth of data questioning central aspects of the current theory. Here, we demonstrate that sex in Midas cichlids is determined by an XY system, and identify and assemble the SDR by combining forward-genetics, long-read sequencing and optical mapping. We show how long-reads aid in the detection of artefacts in genotype-phenotype mapping that arise from incomplete genome assemblies. The male-specific region is restricted to a 100-kb segment on chromosome 4 that harbours transposable elements and a Y-specific duplicate of the anti-Mullerian receptor 2 gene, which has evolved master sex-determining functions repeatedly. Our data suggest that amhr2Y originated by an interchromosomal translocation from chromosome 20 to 4 pre-dating the split of Midas and Flier cichlids. In the latter, it is pseudogenized and translocated to another chromosome. Duplication of anti-Mullerian genes is a common route to establishing new sex determiners, highlighting the role of molecular parallelism in the evolution of sex determination.

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.

Noncanonical DNA structures are drivers of genome evolution

In addition to the canonical right-handed double helix, other DNA structures, termed 'non-B DNA', can form in the genomes across the tree of life. Non-B DNA regulates multiple cellular processes, including replication and transcription, yet its presence is associated with elevated mutagenicity and genome instability. These discordant cellular roles fuel the enormous potential of non-B DNA to drive genomic and phenotypic evolution. Here we discuss recent studies establishing non-B DNA structures as novel functional elements subject to natural selection, affecting evolution of transposable elements (TEs), and specifying centromeres. By highlighting the contributions of non-B DNA to repeated evolution and adaptation to changing environments, we conclude that evolutionary analyses should include a perspective of not only DNA sequence, but also its structure.

Advances in Vertebrate (Cyto)Genomics Shed New Light on Fish Compositional Genome Evolution

Cytogenetic and compositional studies considered fish genomes rather poor in guanine-cytosine content (GC%) because of a putative "sharp increase in genic GC% during the evolution of higher vertebrates". However, the available genomic data have not been exploited to confirm this viewpoint. In contrast, further misunderstandings in GC%, mostly of fish genomes, originated from a misapprehension of the current flood of data. Utilizing public databases, we calculated the GC% in animal genomes of three different, technically well-established fractions: DNA (entire genome), cDNA (complementary DNA), and cds (exons). Our results across chordates help set borders of GC% values that are still incorrect in literature and show: (i) fish in their immense diversity possess comparably GC-rich (or even GC-richer) genomes as higher vertebrates, and fish exons are GC-enriched among vertebrates; (ii) animal genomes generally show a GC-enrichment from the DNA, over cDNA, to the cds level (i.e., not only the higher vertebrates); (iii) fish and invertebrates show a broad(er) inter-quartile range in GC%, while avian and mammalian genomes are more constrained in their GC%. These results indicate no sharp increase in the GC% of genes during the transition to higher vertebrates, as stated and numerously repeated before. We present our results in 2D and 3D space to explore the compositional genome landscape and prepared an online platform to explore the AT/GC compositional genome evolution.

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.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

A DNA Replication Mechanism Can Explain Structural Variation at the Pigeon Recessive Red Locus

For species to adapt to their environment, evolution must act upon genetic variation that is present in the population. Elucidating the molecular mechanisms that give rise to this variation is thus of crucial importance for understanding how organisms evolve. In addition to variation caused by point mutations, structural variation (deletions, duplications, inversions, translocations) is also an important source of variety. Mechanisms involving recombination, transposition and retrotransposition, and replication have been proposed for generating structural variation, and each are capable of explaining certain rearrangements. In this study, we conduct a detailed analysis of two partially overlapping rearrangements (e1 and e2 allele) in domestic rock pigeon (Columba livia) which are both associated with the recessive red phenotype. We find that a replicative mechanism is best able to explain the complex architecture of the e1 allele, and is also compatible with the simpler architecture of the e2 allele as well.

Comparative ontogenetic and transcriptomic analyses shed light on color pattern divergence in cichlid fishes

Stripe patterns are a striking example for a repeatedly evolved color pattern. In the African adaptive radiations of cichlid fishes, stripes evolved several times independently. Previously, it has been suggested that regulatory evolution of a single gene, agouti-related-peptide 2 (agrp2), explains the evolutionary lability of this trait. Here, using a comparative transcriptomic approach, we performed comparisons between (adult) striped and nonstriped cichlid fishes of representatives of Lake Victoria and the two major clades of Lake Malawi (mbuna and non-mbuna lineage). We identify agrp2 to be differentially expressed across all pairwise comparisons, reaffirming its association with stripe pattern divergence. We therefore also provide evidence that agrp2 is associated with the loss of the nonstereotypic oblique stripe of Mylochromis mola. Complementary ontogenetic data give insights into the development of stripe patterns as well as vertical bar patterns that both develop postembryonically. Lastly, using the Lake Victoria species pair Haplochromis sauvagei and Pundamilia nyererei, we investigated the differences between melanic and non-melanic regions to identify additional genes that contribute to the formation of stripes. Expression differences-that most importantly also do not include agrp2-are surprisingly small. This suggests, at least in this species pair, that the stripe phenotype might be caused by a combination of more subtle transcriptomic differences or cellular changes without transcriptional correlates. In summary, our comprehensive analysis highlights the ontogenetic and adult transcriptomic differences between cichlids with different color patterns and serves as a basis for further investigation of the mechanistic underpinnings of their diversification.