Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

From Genes to Pathways: A Curated Gene Approach to Accurate Pathway Reconstruction in Teleost Fish Transcriptomics

Interpreting the vast amounts of data generated by high-throughput sequencing technologies can often present a significant challenge, particularly for non-model organisms. While automated approaches like GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses are widely used, they often lack specificity for non-model organisms. To bridge this gap, we present a manually curated gene list tailored for teleost fish transcriptomics. This resource focuses on key biological processes crucial for understanding teleost fish physiology, development, and adaptation, including hormone signaling, various metabolic pathways, appetite regulation, digestion, gastrointestinal function, vision, ossification, osmoregulation, and pigmentation. Developed through collaborative efforts of specialists in diverse fields, the list prioritizes genes with established roles in teleost physiology, experimental evidence, and conservation across species. This curated list aims to provide researchers with a reliable starting point for transcriptomic analyses, offering a carefully evaluated set of genes relevant to current research priorities. By streamlining the process of gene selection and interpretation, this resource supports the broader teleost fish research community in designing and analyzing studies that investigate molecular responses to developmental and environmental changes. We encourage the scientific community to collaboratively expand and refine this list, ensuring its continued relevance and utility for teleost fish research.

Gene flow throughout the evolutionary history of a colour polymorphic and generalist clownfish

Even seemingly homogeneous on the surface, the oceans display high environmental heterogeneity across space and time. Indeed, different soft barriers structure the marine environment, which offers an appealing opportunity to study various evolutionary processes such as population differentiation and speciation. Here, we focus on Amphiprion clarkii (Actinopterygii; Perciformes), the most widespread of clownfishes that exhibits the highest colour polymorphism. Clownfishes can only disperse during a short pelagic larval phase before their sedentary adult lifestyle, which might limit connectivity among populations, thus facilitating speciation events. Consequently, the taxonomic status of A. clarkii has been under debate. We used whole-genome resequencing data of 67 A. clarkii specimens spread across the Indian and Pacific Oceans to characterize the species' population structure, demographic history and colour polymorphism. We found that A. clarkii spread from the Indo-Pacific Ocean to the Pacific and Indian Oceans following a stepping-stone dispersal and that gene flow was pervasive throughout its demographic history. Interestingly, colour patterns differed noticeably among the Indonesian populations and the two populations at the extreme of the sampling distribution (i.e. Maldives and New Caledonia), which exhibited more comparable colour patterns despite their geographic and genetic distances. Our study emphasizes how whole-genome studies can uncover the intricate evolutionary past of wide-ranging species with diverse phenotypes, shedding light on the complex nature of the species concept paradigm.

The nature and distribution of putative non-functional alleles suggest only two independent events at the origins of Astyanax mexicanus cavefish populations

BACKGROUND

Several studies suggested that cavefish populations of Astyanax mexicanus settled during the Late Pleistocene. This implies that the cavefish's most conspicuous phenotypic changes, blindness and depigmentation, and more cryptic characters important for cave life, evolved rapidly.

RESULTS

Using the published genomes of 47 Astyanax cavefish from la Cueva de El Pachón, El Sótano de la Tinaja, La Cueva Chica and El Sótano de Molino, we searched for putative loss-of-function mutations in previously defined sets of genes, i.e., vision, circadian clock and pigmentation genes. Putative non-functional alleles for four vision genes were identified. Then, we searched genome-wide for putative non-functional alleles in these four cave populations. Among 512 genes with segregating putative non-functional alleles in cavefish that are absent in surface fish, we found an enrichment in visual perception genes. Among cavefish populations, different levels of shared putative non-functional alleles were found. Using a subset of 12 genes for which putative loss-of-function mutations were found, we extend the analysis of shared pseudogenes to 11 cave populations. Using a subset of six genes for which putative loss-of-function mutations were found in the El Sótano del Toro population, where extensive hybridization with surface fish occurs, we found a correlation between the level of eye regression and the amount of putative non-functional alleles.

CONCLUSIONS

