Agouti-related peptide 2 facilitates convergent evolution of stripe patterns across cichlid fish radiations

Uncovering developmental diversity in the field

Many developmental biologists seldom leave the lab for research, relying instead on establishing colonies of traditional and emerging model systems. However, to fully understand the mechanisms and principles of development and evolution, including the role of ecology and the environment, it is important to study a diverse range of organisms in context. In this Perspective, we hear from five research teams from around the world about the importance and challenges of going into the field to study their organisms of interest. We also invite you to share your own fieldwork stories on the Node.

A whole-body micro-CT scan library that captures the skeletal diversity of Lake Malawi cichlid fishes

Here we describe a dataset of freely available, readily processed, whole-body μCT-scans of 56 species (116 specimens) of Lake Malawi cichlid fishes that captures a considerable majority of the morphological variation present in this remarkable adaptive radiation. We contextualise the scanned specimens within a discussion of their respective ecomorphological groupings and suggest possible macroevolutionary studies that could be conducted with these data. In addition, we describe a methodology to efficiently μCT-scan (on average) 23 specimens per hour, limiting scanning time and alleviating the financial cost whilst maintaining high resolution. We demonstrate the utility of this method by reconstructing 3D models of multiple bones from multiple specimens within the dataset. We hope this dataset will enable further morphological study of this fascinating system and permit wider-scale comparisons with other cichlid adaptive radiations.

Two novel genomes of fireflies with different degrees of sexual dimorphism reveal insights into sex-biased gene expression and dosage compensation

Sexual dimorphism arises because of divergent fitness optima between the sexes. Phenotypic divergence between sexes can range from mild to extreme. Fireflies, bioluminescent beetles, present various degrees of sexual dimorphism, with species showing very mild sexual dimorphism to species presenting female-specific neoteny, posing a unique framework to investigate the evolution of sexually dimorphic traits across species. In this work, we present novel assembled genomes of two firefly species, Lamprohiza splendidula and Luciola italica, species with different degrees of sexual dimorphism. We uncover high synteny conservation of the X-chromosome across ~ 180 Mya and find full X-chromosome dosage compensation in our two fireflies, hinting at common mechanism upregulating the single male X-chromosome. Different degrees of sex-biased expressed genes were found across two body parts showing different proportions of expression conservation between species. Interestingly, we do not find X-chromosome enrichment of sex-biased genes, but retrieve autosomal enrichment of sex-biased genes. We further uncover higher nucleotide diversity in the intronic regions of sex-biased genes, hinting at a maintenance of heterozygosity through sexual selection. We identify different levels of sex-biased gene expression divergence including a set of genes showing conserved sex-biased gene expression between species. Divergent and conserved sex-biased genes are good candidates to test their role in the maintenance of sexually dimorphic traits.

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.

Identification of pigmentation genes in skin, muscle and tail of a Thai-flag variety of Siamese fighting fish Betta splendens

Pigmentation genes expressed in skin, body muscle and tail of Thai-flag compared with Blue, White and Red varieties of Siamese fighting fish Betta splendens were identified. In total, 22,919 new unigenes were found. Pearson correlation and PCA analysis revealed that expression profiles of genes in muscle, skin and tail across solid color variety were similar. In contrast, those in skin and red tail part of Thai-flag were closely related but they showed different expression profiles with the white tail part. Moreover, 21,347-64,965 SNPs were identified in exonic regions of identified genes. In total, 28,899 genes were differentially expressed between paired comparisons of libraries where 13,907 genes (48.12 %) were upregulated and 14,992 genes (51.88 %) were downregulated. DEGs between paired libraries were 106-5775 genes relative to the compared libraries (56-2982 and 50-2782 for upregulated and downregulated DEGs). Interestingly, 432 pigmentation genes of B. splendens were found. Of these, 297 DEGs showed differential expression between varieties. Many DEGs in melanogenesis (Bsmcr1r, Bsmcr5r, and Bsslc2a15b), tyrosine metabolism (Bstyr, Bstyrp1b and Bsdct), stripe repressor (BsAsip1 and BsAsip2b), pteridine (Bsgch2) and carotenoid (BsBco2) biosynthesis were downregulated in the Thai-flag compared with solid color varieties. Expression of Bsbco1l, Bsfrem2b, Bskcnj13, Bszic2a and Bspah in skin, muscle and tail of Thai-flag, Blue, Red and White varieties was analyzed by qRT-PCR and revealed differential expression between fish varieties and showed anatomical tissue-preferred expression patterns in the same fish variety. The information could be applied to assist genetic-based development of new B. splendens varieties in the future.

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.

