Closing the genotype-phenotype gap: emerging technologies for evolutionary genetics in ecological model vertebrate systems

Beyond genes-for-behaviour: The potential for genomics to resolve long-standing questions in avian brood parasitism

Behavioural ecology by definition of its founding 'Tinbergian framework' is an integrative field, however, it lags behind in incorporating genomic methods. 'Finding the gene/s for a behaviour' is still rarely feasible or cost-effective in the wild but as we show here, genomic data can be used to address broader questions. Here we use avian brood parasitism, a model system in behavioural ecology as a case study to highlight how behavioural ecologists could use the full potential of state-of-the-art genomic tools. Brood parasite-host interactions are one of the most easily observable and amenable natural laboratories of antagonistic coevolution, and as such have intrigued evolutionary biologists for decades. Using worked examples, we demonstrate how genomic data can be used to study the causes and mechanisms of (co)evolutionary adaptation and answer three key questions for the field: (i) Where and when should brood parasitism evolve?, (ii) When and how should hosts defend?, and (iii) Will coevolution persist with ecological change? In doing so, we discuss how behavioural and molecular ecologists can collaborate to integrate Tinbergen's questions and achieve the coherent science that he promoted to solve the mysteries of nature.

Mechanisms Underlying the Formation and Evolution of Vertebrate Color Patterns

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

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

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

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

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

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.

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.

Host-pathogen protein-nucleic acid interactions: A comprehensive review

Recognition of pathogen-derived nucleic acids by host cells is an effective host strategy to detect pathogenic invasion and trigger immune responses. In the context of pathogen-specific pharmacology, there is a growing interest in mapping the interactions between pathogen-derived nucleic acids and host proteins. Insight into the principles of the structural and immunological mechanisms underlying such interactions and their roles in host defense is necessary to guide therapeutic intervention. Here, we discuss the newest advances in studies of molecular interactions involving pathogen nucleic acids and host factors, including their drug design, molecular structure and specific patterns. We observed that two groups of nucleic acid recognizing molecules, Toll-like receptors (TLRs) and the cytoplasmic retinoic acid-inducible gene (RIG)-I-like receptors (RLRs) form the backbone of host responses to pathogen nucleic acids, with additional support provided by absent in melanoma 2 (AIM2) and DNA-dependent activator of Interferons (IFNs)-regulatory factors (DAI) like cytosolic activity. We review the structural, immunological, and other biological aspects of these representative groups of molecules, especially in terms of their target specificity and affinity and challenges in leveraging host-pathogen protein-nucleic acid interactions (HP-PNI) in drug discovery.

Using phenotypic plasticity to understand the structure and evolution of the genotype-phenotype map

Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.

A data-independent acquisition (DIA) assay library for quantitation of environmental effects on the kidney proteome of Oreochromis niloticus

Interactions of organisms with their environment are complex and environmental regulation at different levels of biological organization is often nonlinear. Therefore, the genotype to phenotype continuum requires study at multiple levels of organization. While studies of transcriptome regulation are now common for many species, quantitative studies of environmental effects on proteomes are needed. Here we report the generation of a data-independent acquisition (DIA) assay library that enables simultaneous targeted proteomics of thousands of Oreochromis niloticus kidney proteins using a label- and gel-free workflow that is well suited for ecologically relevant field samples. We demonstrate the usefulness of this DIA assay library by discerning environmental effects on the kidney proteome of O. niloticus. Moreover, we demonstrate that the DIA assay library approach generates data that are complimentary rather than redundant to transcriptomic data. Transcript and protein abundance differences in kidneys of tilapia acclimated to freshwater and brackish water (25 g/kg) were correlated for 2114 unique genes. A high degree of non-linearity in salinity-dependent regulation of transcriptomes and proteomes was revealed suggesting that the regulation of O. niloticus renal function by environmental salinity relies heavily on post-transcriptional mechanisms. The application of functional enrichment analyses using STRING and KEGG to DIA assay data sets is demonstrated by identifying myo-inositol metabolism, antioxidant and xenobiotic functions, and signalling mechanisms as key elements controlled by salinity in tilapia kidneys. The DIA assay library resource presented here can be adopted for other tissues and other organisms to study proteome dynamics during changing ecological contexts.

Few Fixed Variants between Trophic Specialist Pupfish Species Reveal Candidate Cis-Regulatory Alleles Underlying Rapid Craniofacial Divergence

