The evolution of cichlid fish egg-spots is linked with a cis-regulatory change

The molecular basis of phenotypic evolution: beyond the usual suspects

It has been well documented that mutations in coding DNA or cis-regulatory elements underlie natural phenotypic variation in many organisms. However, the development of sophisticated functional tools in recent years in a wide range of traditionally non-model systems have revealed many 'unusual suspects' in the molecular bases of phenotypic evolution, including upstream open reading frames (uORFs), cryptic splice sites, and small RNAs. Furthermore, large-scale genome sequencing, especially long-read sequencing, has identified a cornucopia of structural variation underlying phenotypic divergence and elucidated the composition of supergenes that control complex multi-trait polymorphisms. In this review article we highlight recent studies that demonstrate this great diversity of molecular mechanisms producing adaptive genetic variation and the panoply of evolutionary paths leading to the 'grandeur of life'.

Keratinocytes drive the epithelial hyperplasia key to sea lice resistance in coho salmon

BACKGROUND

Salmonid species have followed markedly divergent evolutionary trajectories in their interactions with sea lice. While sea lice parasitism poses significant economic, environmental, and animal welfare challenges for Atlantic salmon (Salmo salar) aquaculture, coho salmon (Oncorhynchus kisutch) exhibit near-complete resistance to sea lice, achieved through a potent epithelial hyperplasia response leading to rapid louse detachment. The molecular mechanisms underlying these divergent responses to sea lice are unknown.

RESULTS

We characterized the cellular and molecular responses of Atlantic salmon and coho salmon to sea lice using single-nuclei RNA sequencing. Juvenile fish were exposed to copepodid sea lice (Lepeophtheirus salmonis), and lice-attached pelvic fin and skin samples were collected 12 h, 24 h, 36 h, 48 h, and 60 h after exposure, along with control samples. Comparative analysis of control and treatment samples revealed an immune and wound-healing response that was common to both species, but attenuated in Atlantic salmon, potentially reflecting greater sea louse immunomodulation. Our results revealed unique but complementary roles of three layers of keratinocytes in the epithelial hyperplasia response leading to rapid sea lice rejection in coho salmon. Our results suggest that basal keratinocytes direct the expansion and mobility of intermediate and, especially, superficial keratinocytes, which eventually encapsulate the parasite.

CONCLUSIONS

Our results highlight the key role of keratinocytes in coho salmon's sea lice resistance and the diverged biological response of the two salmonid host species when interacting with this parasite. This study has identified key pathways and candidate genes that could be manipulated using various biotechnological solutions to improve Atlantic salmon sea lice resistance.

Genome-wide analysis of histone modifications can contribute to the identification of candidate cis-regulatory regions in the threespine stickleback fish

BACKGROUND

Cis-regulatory mutations often underlie phenotypic evolution. However, because identifying the locations of promoters and enhancers in non-coding regions is challenging, we have fewer examples of identified causative cis-regulatory mutations that underlie naturally occurring phenotypic variations than of causative amino acid-altering mutations. Because cis-regulatory elements have epigenetic marks of specific histone modifications, we can detect cis-regulatory elements by mapping and analyzing them. Here, we investigated histone modifications and chromatin accessibility with cleavage under targets and tagmentation (CUT&Tag) and assay for transposase-accessible chromatin-sequencing (ATAC-seq).

RESULTS

Using the threespine stickleback (Gasterosteus aculeatus) as a model, we confirmed that the genes for which nearby regions showed active marks, such as H3K4me1, H3K4me3, and high chromatin accessibility, were highly expressed. In contrast, the expression levels of genes for which nearby regions showed repressive marks, such as H3K27me3, were reduced, suggesting that our chromatin analysis protocols overall worked well. Genomic regions with peaks of histone modifications showed higher nucleotide diversity within and between populations. By comparing gene expression in the gills of the marine and stream ecotypes, we identified several insertions and deletions (indels) with transposable element fragments in the candidate cis-regulatory regions.

CONCLUSIONS

Thus, mapping and analyzing histone modifications can help identify cis-regulatory elements and accelerate the identification of causative mutations in the non-coding regions underlying naturally occurring phenotypic variations.

Distinct genetic origins of eumelanin levels and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. To understand the genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

Family resemblance in color-patch size is not affected by stress experience in a cichlid fish