We confirm that very few putative non-functional alleles are present in a large set of vision genes, in accordance with the recent origin of Astyanax mexicanus cavefish. Furthermore, the genome-wide analysis indicates an enrichment of putative loss-of-function alleles in genes with vision-related GO-terms, suggesting that visual perception may be the function chiefly impacted by gene losses related to the shift from a surface to a cave environment. The geographic distribution of putative loss-of-function alleles newly suggests that cave populations from Sierra de Guatemala and Sierra de El Abra share a common origin, albeit followed by independent evolution for a long period. It also supports that populations from the Micos area have an independent origin. In El Sótano del Toro, the troglomorphic phenotype is maintained despite massive introgression of the surface genome.

Anemonefishes: A model system for evolutionary genomics

Anemonefishes are an iconic group of coral reef fish particularly known for their mutualistic relationship with sea anemones. This mutualism is especially intriguing as it likely prompted the rapid diversification of anemonefish. Understanding the genomic architecture underlying this process has indeed become one of the holy grails of evolutionary research in these fishes. Recently, anemonefishes have also been used as a model system to study the molecular basis of highly complex traits such as color patterning, social sex change, larval dispersal and life span. Extensive genomic resources including several high-quality reference genomes, a linkage map, and various genetic tools have indeed enabled the identification of genomic features controlling some of these fascinating attributes, but also provided insights into the molecular mechanisms underlying adaptive responses to changing environments. Here, we review the latest findings and new avenues of research that have led to this group of fish being regarded as a model for evolutionary genomics.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

The chromosome-scale genome assembly of the yellowtail clownfish Amphiprion clarkii provides insights into the melanic pigmentation of anemonefish

Anemonefish are an emerging group of model organisms for studying genetic, ecological, evolutionary, and developmental traits of coral reef fish. The yellowtail clownfish Amphiprion clarkii possesses species-specific characteristics such as inter-species co-habitation, high intra-species color variation, no anemone specificity, and a broad geographic distribution, that can increase our understanding of anemonefish evolutionary history, behavioral strategies, fish-anemone symbiosis, and color pattern evolution. Despite its position as an emerging model species, the genome of A. clarkii is yet to be published. Using PacBio long-read sequencing and Hi-C chromatin capture technology, we generated a high-quality chromosome-scale genome assembly initially comprised of 1,840 contigs with an N50 of 1,203,211 bp. These contigs were successfully anchored into 24 chromosomes of 843,582,782 bp and annotated with 25,050 protein-coding genes encompassing 97.0% of conserved actinopterygian genes, making the quality and completeness of this genome the highest among all published anemonefish genomes to date. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further analyses revealed higher copy numbers of erbb3b (a gene involved in melanocyte development) in A. clarkii compared with other anemonefish, thus suggesting a possible link between erbb3b and the natural melanism polymorphism observed in A. clarkii. The publication of this high-quality genome, along with A. clarkii's many unique traits, position this species as an ideal model organism for addressing scientific questions across a range of disciplines.

Inter-paralog amino acid inversion events in large phylogenies of duplicated proteins

Connecting protein sequence to function is becoming increasingly relevant since high-throughput sequencing studies accumulate large amounts of genomic data. In order to go beyond the existing database annotation, it is fundamental to understand the mechanisms underlying functional inheritance and divergence. If the homology relationship between proteins is known, can we determine whether the function diverged? In this work, we analyze different possibilities of protein sequence evolution after gene duplication and identify “inter-paralog inversions”, i.e., sites where the relationship between the ancestry and the functional signal is decoupled. The amino acids in these sites are masked from being recognized by other prediction tools. Still, they play a role in functional divergence and could indicate a shift in protein function. We develop a method to specifically recognize inter-paralog amino acid inversions in a phylogeny and test it on real and simulated datasets. In a dataset built from the Epidermal Growth Factor Receptor (EGFR) sequences found in 88 fish species, we identify 19 amino acid sites that went through inversion after gene duplication, mostly located at the ligand-binding extracellular domain. Our work uncovers an outcome of protein duplications with direct implications in protein functional annotation and sequence evolution. The developed method is optimized to work with large protein datasets and can be readily included in a targeted protein analysis pipeline.

