Genetics and evidence for balancing selection of a sex-linked colour polymorphism in a songbird

Spread of yellow-bill-color alleles favored by selection in the long-tailed finch hybrid system

Carotenoid pigments produce the yellow and red colors of birds and other vertebrates. Despite their importance in social signaling and sexual selection, our understanding of how carotenoid ornamentation evolves in nature remains limited. Here, we examine the long-tailed finch Poephila acuticauda, an Australian songbird with a yellow-billed western subspecies acuticauda and a red-billed eastern subspecies hecki, which hybridize where their ranges overlap. We found that yellow bills can be explained by the loss of C(4)-oxidation, thus preventing yellow dietary carotenoids from being converted to red. Combining linked-read genomic sequencing and reflectance spectrophotometry measurements of bill color collected from wild-sampled finches and laboratory crosses, we identify four loci that together explain 53% of variance in this trait. The two loci of largest effect contain the genes CYP2J19, an essential enzyme for producing red carotenoids, and TTC39B, an enhancer of carotenoid metabolism. A paucity of protein-coding changes and an enrichment of associated upstream variants suggest that the loss of C(4)-oxidation results from cis-regulatory evolution. Evolutionary genealogy reconstruction indicates that the red-billed phenotype is ancestral and that yellow alleles at CYP2J19 and TTC39B first arose and fixed in acuticauda approximately 100 kya. Yellow alleles subsequently introgressed into hecki less than 5 kya. Across all color loci, acuticauda-derived variants show evidence of selective sweeps, implying that yellow bill coloration has been favored by natural selection. Our study illustrates how evolutionary transitions between yellow and red coloration can be achieved by successive selective events acting on regulatory changes at a few interacting genes.

Common misconceptions of speciation

Speciation is a complex process that can unfold in many different ways. Speciation researchers sometimes simplify core principles in their writing in a way that implies misconceptions about the speciation process. While we think that these misconceptions are usually inadvertently implied (and not actively believed) by the researchers, they nonetheless risk warping how external readers understand speciation. Here we highlight six misconceptions of speciation that are especially widespread. First, species are implied to be clearly and consistently defined entities in nature, whereas in reality species boundaries are often fuzzy and semipermeable. Second, speciation is often implied to be 'good', which is two-fold problematic because it implies both that evolution has a goal and that speciation universally increases the chances of lineage persistence. Third, species-poor clades with species-rich sister clades are considered 'primitive' or 'basal', falsely implying a ladder of progress. Fourth, the evolution of species is assumed to be strictly tree-like, but genomic findings show widespread hybridization more consistent with network-like evolution. Fifth, a lack of association between a trait and elevated speciation rates in macroevolutionary studies is often interpreted as evidence against its relevance in speciation-even if microevolutionary case studies show that it is relevant. Sixth, obvious trait differences between species are sometimes too readily assumed to be (i) barriers to reproduction, (ii) a stepping-stone to inevitable speciation, or (iii) reflective of the species' whole divergence history. In conclusion, we call for caution, particularly when communicating science, because miscommunication of these ideas provides fertile ground for misconceptions to spread.

Supergenes are not necessary to explain the maintenance of complex alternative phenotypes

Evolutionary biology aims to explain the diversity seen in nature. Evolutionary theory provides frameworks to understand how simple polymorphisms or continuous variation are maintained, but phenotypes inherited as discrete suites of quantitative traits are difficult to fit into this framework. Supergenes have been proposed as a solution to this problem-if causal genes are co-located, they can be inherited as if a single gene, thus bridging the gap between simple polymorphisms and continuous traits. We develop models to ask: how are critical supergenes for maintaining phenotypic diversity? In our simplest model, without explicit genetic architectures, three alternative reproductive morphs are maintained in many of the parameter combinations we evaluated. For these same parameter values, models with demographic stochasticity, recombination and mutation (but without explicit genetic architecture) maintained only two of these three morphs, with stochasticity determining which morphs persisted. With explicit genetic architectures, regardless of whether causal loci were co-located in a supergene or distributed randomly, this stochasticity in which morphs are maintained was reduced. Even when phenotypic variation was lost, genetic diversity was maintained. Altogether, categorical traits with polygenic bases exhibited similar evolutionary dynamics to those determined by supergenes. Our work suggests that supergenes are not the only answer to the puzzle of how discrete polygenic phenotypic variation is maintained.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Temporal genomics in Hawaiian crickets reveals compensatory intragenomic coadaptation during adaptive evolution