Animal body coloration is often linked to social dominance and mating success. This is because it can carry information on an animal's body condition and competitive ability by reflecting the genetic quality of individuals or by responding to their current or past living conditions. The present study investigates genetic and environmental effects on a conspicuous color pattern of the cichlid fish Tropheus sp. black "Ikola," in which the size of a carotenoid-based yellow area on the body co-varies with social dominance. To examine environmental plasticity of the color pattern, we tested for effects of early-life stress, induced by reduced feeding of juveniles prior to color pattern formation, as well as effects of a stress treatment administered to fully colored adult fish. None of the stress treatments affected the color pattern as quantified by the width of the yellow bar. However, offspring bar width was correlated to parental values in mid-parent-mid-offspring regression analyses, and animal models estimated significant additive genetic effects on bar width, indicating heritability of the trait. Depending on the random effects structure of the animal models (i.e., whether including or excluding maternal and brood effects), narrow-sense heritability estimates for bar width ranged between 0.2 and 0.8, with the strongest statistical support for the highest estimate. In each of the alternative models, a large proportion of the total variance in bar width was explained by the included random effects, suggesting that bar width is strongly determined by genetic factors or shared maternal and brood environments, with limited scope for environmental influences later in life.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

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.

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.

Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus

Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.

Mechanisms Underlying the Formation and Evolution of Vertebrate Color Patterns

Vertebrates exhibit a wide range of color patterns, which play critical roles in mediating intra- and interspecific communication. Because of their diversity and visual accessibility, color patterns offer a unique and fascinating window into the processes underlying biological organization. In this review, we focus on describing many of the general principles governing the formation and evolution of color patterns in different vertebrate groups. We characterize the types of patterns, review the molecular and developmental mechanisms by which they originate, and discuss their role in constraining or facilitating evolutionary change. Lastly, we outline outstanding questions in the field and discuss different approaches that can be used to address them. Overall, we provide a unifying conceptual framework among vertebrate systems that may guide research into naturally evolved mechanisms underlying color pattern formation and evolution.

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.

Distinct genetic origins of eumelanin intensity and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. In order to understand genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation 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.

Effects of Wolbachia on Transposable Element Expression Vary Between Drosophila melanogaster Host Genotypes

Transposable elements (TEs) are repetitive DNA sequences capable of changing position in host genomes, thereby causing mutations. TE insertions typically have deleterious effects but they can also be beneficial. Increasing evidence of the contribution of TEs to adaptive evolution further raises interest in understanding what factors impact TE activity. Based on previous studies associating the bacterial endosymbiont Wolbachia with changes in the abundance of piRNAs, a mechanism for TE repression, and to transposition of specific TEs, we hypothesized that Wolbachia infection would interfere with TE activity. We tested this hypothesis by studying the expression of 14 TEs in a panel of 25 Drosophila melanogaster host genotypes, naturally infected with Wolbachia and annotated for TE insertions. The host genotypes differed significantly in Wolbachia titers inside individual flies, with broad-sense heritability around 20%, and in the number of TE insertions, which depended greatly on TE identity. By removing Wolbachia from the target host genotypes, we generated a panel of 25 pairs of Wolbachia-positive and Wolbachia-negative lines in which we quantified transcription levels for our target TEs. We found variation in TE expression that was dependent on Wolbachia status, TE identity, and host genotype. Comparing between pairs of Wolbachia-positive and Wolbachia-negative flies, we found that Wolbachia removal affected TE expression in 21.1% of the TE-genotype combinations tested, with up to 2.3 times differences in the median level of transcript. Our data show that Wolbachia can impact TE activity in host genomes, underscoring the importance this endosymbiont can have in the generation of genetic novelty in hosts.

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.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

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.

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 analyzing 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 with 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, for example, 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.

Genetic basis of speciation and adaptation: from loci to causative mutations

Does evolution proceed in small steps or large leaps? How repeatable is evolution? How constrained is the evolutionary process? Answering these long-standing questions in evolutionary biology is indispensable for both understanding how extant biodiversity has evolved and predicting how organisms and ecosystems will respond to changing environments in the future. Understanding the genetic basis of phenotypic diversification and speciation in natural populations is key to properly answering these questions. The leap forward in genome sequencing technologies has made it increasingly easier to not only investigate the genetic architecture but also identify the variant sites underlying adaptation and speciation in natural populations. Furthermore, recent advances in genome editing technologies are making it possible to investigate the functions of each candidate gene in organisms from natural populations. In this article, we discuss how these recent technological advances enable the analysis of causative genes and mutations and how such analysis can help answer long-standing evolutionary biology questions.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