Proteins are critical components of living systems because they facilitate most biological processes like protein synthesis, DNA replication, chemical catalysis, etc. Proteins are encoded in their genes. During evolution, genes accumulate mutations that get translated at the protein level. These mutations can be “neutral” if they do not affect the protein function immediately and directly; otherwise, mutations can be functional if they directly modify protein function. An event that provides an opportunity to study protein function is gene duplication namely, when two copies of a gene encoding the same protein appear. One copy of the protein often retains the same function while the other is free to diverge and specialize to a different function. This work sheds light on an alternative outcome of gene duplication that might be critical to discern between neutral and functional mutations. By looking at 88 fish genomes, we found proteins in which the evolution of their sequences does not follow the expected pattern of divergence after gene duplication. In this case, the protein sequence of a subgroup of species diverges in the copy expected to retain its function, while the sequence is retained in the expectedly divergent one. We called this event “inter-paralog amino acid inversion”. Our data shows that this “inversion” event is correlated to function, and its detection has to be considered for assigning protein functions correctly.

Genetic and correlative light and electron microscopy evidence for the unique differentiation pathway of erythrophores in brown trout skin

Based on their cell ultrastructure, two types of erythrophores in the spotted skin regions of brown trout (Salmo trutta) were previously described. To test the hypothesis regarding the origin of a new cell type following genome duplication, we analysed the gene and paralogue gene expression patterns of erythrophores in brown trout skin. In addition, the ultrastructure of both erythrophore types was precisely examined using transmission electron microscopy (TEM) and correlative light microscopy and electron microscopy (CLEM). Ultrastructural differences between the sizes of erythrophore inclusions were confirmed; however, the overlapping inclusion sizes blur the distinction between erythrophore types, which we have instead defined as cell subtypes. Nevertheless, the red spots of brown trout skin with subtype 2 erythrophores, exhibited unique gene expression patterns. Many of the upregulated genes are involved in melanogenesis or xanthophore differentiation. In addition, sox10, related to progenitor cells, was also upregulated in the red spots. The expressions of paralogues derived from two genome duplication events were also analysed. Multiple paralogues were overexpressed in the red spots compared with other skin regions, suggesting that the duplicated gene copies adopted new functions and contributed to the origin of a new cell subtype that is characteristic for red spot. Possible mechanisms regarding erythrophore origin are proposed and discussed. To the best of our knowledge, this is the first study to evaluate pigment cell types in the black and red spots of brown trout skin using the advanced CLEM approach together with gene expression profiling.

Two Genomic Loci Control Three Eye Colors in the Domestic Pigeon (Columba livia)

The iris of the eye shows striking color variation across vertebrate species, and may play important roles in crypsis and communication. The domestic pigeon (Columba livia) has three common iris colors, orange, pearl (white), and bull (dark brown), segregating in a single species, thereby providing a unique opportunity to identify the genetic basis of iris coloration. We used comparative genomics and genetic mapping in laboratory crosses to identify two candidate genes that control variation in iris color in domestic pigeons. We identified a nonsense mutation in the solute carrier SLC2A11B that is shared among all pigeons with pearl eye color, and a locus associated with bull eye color that includes EDNRB2, a gene involved in neural crest migration and pigment development. However, bull eye is likely controlled by a heterogeneous collection of alleles across pigeon breeds. We also found that the EDNRB2 region is associated with regionalized plumage depigmentation (piebalding). Our study identifies two candidate genes for eye colors variation, and establishes a genetic link between iris and plumage color, two traits that vary widely in the evolution of birds and other vertebrates.

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.

Identification of Antimicrobial Peptide Genes in Black Rockfish Sebastes schlegelii and Their Responsive Mechanisms to Edwardsiella tarda Infection