Towards a comprehensive understanding of the muscle proteome in Schizothorax labiatus: Insights from seasonal variations, metabolic responses, and reproductive signatures in the River Jhelum

Proteomics is a very advanced technique used for defining correlations, compositions and activities of hundreds of proteins from organisms as well as effectively used in identifying particular proteins with varying peptide lengths and amino acid counts. In the present study, an endeavour has been put forth to create muscle proteome expression of snow trout, Schizothorax labiatus. Liquid chromatography-mass spectrometry (LC-MS) using label free quantification (LFQ) technique has extensively been carried out to explore changes in protein metabolism and its composition to discriminate across species, clarify functions and pinpoint protein biomarkers from organisms. In LFQ technique, the abundances of proteins are determined based on the signal intensities of their corresponding peptides in mass spectrometry. The main benefit of using this method is that it doesn't require pre-labelling proteins with isotopic tags, which streamlines the experimental procedure and gets rid of any bias that might have been caused by the labelling process. LFQ techniques frequently offer a wider dynamic range, making it possible to detect and quantify proteins over a broad range of abundances obtained from the complex biological materials including fish muscle. The results of proteomic analysis could provide an insight in understanding about how various proteins are expressed in response to environmental challenges. For proteomic study, two different weight groups of S. labiatus were taken from River Jhelum based on biological, physiological and logistical factors. These groups corresponded to different life stages, such as younger size and adults/brooders in order to capture potential variations in the muscle proteome related to growth and development. The proteomic analysis of S. labiatus depicted that an overall of 220 proteins in male and 228 in female fish of group 1 were noted. However, when male and female S. labiatus were examined based on spectral count and peptide abundance using ProteinLynx Global Software, a total of 10 downregulated and 32 upregulated proteins were found. In group 2 of S. labiatus, a total of 249 proteins in male and 301 in female fish were documented. When the two genders of S. labiatus were likened to one another by LFQ technique, a total of 41 downregulated and 06 upregulated proteins were identified. The variability in the protein numbers between two fish weight groups reflected biological differences, influenced by factors such as age, developmental stages, physiological condition and reproductive activities. During the study, it was observed that S. labiatus exhibited downregulated levels of proteins that were involved in feeding and growth. The contributing factors to this manifestation could be explained by lower feeding and metabolic activity of fish and decreased food availability during winter in River Jhelum. Contrarily, the fish immune response proteins were found to be significantly over-expressed in S. labiatus, indicating that the environment was more likely to undergo increased microbial infection, pollution load and anthropogenic activities. In addition, it was also discovered that there was an upregulated expression of the reproductive proteins in S. labiatus, which could be linked to the fish's pre-spawning time as the fish used in this study was collected in the winter season which is the pre-spawning period of the fish. Therefore, the present study would be useful in obtaining new insights regarding the molecular makeup of species, methods of adaptation and reactions to environmental stresses. This information contributes to our understanding of basic science and may have applications in environmental monitoring, conservation and preservation of fish species.

Genetic basis and expression of ventral colour in polymorphic common lizards

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

recolorize: An R package for flexible colour segmentation of biological images

Colour pattern variation provides biological information in fields ranging from disease ecology to speciation dynamics. Comparing colour pattern geometries across images requires colour segmentation, where pixels in an image are assigned to one of a set of colour classes shared by all images. Manual methods for colour segmentation are slow and subjective, while automated methods can struggle with high technical variation in aggregate image sets. We present recolorize, an R package toolbox for human-subjective colour segmentation with functions for batch-processing low-variation image sets and additional tools for handling images from diverse (high-variation) sources. The package also includes export options for a variety of formats and colour analysis packages. This paper illustrates recolorize for three example datasets, including high variation, batch processing and combining with reflectance spectra, and demonstrates the downstream use of methods that rely on this output.

Multiple Routes to Color Convergence in a Radiation of Neotropical Poison Frogs

Convergent evolution is defined as the independent evolution of similar phenotypes in different lineages. Its existence underscores the importance of external selection pressures in evolutionary history, revealing how functionally similar adaptations can evolve in response to persistent ecological challenges through a diversity of evolutionary routes. However, many examples of convergence, particularly among closely related species, involve parallel changes in the same genes or developmental pathways, raising the possibility that homology at deeper mechanistic levels is an important facilitator of phenotypic convergence. Using the genus Ranitomeya, a young, color-diverse radiation of Neotropical poison frogs, we set out to 1) provide a phylogenetic framework for this group, 2) leverage this framework to determine if color phenotypes are convergent, and 3) to characterize the underlying coloration mechanisms to test whether color convergence occurred through the same or different physical mechanisms. We generated a phylogeny for Ranitomeya using ultraconserved elements and investigated the physical mechanisms underlying bright coloration, focusing on skin pigments. Using phylogenetic comparative methods, we identified several instances of color convergence, involving several gains and losses of carotenoid and pterin pigments. We also found a compelling example of nonparallel convergence, where, in one lineage, red coloration evolved through the red pterin pigment drosopterin, and in another lineage through red ketocarotenoids. Additionally, in another lineage, "reddish" coloration evolved predominantly through structural color mechanisms. Our study demonstrates that, even within a radiation of closely related species, convergent evolution can occur through both parallel and nonparallel mechanisms, challenging the assumption that similar phenotypes among close relatives evolve through the same mechanisms.

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.