Investigating closely related species that rapidly evolved divergent feeding morphology is a powerful approach to identify genetic variation underlying variation in complex traits. This can also lead to the discovery of novel candidate genes influencing natural and clinical variation in human craniofacial phenotypes. We combined whole-genome resequencing of 258 individuals with 50 transcriptomes to identify candidate cis-acting genetic variation underlying rapidly evolving craniofacial phenotypes within an adaptive radiation of Cyprinodon pupfishes. This radiation consists of a dietary generalist species and two derived trophic niche specialists-a molluscivore and a scale-eating species. Despite extensive morphological divergence, these species only diverged 10 kya and produce fertile hybrids in the laboratory. Out of 9.3 million genome-wide SNPs and 80,012 structural variants, we found very few alleles fixed between species-only 157 SNPs and 87 deletions. Comparing gene expression across 38 purebred F1 offspring sampled at three early developmental stages, we identified 17 fixed variants within 10 kb of 12 genes that were highly differentially expressed between species. By measuring allele-specific expression in F1 hybrids from multiple crosses, we found that the majority of expression divergence between species was explained by trans-regulatory mechanisms. We also found strong evidence for two cis-regulatory alleles affecting expression divergence of two genes with putative effects on skeletal development (dync2li1 and pycr3). These results suggest that SNPs and structural variants contribute to the evolution of novel traits and highlight the utility of the San Salvador Island pupfish system as an evolutionary model for craniofacial development.

Different Sources of Allelic Variation Drove Repeated Color Pattern Divergence in Cichlid Fishes

The adaptive radiations of East African cichlid fish in the Great Lakes Victoria, Malawi, and Tanganyika are well known for their diversity and repeatedly evolved phenotypes. Convergent evolution of melanic horizontal stripes has been linked to a single locus harboring the gene agouti-related peptide 2 (agrp2). However, where and when the causal variants underlying this trait evolved and how they drove phenotypic divergence remained unknown. To test the alternative hypotheses of standing genetic variation versus de novo mutations (independently originating in each radiation), we searched for shared signals of genomic divergence at the agrp2 locus. Although we discovered similar signatures of differentiation at the locus level, the haplotypes associated with stripe patterns are surprisingly different. In Lake Malawi, the highest associated alleles are located within and close to the 5' untranslated region of agrp2 and likely evolved through recent de novo mutations. In the younger Lake Victoria radiation, stripes are associated with two intronic regions overlapping with a previously reported cis-regulatory interval. The origin of these segregating haplotypes predates the Lake Victoria radiation because they are also found in more basal riverine and Lake Kivu species. This suggests that both segregating haplotypes were present as standing genetic variation at the onset of the Lake Victoria adaptive radiation with its more than 500 species and drove phenotypic divergence within the species flock. Therefore, both new (Lake Malawi) and ancient (Lake Victoria) allelic variation at the same locus fueled rapid and convergent phenotypic evolution.

De novo assembly of a chromosome-scale reference genome for the northern flicker Colaptes auratus

The northern flicker, Colaptes auratus, is a widely distributed North American woodpecker and a long-standing focal species for the study of ecology, behavior, phenotypic differentiation, and hybridization. We present here a highly contiguous de novo genome assembly of C. auratus, the first such assembly for the species and the first published chromosome-level assembly for woodpeckers (Picidae). The assembly was generated using a combination of short-read Chromium 10× and long-read PacBio sequencing, and further scaffolded with chromatin conformation capture (Hi-C) reads. The resulting genome assembly is 1.378 Gb in size, with a scaffold N50 of 11  and a scaffold L50 of 43.948 Mb. This assembly contains 87.4-91.7% of genes present across four sets of universal single-copy orthologs found in tetrapods and birds. We annotated the assembly both for genes and repetitive content, identifying 18,745 genes and a prevalence of ∼28.0% repetitive elements. Lastly, we used fourfold degenerate sites from neutrally evolving genes to estimate a mutation rate for C. auratus, which we estimated to be 4.007 × 10-9 substitutions/site/year, about 1.5× times faster than an earlier mutation rate estimate of the family. The highly contiguous assembly and annotations we report will serve as a resource for future studies on the genomics of C. auratus and comparative evolution of woodpeckers.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Predicting gene regulatory regions with a convolutional neural network for processing double-strand genome sequence information