The black rockfish, Sebastes schlegelii, is a typical viviparous teleost, which belongs to the family Scorpaenidae. Due to its high economic and ecological values, S. schlegelii has been widely cultured in East Asian countries. With the enlargement of cultivation scale, bacterial and viral diseases have become the main threats to the farming industry of S. schlegelii, which have resulted in significant economic losses. In this study, Illumina shotgun sequencing, single-molecule real-time (SMRT) sequencing, 10× genomics and high-throughput chromosome conformation capture (Hi-C) technologies were collectively applied to assemble the genome of S. schlegelii. Then, we identified the antimicrobial peptide genes (AMPs) in the S. schlegelii genome. In total, 214 AMPs were identified in the S. schlegelii genome, which can be divided into 33 classes according to the annotation and cataloging of the Antimicrobial Peptides Database (APD3). Among these AMPs, thrombin-derived C-terminal peptide (TCP) was the dominant type, followed by RegIIIgamma and chemokine. The amino acid sequences of the TCP, cgUbiquitin, RegIIIalpha, RegIIIgamma, chemokine shared 32.55%, 42.63%, 29.87%, 28.09%, and 32.15% similarities among the same type in S. schlegelii. Meanwhile, the expression patterns of these AMPs in nine healthy tissues and at different infection time points in intestine were investigated. The results showed that the numbers and types of AMPs that responded to Edwardsiella tarda infection gradually increased as the infection progressed. In addition, we analyzed the phylogenetic relationships of hepcidins in teleost. The identification of AMPs based on the whole genome could provide a comprehensive database of potential AMPs, and benefit for the understanding of the molecular mechanisms of immune responses to E. tarda infection in S. schlegelii. This would further offer insights into an accurate and effective design and development of AMP for aquaculture therapy in the future.

Contrasting Gene Decay in Subterranean Vertebrates: Insights from Cavefishes and Fossorial Mammals

Evolution sometimes proceeds by loss, especially when structures and genes become dispensable after an environmental shift relaxes functional constraints. Subterranean vertebrates are outstanding models to analyze this process, and gene decay can serve as a readout. We sought to understand some general principles on the extent and tempo of the decay of genes involved in vision, circadian clock, and pigmentation in cavefishes. The analysis of the genomes of two Cuban species belonging to the genus Lucifuga provided evidence for the largest loss of eye-specific genes and nonvisual opsin genes reported so far in cavefishes. Comparisons with a recently evolved cave population of Astyanax mexicanus and three species belonging to the Chinese tetraploid genus Sinocyclocheilus revealed the combined effects of the level of eye regression, time, and genome ploidy on eye-specific gene pseudogenization. The limited extent of gene decay in all these cavefishes and the very small number of loss-of-function mutations per pseudogene suggest that their eye degeneration may not be very ancient, ranging from early to late Pleistocene. This is in sharp contrast with the identification of several vision genes carrying many loss-of-function mutations in ancient fossorial mammals, further suggesting that blind fishes cannot thrive more than a few million years in cave ecosystems.

Stages of embryonic development in the live-bearing fish, Gambusia holbrooki

BACKGROUND

Divergent morphology and placentation of Poeciliids make them suitable model for investigating how evolutionary selection has altered and conserved the developmental mechanisms. However, there is limited description of their embryonic staging, despite representing a key evolutionary node that shares developmental strategy with placental vertebrates. Here, we describe the embryonic developmental stages of Gambusia holbrooki from zygote to parturition using freshly harvested embryos.

RESULTS

We defined 40 embryonic stages using a numbered (stages 0-39; zygote to parturition, respectively) and named (grouped into seven periods, ie, zygote, cleavage, blastula, gastrula, segmentation, pharyngula, and parturition) staging system. Two sets of quantitative (ie, egg diameter, embryonic total length, otic vesicle closure index, heart rates, the number of caudal fin rays and elements) and qualitative (ie, three-dimensional analysis of images and key morphological criteria) data were acquired and used in combination to describe each stage. All 40 stages are separated by well-defined morphological traits, revealing developmental novelties that are influenced by narrow perivitelline space, placentation, internal gestation, and sex differentiation.

CONCLUSIONS

The principal diagnostic features described are quick, reliable, and easy to apply. This system will benefit researchers investigating molecular ontogeny, particularly sexual differentiation mechanisms in G. holbrooki.

Thyroid hormones regulate the formation and environmental plasticity of white bars in clownfishes