Species differences in hormonally mediated gene expression underlie the evolutionary loss of sexually dimorphic coloration in Sceloporus lizards

Phenotypic sexual dimorphism often involves the hormonal regulation of sex-biased expression for underlying genes. However, it is generally unknown whether the evolution of hormonally mediated sexual dimorphism occurs through upstream changes in tissue sensitivity to hormone signals, downstream changes in responsiveness of target genes, or both. Here, we use comparative transcriptomics to explore these possibilities in 2 species of Sceloporus lizards exhibiting different patterns of sexual dichromatism. Sexually dimorphic S. undulatus develops blue and black ventral coloration in response to testosterone, while sexually monomorphic S. virgatus does not, despite exhibiting similar sex differences in circulating testosterone levels. We administered testosterone implants to juveniles of each species and used RNAseq to quantify gene expression in ventral skin. Transcriptome-wide responses to testosterone were stronger in S. undulatus than in S. virgatus, suggesting species differences in tissue sensitivity to this hormone signal. Species differences in the expression of genes for androgen metabolism and sex hormone-binding globulin were consistent with this idea, but expression of the androgen receptor gene was higher in S. virgatus, complicating this interpretation. Downstream of androgen signaling, we found clear species differences in hormonal responsiveness of genes related to melanin synthesis, which were upregulated by testosterone in S. undulatus, but not in S. virgatus. Collectively, our results indicate that hormonal regulation of melanin synthesis pathways contributes to the development of sexual dimorphism in S. undulatus, and that changes in the hormonal responsiveness of these genes in S. virgatus contribute to the evolutionary loss of ventral coloration.

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.

Colour polymorphism associated with a gene duplication in male wood tiger moths

Colour is often used as an aposematic warning signal, with predator learning expected to lead to a single colour pattern within a population. However, there are many puzzling cases where aposematic signals are also polymorphic. The wood tiger moth, Arctia plantaginis, displays bright hindwing colours associated with unpalatability, and males have discrete colour morphs which vary in frequency between localities. In Finland, both white and yellow morphs can be found, and these colour morphs also differ in behavioural and life-history traits. Here, we show that male colour is linked to an extra copy of a yellow family gene that is only present in the white morphs. This white-specific duplication, which we name valkea, is highly upregulated during wing development. CRISPR targeting valkea resulted in editing of both valkea and its paralog, yellow-e, and led to the production of yellow wings. We also characterise the pigments responsible for yellow, white, and black colouration, showing that yellow is partly produced by pheomelanins, while black is dopamine-derived eumelanin. Our results add to a growing number of studies on the genetic architecture of complex and seemingly paradoxical polymorphisms, and the role of gene duplications and structural variation in adaptive evolution.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

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

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

kcnj13 regulates pigment cell shapes in zebrafish and has diverged by cis-regulatory evolution between Danio species

Teleost fish of the genus Danio are excellent models to study the genetic and cellular bases of pigment pattern variation in vertebrates. The two sister species Danio rerio and Danio aesculapii show divergent patterns of horizontal stripes and vertical bars that are partly caused by the divergence of the potassium channel gene kcnj13. Here, we show that kcnj13 is required only in melanophores for interactions with xanthophores and iridophores, which cause location-specific pigment cell shapes and thereby influence colour pattern and contrast in D. rerio. Cis-regulatory rather than protein coding changes underlie kcnj13 divergence between the two Danio species. Our results suggest that homotypic and heterotypic interactions between the pigment cells and their shapes diverged between species by quantitative changes in kcnj13 expression during pigment pattern diversification.

Neural crest cells as a source of microevolutionary variation