CRISPR knockouts of pmela and pmelb engineered a golden tilapia by regulating relative pigment cell abundance

Premelanosome protein (pmel) is a key gene for melanogenesis. Mutations in this gene are responsible for white plumage in chicken, but its role in pigmentation of fish remains to be demonstrated. In this study we found that most fishes have two pmel genes arising from the teleost-specific whole genome duplication. Both pmela and pmelb were expressed at high levels in the eyes and skin of Nile tilapia. We mutated both genes in tilapia using CRISPR/Cas9. Homozygous mutation of pmela resulted in yellowish body color with weak vertical bars and a hypo-pigmented retinal pigment epithelium (RPE) due to significantly reduced number and size of melanophores. In contrast, we observed an increased number and size of xanthophores in mutants compared to wild-type fish. Homozygous mutation of pmelb resulted in a similar, but milder phenotype than pmela-/- mutants. Double mutation of pmela and pmelb resulted in loss of additional melanophores compared to the pmela-/- mutants, and also an increase in the number and size of xanthophores, producing a golden body color. The RPE pigmentation of pmela-/-;pmelb-/- was similar to pmela-/- mutants, with much less pigmentation than pmelb-/- mutants and wild-type fish. Taken together, our results indicate that, while both pmel genes are important for the formation of body color in tilapia, pmela plays a more important role than pmelb. To our knowledge, this is the first report on mutation of pmelb or both pmela;pmelb in fish. Studies on these mutants suggest new strategies for breeding golden tilapia, and also provide a new model for studies of pmel function in vertebrates.

Oca2 targeting using CRISPR/Cas9 in the Malawi cichlid Astatotilapia calliptera

Identifying genetic loci underlying trait variation provides insights into the mechanisms of diversification, but demonstrating causality and characterizing the role of genetic loci requires testing candidate gene function, often in non-model species. Here we establish CRISPR/Cas9 editing in Astatotilapia calliptera, a generalist cichlid of the remarkably diverse Lake Malawi radiation. By targeting the gene oca2 required for melanin synthesis in other vertebrate species, we show efficient editing and germline transmission. Gene edits include indels in the coding region, probably a result of non-homologous end joining, and a large deletion in the 3' untranslated region due to homology-directed repair. We find that oca2 knock-out A. calliptera lack melanin, which may be useful for developmental imaging in embryos and studying colour pattern formation in adults. As A. calliptera resembles the presumed generalist ancestor of the Lake Malawi cichlids radiation, establishing genome editing in this species will facilitate investigating speciation, adaptation and trait diversification in this textbook radiation.

Genome size evolution in the diverse insect order Trichoptera

Background

Genome size is implicated in the form, function, and ecological success of a species. Two principally different mechanisms are proposed as major drivers of eukaryotic genome evolution and diversity: polyploidy (i.e., whole-genome duplication) or smaller duplication events and bursts in the activity of repetitive elements. Here, we generated de novo genome assemblies of 17 caddisflies covering all major lineages of Trichoptera. Using these and previously sequenced genomes, we use caddisflies as a model for understanding genome size evolution in diverse insect lineages.

Results

We detect a ∼14-fold variation in genome size across the order Trichoptera. We find strong evidence that repetitive element expansions, particularly those of transposable elements (TEs), are important drivers of large caddisfly genome sizes. Using an innovative method to examine TEs associated with universal single-copy orthologs (i.e., BUSCO genes), we find that TE expansions have a major impact on protein-coding gene regions, with TE-gene associations showing a linear relationship with increasing genome size. Intriguingly, we find that expanded genomes preferentially evolved in caddisfly clades with a higher ecological diversity (i.e., various feeding modes, diversification in variable, less stable environments).

Conclusion

Our findings provide a platform to test hypotheses about the potential evolutionary roles of TE activity and TE-gene associations, particularly in groups with high species, ecological, and functional diversities.

Genetic basis of orange spot formation in the guppy (Poecilia reticulata)

Background

To understand the evolutionary significance of female mate choice for colorful male ornamentation, the underlying regulatory mechanisms of such ornamentation must be understood for examining how the ornaments are associated with “male qualities” that increase the fitness or sexual attractiveness of offspring. In the guppy (Poecilia reticulata), an established model system for research on sexual selection, females prefer males possessing larger and more highly saturated orange spots as potential mates. Although previous studies have identified some chromosome regions and genes associated with orange spot formation, the regulation and involvement of these genetic elements in orange spot formation have not been elucidated. In this study, the expression patterns of genes specific to orange spots and certain color developmental stages were investigated using RNA-seq to reveal the genetic basis of orange spot formation.