Theory predicts that compensatory genetic changes reduce negative indirect effects of selected variants during adaptive evolution, but evidence is scarce. Here, we test this in a wild population of Hawaiian crickets using temporal genomics and a high-quality chromosome-level cricket genome. In this population, a mutation, flatwing, silences males and rapidly spread due to an acoustically-orienting parasitoid. Our sampling spanned a social transition during which flatwing fixed and the population went silent. We find long-range linkage disequilibrium around the putative flatwing locus was maintained over time, and hitchhiking genes had functions related to negative flatwing-associated effects. We develop a combinatorial enrichment approach using transcriptome data to test for compensatory, intragenomic coevolution. Temporal changes in genomic selection were distributed genome-wide and functionally associated with the population's transition to silence, particularly behavioural responses to silent environments. Our results demonstrate how 'adaptation begets adaptation'; changes to the sociogenetic environment accompanying rapid trait evolution can generate selection provoking further, compensatory adaptation.

The genomics and evolution of inter-sexual mimicry and female-limited polymorphisms in damselflies

Sex-limited morphs can provide profound insights into the evolution and genomic architecture of complex phenotypes. Inter-sexual mimicry is one particular type of sex-limited polymorphism in which a novel morph resembles the opposite sex. While inter-sexual mimics are known in both sexes and a diverse range of animals, their evolutionary origin is poorly understood. Here, we investigated the genomic basis of female-limited morphs and male mimicry in the common bluetail damselfly. Differential gene expression between morphs has been documented in damselflies, but no causal locus has been previously identified. We found that male mimicry originated in an ancestrally sexually dimorphic lineage in association with multiple structural changes, probably driven by transposable element activity. These changes resulted in ~900 kb of novel genomic content that is partly shared by male mimics in a close relative, indicating that male mimicry is a trans-species polymorphism. More recently, a third morph originated following the translocation of part of the male-mimicry sequence into a genomic position ~3.5 mb apart. We provide evidence of balancing selection maintaining male mimicry, in line with previous field population studies. Our results underscore how structural variants affecting a handful of potentially regulatory genes and morph-specific genes can give rise to novel and complex phenotypic polymorphisms.

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

A complex genomic architecture underlies reproductive isolation in a North American oriole hybrid zone

Natural hybrid zones provide powerful opportunities for identifying the mechanisms that facilitate and inhibit speciation. Documenting the extent of genomic admixture allows us to discern the architecture of reproductive isolation through the identification of isolating barriers. This approach is particularly powerful for characterizing the accumulation of isolating barriers in systems exhibiting varying levels of genomic divergence. Here, we use a hybrid zone between two species-the Baltimore (Icterus galbula) and Bullock's (I. bullockii) orioles-to investigate this architecture of reproductive isolation. We combine whole genome re-sequencing with data from an additional 313 individuals amplityped at ancestry-informative markers to characterize fine-scale patterns of admixture, and to quantify links between genes and the plumage traits. On a genome-wide scale, we document several putative barriers to reproduction, including elevated peaks of divergence above a generally high genomic baseline, a large putative inversion on the Z chromosome, and complex interactions between melanogenesis-pathway candidate genes. Concordant and coincident clines for these different genomic regions further suggest the coupling of pre- and post-mating barriers. Our findings of complex and coupled interactions between pre- and post-mating barriers suggest a relatively rapid accumulation of barriers between these species, and they demonstrate the complexities of the speciation process.

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

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

Genetic divergence among threespine stickleback that differ in nuptial coloration

Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.

A large and diverse autosomal haplotype is associated with sex-linked colour polymorphism in the guppy

Male colour patterns of the Trinidadian guppy (Poecilia reticulata) are typified by extreme variation governed by both natural and sexual selection. Since guppy colour patterns are often inherited faithfully from fathers to sons, it has been hypothesised that many of the colour trait genes must be physically linked to sex determining loci as a ‘supergene’ on the sex chromosome. Here, we phenotype and genotype four guppy ‘Iso-Y lines’, where colour was inherited along the patriline for 40 generations. Using an unbiased phenotyping method, we confirm the breeding design was successful in creating four distinct colour patterns. We find that genetic differentiation among the Iso-Y lines is repeatedly associated with a diverse haplotype on an autosome (LG1), not the sex chromosome (LG12). Moreover, the LG1 haplotype exhibits elevated linkage disequilibrium and evidence of sex-specific diversity in the natural source population. We hypothesise that colour pattern polymorphism is driven by Y-autosome epistasis.

Extreme colour pattern variation in male Trinidadian guppies are influenced by natural selection and sexual selection. Here, the authors phenotype and genotype four guppy lineages finding that colour pattern is associated with a diverse haplotype on an autosome.

