Inter-chromosomal coupling between vision and pigmentation genes during genomic divergence

Beyond cyanogenesis: Temperature gradients drive environmental adaptation in North American white clover (Trifolium repens L.)

Species that repeatedly evolve phenotypic clines across environmental gradients have been highlighted as ideal systems for characterizing the genomic basis of local environmental adaptation. However, few studies have assessed the importance of observed phenotypic clines for local adaptation: conspicuous traits that vary clinally may not necessarily be the most critical in determining local fitness. The present study was designed to fill this gap, using a plant species characterized by repeatedly evolved adaptive phenotypic clines. White clover is naturally polymorphic for its chemical defence cyanogenesis (HCN release with tissue damage); climate-associated cyanogenesis clines have evolved throughout its native and introduced range worldwide. We performed landscape genomic analyses on 415 wild genotypes from 43 locations spanning much of the North American species range to assess the relative importance of cyanogenesis loci vs. other genomic factors in local climatic adaptation. We find clear evidence of local adaptation, with temperature-related climatic variables best describing genome-wide differentiation between sampling locations. The same climatic variables are also strongly correlated with cyanogenesis frequencies and gene copy number variations (CNVs) at cyanogenesis loci. However, landscape genomic analyses indicate no significant contribution of cyanogenesis loci to local adaptation. Instead, several genomic regions containing promising candidate genes for plant response to seasonal cues are identified - some of which are shared with previously identified QTLs for locally adaptive fitness traits in North American white clover. Our findings suggest that local adaptation in white clover is likely determined primarily by genes controlling the timing of growth and flowering in response to local seasonal cues. More generally, this work suggests that caution is warranted when considering the importance of conspicuous phenotypic clines as primary determinants of local adaptation.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

Phenotypic and genomic dissection of colour pattern variation in a reef fish radiation

Coral reefs rank among the most diverse species assemblages on Earth. A particularly striking aspect of coral reef communities is the variety of colour patterns displayed by reef fishes. Colour pattern is known to play a central role in the ecology and evolution of reef fishes through, for example, signalling or camouflage. Nevertheless, colour pattern is a complex trait in reef fishes-actually a collection of traits-that is difficult to analyse in a quantitative and standardized way. This is the challenge that we address in this study using the hamlets (Hypoplectrus spp., Serranidae) as a model system. Our approach involves a custom underwater camera system to take orientation- and size-standardized photographs in situ, colour correction, alignment of the fish images with a combination of landmarks and Bézier curves, and principal component analysis on the colour value of each pixel of each aligned fish. This approach identifies the major colour pattern elements that contribute to phenotypic variation in the group. Furthermore, we complement the image analysis with whole-genome sequencing to run a multivariate genome-wide association study for colour pattern variation. This second layer of analysis reveals sharp association peaks along the hamlet genome for each colour pattern element and allows to characterize the phenotypic effect of the single nucleotide polymorphisms that are most strongly associated with colour pattern variation at each association peak. Our results suggest that the diversity of colour patterns displayed by the hamlets is generated by a modular genomic and phenotypic architecture.

Genetics and the Evolution of Prezygotic Isolation

The significance of prezygotic isolation for speciation has been recognized at least since the Modern Synthesis. However, fundamental questions remain. For example, how are genetic associations between traits that contribute to prezygotic isolation maintained? What is the source of genetic variation underlying the evolution of these traits? And how do prezygotic barriers affect patterns of gene flow? We address these questions by reviewing genetic features shared across plants and animals that influence prezygotic isolation. Emerging technologies increasingly enable the identification and functional characterization of the genes involved, allowing us to test established theoretical expectations. Embedding these genes in their developmental context will allow further predictions about what constrains the evolution of prezygotic isolation. Ongoing improvements in statistical and computational tools will reveal how pre- and postzygotic isolation may differ in how they influence gene flow across the genome. Finally, we highlight opportunities for progress by combining theory with appropriate data.

A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics

The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.

A few genetic variants go a long way in differentiating Penstemon species

The integrity of hybridizing species is usually maintained by genome-wide selection or by selection on a few genomic regions. A study published in PLOS Biology finds a different pattern-60 SNPs spread across the genome differentiate a Penstemon species pair.

Deep genome-wide divergences among species in White Cloud Mountain minnow Tanichthys albonubes (Cypriniformes: Tanichthyidae) complex: Conservation and species management implications