Vertebrates have some of the most complex and diverse features in animals, from varied craniofacial morphologies to colorful pigmentation patterns and elaborate social behaviors. All of these traits have their developmental origins in a multipotent embryonic lineage of neural crest cells. This "fourth germ layer" is a vertebrate innovation and the source of a wide range of adult cell types. While others have discussed the role of neural crest cells in human disease and animal domestication, less is known about their role in contributing to adaptive changes in wild populations. Here, we review how variation in the development of neural crest cells and their derivatives generates considerable phenotypic diversity in nature. We focus on the broad span of traits under natural and sexual selection whose variation may originate in the neural crest, with emphasis on behavioral factors such as intraspecies communication that are often overlooked. In all, we encourage the integration of evolutionary ecology with developmental biology and molecular genetics to gain a more complete understanding of the role of this single cell type in trait covariation, evolutionary trajectories, and vertebrate diversity.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

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

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.

High level of novelty under the hood of convergent evolution

Little is known about the extent to which species use homologous regulatory architectures to achieve phenotypic convergence. By characterizing chromatin accessibility and gene expression in developing wing tissues, we compared the regulatory architecture of convergence between a pair of mimetic butterfly species. Although a handful of color pattern genes are known to be involved in their convergence, our data suggest that different mutational paths underlie the integration of these genes into wing pattern development. This is supported by a large fraction of accessible chromatin being exclusive to each species, including the de novo lineage-specific evolution of a modular optix enhancer. These findings may be explained by a high level of developmental drift and evolutionary contingency that occurs during the independent evolution of mimicry.

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

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

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

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

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.

Evolutionary convergence in body shape obscures taxonomic diversity in species of the African Labeo forskalii group: Case study of L. parvus Boulenger 1902 and L. ogunensis Boulenger 1910

Labeo is the third most diverse genus of African cyprinids and is widely distributed across the continent. Labeo parvus, a small species originally described from the Congo basin, has been considered the only species of the L. forskalii group distributed across five African ichthyofaunal provinces (Nilo-Sudan, Congo, Cuanza, and Upper and Lower Guinea). However, morphological similarity between L. parvus and numerous congeners remains a central cause of taxonomic confusion within the genus. Here we employed a phylogenetic comparative approach to assess phenotypic convergence among species of the L. forskalii group, investigate the taxonomic status of L. parvus sensu lato (sl) in west Africa, and reevaluate the composition and distribution of L. parvus sensu stricto (ss). Our phylogenetic analysis provides no support for a sister relationship between L. parvus ss and any of the west African Labeo parvus-like species. Geometric morphometric and molecular phylogenetic data indicate that L. parvus ss is a Congo basin endemic, and seemingly ecologically equivalent species found in west Africa are L. ogunensis, L. obscurus and other undescribed or previously synonymized species. We discuss our findings in terms of convergent evolution using phylomorphospace and tests for phylogenetic signal.

Transcriptome comparison for identification of pigmentation-related genes in different color varieties of Siamese fighting fish Betta splendens

Transcriptome comparison was performed to identify genes expressed in skin, muscle and tails of mono-color (Red, Blue, Black, White and Yellow), bi-color (Cambodian) and multi-color (Marble) varieties of Siamese fighting fish Betta splendens. In total, 163,140 unigenes covering 26.348 Gb were found. Of these, 93,899 (57.55 %) unigenes significantly matched at least one database. In total, 5039 differentially expressed genes (DEGs) were found where 2415 genes (47.93 %) showed higher expression and 2624 genes (52.07 %) showed lower expression for all pairwise comparisons. DEGs between paired color varieties were 133-443. Of these, 38-220 genes were more highly expressed while 37-280 genes were more lowly expressed relative to the compared varieties. A total of 897 sequences (148 genes) significantly matched pigmentation-related genes of Danio rerio (E-value < 1e-06). Of these, 19 DEGs were identified. Examples are tyrosinase-related protein 1a (BsTyrp1a), epidermal growth factor receptor (BsEgfr) and neurofibronin 1a (BsNf1a). Moreover, 711,123 SNPs were identified and 1365 of these were located in pigmentation-related genes. Interestingly, an A > C₄₇₄ SNP in the gene BsTrpm7 and an indel (position 3571) in the BsItgb1a gene were found only in Cambodian. A C > T₂₅₂₀ SNP in BsFzd4 and 10 of 11 SNPs in BsTyrp1a were found only in Black. Different expression levels (P < 0.05) were found for tyrosinase (BsTyr), BsTyrp1a, BsNf1a and BsEgf1 among skin, body muscle and tails of the same variety and among the same tissues of different varieties (Red, Green, Blue, Black, Cambodian and Multi-colors, N = 5 each).

The evolution and diversification of oakleaf butterflies

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

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.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Speciation and adaptation research meets genome editing