With advances in sequencing technology, a vast amount of genomic sequence information has become available. However, annotating biological functions particularly of non-protein-coding regions in genome sequences without experiments is still a challenging task. Recently deep learning-based methods were shown to have the ability to predict gene regulatory regions from genome sequences, promising to aid the interpretation of genomic sequence data. Here, we report an improvement of the prediction accuracy for gene regulatory regions by using the design of convolution layers that efficiently process genomic sequence information, and developed a software, DeepGMAP, to train and compare different deep learning-based models (https://github.com/koonimaru/DeepGMAP). First, we demonstrate that our convolution layers, termed forward- and reverse-sequence scan (FRSS) layers, integrate both forward and reverse strand information, and enhance the power to predict gene regulatory regions. Second, we assessed previous studies and identified problems associated with data structures that caused overfitting. Finally, we introduce visualization methods to examine what the program learned. Together, our FRSS layers improve the prediction accuracy for gene regulatory regions.

Prospects for sociogenomics in avian cooperative breeding and parental care

For the last 40 years, the study of cooperative breeding (CB) in birds has proceeded primarily in the context of discovering the ecological, geographical, and behavioral drivers of helping. The advent of molecular tools in the early 1990s assisted in clarifying the relatedness of helpers to those helped, in some cases, confirming predictions of kin selection theory. Methods for genome-wide analysis of sequence variation, gene expression, and epigenetics promise to add new dimensions to our understanding of avian CB, primarily in the area of molecular and developmental correlates of delayed breeding and dispersal, as well as the ontogeny of achieving parental status in nature. Here, we outline key ways in which modern -omics approaches, in particular genome sequencing, transcriptomics, and epigenetic profiling such as ATAC-seq, can be used to add a new level of analysis of avian CB. Building on recent and ongoing studies of avian social behavior and sociogenomics, we review how high-throughput sequencing of a focal species or clade can provide a robust foundation for downstream, context-dependent destructive and non-destructive sampling of specific tissues or physiological states in the field for analysis of gene expression and epigenetics. -Omics approaches have the potential to inform not only studies of the diversification of CB over evolutionary time, but real-time analyses of behavioral interactions in the field or lab. Sociogenomics of birds represents a new branch in the network of methods used to study CB, and can help clarify ways in which the different levels of analysis of CB ultimately interact in novel and unexpected ways.

Hierarchically embedded interaction networks represent a missing link in the study of behavioral and community ecology

Because genes and phenotypes are embedded within individuals, and individuals within populations, interactions within one level of biological organization are inherently linked to interactors at others. Here, we expand the network paradigm to consider that nodes can be embedded within other nodes, and connections (edges) between nodes at one level of organization form "bridges" for connections between nodes embedded within them. Such hierarchically embedded networks highlight two central properties of biological systems: 1) processes occurring across multiple levels of organization shape connections among biological units at any given level of organization and 2) ecological effects occurring at a given level of organization can propagate up or down to additional levels. Explicitly considering the embedded structure of evolutionary and ecological networks can capture otherwise hidden feedbacks and generate new insights into key biological phenomena, ultimately promoting a broader understanding of interactions in evolutionary theory.

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

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

De Novo Assembly of a High-Quality Reference Genome for the Horned Lark (Eremophila alpestris)

The Horned Lark (Eremophila alpestris) is a small songbird that exhibits remarkable geographic variation in appearance and habitat across an expansive distribution. While E. alpestris has been the focus of many ecological and evolutionary studies, we still lack a highly contiguous genome assembly for the Horned Lark and related taxa (Alaudidae). Here, we present CLO_EAlp_1.0, a highly contiguous assembly for E. alpestris generated from a blood sample of a wild, male bird captured in the Altiplano Cundiboyacense of Colombia. By combining short-insert and mate-pair libraries with the ALLPATHS-LG genome assembly pipeline, we generated a 1.04 Gb assembly comprised of 2713 scaffolds, with a largest scaffold size of 31.81 Mb, a scaffold N50 of 9.42 Mb, and a scaffold L50 of 30. These scaffolds were assembled from 23685 contigs, with a largest contig size of 1.69 Mb, a contig N50 of 193.81 kb, and a contig L50 of 1429. Our assembly pipeline also produced a single mitochondrial DNA contig of 14.00 kb. After polishing the genome, we identified 94.5% of single-copy gene orthologs from an Aves data set and 97.7% of single-copy gene orthologs from a vertebrata data set, which further demonstrates the high quality of our assembly. We anticipate that this genomic resource will be useful to the broader ornithological community and those interested in studying the evolutionary history and ecological interactions of larks, which comprise a widespread, yet understudied lineage of songbirds.