Determining how plasticity of developmental traits responds to environmental conditions is a challenge that must combine evolutionary sciences, ecology, and developmental biology. During metamorphosis, fish alter their morphology and color pattern according to environmental cues. We observed that juvenile clownfish (Amphiprion percula) modulate the developmental timing of their adult white bar formation during metamorphosis depending on the sea anemone species in which they are recruited. We observed an earlier formation of white bars when clownfish developed with Stichodactyla gigantea (Sg) than with Heteractis magnifica (Hm). As these bars, composed of iridophores, form during metamorphosis, we hypothesized that timing of their development may be thyroid hormone (TH) dependent. We treated clownfish larvae with TH and found that white bars developed earlier than in control fish. We further observed higher TH levels, associated with rapid white bar formation, in juveniles recruited in Sg than in Hm, explaining the faster white bar formation. Transcriptomic analysis of Sg recruits revealed higher expression of duox, a dual oxidase implicated in TH production as compared to Hm recruits. Finally, we showed that duox is an essential regulator of iridophore pattern timing in zebrafish. Taken together, our results suggest that TH controls the timing of adult color pattern formation and that shifts in duox expression and TH levels are associated with ecological differences resulting in divergent ontogenetic trajectories in color pattern development.

Phylogenetic and developmental analyses indicate complex functions of calcium-activated potassium channels in zebrafish embryonic development

BACKGROUND

Calcium-activated potassium channels (KCa) are a specific type of potassium channel activated by intracellular calcium concentration changes. This group of potassium channels plays fundamental roles ranging from regulating neuronal excitability to immune cell activation. Many human diseases such as schizophrenia, hypertension, epilepsy, and cancers have been linked to mutations in this group of potassium channels. Although the KCa channels have been extensively studied electrophysiologically and pharmacologically, their spatiotemporal gene expression during embryogenesis remains mostly unknown.

RESULTS

Using zebrafish as a model, we identified and renamed 14 KCa genes. We further performed phylogenetic and syntenic analyses on vertebrate KCa genes. Our data revealed that the number of KCa genes in zebrafish was increased, most likely due to teleost-specific whole-genome duplication. Moreover, we examined zebrafish KCa gene expression during early embryogenesis. The duplicated ohnologous genes show distinct and overlapped gene expression. Furthermore, we found that zebrafish KCa genes are expressed in various tissues and organs (somites, fins, olfactory regions, eye, kidney, and so on) and neuronal tissues, suggesting that they may play important roles during zebrafish embryogenesis.

CONCLUSIONS

Our phylogenetic and developmental analyses shed light on the potential functions of the KCa genes during embryogenesis related to congenital diseases and human channelopathies.

Overview of the betta fish genome regarding species radiation, parental care, behavioral aggression, and pigmentation model relevant to humans

BACKGROUND

The Siamese fighting fish (Betta splendens, also known as the betta) is well known in aquarium markets, and also presents an exciting new research model for studying parental care, aggressive behavior, and cryptically diverse pigmentation. However, concentrated efforts are required, both in the context of conservation biology and in its genetics, to address the problems of ongoing outbreeding depression, loss of biodiversity, and lack of scientific biological information.

OBJECTIVE

The evolutionary dynamics of the betta must be better understood at the genomic scale in order to resolve the phylogenetic status of unrecognized species, develop molecular markers to study variation in traits, and identify interesting sets of genes encoding various bioresource functions.

METHODS

The recent revolution in multi-omics approaches such as genomics, transcriptomics, epigenomics, and proteomics has uncovered genetic diversity and gained insights into many aspects of betta bioresources.

RESULTS

Here, we present current research and future plans in an ongoing megaproject to characterize the betta genome as de novo assemblies, genes and repeat annotations, generating data to study diverse biological phenomena. We highlight key questions that require answers and propose new directions and recommendations to develop bioresource management to protect and enhance the betta genus.

CONCLUSION

Successful accomplishment of these plans will allow the creation of a reference annotated genome and provide valuable information at the molecular level that can be utilized to sustain biodiversity and eco-management of the betta to improve breeding programs for future biomedical research.

Global tissue transcriptomic analysis to improve genome annotation and unravel skin pigmentation in goldfish

Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.

Evolution of the potassium channel gene Kcnj13 underlies colour pattern diversification in Danio fish