Understanding the genetic basis of reproductive isolation and adaptive traits in natural populations is one of the fundamental goals in evolutionary biology. Genome editing technologies based on CRISPR-Cas systems and site-specific recombinases have enabled us to modify a targeted genomic region as desired and thus to conduct functional analyses of target loci, genes and mutations even in non-conventional model organisms. Here, we review the technical properties of genome editing techniques by classifying them into the following applications: targeted gene knock-out for investigating causative gene functions, targeted gene knock-in of marker genes for visualizing expression patterns and protein functions, precise gene replacement for identifying causative alleles and mutations, and targeted chromosomal rearrangement for investigating the functional roles of chromosomal structural variations. We describe examples of their application to demonstrate functional analysis of naturally occurring genetic variations and discuss how these technologies can be applied to speciation and adaptation research. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

An enhancer of Agouti contributes to parallel evolution of cryptically colored beach mice

Identifying the genetic basis of repeatedly evolved traits provides a way to reconstruct their evolutionary history and ultimately investigate the predictability of evolution. Here, we focus on the oldfield mouse (Peromyscus polionotus), which occurs in the southeastern United States, where it exhibits considerable color variation. Dorsal coats range from dark brown in mainland mice to near white in mice inhabiting sandy beaches; this light pelage has evolved independently on Florida's Gulf and Atlantic coasts as camouflage from predators. To facilitate genomic analyses, we first generated a chromosome-level genome assembly of Peromyscus polionotus subgriseus. Next, in a uniquely variable mainland population (Peromyscus polionotus albifrons), we scored 23 pigment traits and performed targeted resequencing in 168 mice. We find that pigment variation is strongly associated with an ∼2-kb region ∼5 kb upstream of the Agouti signaling protein coding region. Using a reporter-gene assay, we demonstrate that this regulatory region contains an enhancer that drives expression in the dermis of mouse embryos during the establishment of pigment prepatterns. Moreover, extended tracts of homozygosity in this Agouti region indicate that the light allele experienced recent and strong positive selection. Notably, this same light allele appears fixed in both Gulf and Atlantic coast beach mice, despite these populations being separated by >1,000 km. Together, our results suggest that this identified Agouti enhancer allele has been maintained in mainland populations as standing genetic variation and from there, has spread to and been selected in two independent beach mouse lineages, thereby facilitating their rapid and parallel evolution.

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.

Genomic changes underlying repeated niche shifts in an adaptive radiation

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space.

Parallel evolutionary paths of rove beetle myrmecophiles: replaying a deep-time tape of life

The rise of ants over the past ~100 million years reshaped the biosphere, presenting ecological challenges for many organisms, but also opportunities. No insect group has been so adept at exploiting niches inside ant colonies as the rove beetles (Staphylinidae) - a global clade of>64,000 predominantly free-living predators from which numerous socially parasitic 'myrmecophile' lineages have emerged. Myrmecophilous staphylinids are specialized for colony life through changes in behavior, chemistry, anatomy, and life history that are often strikingly convergent, and hence potentially adaptive for this symbiotic way of life. Here, we examine how the interplay between ecological pressures and molecular, cellular, and neurobiological mechanisms shape the evolutionary trajectories of symbiotic lineages in this ancient, convergent system.

Phylogenomic analyses show repeated evolution of hypertrophied lips among Lake Malawi cichlid fishes

Cichlid fishes have repeatedly evolved an astounding diversity of trophic morphologies. For example, hypertrophied lips have evolved multiple times in both African and Neotropical cichlids and could have even evolved convergently within single species assemblages such as African Lake Malawi cichlids. However, the extremely high diversification rate in Lake Malawi cichlids and extensive potential for hybridization has cast doubt on whether even genome-level phylogenetic reconstructions could delineate if these types of adaptations have evolved once or multiple times. To examine the evolution of this iconic trait using protein-coding and noncoding single nucleotide polymorphisms (SNPs), we analyzed the genomes of 86 Lake Malawi cichlid species, including 33 de novo resequenced genomes. Surprisingly, genome-wide protein-coding SNPs exhibited enough phylogenetic informativeness to reconstruct interspecific and intraspecific relationships of hypertrophied lip cichlids, although noncoding SNPs provided better support. However, thinning of noncoding SNPs indicated most discrepancies come from the relatively smaller number of protein-coding sites and not from fundamental differences in their phylogenetic informativeness. Both coding and noncoding reconstructions showed that several “sand-dwelling” hypertrophied lip species, sampled intraspecifically, form a clade interspersed with a few other nonhypertrophied lip lineages. We also recovered Abactochromis labrosus within the rock-dwelling “mbuna” lineage, starkly contrasting with the affinities of other hypertrophied lip taxa found in the largely sand-dwelling “nonmbuna” component of this radiation. Comparative analyses coupled with tests for introgression indicate there is no widespread introgression between the hypertrophied lip lineages and taken together suggest this trophic phenotype has likely evolved at least twice independently within-lake Malawi.

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.