Identifying cryptic species is important for the assessments of biodiversity. Further, untangling mechanisms underlying the origins of cryptic species can facilitate our understanding of evolutionary processes. Advancements in genomic approaches for non-model systems have offered unprecedented opportunities to investigate these areas. The White Cloud Mountain minnow (Tanichthys albonubes) is a popular freshwater pet fish worldwide but its wild populations in China are critically endangered. Recent research based on a few molecular markers suggested that this species in fact comprised seven cryptic species, of which six were previously unknown. Here, we tested six of these cryptic species and quantified genomic interspecific divergences between species in the T. albonubes complex by analyzing genome-wide restriction site-associated DNA sequencing (RADseq) data generated from 189 individuals sampled from seven populations (including an outgroup congeneric species, T. micagemmae). We found that six cryptic species previously suggested were well supported by RADseq data. The genetic diversity of each species in the T. albonubes complex was low compared with T. micagemmae and the contemporary effective population sizes (Ne) of each cryptic species were small. Phylogenetic analysis showed seven clades with high support values confirmed with Neighbor-Net trees. The pairwise divergences between species in T. albonubes complex were deep and the highly differentiated loci were evenly distributed across the genome. We proposed that the divergence level of T. albonubes complex is at a late stage of cryptic speciation and lacking gene flow. Our findings provide new insights into cryptic speciation and have important implications for conservation and species management of T. albonubes complex.

Interchromosomal linkage disequilibrium and linked fitness cost loci associated with selection for herbicide resistance

The adaptation of weeds to herbicide is both a significant problem in agriculture and a model of rapid adaptation. However, significant gaps remain in our knowledge of resistance controlled by many loci and the evolutionary factors that influence the maintenance of resistance. Here, using herbicide-resistant populations of the common morning glory (Ipomoea purpurea), we perform a multilevel analysis of the genome and transcriptome to uncover putative loci involved in nontarget-site herbicide resistance (NTSR) and to examine evolutionary forces underlying the maintenance of resistance in natural populations. We found loci involved in herbicide detoxification and stress sensing to be under selection and confirmed that detoxification is responsible for glyphosate (RoundUp) resistance using a functional assay. We identified interchromosomal linkage disequilibrium (ILD) among loci under selection reflecting either historical processes or additive effects leading to the resistance phenotype. We further identified potential fitness cost loci that were strongly linked to resistance alleles, indicating the role of genetic hitchhiking in maintaining the cost. Overall, our work suggests that NTSR glyphosate resistance in I. purpurea is conferred by multiple genes which are potentially maintained through generations via ILD, and that the fitness cost associated with resistance in this species is likely a by-product of genetic hitchhiking.

Meeting Zn Needs during Medaka Eye Development: Nanoscale Visualization of Retina by Expansion Microscopy

Fish eyes require high Zn levels to support their early development. Although numerous studies have been conducted on the nutritional and toxic effects of Zn on the eye, the Zn requirement for retinal cell development is still debatable. Moreover, due to the complexity of the retinal structure, it is difficult to clearly visualize each retinal layer and accurately separate cell morphology in vivo by conventional methods. In the present study, we for the first time have achieved nanoscale imaging of retinal anatomy affected by dietary and waterborne Zn exposure by novel expansion microscopy. We demonstrated that the fish retina showed different developmental strategies in response to dietary and aqueous Zn exposures. Excess dietary Zn produced toxicity to retinal photoreceptor cells, resulting in a reduction in cell number and cell area, and this toxicity became severe with biological development. In contrast, waterborne Zn in the natural environment probably failed to meet the Zn requirements of retinal development. Overall, our results indicated that during early development, the Zn requirement of the fish eyes was sensitive, and oversupplementation led to impaired photoreceptor cell development. Our study has provided new perspectives using the powerful and novel expansion microscopy technique in toxicity assessment, enabling ultra-clear visualization of small but complex organ development.

Evolution of physical linkage between loci controlling ecological traits and mating preferences

Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.

Allorecognition genes drive reproductive isolation in Podospora anserina

Allorecognition, the capacity to discriminate self from conspecific non-self, is a ubiquitous organismal feature typically governed by genes evolving under balancing selection. Here, we show that in the fungus Podospora anserina, allorecognition loci controlling vegetative incompatibility (het genes), define two reproductively isolated groups through pleiotropic effects on sexual compatibility. These two groups emerge from the antagonistic interactions of the unlinked loci het-r (encoding a NOD-like receptor) and het-v (encoding a methyltransferase and an MLKL/HeLo domain protein). Using a combination of genetic and ecological data, supported by simulations, we provide a concrete and molecularly defined example whereby the origin and coexistence of reproductively isolated groups in sympatry is driven by pleiotropic genes under balancing selection.