Neo-sex chromosome evolution and phenotypic differentiation across an elevational gradient in horned larks (Eremophila Alpestris)

Genetic structure and phenotypic variation among populations is affected by both geographic distance and environmental variation across species' distributions. Understanding the relative contributions of isolation by distance (IBD) and isolation by environment (IBE) is important for elucidating population dynamics across habitats and ecological gradients. In this study, we compared phenotypic and genetic variation among Horned Lark (Eremophila alpestris) populations from 10 sites encompassing an elevational gradient from low-elevation desert scrub in Death Valley (285 a.s.l.) to high-elevation meadows in the White Mountains of the Sierra Nevada of California (greater than 3000 m a.s.l.). Using a ddRAD dataset of 28,474 SNPs aligned to a high-quality reference genome, we compared genetic structure with elevational, environmental, and spatial distance to quantify how different aspects of the landscape drive genomic and phenotypic differentiation in Horned Larks. We found larger-bodied birds were associated with sites that had less seasonality and higher annual precipitation, and longer spurs occurred in soils with more clay and silt content, less sand, and finer fragments. Larks have large neo-sex chromosomes, and we found that associations with elevation and environmental variation were much stronger among neo-sex chromosomes compared to autosomes. Furthermore, we found that putative chromosomal translocations, fusions, and inversions were associated with elevation and may underlie local adaptation across an elevational gradient in Horned Larks. Our results suggest that genetic variation in Horned Larks is affected more by IBD than IBE, but specific phenotypes and genomic regions-particually on neo-sex chromosomes-bear stronger associations with the environment.

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

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

Comparative Genomics and Evolution of Avian Specialized Traits

Genomic data are important for understanding the origin and evolution of traits. Under the context of rapidly developing of sequencing technologies and more widely available genome sequences, researchers are able to study evolutionary mechanisms of traits via comparative genomic methods. Compared with other vertebrates, bird genomes are relatively small and exhibit conserved synteny with few repetitive elements, which makes them suitable for evolutionary studies. Increasing genomic progress has been reported on the evolution of powered flight, body size variation, beak morphology, plumage colouration, high-elevation colonization, migration, and vocalization. By summarizing previous studies, we demonstrate the genetic bases of trait evolution, highlighting the roles of small-scale sequence variation, genomic structural variation, and changes in gene interaction networks. We suggest that future studies should focus on improving the quality of reference genomes, exploring the evolution of regulatory elements and networks, and combining genomic data with morphological, ecological, behavioural, and developmental biology data.

Individual Variation in Dietary Wariness Is Predicted by Head Color in a Specialist Feeder, the Gouldian Finch

Shifts in resource availability due to environmental change are increasingly confronting animals with unfamiliar food types. Species that can rapidly accept new food types may be better adapted to ecological change. Intuitively, dietary generalists are expected to accept new food types when resources change, while dietary specialists would be more averse to adopting novel food. However, most studies investigating changes in dietary breadth focus on generalist species and do not delve into potential individual predictors of dietary wariness and the social factors modulating these responses. We investigated dietary wariness in the Gouldian finch, a dietary specialist, that is expected to avoid novel food. This species occurs in two main head colors (red, black), which signal personality in other contexts. We measured their initial neophobic responses (approach attempts before first feed and latency to first feed) and willingness to incorporate novel food into their diet (frequency of feeding on novel food after first feed). Birds were tested in same-sex pairs in same and different head color pairings balanced across experiments 1 and 2. Familiar and novel food (familiar food dyed) were presented simultaneously across 5 days for 3 h, each. Gouldian finches fed on the familiar food first demonstrating food neophobia, and these latencies were repeatable. Birds made more approach attempts before feeding on novel than familiar food, particularly red-headed birds in experiment 1 and when partnered with a black-headed bird. Individuals consistently differed in their rate of incorporation of novel food, with clear differences between head colors; red-headed birds increased their feeding visits to novel food across experimentation equaling their familiar food intake by day five, while black-headed birds continually favored familiar food. Results suggest consistent among individual differences in response to novel food with red-headed birds being adventurous consumers and black-headed birds dietary conservatives. The differences in food acceptance aligned with responses to novel environments on the individual level (found in an earlier study) providing individuals with an adaptive combination of novelty responses across contexts in line with potential differences in movement patterns. Taken together, these novelty responses could aid in population persistence when faced with environmental changes.

A multispecies BCO2 beak color polymorphism in the Darwin's finch radiation