Utilizing the blind cavefish Astyanax mexicanus to understand the genetic basis of behavioral evolution

Colonization of novel habitats often results in the evolution of diverse behaviors. Comparisons between individuals from closely related populations that have evolved divergent behaviors in different environments can be used to investigate behavioral evolution. However, until recently, functionally connecting genotypes to behavioral phenotypes in these evolutionarily relevant organisms has been difficult. The development of gene editing tools will facilitate functional genetic analysis of genotype–phenotype connections in virtually any organism, and has the potential to significantly transform the field of behavioral genetics when applied to ecologically and evolutionarily relevant organisms. The blind cavefish Astyanax mexicanus provides a remarkable example of evolution associated with colonization of a novel habitat. These fish consist of a single species that includes sighted surface fish that inhabit the rivers of Mexico and southern Texas and at least 29 populations of blind cavefish from the Sierra Del Abra and Sierra de Guatemala regions of Northeast Mexico. Although eye loss and albinism have been studied extensively in A. mexicanus, derived behavioral traits including sleep loss, alterations in foraging and reduction in social behaviors are now also being investigated in this species to understand the genetic and neural basis of behavioral evolution. Astyanax mexicanus has emerged as a powerful model system for genotype–phenotype mapping because surface and cavefish are interfertile. Further, the molecular basis of repeated trait evolution can be examined in this species, as multiple cave populations have independently evolved the same traits. A sequenced genome and the implementation of gene editing in A. mexicanus provides a platform for gene discovery and identification of the contributions of naturally occurring variation to behaviors. This review describes the current knowledge of behavioral evolution in A. mexicanus with an emphasis on the molecular and genetic underpinnings of evolved behaviors. Multiple avenues of new research that can be pursued using gene editing tools are identified, and how these will enhance our understanding of behavioral evolution is discussed.

Evolution of cellular stress response mechanisms

The cellular stress response (CSR) is pervasive to all domains of life. It has shaped the interaction between organisms and their environment since the origin of the first cell. Although the CSR has been subject to a myriad of nuanced modifications in the various branches of life present today, its core features remain preserved. The scientific literature covering the CSR is enormous and the broad scope of this brief overview was challenging. However, it is critical to conceptually understand how cells respond to stress in a holistic sense and to point out how fundamental aspects of the CSR framework are integrated. It was necessary to be extremely selective and not feasible to even mention many interesting and important developments in this expansive field. The purpose of this overview is to sketch out general and emerging CSR concepts with an emphasis on the initial cellular strain resulting from stress (macromolecular damage) and the evolutionarily most highly conserved elements of the CSR. Examples emphasize fish and aquatic invertebrates to highlight what is known in organisms beyond mammals, yeast, and other common models. Nonetheless, select pioneering studies using canonical models are also considered and the concepts discussed are applicable to all cells. More detail on important aspects of the CSR in aquatic animals is provided in the accompanying articles of this special issue.

The Cre/Lox System to Assess the Development of the Mouse Brain

Cre-mediated recombination has become a powerful tool to confine gene deletions (conditional knockouts) or overexpression of genes (conditional knockin/overexpression). By spatiotemporal restriction of genetic manipulations, major problems of classical knockouts such as embryonic lethality or pleiotropy can be circumvented. Furthermore, Cre-mediated recombination has broad applications in the analysis of the cellular behavior of subpopulations and cell types as well as for genetic fate mapping. This chapter gives an overview about applications for the Cre/LoxP system and their execution.

Evolutionary Dynamics of Structural Variation at a Key Locus for Color Pattern Diversification in Cichlid Fishes