Modeling allelic diversity of multiparent mapping populations affects detection of quantitative trait loci

The search for quantitative trait loci that explain complex traits such as yield and drought tolerance has been ongoing in all crops. Methods such as biparental quantitative trait loci mapping and genome-wide association studies each have their own advantages and limitations. Multiparent advanced generation intercross populations contain more recombination events and genetic diversity than biparental mapping populations and are better able to estimate effect sizes of rare alleles than association mapping populations. Here, we discuss the results of using a multiparent advanced generation intercross population of doubled haploid maize lines created from 16 diverse founders to perform quantitative trait loci mapping. We compare 3 models that assume bi-allelic, founder, and ancestral haplotype allelic states for quantitative trait loci. The 3 methods have differing power to detect quantitative trait loci for a variety of agronomic traits. Although the founder approach finds the most quantitative trait loci, all methods are able to find unique quantitative trait loci, suggesting that each model has advantages for traits with different genetic architectures. A closer look at a well-characterized flowering time quantitative trait loci, qDTA8, which contains vgt1, highlights the strengths and weaknesses of each method and suggests a potential epistatic interaction. Overall, our results reinforce the importance of considering different approaches to analyzing genotypic datasets, and shows the limitations of binary SNP data for identifying multiallelic quantitative trait loci.

Molecular parallelism in signaling function across different sexually selected ornaments in a warbler

Extravagant ornaments are thought to signal male quality to females choosing mates, but the evidence linking ornament size to male quality is controversial, particularly in cases in which females prefer different ornaments in different populations. Here, we use whole-genome sequencing and transcriptomics to determine the genetic basis of ornament size in two populations of a widespread warbler, the common yellowthroat (Geothlypis trichas). Within a single subspecies, females in a Wisconsin population prefer males with larger black masks as mates, while females in a New York population prefer males with larger yellow bibs. Despite being produced by different pigments in different patches on the body, the size of the ornament preferred by females in each population was linked to numerous genes that function in many of the same core aspects of male quality (e.g., immunity and oxidative balance). These relationships confirm recent hypotheses linking the signaling function of ornaments to male quality. Furthermore, the parallelism in signaling function provides the flexibility for different types of ornaments to be used as signals of similar aspects of male quality. This could facilitate switches in female preference for different ornaments, a potentially important step in the early stages of divergence among populations.

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.

Rapid radiation in a highly diverse marine environment

Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.

Assortative mating and epistatic mating-trait architecture induce complex movement of the crow hybrid zone

Hybrid zones provide a window into the evolutionary processes governing species divergence. Yet, the contribution of mate choice to the temporal and spatial stability of hybrid zones remains poorly explored. Here, we investigate the effects of assortative mating on hybrid-zone dynamics by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and gray-coated hooded crows. In the best-fit model, narrow clines of the two mating-trait loci were maintained by a moderate degree of assortative mating inducing pre- and postzygotic isolation via positive frequency-dependent selection. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favoring gray-coated phenotypes ∼ 1 200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of the genetic mating-trait architecture for spatio-temporal hybrid-zone dynamics.

A widely diverged locus involved in locomotor adaptation in Heliconius butterflies

Heliconius butterflies have undergone adaptive radiation and therefore serve as an excellent system for exploring the continuum of speciation and adaptive evolution. However, there is a long-lasting paradox between their convergent mimetic wing patterns and rapid divergence in speciation. Here, we characterize a locus that consistently displays high divergence among Heliconius butterflies and acts as an introgression hotspot. We further show that this locus contains multiple genes related to locomotion and conserved in Lepidoptera. In light of these findings, we consider that locomotion traits may be under selection, and if these are heritable traits that are selected for, then they might act as species barriers.

Sexual Dichromatism Is Decoupled from Diversification over Deep Time in Fishes