Carotenoid-based polymorphisms are widespread in populations of birds, fish, and reptiles,¹ but generally little is known about the factors affecting their maintenance in populations.² We report a combined field and molecular-genetic investigation of a nestling beak color polymorphism in Darwin's finches. Beaks are pink or yellow, and yellow is recessive.³ Here we show that the polymorphism arose in the Galápagos half a million years ago through a mutation associated with regulatory change in the BCO2 gene and is shared by 14 descendant species. The polymorphism is probably a balanced polymorphism, maintained by ecological selection associated with survival and diet. In cactus finches, the frequency of the yellow genotype is correlated with cactus fruit abundance and greater hatching success and may be altered by introgressive hybridization. Polymorphisms that are hidden as adults, as here, may be far more common than is currently recognized, and contribute to diversification in ways that are yet to be discovered.

Sex and morph differences in age-dependent trait changes in a polymorphic songbird

There is growing evidence of important variation in how animals age, in particular in how the expression of traits changes with age among different species and populations. However, less is known about variation within populations, which may include variation in ageing patterns between different types of individuals (e.g. sexes or distinct polymorphisms) and between different types of traits (e.g. general traits versus those used in social signalling contexts). We used 6 years of longitudinal data to examine age-related changes in trait expression in a captive population of Gouldian finches (Erythrura gouldiae), a socially monogamous songbird with genetically determined colour morphs that differ in behaviour and physiology. We contrasted ageing patterns of different types of traits (social signalling vs. size-related) in both sexes and in two colour morphs, using a mixed model approach to account for both within- and between-individual effects. We found pronounced sex differences in how social signalling traits change with age, showing a quadratic pattern in males, but not changing with age in females. In contrast, we observed no sex-specific ageing patterns in size traits. We also found subtle morph differences in how size-related traits changed with age, with black morphs stable or increasing with age while red morphs showing a decline with age. Finally, we found an interesting sex by morph interaction in one important social signal (headband width). These results highlight the importance of using within-individual approaches to understand ageing patterns across types of individuals (sex, morph, etc.) and the need for further research on the ageing patterns of traits that may experience different selective pressures.

Rare morph Lake Malawi mbuna cichlids benefit from reduced aggression from con- and hetero-specifics

Balancing selection is important for the maintenance of polymorphism as it can prevent either fixation of one morph through directional selection or genetic drift, or speciation by disruptive selection. Polychromatism, the presence of multiple genetically determined colour phenotypes, can be maintained if the fitness of alternative morphs depends on the relative frequency in a population. In aggressive species, negative frequency-dependent antagonism can prevent an increase in the frequency of rare morphs as they would only benefit from increased fitness while they are rare. Heterospecific aggression is common in nature and has the potential to contribute to rare morph advantage. Here we carry out field observations and laboratory aggression experiments with mbuna cichlids from Lake Malawi, to investigate the role of con- and heterospecific aggression in the maintenance of polychromatism and identify benefits to rare morphs which are likely to result from reduced aggression. We hypothesize that rare morph individuals receive less aggression than common morph individuals and therefore have an ecological advantage. Within species we found that males and females bias aggression towards their own morph, adding to the evidence that inherent own-morph aggression biases can contribute to balancing selection. Over-representation of rare morph territory owners may be influenced by two factors; higher tolerance of different morph individuals as neighbours, and the ability of rare morphs to spend more time feeding. Reduced aggression to rare morph individuals by heterospecifics may also contribute to rare morph advantage.

Avian Coloration Genetics: Recent Advances and Emerging Questions

The colorful phenotypes of birds have long provided rich source material for evolutionary biologists. Avian plumage, beaks, skin, and eggs-which exhibit a stunning range of cryptic and conspicuous forms-inspired early work on adaptive coloration. More recently, avian color has fueled discoveries on the physiological, developmental, and-increasingly-genetic mechanisms responsible for phenotypic variation. The relative ease with which avian color traits can be quantified has made birds an attractive system for uncovering links between phenotype and genotype. Accordingly, the field of avian coloration genetics is burgeoning. In this review, we highlight recent advances and emerging questions associated with the genetic underpinnings of bird color. We start by describing breakthroughs related to 2 pigment classes: carotenoids that produce red, yellow, and orange in most birds and psittacofulvins that produce similar colors in parrots. We then discuss structural colors, which are produced by the interaction of light with nanoscale materials and greatly extend the plumage palette. Structural color genetics remain understudied-but this paradigm is changing. We next explore how colors that arise from interactions among pigmentary and structural mechanisms may be controlled by genes that are co-expressed or co-regulated. We also identify opportunities to investigate genes mediating within-feather micropatterning and the coloration of bare parts and eggs. We conclude by spotlighting 2 research areas-mechanistic links between color vision and color production, and speciation-that have been invigorated by genetic insights, a trend likely to continue as new genomic approaches are applied to non-model species.