Results

Comparing the gene expression levels of male guppy skin with orange spots (orange skin) with those without any color spots (dull skin) from the same individuals identified 1102 differentially expressed genes (DEGs), including 630 upregulated genes and 472 downregulated genes in the orange skin. Additionally, the gene expression levels of the whole trunk skin were compared among the three developmental stages and 2247 genes were identified as DEGs according to color development. These analyses indicated that secondary differentiation of xanthophores may affect orange spot formation.

Conclusions

The results suggested that orange spots might be formed by secondary differentiation, rather than de novo generation, of xanthophores, which is induced by Csf1 and thyroid hormone signaling pathways. Furthermore, we suggested candidate genes associated with the areas and saturation levels of orange spots, which are both believed to be important for female mate choice and independently regulated. This study provides insights into the genetic and cellular regulatory mechanisms underlying orange spot formation, which would help to elucidate how these processes are evolutionarily maintained as ornamental traits relevant to sexual selection.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01942-2.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Mapping epigenetic divergence in the massive radiation of Lake Malawi cichlid fishes

Epigenetic variation modulates gene expression and can be heritable. However, knowledge of the contribution of epigenetic divergence to adaptive diversification in nature remains limited. The massive evolutionary radiation of Lake Malawi cichlid fishes displaying extensive phenotypic diversity despite extremely low sequence divergence is an excellent system to study the epigenomic contribution to adaptation. Here, we present a comparative genome-wide methylome and transcriptome study, focussing on liver and muscle tissues in phenotypically divergent cichlid species. In both tissues we find substantial methylome divergence among species. Differentially methylated regions (DMR), enriched in evolutionary young transposons, are associated with transcription changes of ecologically-relevant genes related to energy expenditure and lipid metabolism, pointing to a link between dietary ecology and methylome divergence. Unexpectedly, half of all species-specific DMRs are shared across tissues and are enriched in developmental genes, likely reflecting distinct epigenetic developmental programmes. Our study reveals substantial methylome divergence in closely-related cichlid fishes and represents a resource to study the role of epigenetics in species diversification.

Transposable elements and introgression introduce genetic variation in the invasive ant Cardiocondyla obscurior

Introduced populations of invasive organisms have to cope with novel environmental challenges, while having reduced genetic variation caused by founder effects. The mechanisms associated with this "genetic paradox of invasive species" has received considerable attention, yet few studies have examined the genomic architecture of invasive species. Populations of the heart node ant Cardiocondyla obscurior belong to two distinct lineages, a New World lineage so far only found in Latin America and a more globally distributed Old World lineage. In the present study, we use population genomic approaches to compare populations of the two lineages with apparent divergent invasive potential. We find that the strong genetic differentiation of the two lineages began at least 40,000 generations ago and that activity of transposable elements (TEs) has contributed significantly to the divergence of both lineages, possibly linked to the very unusual genomic distribution of TEs in this species. Furthermore, we show that introgression from the Old World lineage is a dominant source of genetic diversity in the New World lineage, despite the lineages' strong genetic differentiation. Our study uncovers mechanisms underlying novel genetic variation in introduced populations of C. obscurior that could contribute to the species' adaptive potential.

Transcriptomics unravels molecular players shaping dorsal lip hypertrophy in the vacuum cleaner cichlid, Gnathochromis permaxillaris

BACKGROUND

Teleosts display a spectacular diversity of craniofacial adaptations that often mediates ecological specializations. A considerable amount of research has revealed molecular players underlying skeletal craniofacial morphologies, but less is known about soft craniofacial phenotypes. Here we focus on an example of lip hypertrophy in the benthivorous Lake Tangnayika cichlid, Gnathochromis permaxillaris, considered to be a morphological adaptation to extract invertebrates out of the uppermost layer of mud bottom. We investigate the molecular and regulatory basis of lip hypertrophy in G. permaxillaris using a comparative transcriptomic approach.

RESULTS