The genetic basis of morphological variation provides a major topic in evolutionary developmental biology. Fish of the genus Danio display colour patterns ranging from horizontal stripes, to vertical bars or spots. Stripe formation in zebrafish, Danio rerio, is a self-organizing process based on cell-contact mediated interactions between three types of chromatophores with a leading role of iridophores. Here we investigate genes known to regulate chromatophore interactions in zebrafish that might have evolved to produce a pattern of vertical bars in its sibling species, Danio aesculapii. Mutant D. aesculapii indicate a lower complexity in chromatophore interactions and a minor role of iridophores in patterning. Reciprocal hemizygosity tests identify the potassium channel gene obelix/Kcnj13 as evolved between the two species. Complementation tests suggest evolutionary change through divergence in Kcnj13 function in two additional Danio species. Thus, our results point towards repeated and independent evolution of this gene during colour pattern diversification.

Melanin-concentrating hormone like and somatolactin. A teleost-specific hypothalamic-hypophyseal axis system linking physiological and morphological pigmentation

Plastic adaptation to match the skin colour to the surrounding is key to survival. Two biological responses in skin colour are associated with background adaptation. A fast "physiological response" that aggregates/disperses the pigment organelles of skin chromatophores, and a slow "morphological response" that alters the type and/or density of pigment cells in the skin. Both responses are linked by unknown mechanisms. In this review, we discuss the role in skin colour regulation of two molecules that form part of a hypothalamic-hypophyseal pathway unique to teleosts, melanin-concentrating hormone "like" (MCHL) (previously known as MCH), and somatolactin. MCHL neurons project to the neurohypophysis and to the pars intermedia pituitary, where they interact with somatolactin-expressing cells. With a white background MCHL is released into the circulation to induce rapid melanosome aggregation and skin lightening. Somatolactin is also a fish-specific peptide whose expression and secretion are altered in organisms adapted chronically to white/black backgrounds, and that regulates morphological pigmentation. We discuss the evidence for a model whereby in teleosts, MCHL and somatolactin provide the previously unknown link between physiological and morphological pigmentation.

Synteny-Guided Resolution of Gene Trees Clarifies the Functional Impact of Whole-Genome Duplications

Whole-genome duplications (WGDs) have major impacts on the evolution of species, as they produce new gene copies contributing substantially to adaptation, isolation, phenotypic robustness, and evolvability. They result in large, complex gene families with recurrent gene losses in descendant species that sequence-based phylogenetic methods fail to reconstruct accurately. As a result, orthologs and paralogs are difficult to identify reliably in WGD-descended species, which hinders the exploration of functional consequences of WGDs. Here, we present Synteny-guided CORrection of Paralogies and Orthologies (SCORPiOs), a novel method to reconstruct gene phylogenies in the context of a known WGD event. WGDs generate large duplicated syntenic regions, which SCORPiOs systematically leverages as a complement to sequence evolution to infer the evolutionary history of genes. We applied SCORPiOs to the 320-My-old WGD at the origin of teleost fish. We find that almost one in four teleost gene phylogenies in the Ensembl database (3,394) are inconsistent with their syntenic contexts. For 70% of these gene families (2,387), we were able to propose an improved phylogenetic tree consistent with both the molecular substitution distances and the local syntenic information. We show that these synteny-guided phylogenies are more congruent with the species tree, with sequence evolution and with expected expression conservation patterns than those produced by state-of-the-art methods. Finally, we show that synteny-guided gene trees emphasize contributions of WGD paralogs to evolutionary innovations in the teleost clade.

The genome of Mekong tiger perch (Datnioides undecimradiatus) provides insights into the phylogenetic position of Lobotiformes and biological conservation