Metabolic rates of aggressive and submissive phenotypes are colour blind in the polymorphic Gouldian finch

Evidence from a number of species suggests that behaviours associated with social rank are positively correlated with metabolic rate. These studies, however, are based on metabolic measurements of isolated individuals, thereby ignoring potential effects of social interactions on metabolic rates. Here, we characterised three pertinent metabolic indices in the two predominant genetic colour morphs of the Gouldian finch (Erythrura gouldiae): diurnal resting metabolic rate (RMR), nocturnal basal metabolic rate (BMR) and exercise-induced maximal metabolic rate (MMR). Research reveals that red-headed morphs consistently dominate the less aggressive black-headed morphs and that the two morphs differ in other behavioural and physiological traits. We measured daytime RMR of intermorph naïve birds (first-year virgin males maintained in total isolation from opposite colour morphs) and their metabolic responses to viewing a socially unfamiliar bird of each colour. Subsequently, each bird was placed in a home cage with an opposite colour morph (intermorph exposed) and the series of measurements was repeated. Daytime RMR was indistinguishable between the two morphs, regardless of whether they were intermorph naïve or intermorph exposed. However, both red- and black-headed birds showed a greater short-term increase in metabolic rate when viewing an unfamiliar red-headed bird than when seeing a black-headed bird, but only when intermorph naïve. Measurements of BMR and exercise-induced MMR did not differ between the two morphs, and consequently, aerobic scope was indistinguishable between them. We propose that the behavioural differences between these two sympatric morphs are functionally complementary and represent evolutionary stable strategies permitting establishment of dominance status in the absence of metabolic costs.

Genomic investigation of colour polymorphism and phylogeographic variation among populations of black-headed bulbul (Brachypodius atriceps) in insular southeast Asia

Intraspecific polymorphism in birds, especially plumage colour polymorphism, and the mechanisms that control it are an area of active research in evolutionary biology. The black-headed bulbul (Brachypodius atriceps) is a polymorphic species with two distinct morphs, yellow and grey. This species inhabits the mainland and virtually all continental islands of Southeast Asia where yellow morphs predominate, but on two islands in the Sunda region, Bawean and Maratua, grey morphs are common or exclusive. Here, we generated a high-quality reference genome of a yellow individual and resequenced genomes of multiple individuals of both yellow and grey morphs to study the genetic basis of coloration and population history of the species. Using PCA and STRUCTURE analysis, we found the Maratua Island population (which is exclusively grey) to be distinct from all other B. atriceps populations, having been isolated c. 1.9 million years ago (Ma). In contrast, Bawean grey individuals (a subset of yellow and grey individuals on that island) are embedded within an almost panmictic Sundaic clade of yellow birds. Using F_ST and d_xy to compare variable genomic segments between Maratua and yellow individuals, we located peaks of divergence and identified candidate loci involved in the colour polymorphism. Tests of selection among coding-proteins in high F_ST regions, however, did not indicate selection on the candidate genes. Overall, we report on some loci that are potentially responsible for the grey/yellow polymorphism in a species that otherwise shows little genetic diversification across most of its range.

Transcriptomic Analysis of Skin Color in Anole Lizards

Color and color pattern are critical for animal camouflage, reproduction, and defense. Few studies, however, have attempted to identify candidate genes for color and color pattern in squamate reptiles, a colorful group with over 10,000 species. We used comparative transcriptomic analyses between white, orange, and yellow skin in a color-polymorphic species of anole lizard to 1) identify candidate color and color-pattern genes in squamates and 2) assess if squamates share an underlying genetic basis for color and color pattern variation with other vertebrates. Squamates have three types of chromatophores that determine color pattern: guanine-filled iridophores, carotenoid- or pteridine-filled xanthophores/erythrophores, and melanin-filled melanophores. We identified 13 best candidate squamate color and color-pattern genes shared with other vertebrates: six genes linked to pigment synthesis pathways, and seven genes linked to chromatophore development and maintenance. In comparisons of expression profiles between pigment-rich and white skin, pigment-rich skin upregulated the pteridine pathway as well as xanthophore/erythrophore development and maintenance genes; in comparisons between orange and yellow skin, orange skin upregulated the pteridine and carotenoid pathways as well as melanophore maintenance genes. Our results corroborate the predictions that squamates can produce similar colors using distinct color-reflecting molecules, and that both color and color-pattern genes are likely conserved across vertebrates. Furthermore, this study provides a concise list of candidate genes for future functional verification, representing a first step in determining the genetic basis of color and color pattern in anoles.