We identified a gene regulatory network involved in tissue overgrowth and cellular hypertrophy, potentially associated with the formation of a locally restricted hypertrophic lip in a teleost fish species. Of particular interest were the increased expression level of apoda and fhl2, as well as reduced expression of cyp1a, gimap8, lama5 and rasal3, in the hypertrophic lip region which have been implicated in lip formation in other vertebrates. Among the predicted upstream transcription factors, we found reduced expression of foxp1 in the hypertrophic lip region, which is known to act as repressor of cell growth and proliferation, and its function has been associated with hypertrophy of upper lip in human.

CONCLUSION

Our results provide a genetic foundation for future studies of molecular players shaping soft and exaggerated, but locally restricted, craniofacial morphological changes in fish and perhaps across vertebrates. In the future, we advocate integrating gene regulatory networks of various craniofacial phenotypes to understand how they collectively govern trophic and behavioural adaptations.

Nile Tilapia: A Model for Studying Teleost Color Patterns

The diverse color patterns of cichlid fishes play an important role in mate choice and speciation. Here we develop the Nile tilapia (Oreochromis niloticus) as a model system for studying the developmental genetics of cichlid color patterns. We identified 4 types of pigment cells: melanophores, xanthophores, iridophores and erythrophores, and characterized their first appearance in wild-type fish. We mutated 25 genes involved in melanogenesis, pteridine metabolism, and the carotenoid absorption and cleavage pathways. Among the 25 mutated genes, 13 genes had a phenotype in both the F0 and F2 generations. None of F1 heterozygotes had phenotype. By comparing the color pattern of our mutants with that of red tilapia (Oreochromis spp), a natural mutant produced during hybridization of tilapia species, we found that the pigmentation of the body and eye is controlled by different genes. Previously studied genes like mitf, kita/kitlga, pmel, tyrb, hps4, gch2, csf1ra, pax7b, and bco2b were proved to be of great significance for color patterning in tilapia. Our results suggested that tilapia, a fish with 4 types of pigment cells and a vertically barred wild-type color pattern, together with various natural and artificially induced color gene mutants, can serve as an excellent model system for study color patterning in vertebrates.

Evolution of pigment cells and patterns: recent insights from teleost fishes

Skin pigment patterns of vertebrates are stunningly diverse, and nowhere more so than in teleost fishes. Several species, including relatives of zebrafish, recently evolved cichlid fishes of East Africa, clownfishes, deep sea fishes, and others are providing insights into pigment pattern evolution. This overview describes recent advances in understanding periodic patterns, like stripes and spots, the loss of patterns, and the role of cell-type diversification in generating pigmentation phenotypes. Advances in this area are being facilitated by the application of modern methods of gene editing, genomics, computational analysis, and other approaches to non-traditional model organisms having interesting pigmentary phenotypes. Several topics worthy of future attention are outlined as well.

Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes

Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.

Movement of transposable elements contributes to cichlid diversity

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

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

Ancient Lake Tanganyika in East Africa houses the world's ecologically and morphologically most diverse assemblage of cichlid fishes, and the third most species-rich after lakes Malawi and Victoria. Despite long-lasting scientific interest in the cichlid species flocks of the East African Great Lakes, for example in the context of adaptive radiation and explosive diversification, their taxonomy and systematics are only partially explored; and many cichlid species still await their formal description. Here, we provide a current inventory of the cichlid fish fauna of Lake Tanganyika, providing a complete list of all valid 208 Tanganyikan cichlid species, and discuss the taxonomic status of more than 50 undescribed taxa on the basis of the available literature as well as our own observations and collections around the lake. This leads us to conclude that there are at least 241 cichlid species present in Lake Tanganyika, all but two are endemic to the basin. We finally summarize some of the major taxonomic challenges regarding Lake Tanganyika's cichlid fauna. The taxonomic inventory of the cichlid fauna of Lake Tanganyika presented here will facilitate future research on the taxonomy and systematics and the ecology and evolution of the species flock, as well as its conservation.

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

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

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

The Genomic Substrate for Adaptive Radiation: Copy Number Variation across 12 Tribes of African Cichlid Species