Mekong tiger perch (Datnioides undecimradiatus) is an ornamental and vulnerable freshwater fish native to the Mekong basin in Indochina, belonging to the order Lobotiformes. Here, we generated 121X stLFR co-barcode clean reads and 18X Oxford Nanopore MinION reads and obtained a 595 Mb Mekong tiger perch genome, which is the first whole genome sequence in the order Lobotiformes. Based on this genome, the phylogenetic tree analysis suggested that Lobotiformes is more closely related to Sciaenidae than to Tetraodontiformes, resolving a long-time dispute. We depicted the genes involved in pigment development in Mekong tiger perch and results confirmed that the four rate-limiting genes of pigment synthesis had been retained after fish-specific genome duplication. We also estimated the demographic history of Mekong tiger perch, which showed that the effective population size suffered a continuous reduction possibly related to the contraction of immune-related genes. Our study provided a reference genome resource for the Lobotiformes, as well as insights into the phylogenetic position of Lobotiformes and biological conservation.

Unique and non-redundant function of csf1r paralogues in regulation and evolution of post-embryonic development of the zebrafish

Evolution is replete with reuse of genes in different contexts, leading to multifunctional roles of signaling factors during development. Here, we explore osteoclast regulation during skeletal development through analysis of colony-stimulating factor 1 receptor (csf1r) function in the zebrafish. A primary role of Csf1r signaling is to regulate the proliferation, differentiation and function of myelomonocytic cells, including osteoclasts. We demonstrate the retention of two functional paralogues of csf1r in zebrafish. Mutant analysis indicates that the paralogues have shared, non-redundant roles in regulating osteoclast activity during the formation of the adult skeleton. csf1ra, however, has adopted unique roles in pigment cell patterning not seen in the second paralogue. We identify a unique noncoding element within csf1ra of fishes that is sufficient for controlling gene expression in pigment cells during development. As a role for Csf1r signaling in pigmentation is not observed in mammals or birds, it is likely that the overlapping roles of the two paralogues released functional constraints on csf1ra, allowing the signaling capacity of Csf1r to serve a novel function in the evolution of pigment pattern in fishes.

The identification of genes involved in the evolution of color patterns in fish

The genetic basis of morphological variation, both within and between species, provides a major topic in evolutionary biology. Teleost fish produce most elaborate color patterns, and among the more than 20000 species a number have been chosen for more detailed analyses because they are suitable to study particular aspects of color pattern evolution. In several fish species, color variants and pattern variants have been collected, transcriptome analyses have been carried out, and the recent advent of gene editing tools, such as CRISPR/Cas9, has allowed the production of mutants. Covering mostly the literature from the last three years, we discuss the cellular basis of coloration and the identification of loci involved in color pattern differences between sister species in cichlids and Danio species, in which cis-regulatory changes seem to prevail.

Plasticity for colour adaptation in vertebrates explained by the evolution of the genes pomc, pmch and pmchl

Different camouflages work best with some background matching colour. Our understanding of the evolution of skin colour is based mainly on the genetics of pigmentation ("background matching"), with little known about the evolution of the neuroendocrine systems that facilitate "background adaptation" through colour phenotypic plasticity. To address the latter, we studied the evolution in vertebrates of three genes, pomc, pmch and pmchl, that code for α-MSH and two melanin-concentrating hormones (MCH and MCHL). These hormones induce either dispersion/aggregation or the synthesis of pigments. We find that α-MSH is highly conserved during evolution, as is its role in dispersing/synthesizing pigments. Also conserved is the three-exon pmch gene that encodes MCH, which participates in feeding behaviours. In contrast, pmchl (known previously as pmch), is a teleost-specific intron-less gene. Our data indicate that in zebrafish, pmchl-expressing neurons extend axons to the pituitary, supportive of an MCHL hormonal role, whereas zebrafish and Xenopus pmch+ neurons send axons dorsally in the brain. The evolution of these genes and acquisition of hormonal status for MCHL explain different mechanisms used by vertebrates to background-adapt.

Insights into the Genomics of Clownfish Adaptive Radiation: Genetic Basis of the Mutualism with Sea Anemones