The repeated evolution of stripe patterns is correlated with body morphology in the adaptive radiations of East African cichlid fishes

Color patterns are often linked to the behavioral and morphological characteristics of an animal, contributing to the effectiveness of such patterns as antipredatory strategies. Species-rich adaptive radiations, such as the freshwater fish family Cichlidae, provide an exciting opportunity to study trait correlations at a macroevolutionary scale. Cichlids are also well known for their diversity and repeated evolution of color patterns and body morphology. To study the evolutionary dynamics between color patterns and body morphology, we used an extensive dataset of 461 species. A phylogenetic supertree of these species shows that stripe patterns evolved ~70 times independently and were lost again ~30 times. Moreover, stripe patterns show strong signs of correlated evolution with body elongation, suggesting that the stripes' effectiveness as antipredatory strategy might differ depending on the body shape. Using pedigree-based analyses, we show that stripes and body elongation segregate independently, indicating that the two traits are not genetically linked. Their correlation in nature is therefore likely maintained by correlational selection. Lastly, by performing a mate preference assay using a striped CRISPR-Cas9 mutant of a nonstriped species, we show that females do not differentiate between striped CRISPR mutant males and nonstriped wild-type males, suggesting that these patterns might be less important for species recognition and mate choice. In summary, our study suggests that the massive rates of repeated evolution of stripe patterns are shaped by correlational selection with body elongation, but not by sexual selection.

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

Polymorphisms have fascinated biologists for a long time, but their genetic underpinnings often remain elusive. Here, we aim to uncover the genetic basis of the gold/dark polymorphism that is eponymous of Midas cichlid fish (Amphilophus spp.) adaptive radiations in Nicaraguan crater lakes. While most Midas cichlids are of the melanic "dark morph", about 10% of individuals lose their melanic pigmentation during their ontogeny and transition into a conspicuous "gold morph". Using a new haplotype-resolved long-read assembly we discover an 8.2 kb, transposon-derived inverted repeat in an intron of an undescribed gene, which we term goldentouch in reference to the Greek myth of King Midas. The gene goldentouch is differentially expressed between morphs, presumably due to structural implications of inverted repeats in both DNA and/or RNA (cruciform and hairpin formation). The near-perfect association of the insertion with the phenotype across independent populations suggests that it likely underlies this trans-specific, stable polymorphism.

Concerted variation in melanogenesis genes underlies emergent patterning of plumage in capuchino seedeaters

Coloration traits are central to animal communication; they often govern mate choice, promote reproductive isolation and catalyse speciation. Specific genetic changes can cause variation in coloration, yet far less is known about how overall coloration patterns-which involve combinations of multiple colour patches across the body-can arise and are genomically controlled. We performed genome-wide association analyses to link genomic changes to variation in melanin (eumelanin and pheomelanin) concentration in feathers from different body parts in the capuchino seedeaters, an avian radiation with diverse colour patterns despite remarkably low genetic differentiation across species. Cross-species colour variation in each plumage patch is associated with unique combinations of variants at a few genomic regions, which include mostly non-coding (presumably regulatory) areas close to known pigmentation genes. Genotype-phenotype associations can vary depending on patch colour and are stronger for eumelanin pigmentation, suggesting eumelanin production is tightly regulated. Although some genes are involved in colour variation in multiple patches, in some cases, the SNPs associated with colour changes in different patches segregate spatially. These results suggest that coloration patterning in capuchinos is generated by the modular combination of variants that regulate multiple melanogenesis genes, a mechanism that may have promoted this rapid radiation.

OUP accepted manuscript

Whole genome sequences are beginning to revolutionize our understanding of phylogenetic relationships. Yet, even whole genome sequences can fail to resolve the evolutionary history of the most rapidly radiating lineages, where incomplete lineage sorting, standing genetic variation, introgression, and other factors obscure the phylogenetic history of the group. To overcome such challenges, one emerging strategy is to integrate results across different methods. Most such approaches have been implemented on reduced representation genomic data sets, but whole genomes should provide the maximum possible evidence approach. Here, we test the ability of single nucleotide polymorphisms extracted from whole genome resequencing data, implemented in an integrative genomic approach, to resolve key nodes in the phylogeny of the mbuna, rock-dwelling cichlid fishes of Lake Malaŵi, which epitomize the phylogenetic intractability that often accompanies explosive lineage diversification. This monophyletic radiation has diversified at an unparalleled rate into several hundred species in less than 2 million years. Using an array of phylogenomic methods, we consistently recovered four major clades of mbuna, but a large basal polytomy among them. Although introgression between clades apparently contributed to the challenge of phylogenetic reconstruction, reduction of the data set to nonintrogressed sites still did not help to resolve the basal polytomy. On the other hand, relationships among six congeneric species pairs were resolved without ambiguity, even in one case where existing data led us to predict that resolution would be difficult. We conclude that the bursts of diversification at the earliest stages of the mbuna radiation may be phylogenetically unresolvable, but other regions of the tree are phylogenetically clearly supported. Integration of multiple phylogenomic approaches will continue to increase confidence in relationships inferred from these and other whole-genome data sets. [Incomplete lineage sorting; introgression; linkage disequilibrium; multispecies coalescence; rapid radiation; soft polytomy.]