AbstractSexually selected traits have long been thought to drive diversification, but support for this hypothesis has been persistently controversial. In fishes, sexually dimorphic coloration is associated with assortative mating and speciation among closely related species, as shown in classic studies. However, it is unclear whether these results can generalize to explain diversity patterns across ray-finned fishes, which contain the majority of vertebrate species and 96% of fishes. Here, we use phylogenetic approaches to test for an association between sexual dichromatism and diversification rates (speciation minus extinction) in ray-finned fishes. We assembled dichromatism data for 10,898 species, a data set of unprecedented size. We found no difference in diversification rates between monochromatic and dichromatic species when including all ray-finned fishes. However, at lower phylogenetic scales (within orders and families), some intermediate-sized clades did show an effect of dichromatism on diversification. Surprisingly, dichromatism could significantly increase or decrease diversification rates. Moreover, we found no effect in many of the clades initially used to link dichromatism to speciation in fishes (e.g., cichlids) or an effect only at shallow scales (within subclades). Overall, we show how the effects of dichromatism on diversification are highly variable in direction and restricted to certain clades and phylogenetic scales.

Correlational selection in the age of genomics

Ecologists and evolutionary biologists are well aware that natural and sexual selection do not operate on traits in isolation, but instead act on combinations of traits. This long-recognized and pervasive phenomenon is known as multivariate selection, or-in the particular case where it favours correlations between interacting traits-correlational selection. Despite broad acknowledgement of correlational selection, the relevant theory has often been overlooked in genomic research. Here, we discuss theory and empirical findings from ecological, quantitative genetic and genomic research, linking key insights from different fields. Correlational selection can operate on both discrete trait combinations and quantitative characters, with profound implications for genomic architecture, linkage, pleiotropy, evolvability, modularity, phenotypic integration and phenotypic plasticity. We synthesize current knowledge and discuss promising research approaches that will enable us to understand how correlational selection shapes genomic architecture, thereby linking quantitative genetic approaches with emerging genomic methods. We suggest that research on correlational selection has great potential to integrate multiple fields in evolutionary biology, including developmental and functional biology, ecology, quantitative genetics, phenotypic polymorphisms, hybrid zones and speciation processes.

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

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

Drivers of eyespot evolution in coral reef fishes

Evolution via natural selection has continually shaped the coloration of numerous organisms. One coloration of particular importance is the eyespot: a phylogenetically widespread, conspicuous marking that has been shown to effectively reduce predation, often through its resemblance to the eye. Although widely studied, most research has been experimental in nature. We approach eyespots using a comparative phylogenetic framework that is global in scope. Herein, we identify the potential drivers of eyespot evolution in coral reef fishes; essentially the rules that govern their appearance in this group of organisms. We surveyed 2664 reef fish species (42% of all described reef fish species) and found that eyespots are present in approximately one in every 10 species. Most eyespots occur in closely related species and have been present in some families for over 50 million years. Focusing on damselfishes (family: Pomacentridae) as a study group, we reveal that eyespots are rare in planktivorous species, which is likely driven by the predation risk associated with their feeding location. Using a heatmapping technique, we also show that the location of eyespots is fundamentally different in active fishes that swim above the benthos vs. cryptobenthic fishes that rest on the benthos. These location differences may reflect different functions of eyespots among reef fish species.

Phylogenomic Analysis of Concatenated Ultraconserved Elements Reveals the Recent Evolutionary Radiation of the Fairy Wrasses (Teleostei: Labridae: Cirrhilabrus)

The fairy wrasses (genus Cirrhilabrus) are among the most successful of the extant wrasse lineages (Teleostei: Labridae), with their 61 species accounting for nearly 10% of the family. Although species complexes within the genus have been diagnosed on the basis of coloration patterns and synapomorphies, attempts to resolve evolutionary relationships among these groups using molecular and morphological data have largely been unsuccessful. Here we use a phylogenomic approach with a data set comprising 991 ultraconserved elements (UCEs) and mitochondrial COI to uncover the evolutionary history and patterns of temporal and spatial diversification of the fairy wrasses. Our analyses of phylogenetic signal suggest that most gene-tree incongruence is caused by estimation error, leading to poor resolution in a summary-coalescent analysis of the data. In contrast, analyses of concatenated sequences are able to resolve the major relationships of Cirrhilabrus. We determine the placements of species that were previously regarded as incertae sedis and find evidence for the nesting of Conniella, an unusual, monotypic genus, within Cirrhilabrus. Our relaxed-clock dating analysis indicates that the major divergences within the genus occurred around the Miocene-Pliocene boundary, followed by extensive cladogenesis of species complexes in the Pliocene-Pleistocene. Biogeographic reconstruction suggests that the fairy wrasses emerged within the Coral Triangle, with episodic fluctuations of sea levels during glacial cycles coinciding with shallow divergence events but providing few opportunities for more widespread dispersal. Our study demonstrates both the resolving power and limitations of UCEs across shallow timescales where there is substantial estimation error in individual gene trees.