The initial sequencing of five cichlid genomes revealed an accumulation of genetic variation, including extensive copy number variation in cichlid lineages particularly those that have undergone dramatic evolutionary radiation. Gene duplication has the potential to generate substantial molecular substrate for the origin of evolutionary novelty. We use array-based comparative heterologous genomic hybridization to identify copy number variation events (CNVEs) for 168 samples representing 53 cichlid species including the 5 species for which full genome sequence is available. We identify an average of 50-100 CNVEs per individual. For those species represented by multiple samples, we identify 150-200 total CNVEs suggesting a substantial amount of intraspecific variation. For these species, only ∼10% of the detected CNVEs are fixed. Hierarchical clustering of species according to CNVE data recapitulates phylogenetic relationships fairly well at both the tribe and radiation level. Although CNVEs are detected on all linkage groups, they tend to cluster in "hotspots" and are likely to contain and be flanked by transposable elements. Furthermore, we show that CNVEs impact functional categories of genes with potential roles in adaptive phenotypes that could reasonably promote divergence and speciation in the cichlid clade. These data contribute to a more complete understanding of the molecular basis for adaptive natural selection, speciation, and evolutionary radiation.

Developmental and Cellular Basis of Vertical Bar Color Patterns in the East African Cichlid Fish Haplochromis latifasciatus

The East African adaptive radiations of cichlid fishes are renowned for their diversity in coloration. Yet, the developmental basis of pigment pattern formation remains largely unknown. One of the most common melanic patterns in cichlid fishes are vertical bar patterns. Here we describe the ontogeny of this conspicuous pattern in the Lake Kyoga species Haplochromis latifasciatus. Beginning with the larval stages we tracked the formation of this stereotypic color pattern and discovered that its macroscopic appearance is largely explained by an increase in melanophore density and accumulation of melanin during the first 3 weeks post-fertilization. The embryonal analysis is complemented with cytological quantifications of pigment cells in adult scales and the dermis beneath the scales. In adults, melanic bars are characterized by a two to threefold higher density of melanophores than in the intervening yellow interbars. We found no strong support for differences in other pigment cell types such as xanthophores. Quantitative PCRs for twelve known pigmentation genes showed that expression of melanin synthesis genes tyr and tyrp1a is increased five to sixfold in melanic bars, while xanthophore and iridophore marker genes are not differentially expressed. In summary, we provide novel insights on how vertical bars, one of the most widespread vertebrate color patterns, are formed through dynamic control of melanophore density, melanin synthesis and melanosome dispersal.

Zebrafish Pigment Pattern Formation: Insights into the Development and Evolution of Adult Form

Vertebrate pigment patterns are diverse and fascinating adult traits that allow animals to recognize conspecifics, attract mates, and avoid predators. Pigment patterns in fish are among the most amenable traits for studying the cellular basis of adult form, as the cells that produce diverse patterns are readily visible in the skin during development. The genetic basis of pigment pattern development has been most studied in the zebrafish, Danio rerio. Zebrafish adults have alternating dark and light horizontal stripes, resulting from the precise arrangement of three main classes of pigment cells: black melanophores, yellow xanthophores, and iridescent iridophores. The coordination of adult pigment cell lineage specification and differentiation with specific cellular interactions and morphogenetic behaviors is necessary for stripe development. Besides providing a nice example of pattern formation responsible for an adult trait of zebrafish, stripe-forming mechanisms also provide a conceptual framework for posing testable hypotheses about pattern diversification more broadly. Here, we summarize what is known about lineages and molecular interactions required for pattern formation in zebrafish, we review some of what is known about pattern diversification in Danio, and we speculate on how patterns in more distant teleosts may have evolved to produce a stunningly diverse array of patterns in nature.

Conservation and novelty in the microRNA genomic landscape of hyperdiverse cichlid fishes

MicroRNAs (miRNAs) play crucial roles in the post-transcriptional control of messenger RNA (mRNA). These miRNA-mRNA regulatory networks are present in nearly all organisms and contribute to development, phenotypic divergence, and speciation. To examine the miRNA landscape of cichlid fishes, one of the most species-rich families of vertebrates, we profiled the expression of both miRNA and mRNA in a diverse set of cichlid lineages. Among these, we found that conserved miRNAs differ from recently arisen miRNAs (i.e. lineage specific) in average expression levels, number of target sites, sequence variability, and physical clustering patterns in the genome. Furthermore, conserved miRNA target sites tend to be enriched at the 5' end of protein-coding gene 3' UTRs. Consistent with the presumed regulatory role of miRNAs, we detected more negative correlations between the expression of miRNA-mRNA functional pairs than in random pairings. Finally, we provide evidence that novel miRNA targets sites are enriched in genes involved in protein synthesis pathways. Our results show how conserved and evolutionarily novel miRNAs differ in their contribution to the genomic landscape and highlight their particular evolutionary roles in the adaptive diversification of cichlids.

Evolution of a soft-tissue foraging adaptation in African cichlids: Roles for novelty, convergence, and constraint