Studies of Turing pattern formation in zebrafish skin

Skin patterns are the first example of the existence of Turing patterns in living organisms. Extensive research on zebrafish, a model organism with stripes on its skin, has revealed the principles of pattern formation at the molecular and cellular levels. Surprisingly, although the networks of cell-cell interactions have been observed to satisfy the 'short-range activation and long-range inhibition' prerequisites for Turing pattern formation, numerous individual reactions were not envisioned based on the classical reaction-diffusion model. For example, in real skin, it is not an alteration in concentrations of chemicals, but autonomous migration and proliferation of pigment cells that establish patterns, and cell-cell interactions are mediated via direct contact through cell protrusions. Therefore, the classical reaction-diffusion mechanism cannot be used as it is for modelling skin pattern formation. Various studies are underway to adapt mathematical models to the experimental findings on research into skin patterns, and the purpose of this review is to organize and present them. These novel theoretical methods could be applied to autonomous pattern formation phenomena other than skin patterns. This article is part of the theme issue 'Recent progress and open frontiers in Turing's theory of morphogenesis'.

Dual function and associated costs of a highly exaggerated trait in a cichlid fish

Exaggerated secondary sexual characteristics are apparently costly and seem to defy natural selection. This conundrum promoted the theory of sexual selection. Accordingly, exaggerated secondary sexual characteristics might be ornaments on which female choice is based and/or armaments used during male-male competition. Males of many cichlid fish species, including the adaptive radiation of Nicaraguan Midas cichlids, develop a highly exaggerated nuchal hump, which is thought to be a sexually selected trait. To test this hypothesis, we conducted a series of behavioral assays in F2 hybrids obtained from crossing a species with a relatively small hump and one with an exaggerated hump. Mate-choice experiments showed a clear female preference for males with large humps. In an open-choice experiment with limited territories, couples including large humped males were more successful in acquiring these territories. Therefore, nuchal humps appear to serve dual functions as an ornament for attracting mates and as an armament for direct contest with rivals. Although being beneficial in terms of sexual selection, this trait also imposes fitness costs on males possessing disproportionally large nuchal humps since they exhibit decreased endurance and increased energetic costs when swimming. We conclude that these costs illustrate trade-offs associated with large hump size between sexual and natural selection, which causes the latter to limit further exaggeration of this spectacular male trait.

Can gene-inactivating mutations lead to evolutionary novelty?

Evolutionary novelty is difficult to define. It typically involves shifts in organismal or biochemical phenotypes that can be seen as qualitative as well as quantitative changes. In laboratory-based experimental evolution of novel phenotypes and the human domestication of crops, the majority of the mutations that lead to adaptation are loss-of-function mutations that impair or eliminate the function of genes rather than gain-of-function mutations that increase or qualitatively alter the function of proteins. Here, I speculate that easier access to loss-of-function mutations has led them to play a major role in the adaptive radiations that occur when populations have access to many unoccupied ecological niches. I discuss five possible objections to this claim: that genes can only survive if they confer benefits to the organisms that bear them, antagonistic pleiotropy, the importance of pre-existing genetic variation in populations, the danger that adaptation by breaking genes will, over long times, cause organisms to run out of genes, and the recessive nature of most loss-of-function mutations.

The application of genome editing technology in fish

The advent and development of genome editing technology has opened up the possibility of directly targeting and modifying genomic sequences in the field of life sciences with rapid developments occurring in the last decade. As a powerful tool to decipher genome data at the molecular biology level, genome editing technology has made important contributions to elucidating many biological problems. Currently, the three most widely used genome editing technologies include: zinc finger nucleases (ZFN), transcription activator like effector nucleases (TALEN), and clustered regularly interspaced short palindromic repeats (CRISPR). Researchers are still striving to create simpler, more efficient, and accurate techniques, such as engineered base editors and new CRISPR/Cas systems, to improve editing efficiency and reduce off-target rate, as well as a near-PAMless SpCas9 variants to expand the scope of genome editing. As one of the important animal protein sources, fish has significant economic value in aquaculture. In addition, fish is indispensable for research as it serves as the evolutionary link between invertebrates and higher vertebrates. Consequently, genome editing technologies were applied extensively in various fish species for basic functional studies as well as applied research in aquaculture. In this review, we focus on the application of genome editing technologies in fish species detailing growth, gender, and pigmentation traits. In addition, we have focused on the construction of a zebrafish (Danio rerio) disease model and high-throughput screening of functional genes. Finally, we provide some of the future perspectives of this technology.