Birds of different morphs use slightly different strategies to achieve similar reproductive performance following heatwave exposure

Responses to extreme climatic events may differ between individuals of distinct morphs which differ in life-history strategies, resulting in climate change 'winners' and 'losers' within species. We examined the reproductive performance and carry-over effects on offspring of black- and red-headed Gouldian finches Erythrura gouldiae after exposure to simulated heatwaves of moderate or severe intensity. We expected black-headed pairs' reproductive performance to decline after the severe heatwave because only the condition of black-headed females deteriorates during such a heatwave. Supporting the fact that Gouldian finches of different morphs use alternative reproductive strategies, we found that black-headed females initiated egg-laying a month earlier than red-headed females after experiencing a severe heatwave. We also found that this severe heatwave resulted in shorter spermatozoa in males irrespective of their morph. Despite these effects associated with heatwave intensity, the overall reproductive performance of both morphs was not affected by this factor, which was possibly due to an increased nestling provisioning rate by parents after exposure to the severe heatwave. However, offspring still bore the cost of parental exposure to the severe heatwave, as they showed a reduced condition (lower plasma antioxidant capacity and transient lower breathing rate) and higher oxidative damage (at least in fledglings with black-headed parents). These results suggest that inter-morph phenotypic variability in the Gouldian finch does not result in clear differences in reproductive performance following heatwave exposure, despite basal phenotypic differences between morphs. Whether animals using alternative reproductive strategies are, in the end, differently affected by climate changes will likely depend on the capacity of their offspring to recover from altered developmental conditions.

Experimental tests of selection against heterospecific aggression as a driver of avian colour pattern divergence

Signal divergence is thought to reduce the costs of co-occurrence for closely related species and may thereby be important in the generation and maintenance of new biodiversity. In birds, closely related, sympatric species are more divergent in their colour patterns than those that live apart, but the selective pressures driving sympatric divergence in colour pattern are not well-understood. Here, we conducted field experiments on naïve birds using spectrometer-matched, painted, 3D-printed models to test whether selection against heterospecific aggression might drive colour pattern divergence in the genus Poecile. We found that territorial male black-capped chickadees (P. atricapillus) are equally likely to attack sympatric and allopatric congeners, and wintering flocks are equally likely to visit feeders occupied by sympatric and allopatric congeners, despite sympatric congeners being more divergent in colour pattern. These results suggest that either the concerted evolution of additional traits (e.g. discrimination), or interactions in sympatry that promote learning, is required if colour pattern divergence among sympatric species is to reduce heterospecific aggression. Alternatively, colour pattern divergence among sympatric species may be caused by other selective pressures, such as selection against hybridization or habitat partitioning and secondary signal adaptation.

Extensive historical and contemporary hybridization suggests premating isolation in Vermivora warblers is not strong: A reply to Confer et al

We present comments on an article published by Confer et al. (Ecology and Evolution, 10, 2020). Confer et al. (2020) aggregate data from multiple studies of social pairing between Vermivora chrysoptera and V. cyanoptera, two wood warblers in the family Parulidae that hybridize extensively where they co‐occur. From analysis of these data, they conclude there is near‐complete reproductive isolation between these two species. In our reply, we show that this finding is not supported by other lines of evidence, and significant drawbacks of their study design preclude such strong conclusions. In our critique, we show that (a) coarse‐scale plumage classifications cannot be used to accurately estimate hybrid ancestry in Vermivora; (b) extra‐pair paternity is very high in Vermivora and is likely facilitating hybridization, yet was not considered by Confer et al. (2020), and we suggest this will have a substantial influence on the interpretation of reproductive isolation in the system; and (c) the central finding of strong total reproductive isolation is not compatible with the results of other long‐term studies, which demonstrate low isolation and high gene flow. We conclude with a more comprehensive interpretation of hybridization and reproductive isolation in Vermivora warblers.

Studying models of balancing selection using phase-type theory