Color patterns in African cichlid fishes vary spectacularly. Although phylogenetic analysis showed already 30 years ago that many color patterns evolved repeatedly in these adaptive radiations, only recently have we begun to understand the genomic basis of color variation. Horizontal stripe patterns evolved and were lost several times independently across the adaptive radiations of Lake Victoria, Malawi, and Tanganyika and regulatory evolution of agouti-related peptide 2 (agrp2/asip2b) has been linked to this phenotypically labile trait. Here, we asked whether the agrp2 locus exhibits particular characteristics that facilitate divergence in color patterns. Based on comparative genomic analyses, we discovered several recent duplications, insertions, and deletions. Interestingly, one of these events resulted in a tandem duplication of the last exon of agrp2. The duplication likely precedes the East African radiations that started 8-12 Ma, is not fixed within any of the radiations, and is found to vary even within some species. Moreover, we also observed variation in copy number (two to five copies) and secondary loss of the duplication, illustrating a surprising dynamic at this locus that possibly promoted functional divergence of agrp2. Our work suggests that such instances of exon duplications are a neglected mechanism potentially involved in the repeated evolution and diversification that deserves more attention.

Genomics of host-pathogen interactions: challenges and opportunities across ecological and spatiotemporal scales

Evolutionary genomics has recently entered a new era in the study of host-pathogen interactions. A variety of novel genomic techniques has transformed the identification, detection and classification of both hosts and pathogens, allowing a greater resolution that helps decipher their underlying dynamics and provides novel insights into their environmental context. Nevertheless, many challenges to a general understanding of host-pathogen interactions remain, in particular in the synthesis and integration of concepts and findings across a variety of systems and different spatiotemporal and ecological scales. In this perspective we aim to highlight some of the commonalities and complexities across diverse studies of host-pathogen interactions, with a focus on ecological, spatiotemporal variation, and the choice of genomic methods used. We performed a quantitative review of recent literature to investigate links, patterns and potential tradeoffs between the complexity of genomic, ecological and spatiotemporal scales undertaken in individual host-pathogen studies. We found that the majority of studies used whole genome resolution to address their research objectives across a broad range of ecological scales, especially when focusing on the pathogen side of the interaction. Nevertheless, genomic studies conducted in a complex spatiotemporal context are currently rare in the literature. Because processes of host-pathogen interactions can be understood at multiple scales, from molecular-, cellular-, and physiological-scales to the levels of populations and ecosystems, we conclude that a major obstacle for synthesis across diverse host-pathogen systems is that data are collected on widely diverging scales with different degrees of resolution. This disparity not only hampers effective infrastructural organization of the data but also data granularity and accessibility. Comprehensive metadata deposited in association with genomic data in easily accessible databases will allow greater inference across systems in the future, especially when combined with open data standards and practices. The standardization and comparability of such data will facilitate early detection of emerging infectious diseases as well as studies of the impact of anthropogenic stressors, such as climate change, on disease dynamics in humans and wildlife.

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

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

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

The color patterns of African cichlid fishes provide notable examples of phenotypic convergence. Across the more than 1200 East African rift lake species, melanic horizontal stripes have evolved numerous times. We discovered that regulatory changes of the gene agouti-related peptide 2 (agrp2) act as molecular switches controlling this evolutionarily labile phenotype. Reduced agrp2 expression is convergently associated with the presence of stripe patterns across species flocks. However, cis-regulatory mutations are not predictive of stripes across radiations, suggesting independent regulatory mechanisms. Genetic mapping confirms the link between the agrp2 locus and stripe patterns. The crucial role of agrp2 is further supported by a CRISPR-Cas9 knockout that reconstitutes stripes in a nonstriped cichlid. Thus, we unveil how a single gene affects the convergent evolution of a complex color pattern.

Keystone Genes

The keystone species concept is used in ecology to describe individual species with disproportionately large effects on their communities. We extend this idea to the level of genes with disproportionately large effects on ecological processes. Such 'keystone genes' (KGs) would underlie traits involved in species interactions or causing critical biotic and/or abiotic changes that influence emergent community and ecosystem properties. We propose a general framework for how KGs could be identified, while keeping KGs under the umbrella of 'ecologically important genes' (EIGs) that also include categories such as 'foundation genes', 'ecosystem engineering genes', and more. Although likely rare, KGs and other EIGs could dominate certain ecological processes; thus, their discovery and study are relevant for understanding eco-evolutionary dynamics.

Early developmental and allometric patterns in the electric yellow cichlid Labidochromis caeruleus