Understanding the origins of biodiversity demands consideration of both extrinsic (e.g., ecological opportunity) and intrinsic (e.g., developmental constraint) factors. Here, we use a combination of phylogenetic and genetic tools to address the origin of novelty in African cichlids. In particular, we focus on an extreme hypertrophied snout that is structurally integrated with the upper jaw. We show that this bizarre trait has evolved independently in at least two distinct and ecologically successful cichlid clades. We find that snout dimensions are decoupled both phenotypically and genetically, which has enabled it to evolve independently in multiple directions. Further, patterns of variation among species and within a genetic mapping pedigree suggest that relative to snout length, depth is under greater genetic and/or developmental constraint. Models of evolution suggest that snout shape is under selection for feeding behavior, with snout depth being important for algae scraping and snout length for sand sifting. Indeed, the deep snout of some algivores is achieved via an expansion of the intermaxillary ligament, which is important for jaw stability and may increase feeding performance. Overall, our data imply that the evolution of exaggerated snout depth required overcoming a genetic/developmental constraint, which led to expanded ecological opportunity via foraging adaptation.

Historical contingency shapes adaptive radiation in Antarctic fishes

Adaptive radiation illustrates links between ecological opportunity, natural selection and the generation of biodiversity. Central to adaptive radiation is the association between a diversifying lineage and the evolution of phenotypic variation that facilitates the use of new environments or resources. However, is not clear whether adaptive evolution or historical contingency is more important for the origin of key phenotypic traits in adaptive radiation. Here we use targeted sequencing of >250,000 loci across 46 species to examine hypotheses concerning the origin and diversification of key traits in the adaptive radiation of Antarctic notothenioid fishes. Contrary to expectations of adaptive evolution, we show that notothenioids experienced a punctuated burst of genomic diversification and evolved key skeletal modifications before the onset of polar conditions in the Southern Ocean. We show that diversifying selection in pathways associated with human skeletal dysplasias facilitates ecologically important variation in buoyancy among Antarctic notothenioid species, and demonstrate the sufficiency of altered trip11, col1a2 and col1a1a function in zebrafish (Danio rerio) to phenocopy skeletal reduction in Antarctic notothenioids. Rather than adaptation being driven by the cooling of the Antarctic, our results highlight the role of historical contingency in shaping the adaptive radiation of notothenioids. Understanding the historical and environmental context for the origin of key traits in adaptive radiations extends beyond reconstructing events that result in evolutionary innovation, as it also provides a context in forecasting the effects of climate change on the stability and evolvability of natural populations.

Genome of the Malawi golden cichlid fish (Melanochromis auratus) reveals exon loss of oca2 in an amelanistic morph

The tropical freshwater fish family Cichlidae is famous for its record-breaking rates of speciation and diversity in colors and color patterns. Here, we sequenced the genome of the Lake Malawi cichlid Melanochromis auratus to study the genetic basis of an amelanistic morph of this species that lacks the typical melanic stripes and markings. Genome sequencing of the amelanistic and wild-type morph revealed the loss of the second exon of the known pigmentation gene oculocutaneous albinism II (oca2), also known as p(ink-eyed dilution) gene or melanocyte-specific transporter gene. Additional genotyping confirms the complete association with this recessive Mendelian phenotype. The deletion results in a shorter transcript, lacking an acidic di-leucine domain that is crucial for trafficking of the Oca2 protein to melanosomes. The fact that oca2 is involved in a wide range of amelanistic morphs across vertebrates demonstrates its highly conserved function.

East African cichlid lineages (Teleostei: Cichlidae) might be older than their ancient host lakes: new divergence estimates for the east African cichlid radiation

BACKGROUND

Cichlids are a prime model system in evolutionary research and several of the most prominent examples of adaptive radiations are found in the East African Lakes Tanganyika, Malawi and Victoria, all part of the East African cichlid radiation (EAR). In the past, great effort has been invested in reconstructing the evolutionary and biogeographic history of cichlids (Teleostei: Cichlidae). In this study, we present new divergence age estimates for the major cichlid lineages with the main focus on the EAR based on a dataset encompassing representative taxa of almost all recognized cichlid tribes and ten mitochondrial protein genes. We have thoroughly re-evaluated both fossil and geological calibration points, and we included the recently described fossil †Tugenchromis pickfordi in the cichlid divergence age estimates.

RESULTS