Of bars and stripes: A Malawi cichlid hybrid cross provides insights into genetic modularity and evolution of modifier loci underlying colour pattern diversification

Understanding the origins of phenotypic diversity among closely related species remains an important largely unsolved question in evolutionary biology. With over 800 species, Lake Malawi haplochromine cichlid fishes are a prominent example of extremely fast evolution of diversity including variation in colouration. Previously, a single major effect gene, agrp2 (asip2b), has been linked to evolutionary losses and gains of horizontal stripe patterns in cichlids, but it remains unknown what causes more fine-scale variation in the number and continuity of the stripes. Also, the genetic basis of the most common colour pattern in African cichlids, vertical bars, and potential interactions between the two colour patterns remain unknown. Based on a hybrid cross of the horizontally striped Lake Malawi cichlid Pseudotropheus cyaneorhabdos and the vertically barred species Chindongo demasoni we investigated the genetic basis of both colour patterns. The distribution of phenotypes in the F₂  generation of the cross indicates that horizontal stripes and vertical bars are independently inherited patterns that are caused by two sets of genetic modules. While horizontal stripes are largely controlled by few major effect loci, vertical bars are a highly polygenic trait. Horizontal stripes show substantial variation in the F₂  generation that, interestingly, resemble naturally occurring phenotypes found in other Lake Malawi cichlid species. Quantitative trait loci (QTL) mapping of this cross reveals known (agrp2) and unknown loci underlying horizontal stripe patterns. These findings provide novel insights into the incremental fine-tuning of an adaptive trait that diversified through the evolution of additional modifier loci.

Manipulation of the Tyrosinase gene permits improved CRISPR/Cas editing and neural imaging in cichlid fish

Direct tests of gene function have historically been performed in a limited number of model organisms. The CRISPR/Cas system is species-agnostic, offering the ability to manipulate genes in a range of models, enabling insights into evolution, development, and physiology. Astatotilapia burtoni, a cichlid fish from the rivers and shoreline around Lake Tanganyika, has been extensively studied in the laboratory to understand evolution and the neural control of behavior. Here we develop protocols for the creation of CRISPR-edited cichlids and create a broadly useful mutant line. By manipulating the Tyrosinase gene, which is necessary for eumelanin pigment production, we describe a fast and reliable approach to quantify and optimize gene editing efficiency. Tyrosinase mutants also remove a major obstruction to imaging, enabling visualization of subdermal structures and fluorophores in situ. These protocols will facilitate broad application of CRISPR/Cas9 to studies of cichlids as well as other non-traditional model aquatic species.

Sensory-based quantification of male colour patterns in Trinidadian guppies reveals no support for parallel phenotypic evolution in multivariate trait space

Parallel evolution, in which independent populations evolve along similar phenotypic trajectories, offers insights into the repeatability of adaptive evolution. Here, we revisit a classic example of parallelism, that of repeated evolution of brighter males in the Trinidadian guppy (Poecilia reticulata). In guppies, colonisation of low predation habitats is associated with emergence of 'more colourful' phenotypes since predator-induced viability selection for crypsis weakens while sexual selection by female preference for conspicuousness remains strong. Our study differs from previous investigations in three respects. First, we adopted a multivariate phenotyping approach to characterise parallelism in multitrait space. Second, we used ecologically-relevant colour traits defined by the visual systems of the two selective agents (i.e., guppy, predatory cichlid). Third, we estimated population genetic structure to test for adaptive (parallel) evolution against a model of neutral phenotypic divergence. We find strong phenotypic differentiation that is inconsistent with a neutral model but very limited support for the predicted pattern of greater conspicuousness at low predation. Effects of predation regime on each trait were in the expected direction, but weak, largely nonsignificant, and explained little among-population variation. In multitrait space, phenotypic trajectories of lineages colonising low from high predation regimes were not parallel. Our results are consistent with reduced predation risk facilitating adaptive differentiation, potentially by female choice, but suggest that this proceeds in independent directions of multitrait space across lineages. Pool-sequencing data also revealed SNPs showing greater differentiation than expected under neutrality, among which some are found in genes contributing to colour pattern variation, presenting opportunities for future genetic study.