The electric yellow cichlid Labidochromis caeruleus is a mouth-brooding haplochromine cichlid from Lake Malawi and one of the most popular cichlids in the ornamental fish industry. To investigate the early development of L. caeruleus from hatching until the juvenile stage, we studied its morphological development and allometric growth patterns. In newly-hatched larvae, most organs and body parts were not yet differentiated and continued to develop until 15 days post hatching (dph). The yolk sac was depleted at 13 dph. There was allometric growth, primarily in the anterior and posterior regions of the body, and inflection points when trajectories of allometric growth changed. Head and tail growth was prioritized, suggesting that body parts linked to feeding and swimming behaviour mature earlier than the rest of the body. Additionally, growth patterns revealed that development of organs related to vital functions such as branchial respiration, sensation, exogenous feeding and swimming was prioritized. Comparisons with other African and Neotropical cichlids revealed differences in ontogenetic processes and allometric growth along the anterior-posterior axis as well as variation in developmental timing. These results indicate how early morphological development and ontogenic processes might respond to the distinctive parental care observed in mouth-brooding cichlids.

Tol2 transposon-mediated transgenesis in the Midas cichlid (Amphilophus citrinellus) - towards understanding gene function and regulatory evolution in an ecological model system for rapid phenotypic diversification

Background

The Midas cichlid species complex (Amphilophus spp.) is widely known among evolutionary biologists as a model system for sympatric speciation and adaptive phenotypic divergence within extremely short periods of time (a few hundred generations). The repeated parallel evolution of adaptive phenotypes in this radiation, combined with their near genetic identity, makes them an excellent model for studying phenotypic diversification. While many ecological and evolutionary studies have been performed on Midas cichlids, the molecular basis of specific phenotypes, particularly adaptations, and their underlying coding and cis-regulatory changes have not yet been studied thoroughly.

Results

For the first time in any New World cichlid, we use Tol2 transposon-mediated transgenesis in the Midas cichlid (Amphilophus citrinellus). By adapting existing microinjection protocols, we established an effective protocol for transgenesis in Midas cichlids. Embryos were injected with a Tol2 plasmid construct that drives enhanced green fluorescent protein (eGFP) expression under the control of the ubiquitin promoter. The transgene was successfully integrated into the germline, driving strong ubiquitous expression of eGFP in the first transgenic Midas cichlid line. Additionally, we show transient expression of two further transgenic constructs, ubiquitin::tdTomato and mitfa::eGFP. Transgenesis in Midas cichlids will facilitate further investigation of the genetic basis of species-specific traits, many of which are adaptations.

Conclusion

Transgenesis is a versatile tool not only for studying regulatory elements such as promoters and enhancers, but also for testing gene function through overexpression of allelic gene variants. As such, it is an important first step in establishing the Midas cichlid as a powerful model for studying adaptive coding and non-coding changes in an ecological and evolutionary context.

The Long Journey of Pontine Nuclei Neurons: From Rhombic Lip to Cortico-Ponto-Cerebellar Circuitry

The pontine nuclei (PN) are the largest of the precerebellar nuclei, neuronal assemblies in the hindbrain providing principal input to the cerebellum. The PN are predominantly innervated by the cerebral cortex and project as mossy fibers to the cerebellar hemispheres. Here, we comprehensively review the development of the PN from specification to migration, nucleogenesis and circuit formation. PN neurons originate at the posterior rhombic lip and migrate tangentially crossing several rhombomere derived territories to reach their final position in ventral part of the pons. The developing PN provide a classical example of tangential neuronal migration and a study system for understanding its molecular underpinnings. We anticipate that understanding the mechanisms of PN migration and assembly will also permit a deeper understanding of the molecular and cellular basis of cortico-cerebellar circuit formation and function.

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

Osmolality/salinity-responsive enhancers (OSREs) control induction of osmoprotective genes in euryhaline fish