Clownfishes are an iconic group of coral reef fishes, especially known for their mutualism with sea anemones. This mutualism is particularly interesting as it likely acted as the key innovation that triggered clownfish adaptive radiation. Indeed, after the acquisition of the mutualism, clownfishes diversified into multiple ecological niches linked with host and habitat use. However, despite the importance of this mutualism, the genetic mechanisms allowing clownfishes to interact with sea anemones are still unclear. Here, we used a comparative genomics and molecular evolutionary analyses to investigate the genetic basis of clownfish mutualism with sea anemones. We assembled and annotated the genome of nine clownfish species and one closely related outgroup. Orthologous genes inferred between these species and additional publicly available teleost genomes resulted in almost 16,000 genes that were tested for positively selected substitutions potentially involved in the adaptation of clownfishes to live in sea anemones. We identified 17 genes with a signal of positive selection at the origin of clownfish radiation. Two of them (Versican core protein and Protein O-GlcNAse) show particularly interesting functions associated with N-acetylated sugars, which are known to be involved in sea anemone discharge of toxins. This study provides the first insights into the genetic mechanisms of clownfish mutualism with sea anemones. Indeed, we identified the first candidate genes likely to be associated with clownfish protection form sea anemones, and thus the evolution of their mutualism. Additionally, the genomic resources acquired represent a valuable resource for further investigation of the genomic basis of clownfish adaptive radiation.

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.

Evolution of Endothelin signaling and diversification of adult pigment pattern in Danio fishes

Fishes of the genus Danio exhibit diverse pigment patterns that serve as useful models for understanding the genes and cell behaviors underlying the evolution of adult form. Among these species, zebrafish D. rerio exhibit several dark stripes of melanophores with sparse iridophores that alternate with light interstripes of dense iridophores and xanthophores. By contrast, the closely related species D. nigrofasciatus has an attenuated pattern with fewer melanophores, stripes and interstripes. Here we demonstrate species differences in iridophore development that presage the fully formed patterns. Using genetic and transgenic approaches we identify the secreted peptide Endothelin-3 (Edn3)—a known melanogenic factor of tetrapods—as contributing to reduced iridophore proliferation and fewer stripes and interstripes in D. nigrofasciatus. We further show the locus encoding this factor is expressed at lower levels in D. nigrofasciatus owing to cis-regulatory differences between species. Finally, we show that functions of two paralogous loci encoding Edn3 have been partitioned between skin and non-skin iridophores. Our findings reveal genetic and cellular mechanisms contributing to pattern differences between these species and suggest a model for evolutionary changes in Edn3 requirements for pigment patterning and its diversification across vertebrates.

Neural crest derived pigment cells generate the spectacular variation in skin pigment patterns among vertebrates. Mammals and birds have just a single skin pigment cell, the melanocyte, whereas ectothermic vertebrates have several pigment cells including melanophores, iridophores and xanthophores, that together organize into a diverse array of patterns. In the teleost zebrafish, Danio rerio, an adult pattern of stripes depends on interactions between pigment cell classes and between pigment cells and their tissue environment. The close relative D. nigrofasciatus has fewer stripes and prior analyses suggested a difference between these species that lies extrinsic to the pigment cells themselves. A candidate for mediating this difference is Endothelin-3 (Edn3), essential for melanocyte development in warm-blooded animals, and required by all three classes of pigment cells in an amphibian. We show that Edn3 specifically promotes iridophore development in Danio, and that differences in Edn3 expression contribute to differences in iridophore complements, and striping, between D. rerio and D. nigrofasciatus. Our study reveals a novel function for Edn3 and provides new insights into how changes in gene expression yield morphogenetic outcomes to effect diversification of adult form.

Colour pattern component phenotypic divergence can be predicted by the light environment

The sensory drive hypothesis predicts that across different light environments sexually selected colour patterns will change to increase an animal's visual communication efficiency within different habitats. This is because individuals with more efficient signal components are likely to have more successful matings and hence produce more offspring. However, how colour pattern signals change over multiple generations under different light environmental conditions has not been tested experimentally. Here, we manipulated colour pattern signal efficiency by providing different ambient light environments over multiple generations to examine whether male colour pattern components change within large replicated populations of guppies (Poecilia reticulata). We report that colour patches change within populations over time and are phenotypically different among our three different light environments. Visual modelling suggests that the majority of these changes can be understood by considering the chroma, hue and luminance of each colour patch as seen by female guppies under each light environment. Taken together, our results support the hypothesis that different environmental conditions during signal reception can directly or indirectly drive the phenotypic diversification of visual signals within species.