Genomic landscape of geographically structured colour polymorphism in a temperate marine fish

The study of phenotypic variation patterns among populations is fundamental to elucidate the drivers of evolutionary processes. Empirical evidence that supports ongoing genetic divergence associated with phenotypic variation remains very limited for marine species where larval dispersal is a common homogenizing force. We present a genome‐wide analysis of a marine fish, Labrus bergylta, comprising 144 samples distributed from Norway to Spain, a large geographical area that harbours a gradient of phenotypic differentiation. We analysed 39,602 biallelic single nucleotide polymorphisms and found a clear latitudinal gradient of genomic differentiation strongly correlated with the variation in phenotypic morph frequencies observed across the North Atlantic. We also detected a strong association between the latitude and the number of loci that appear to be under divergent selection, which increased with differences in coloration but not with overall genetic differentiation. Our results demonstrate that strong reproductive isolation is occurring between sympatric colour morphs of L. bergylta found at the southern areas and provide important new insights into the genomic changes shaping early stages of differentiation that might precede speciation with gene flow.

Ecology shapes epistasis in a genotype-phenotype-fitness map for stick insect colour

Genetic interactions such as epistasis are widespread in nature and can shape evolutionary dynamics. Epistasis occurs due to nonlinearity in biological systems, which can arise via cellular processes that convert genotype to phenotype and via selective processes that connect phenotype to fitness. Few studies in nature have connected genotype to phenotype to fitness for multiple potentially interacting genetic variants. Thus, the causes of epistasis in the wild remain poorly understood. Here, we show that epistasis for fitness is an emergent and predictable property of nonlinear selective processes. We do so by measuring the genetic basis of cryptic colouration and survival in a field experiment with stick insects. We find that colouration shows a largely additive genetic basis but with some effects of epistasis that enhance differentiation between colour morphs. In terms of fitness, different combinations of loci affecting colouration confer high survival in one host-plant treatment. Specifically, nonlinear correlational selection for specific combinations of colour traits in this treatment drives the emergence of pairwise and higher-order epistasis for fitness at loci underlying colour. In turn, this results in a rugged fitness landscape for genotypes. In contrast, fitness epistasis was dampened in another treatment, where selection was weaker. Patterns of epistasis that are shaped by ecologically based selection could be common and central to understanding fitness landscapes, the dynamics of evolution and potentially other complex systems.

Aggressive mimicry in a coral reef fish: The prey's view

Since all forms of mimicry are based on perceptual deception, the sensory ecology of the intended receiver is of paramount importance to test the necessary precondition for mimicry to occur, that is, model-mimic misidentification, and to gain insight in the origin and evolutionary trajectory of the signals. Here we test the potential for aggressive mimicry by a group of coral reef fishes, the color polymorphic Hypoplectrus hamlets, from the point of view of their most common prey, small epibenthic gobies and mysid shrimp. We build visual models based on the visual pigments and spatial resolution of the prey, the underwater light spectrum and color reflectances of putative models and their hamlet mimics. Our results are consistent with one mimic-model relationship between the butter hamlet H. unicolor and its model the butterflyfish Chaetodon capistratus but do not support a second proposed mimic-model pair between the black hamlet H. nigricans and the dusky damselfish Stegastes adustus. We discuss our results in the context of color morphs divergence in the Hypoplectrus species radiation and suggest that aggressive mimicry in H. unicolor might have originated in the context of protective (Batesian) mimicry by the hamlet from its fish predators rather than aggressive mimicry driven by its prey.

Twelve quick steps for genome assembly and annotation in the classroom

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Visual mate preference evolution during butterfly speciation is linked to neural processing genes

Many animal species remain separate not because their individuals fail to produce viable hybrids but because they "choose" not to mate. However, we still know very little of the genetic mechanisms underlying changes in these mate preference behaviours. Heliconius butterflies display bright warning patterns, which they also use to recognize conspecifics. Here, we couple QTL for divergence in visual preference behaviours with population genomic and gene expression analyses of neural tissue (central brain, optic lobes and ommatidia) across development in two sympatric Heliconius species. Within a region containing 200 genes, we identify five genes that are strongly associated with divergent visual preferences. Three of these have previously been implicated in key components of neural signalling (specifically an ionotropic glutamate receptor and two regucalcins), and overall our candidates suggest shifts in behaviour involve changes in visual integration or processing. This would allow preference evolution without altering perception of the wider environment.