Fish respond to salinity stress by transcriptional induction of many genes, but the mechanism of their osmotic regulation is unknown. We developed a reporter assay using cells derived from the brain of the tilapia Oreochromis mossambicus (OmB cells) to identify osmolality/salinity-responsive enhancers (OSREs) in the genes of Omossambicus Genomic DNA comprising the regulatory regions of two strongly salinity-induced genes, inositol monophosphatase 1 (IMPA1.1) and myo-inositol phosphate synthase (MIPS), was isolated and analyzed with dual luciferase enhancer trap reporter assays. We identified five sequences (two in IMPA1.1 and three in MIPS) that share a common consensus element (DDKGGAAWWDWWYDNRB), which we named "OSRE1." Additional OSREs that were less effective in conferring salinity-induced trans-activation and do not match the OSRE1 consensus also were identified in both MIPS and IMPA1.1 Although OSRE1 shares homology with the mammalian osmotic-response element/tonicity-responsive enhancer (ORE/TonE) enhancer, the latter is insufficient to confer osmotic induction in fish. Like other enhancers, OSRE1 trans-activates genes independent of orientation. We conclude that OSRE1 is a cis-regulatory element (CRE) that enhances the hyperosmotic induction of osmoregulated genes in fish. Our study also shows that tailored reporter assays developed for OmB cells facilitate the identification of CREs in fish genomes. Knowledge of the OSRE1 motif allows affinity-purification of the corresponding transcription factor and computational approaches for enhancer screening of fish genomes. Moreover, our study enables targeted inactivation of OSRE1 enhancers, a method superior to gene knockout for functional characterization because it confines impairment of gene function to a specific context (salinity stress) and eliminates pitfalls of constitutive gene knockouts (embryonic lethality, developmental compensation).

Evolutionary rescue and local adaptation under different rates of temperature increase: a combined analysis of changes in phenotype expression and genotype frequency in Paramecium microcosms

Evolutionary rescue (ER) occurs when populations, which have declined due to rapid environmental change, recover through genetic adaptation. The success of this process and the evolutionary trajectory of the population strongly depend on the rate of environmental change. Here we investigated how different rates of temperature increase (from 23 to 32 °C) affect population persistence and evolutionary change in experimental microcosms of the protozoan Paramecium caudatum. Consistent with theory on ER, we found that those populations experiencing the slowest rate of temperature increase were the least likely to become extinct and tended to be the best adapted to the new temperature environment. All high-temperature populations were more tolerant to severe heat stress (35, 37 °C), indicating a common mechanism of heat protection. High-temperature populations also had superior growth rates at optimum temperatures, leading to the absence of a pattern of local adaptation to control (23 °C) and high-temperature (32 °C) environments. However, high-temperature populations had reduced growth at low temperatures (5-9 °C), causing a shift in the temperature niche. In part, the observed evolutionary change can be explained by selection from standing variation. Using mitochondrial markers, we found complete divergence between control and high-temperature populations in the frequencies of six initial founder genotypes. Our results confirm basic predictions of ER and illustrate how adaptation to an extreme local environment can produce positive as well as negative correlated responses to selection over the entire range of the ecological niche.

Signalling pathways in trophic skeletal development and morphogenesis: Insights from studies on teleost fish

During the development of the vertebrate feeding apparatus, a variety of complicated cellular and molecular processes participate in the formation and integration of individual skeletal elements. The molecular mechanisms regulating the formation of skeletal primordia and their development into specific morphological structures are tightly controlled by a set of interconnected signalling pathways. Some of these pathways, such as Bmp, Hedgehog, Notch and Wnt, are long known for their pivotal roles in craniofacial skeletogenesis. Studies addressing the functional details of their components and downstream targets, the mechanisms of their interactions with other signals as well as their potential roles in adaptive morphological divergence, are currently attracting considerable attention. An increasing number of signalling pathways that had previously been described in different biological contexts have been shown to be important in the regulation of jaw skeletal development and morphogenesis. In this review, I provide an overview of signalling pathways involved in trophic skeletogenesis emphasizing studies of the most species-rich group of vertebrates, the teleost fish, which through their evolutionary history have undergone repeated episodes of spectacular trophic diversification.

Evolution: tinkering within gene regulatory landscapes

Recently evolved enhancers dominate mammalian gene regulatory landscapes. Mostly exapted from ancestral DNA sequences, many are linked to genes under positive selection. Just as RNA-seq some years ago, unbiased enhancer mapping is on the verge of changing evolutionary research.

Molecular Evolution of the Neural Crest Regulatory Network in Ray-Finned Fish