Our results estimate the origin of the EAR to Late Eocene/Early Oligocene (28.71 Ma; 95% HPD: 24.43-33.15 Ma). More importantly divergence ages of the most recent common ancestor (MRCA) of several Tanganyika cichlid tribes were estimated to be substantially older than the oldest estimated maximum age of the Lake Tanganyika: Trematocarini (16.13 Ma, 95% HPD: 11.89-20.46 Ma), Bathybatini (20.62 Ma, 95% HPD: 16.88-25.34 Ma), Lamprologini (15.27 Ma; 95% HPD: 12.23-18.49 Ma). The divergence age of the crown haplochromine H-lineage is estimated to 22.8 Ma (95% HPD: 14.40-26.32 Ma) and of the Lake Malawi radiation to 4.07 Ma (95% HDP: 2.93-5.26 Ma). In addition, we recovered a novel lineage within the Lamprologini tribe encompassing only Lamprologus of the lower and central Congo drainage with its divergence estimated to the Late Miocene or early Pliocene. Furthermore we recovered two novel mitochondrial haplotype lineages within the Haplochromini tribe: 'Orthochromis' indermauri and 'Haplochormis' vanheusdeni.

CONCLUSIONS

Divergence time estimates of the MRCA of several Tanganyika cichlid tribes predate the age of the extant Lake Tanganyika basin, and hence are in line with the recently formulated "Melting-Pot Tanganyika" hypothesis. The radiation of the 'Lower Congo Lamprologus clade' might be linked with the Pliocene origin of the modern lower Congo rapids as has been shown for other Lower Congo cichlid assemblages. Finally, the age of origin of the Lake Malawi cichlid flock agrees well with the oldest age estimate for lacustrine conditions in Lake Malawi.

Chromosome-scale assemblies reveal the structural evolution of African cichlid genomes

BACKGROUND

African cichlid fishes are well known for their rapid radiations and are a model system for studying evolutionary processes. Here we compare multiple, high-quality, chromosome-scale genome assemblies to elucidate the genetic mechanisms underlying cichlid diversification and study how genome structure evolves in rapidly radiating lineages.

RESULTS

We re-anchored our recent assembly of the Nile tilapia (Oreochromis niloticus) genome using a new high-density genetic map. We also developed a new de novo genome assembly of the Lake Malawi cichlid, Metriaclima zebra, using high-coverage Pacific Biosciences sequencing, and anchored contigs to linkage groups (LGs) using 4 different genetic maps. These new anchored assemblies allow the first chromosome-scale comparisons of African cichlid genomes. Large intra-chromosomal structural differences (∼2-28 megabase pairs) among species are common, while inter-chromosomal differences are rare (<10 megabase pairs total). Placement of the centromeres within the chromosome-scale assemblies identifies large structural differences that explain many of the karyotype differences among species. Structural differences are also associated with unique patterns of recombination on sex chromosomes. Structural differences on LG9, LG11, and LG20 are associated with reduced recombination, indicative of inversions between the rock- and sand-dwelling clades of Lake Malawi cichlids. M. zebra has a larger number of recent transposable element insertions compared with O. niloticus, suggesting that several transposable element families have a higher rate of insertion in the haplochromine cichlid lineage.

CONCLUSION

This study identifies novel structural variation among East African cichlid genomes and provides a new set of genomic resources to support research on the mechanisms driving cichlid adaptation and speciation.

Developmental and comparative transcriptomic identification of iridophore contribution to white barring in clownfish

Actinopterygian fishes harbor at least eight distinct pigment cell types, leading to a fascinating diversity of colors. Among this diversity, the cellular origin of the white color appears to be linked to several pigment cell types such as iridophores or leucophores. We used the clownfish Amphiprion ocellaris, which has a color pattern consisting of white bars over a darker body, to characterize the pigment cells that underlie the white hue. We observe by electron microscopy that cells in white bars are similar to iridophores. In addition, the transcriptomic signature of clownfish white bars exhibits similarities with that of zebrafish iridophores. We further show by pharmacological treatments that these cells are necessary for the white color. Among the top differentially expressed genes in white skin, we identified several genes (fhl2a, fhl2b, saiyan, gpnmb, and apoD1a) and show that three of them are expressed in iridophores. Finally, we show by CRISPR/Cas9 mutagenesis that these genes are critical for iridophore development in zebrafish. Our analyses provide clues to the genomic underpinning of color diversity and allow identification of new iridophore genes in fish.