Genetic Evidence for a Mixed Composition of the Genus Myoxocephalus (Cottoidei: Cottidae) Necessitates Generic Realignment

Sculpin fishes belonging to the family Cottidae represent a large and complex group, inhabiting a wide range of freshwater, brackish-water, and marine environments. Numerous studies based on analysis of their morphology and genetic makeup frequently provided controversial results. In the present work, we sequenced complete mitochondrial (mt) genomes and fragments of nuclear ribosomal DNA (rDNA) of the fourhorn sculpin Myoxocephalus quadricornis and some related cottids to increase the power of phylogenetic and taxonomic analyses of this complex fish group. A comparison of the My. quadricornis mt genomes obtained by us with other complete mt genomes available in GenBank has revealed a surprisingly low divergence (3.06 ± 0.12%) with Megalocottus platycephalus and, at the same time, a significantly higher divergence (7.89 ± 0.16%) with the species of the genus Myoxocephalus. Correspondingly, phylogenetic analyses have shown that My. quadricornis is clustered with Me. platycephalus but not with the Myoxocephalus species. Completely consistent patterns of divergence and tree topologies have been obtained based on nuclear rDNA. Thus, the multi-gene data in the present work indicates obvious contradictions in the relationships between the Myoxocephalus and Megalocottus species studied. An extensive phylogenetic analysis has provided evidence for a closer affinity of My. quadricornis with the species of the genus Megalocottus than with the species of the genus Myoxocephalus. A recombination analysis, along with the additional GenBank data, excludes introgression and/or incorrect taxonomic identification as the possible causative factors responsible for the observed closer affinity between the two species from different genera. The above facts necessitate realignment of the genera Myoxocephalus and Megalocottus. The genetic data supports the two recognized genera, Myoxocephalus and Megalocottus, but suggests changing their compositions through transferring My. quadricornis to the genus Megalocottus. The results of the present study resolve the relationships within a complex group of sculpin fishes and show a promising approach to phylogenetic systematics (as a key organizing principle in biodiversity research) for a better understanding of the taxonomy and evolution of fishes and for supplying relevant information to address various fish biodiversity conservation and management issues.

Chemical signals act as the main reproductive barrier between sister and mimetic Heliconius butterflies

Colour pattern is the main trait that drives mate recognition between Heliconius species that are phylogenetically close. However, when this cue is compromised such as in cases of mimetic, sympatric and closely related species, alternative mating signals must evolve to ensure reproductive isolation and species integrity. The closely related species Heliconius melpomene malleti and H. timareta florencia occur in the same geographical region, and despite being co-mimics, they display strong reproductive isolation. In order to test which cues differ between species, and potentially contribute to reproductive isolation, we quantified differences in the wing phenotype and the male chemical profile. As expected, the wing colour pattern was indistinguishable between the two species, while the chemical profile of the androconial and genital males' extracts showed marked differences. We then conducted behavioural experiments to study the importance of these signals in mate recognition by females. In agreement with our previous results, we found that chemical blends and not wing colour pattern drive the preference of females for conspecific males. Also, experiments with hybrid males and females suggested an important genetic component for both chemical production and preference. Altogether, these results suggest that chemicals are the major reproductive barrier opposing gene flow between these two sister and co-mimic species.

Stable species boundaries despite ten million years of hybridization in tropical eels

Genomic evidence is increasingly underpinning that hybridization between taxa is commonplace, challenging our views on the mechanisms that maintain their boundaries. Here, we focus on seven catadromous eel species (genus Anguilla) and use genome-wide sequence data from more than 450 individuals sampled across the tropical Indo-Pacific, morphological information, and three newly assembled draft genomes to compare contemporary patterns of hybridization with signatures of past introgression across a time-calibrated phylogeny. We show that the seven species have remained distinct for up to 10 million years and find that the current frequencies of hybridization across species pairs contrast with genomic signatures of past introgression. Based on near-complete asymmetry in the directionality of hybridization and decreasing frequencies of later-generation hybrids, we suggest cytonuclear incompatibilities, hybrid breakdown, and purifying selection as mechanisms that can support species cohesion even when hybridization has been pervasive throughout the evolutionary history of clades.

Evolution: Genomic Signatures of Mimicry and Mimicry of Genomic Signatures

How new species form in the ocean, and thus what determines the diversity of fish in the sea, is not well understood. A study in Caribbean coral-reef fishes sheds light on the genomic underpinnings of diversification in the marine realm.