Gene regulatory networks (GRN) are central to developmental processes. They are composed of transcription factors and signaling molecules orchestrating gene expression modules that tightly regulate the development of organisms. The neural crest (NC) is a multipotent cell population that is considered a key innovation of vertebrates. Its derivatives contribute to shaping the astounding morphological diversity of jaws, teeth, head skeleton, or pigmentation. Here, we study the molecular evolution of the NC GRN by analyzing patterns of molecular divergence for a total of 36 genes in 16 species of bony fishes. Analyses of nonsynonymous to synonymous substitution rate ratios (dN/dS) support patterns of variable selective pressures among genes deployed at different stages of NC development, consistent with the developmental hourglass model. Model-based clustering techniques of sequence features support the notion of extreme conservation of NC-genes across the entire network. Our data show that most genes are under strong purifying selection that is maintained throughout ray-finned fish evolution. Late NC development genes reveal a pattern of increased constraints in more recent lineages. Additionally, seven of the NC-genes showed signs of relaxation of purifying selection in the famously species-rich lineage of cichlid fishes. This suggests that NC genes might have played a role in the adaptive radiation of cichlids by granting flexibility in the development of NC-derived traits-suggesting an important role for NC network architecture during the diversification in vertebrates.

Embryonic and larval development in the Midas cichlid fish species flock (Amphilophus spp.): a new evo-devo model for the investigation of adaptive novelties and species differences

Background

Central American crater lake cichlid fish of the Midas species complex (Amphilophus spp.) are a model system for sympatric speciation and fast ecological diversification and specialization. Midas cichlids have been intensively analyzed from an ecological and morphological perspective. Genomic resources such as transcriptomic and genomic data sets, and a high-quality draft genome are available now. Many ecologically relevant species-specific traits and differences such as pigmentation and cranial morphology arise during development. Detailed descriptions of the early development of the Midas cichlid in particular, will help to investigate the ontogeny of species differences and adaptations.

Results

We describe the embryonic and larval development of the crater lake cichlid, Amphilophus xiloaensis, until seven days after fertilization. Similar to previous studies on teleost development, we describe six periods of embryogenesis - the zygote, cleavage, blastula, gastrula, segmentation, and post-hatching period. Furthermore, we define homologous stages to well-described teleost models such as medaka and zebrafish, as well as other cichlid species such as the Nile tilapia and the South American cichlid Cichlasoma dimerus. Key morphological differences between the embryos of Midas cichlids and other teleosts are highlighted and discussed, including the presence of adhesive glands and different early chromatophore patterns, as well as variation in developmental timing.

Conclusions

The developmental staging of the Midas cichlid will aid researchers in the comparative investigation of teleost ontogenies. It will facilitate comparative developmental biological studies of Neotropical and African cichlid fish in particular. In the past, the species flocks of the African Great Lakes have received the most attention from researchers, but some lineages of the 300–400 species of Central American lakes are fascinating model systems for adaptive radiation and rapid phenotypic evolution. The availability of genetic resources, their status as a model system for evolutionary research, and the possibility to perform functional experiments including transgenesis makes the Midas cichlid complex a very attractive model for evolutionary-developmental research.

Transcriptomics of two evolutionary novelties: how to make a sperm-transfer organ out of an anal fin and a sexually selected "sword" out of a caudal fin

Swords are exaggerated male ornaments of swordtail fishes that have been of great interest to evolutionary biologists ever since Darwin described them in the Descent of Man (1871). They are a novel sexually selected trait derived from modified ventral caudal fin rays and are only found in the genus Xiphophorus. Another phylogenetically more widespread and older male trait is the gonopodium, an intromittent organ found in all poeciliid fishes, that is derived from a modified anal fin. Despite many evolutionary and behavioral studies on both traits, little is known so far about the molecular mechanisms underlying their development. By investigating transcriptomic changes (utilizing a RNA-Seq approach) in response to testosterone treatment in the swordtail fish, Xiphophorus hellerii, we aimed to better understand the architecture of the gene regulatory networks underpinning the development of these two evolutionary novelties. Large numbers of genes with tissue-specific expression patterns were identified. Among the "sword genes" those involved in embryonic organ development, sexual character development and coloration were highly expressed, while in the gonopodium rather more morphogenesis-related genes were found. Interestingly, many genes and genetic pathways are shared between both developing novel traits derived from median fins: the sword and the gonopodium. Our analyses show that a larger set of gene networks was co-opted during the development and evolution of the "older" gonopodium than in the "younger," and morphologically less complex trait, the sword. We provide a catalog of candidate genes for future efforts to dissect the development of those sexually selected exaggerated male traits in swordtails.