Balancing selection (BLS) is the evolutionary force that maintains high levels of genetic variability in many important genes. To further our understanding of its evolutionary significance, we analyze models with BLS acting on a biallelic locus: an equilibrium model with long-term BLS, a model with long-term BLS and recent changes in population size, and a model of recent BLS. Using phase-type theory, a mathematical tool for analyzing continuous time Markov chains with an absorbing state, we examine how BLS affects polymorphism patterns in linked neutral regions, as summarized by nucleotide diversity, the expected number of segregating sites, the site frequency spectrum, and the level of linkage disequilibrium (LD). Long-term BLS affects polymorphism patterns in a relatively small genomic neighborhood, and such selection targets are easier to detect when the equilibrium frequencies of the selected variants are close to 50%, or when there has been a population size reduction. For a new mutation subject to BLS, its initial increase in frequency in the population causes linked neutral regions to have reduced diversity, an excess of both high and low frequency derived variants, and elevated LD with the selected locus. These patterns are similar to those produced by selective sweeps, but the effects of recent BLS are weaker. Nonetheless, compared to selective sweeps, nonequilibrium polymorphism and LD patterns persist for a much longer period under recent BLS, which may increase the chance of detecting such selection targets. An R package for analyzing these models, among others (e.g., isolation with migration), is available.

The Genetic Architecture of Variation in the Sexually Selected Sword Ornament and Its Evolution in Hybrid Populations

Biologists since Darwin have been fascinated by the evolution of sexually selected ornaments, particularly those that reduce viability. Uncovering the genetic architecture of these traits is key to understanding how they evolve and are maintained. Here, we investigate the genetic architecture and evolutionary loss of a sexually selected ornament, the "sword" fin extension that characterizes many species of swordtail fish (Xiphophorus). Using sworded and swordless sister species of Xiphophorus, we generated a mapping population and show that the sword ornament is polygenic-with ancestry across the genome explaining substantial variation in the trait. After accounting for the impacts of genome-wide ancestry, we identify one major-effect quantitative trait locus (QTL) that explains ~5% of the overall variation in the trait. Using a series of approaches, we narrow this large QTL interval to several likely candidate genes, including genes involved in fin regeneration and growth. Furthermore, we find evidence of selection on ancestry at one of these candidates in four natural hybrid populations, consistent with selection against the sword in these populations.

Extensive hybridization reveals multiple coloration genes underlying a complex plumage phenotype

Coloration is an important target of both natural and sexual selection. Discovering the genetic basis of colour differences can help us to understand how this visually striking phenotype evolves. Hybridizing taxa with both clear colour differences and shallow genomic divergences are unusually tractable for associating coloration phenotypes with their causal genotypes. Here, we leverage the extensive admixture between two common North American woodpeckers-yellow-shafted and red-shafted flickers-to identify the genomic bases of six distinct plumage patches involving both melanin and carotenoid pigments. Comparisons between flickers across approximately 7.25 million genome-wide SNPs show that these two forms differ at only a small proportion of the genome (mean F_ST = 0.008). Within the few highly differentiated genomic regions, we identify 368 SNPs significantly associated with four of the six plumage patches. These SNPs are linked to multiple genes known to be involved in melanin and carotenoid pigmentation. For example, a gene (CYP2J19) known to cause yellow to red colour transitions in other birds is strongly associated with the yellow versus red differences in the wing and tail feathers of these flickers. Additionally, our analyses suggest novel links between known melanin genes and carotenoid coloration. Our finding of patch-specific control of plumage coloration adds to the growing body of literature suggesting colour diversity in animals could be created through selection acting on novel combinations of coloration genes.

Rapid parallel adaptation despite gene flow in silent crickets

Gene flow is predicted to impede parallel adaptation via de novo mutation, because it can introduce pre-existing adaptive alleles from population to population. We test this using Hawaiian crickets (Teleogryllus oceanicus) in which 'flatwing' males that lack sound-producing wing structures recently arose and spread under selection from an acoustically-orienting parasitoid. Morphometric and genetic comparisons identify distinct flatwing phenotypes in populations on three islands, localized to different loci. Nevertheless, we detect strong, recent and ongoing gene flow among the populations. Using genome scans and gene expression analysis we find that parallel evolution of flatwing on different islands is associated with shared genomic hotspots of adaptation that contain the gene doublesex, but the form of selection differs among islands and corresponds to known flatwing demographics in the wild. We thus show how parallel adaptation can occur on contemporary timescales despite gene flow, indicating that it could be less constrained than previously appreciated.

Dog colour patterns explained by modular promoters of ancient canid origin

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Trait differences among discrete morphs of a color polymorphic lizard, Podarcis erhardii

Color polymorphism defies evolutionary expectations as striking phenotypic variation is maintained within a single species. Color and other traits mediate social interactions, and stable polymorphism within a population is hypothesized to be related to correlational selection of other phenotypic traits among color morphs. Here, we report on a previously unknown throat color polymorphism in the Aegean Wall Lizard (Podarcis erhardii) and examine morph-correlated differences in traits important to social behavior and communication: maximum bite force capacity and chemical signal profile. We find that both sexes of P. erhardii have three color morphs: orange, yellow, and white. Moreover, orange males are significantly larger and tend to bite harder than yellow and white males. Although the established color polymorphism only partially matches the observed intraspecific variation in chemical signal signatures, the chemical profile of the secretions of orange males is significantly divergent from that of white males. Our findings suggest that morph colors are related to differences in traits that are crucial for social interactions and competitive ability, illustrating the need to look beyond color when studying polymorphism evolution.

CYP2J19 mediates carotenoid colour introgression across a natural avian hybrid zone

It has long been of interest to identify the phenotypic traits that mediate reproductive isolation between related species, and more recently, the genes that underpin them. Much work has focused on identifying genes associated with animal colour, with the candidate gene CYP2J19 identified in laboratory studies as the ketolase converting yellow dietary carotenoids to red ketocarotenoids in birds with red pigments. However, evidence that CYP2J19 explains variation between red and yellow feather coloration in wild populations of birds is lacking. Hybrid zones provide the opportunity to identify genes associated with specific traits. Here we investigate genomic regions associated with colour in red-fronted and yellow-fronted tinkerbirds across a hybrid zone in southern Africa. We sampled 85 individuals, measuring spectral reflectance of forecrown feathers and scoring colours from photographs, while testing for carotenoid presence with Raman spectroscopy. We performed a genome-wide association study to identify associations with carotenoid-based coloration, using double-digest RAD sequencing aligned to a short-read whole genome of a Pogoniulus tinkerbird. Admixture mapping using 104,933 single nucleotide polymorphisms (SNPs) identified a region of chromosome 8 that includes CYP2J19 as the only locus with more than two SNPs significantly associated with both crown hue and crown score, while Raman spectra provided evidence of ketocarotenoids in red feathers. Asymmetric backcrossing in the hybrid zone suggests that yellow-fronted females mate more often with red-fronted males than vice versa. Female red-fronted tinkerbirds mating assortatively with red-crowned males is consistent with the hypothesis that converted carotenoids are an honest signal of quality.

LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60 000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

De Novo Assembly of the Northern Cardinal (Cardinalis cardinalis) Genome Reveals Candidate Regulatory Regions for Sexually Dichromatic Red Plumage Coloration

Northern cardinals (Cardinalis cardinalis) are common, mid-sized passerines widely distributed in North America. As an iconic species with strong sexual dichromatism, it has been the focus of extensive ecological and evolutionary research, yet genomic studies investigating the evolution of genotype-phenotype association of plumage coloration and dichromatism are lacking. Here we present a new, highly-contiguous assembly for C. cardinalis We generated a 1.1 Gb assembly comprised of 4,762 scaffolds, with a scaffold N50 of 3.6 Mb, a contig N50 of 114.4 kb and a longest scaffold of 19.7 Mb. We identified 93.5% complete and single-copy orthologs from an Aves dataset using BUSCO, demonstrating high completeness of the genome assembly. We annotated the genomic region comprising the CYP2J19 gene, which plays a pivotal role in the red coloration in birds. Comparative analyses demonstrated non-exonic regions unique to the CYP2J19 gene in passerines and a long insertion upstream of the gene in C. cardinalis Transcription factor binding motifs discovered in the unique insertion region in C. cardinalis suggest potential androgen-regulated mechanisms underlying sexual dichromatism. Pairwise Sequential Markovian Coalescent (PSMC) analysis of the genome reveals fluctuations in historic effective population size between 100,000-250,000 in the last 2 millions years, with declines concordant with the beginning of the Pleistocene epoch and Last Glacial Period. This draft genome of C. cardinalis provides an important resource for future studies of ecological, evolutionary, and functional genomics in cardinals and other birds.

The genomics of coloration provides insights into adaptive evolution

Coloration is an easily quantifiable visual trait that has proven to be a highly tractable system for genetic analysis and for studying adaptive evolution. The application of genomic approaches to evolutionary studies of coloration is providing new insight into the genetic architectures underlying colour traits, including the importance of large-effect mutations and supergenes, the role of development in shaping genetic variation and the origins of adaptive variation, which often involves adaptive introgression. Improved knowledge of the genetic basis of traits can facilitate field studies of natural selection and sexual selection, making it possible for strong selection and its influence on the genome to be demonstrated in wild populations.

More than the eye can see: Genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean

The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon - the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters.