Admixture mapping in a hybrid zone reveals loci associated with avian feather coloration

Museum genomics provide insight into the extinction of a specialist North American warbler species

Museum genomics provide an opportunity to investigate population demographics of extinct species, especially valuable when research prior to extinction was minimal. The Bachman's warbler (Vermivora bachmanii) is hypothesized to have gone extinct due to loss of its specialized habitat. However, little is known about other potential contributing factors such as natural rarity or changes to connectivity following habitat fragmentation. We examined mitochondrial DNA (mtDNA) and genome-wide SNPs using specimens collected from breeding and migration sites across the range of the Bachman's warbler. We found no signals of strong population structuring across the breeding range of Bachman's warblers in both mtDNA and genome-wide SNPs. Thus, long-term population isolation did not appear to be a significant contributor to the extinction of the Bachman's warbler. Instead, our findings support the theory that Bachman's warblers underwent a rapid decline likely driven by habitat destruction, which may have been exacerbated by the natural rarity, habitat specificity and low genetic diversity of the species.

Waves of Colonization and Gene Flow in a Great Speciator

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba)

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

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.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.

Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas

In the Himalayas, a number of secondary contact zones have been described for vicariant vertebrate taxa. However, analyses of genetic divergence and admixture are missing for most of these examples. In this study, we provide a population genetic analysis for the coal tit (Periparus ater) hybrid zone in Nepal. Intermediate phenotypes between the distinctive western "spot-winged tit" (P. a. melanolophus) and Eastern Himalayan coal tits (P. a. aemodius) occur across a narrow range of <100 km in western Nepal. As a peculiarity, another distinctive cinnamon-bellied form is known from a single population so far. Genetic admixture of western and eastern mitochondrial lineages was restricted to the narrow zone of phenotypically intermediate populations. The cline width was estimated 46 km only with a center close to the population of the cinnamon-bellied phenotype. In contrast, allelic introgression of microsatellite loci was asymmetrical from eastern P. a. aemodius into far western populations of phenotypic P. a. melanolophus but not vice versa. Accordingly, the microsatellite cline was about 3.7 times wider than the mitochondrial one.

A supergene underlies linked variation in color and morphology in a Holarctic songbird

The genetic architecture of a phenotype can have considerable effects on the evolution of a trait or species. Characterizing genetic architecture provides insight into the complexity of a given phenotype and, potentially, the role of the phenotype in evolutionary processes like speciation. We use genome sequences to investigate the genetic basis of phenotypic variation in redpoll finches (Acanthis spp.). We demonstrate that variation in redpoll phenotype is broadly controlled by a ~55-Mb chromosomal inversion. Within this inversion, we find multiple candidate genes related to melanogenesis, carotenoid coloration, and bill shape, suggesting the inversion acts as a supergene controlling multiple linked traits. A latitudinal gradient in ecotype distribution suggests supergene driven variation in color and bill morphology are likely under environmental selection, maintaining supergene haplotypes as a balanced polymorphism. Our results provide a mechanism for the maintenance of ecotype variation in redpolls despite a genome largely homogenized by gene flow.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Admixture mapping reveals loci for carcass mass in red deer x sika hybrids in Kintyre, Scotland

We deployed admixture mapping on a sample of 386 deer from a hybrid swarm between native red deer (Cervus elaphus) and introduced Japanese sika (Cervus nippon) sampled in Kintyre, Scotland to search for quantitative trait loci (QTLs) underpinning phenotypic differences between the species. These two species are highly diverged genetically [Fst between pure species, based on 50K single nucleotide polymorphism (SNPs) = 0.532] and phenotypically: pure red have on average twice the carcass mass of pure sika in our sample (38.7 kg vs 19.1 kg). After controlling for sex, age, and population genetic structure, we found 10 autosomal genomic locations with QTL for carcass mass. Effect sizes ranged from 0.191 to 1.839 kg and as expected, in all cases the allele derived from sika conferred lower carcass mass. The sika population was fixed for all small carcass mass alleles, whereas the red deer population was typically polymorphic. GO term analysis of genes lying in the QTL regions are associated with oxygen transport. Although body mass is a likely target of selection, none of the SNPs marking QTL are introgressing faster or slower than expected in either direction.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F_ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F_ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F_ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species.

Sexual selection and introgression in avian hybrid zones: Spotlight on Manacus

Hybrid zones offer a window into the processes and outcomes of evolution, from species formation or fusion to genomic underpinnings of specific traits and isolating mechanisms. Sexual selection is believed to be an important factor in speciation processes, and hybrid zones present special opportunities to probe its impact. The manakins (Aves, Pipridae) are a promising group in which to study the interplay of sexual selection and natural hybridization: they show substantial variation across the family in the strength of sexual selection they experience, they readily hybridize within and between genera, and they appear to have formed hybrid species, a rare event in birds. A hybrid zone between two manakins in the genus Manacus is unusual in that plumage and behavioral traits of one species have introgressed asymmetrically into populations of the second species through positive sexual selection, then apparently stalled at a river barrier. This is one of a handful of documented examples of asymmetric sexual trait introgression with a known selective mechanism. It offers opportunities to examine reproductive isolation, introgression, plumage color evolution, and natural factors enhancing or constraining the effects of sexual selection in real time. Here, we review previous work in this system, propose new hypotheses for observed patterns, and recommend approaches to test them.

Processes underlying complex patterns of song trait evolution in a Setophaga hybrid zone

During secondary contact between two species when hybrids are less fit than parents, mating signals are expected to diverge, while aggressive signals are expected to converge. If a single signal trait is used in both mating and aggression, then the dynamics between these two forces could influence the evolutionary trajectory of that trait. We studied such a situation in an avian hybrid zone between two Setophaga species, where birdsong is used in both mate attraction and territory defense. We hypothesized that song modules of the two species will show separate and distinct geographic patterns due to the influence of selective pressures for effective territorial aggression and for effective mate attraction. We conducted geographic cline analyses and playback experiments across this hybrid zone. We found an unexpected geographic pattern of asymmetric introgression of song rhythm, which may be explained by results of the playback experiments that suggest that differences in song rhythm serve a greater role in mate attraction than in territory defense. In contrast, differences in syllable morphology show little evidence of importance in mate attraction or territorial defense. Song features converge in the hybrid zone, yet patterns of trait change suggest that the song production modules may vary in their modes of development and inheritance. Syringeal motor gesturing, which gives rise to syllable morphology, shows a nonclinal mosaic pattern, suggesting that this trait may be predominantly learned. In contrast, respiratory patterning, which forms song rhythm, shows a clinal geographic transition, suggesting that this trait could be more innate. The results indicate that opposing forces act independently on song via distinct modules of the song production mechanism, driving complex patterns of song trait evolution.

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.

Selection on a small genomic region underpins differentiation in multiple color traits between two warbler species

Speciation is one of the most important processes in biology, yet the study of the genomic changes underlying this process is in its infancy. North American warbler species Setophaga townsendi and Setophaga occidentalis hybridize in a stable hybrid zone, following a period of geographic separation. Genomic differentiation accumulated during geographic isolation can be homogenized by introgression at secondary contact, whereas genetic regions that cause low hybrid fitness can be shielded from such introgression. Here, we examined the genomic underpinning of speciation by investigating (1) the genetic basis of divergent pigmentation traits between species, (2) variation in differentiation across the genome, and (3) the evidence for selection maintaining differentiation in the pigmentation genes. Using tens of thousands of single nucleotide polymorphisms (SNPs) genotyped in hundreds of individuals within and near the hybrid zone, genome-wide association mapping revealed a single SNP associated with cheek, crown, breast coloration, and flank streaking, reflecting pleiotropy (one gene affecting multiple traits) or close physical linkage of different genes affecting different traits. This SNP is within an intron of the RALY gene, hence we refer to it as the RALY SNP. We then examined between-species genomic differentiation, using both genotyping-by-sequencing and whole genome sequencing. We found that the RALY SNP is within one of the highest peaks of differentiation, which contains three genes known to influence pigmentation: ASIP, EIF2S2, and RALY (the ASIP-RALY gene block). Heterozygotes at this gene block are likely of reduced fitness, as the geographic cline of the RALY SNP has been narrow over two decades. Together, these results reflect at least one barrier to gene flow within this narrow (∼200 kb) genomic region that modulates plumage difference between species. Despite extensive gene flow between species across the genome, this study provides evidence that selection on a phenotype-associated genomic region maintains a stable species boundary.

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.

Testosterone regulates CYP2J19-linked carotenoid signal expression in male red-backed fairywrens (Malurus melanocephalus)

Carotenoid pigments produce most red, orange and yellow colours in vertebrates. This coloration can serve as an honest signal of quality that mediates social and mating interactions, but our understanding of the underlying mechanisms that control carotenoid signal production, including how different physiological pathways interact to shape and maintain these signals, remains incomplete. We investigated the role of testosterone in mediating gene expression associated with a red plumage sexual signal in red-backed fairywrens (Malurus melanocephalus). In this species, males within a single population can flexibly produce either red/black nuptial plumage or female-like brown plumage. Combining correlational analyses with a field-based testosterone implant experiment and quantitative polymerase chain reaction, we show that testosterone mediates expression of carotenoid-based plumage in part by regulating expression of CYP2J19, a ketolase gene associated with ketocarotenoid metabolism and pigmentation in birds. This is, to our knowledge, the first time that hormonal regulation of a specific genetic locus has been linked to carotenoid production in a natural context, revealing how endocrine mechanisms produce sexual signals that shape reproductive success.

Genomic identification of intergeneric hybrids in New World wood-warblers (Aves: Parulidae)

The documentation of hybrids between distantly related taxa can illustrate an initial step to explain how genes might move between species that do not exhibit complete reproductive isolation. In birds, some of the most phylogenetically distant hybrid combinations occur between genera. Traditionally, morphological and plumage characters have been used to assign the identity of the parental species of a putative hybrid, although recently, nuclear introns also have been used. Here, we demonstrate how high-throughput short-read DNA sequence data can be used to identify the parentage of a putative intergeneric hybrid, in this case between a blue-winged warbler (Vermivora cyanoptera) and a cerulean warbler (Setophaga cerulea). This hybrid had mitochondrial DNA of a cerulean warbler, indicating the maternal parent. For hundreds of single nucleotide polymorphisms within six regions of the nuclear genome that differentiate blue-winged warblers and golden-winged warblers (Vermivora chrysoptera), the hybrid had roughly equal ancestry assignment to blue-winged and cerulean warblers, suggesting a blue-winged warbler as the paternal parent species and demonstrating that this was a first generation (F1) hybrid between these species. Unlike other recently characterized intergeneric warbler hybrids, this individual hybrid learned to song match its maternal parent species, suggesting that it might have been the result of an extra-pair mating and raised in a cerulean warbler nest.

Evolution of strong reproductive isolation in plants: broad-scale patterns and lessons from a perennial model group

Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the 'escape' of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.

Multiple Loci Control Eyespot Number Variation on the Hindwings of Bicyclus anynana Butterflies

The underlying genetic changes that regulate the appearance and disappearance of repeated traits, or serial homologs, remain poorly understood. One hypothesis is that variation in genomic regions flanking master regulatory genes, also known as input-output genes, controls variation in trait number, making the locus of evolution almost predictable. Another hypothesis implicates genetic variation in up- or downstream loci of master control genes. Here, we use the butterfly Bicyclus anynana, a species that exhibits natural variation in eyespot number on the dorsal hindwing, to test these two hypotheses. We first estimated the heritability of dorsal hindwing eyespot number by breeding multiple butterfly families differing in eyespot number and regressing eyespot numbers of offspring on midparent values. We then estimated the number and identity of independent genetic loci contributing to eyespot number variation by performing a genome-wide association study with restriction site-associated DNA sequencing from multiple individuals varying in number of eyespots sampled across a freely breeding laboratory population. We found that dorsal hindwing eyespot number has a moderately high heritability of ∼0.50 and is characterized by a polygenic architecture. Previously identified genomic regions involved in eyespot development, and novel ones, display high association with dorsal hindwing eyespot number, suggesting that homolog number variation is likely determined by regulatory changes at multiple loci that build the trait, and not by variation at single master regulators or input-output genes.

Admixture mapping in interspecific Populus hybrids identifies classes of genomic architectures for phytochemical, morphological and growth traits

The genomic architecture of functionally important traits is key to understanding the maintenance of reproductive barriers and trait differences when divergent populations or species hybridize. We conducted a genome-wide association study (GWAS) to study trait architecture in natural hybrids of two ecologically divergent Populus species. We genotyped 472 seedlings from a natural hybrid zone of Populus alba and Populus tremula for genome-wide markers from reduced representation sequencing, phenotyped the plants in common gardens for 46 phytochemical (phenylpropanoid), morphological and growth traits, and used a Bayesian polygenic model for mapping. We detected three classes of genomic architectures: traits with finite, detectable associations of genetic loci with phenotypic variation in addition to highly polygenic heritability; traits with indications for polygenic heritability only; and traits with no detectable heritability. For the first class, we identified genome regions with plausible candidate genes for phenylpropanoid biosynthesis or its regulation, including MYB transcription factors and glycosyl transferases. GWAS in natural, recombinant hybrids represent a promising step towards resolving the genomic architecture of phenotypic traits in long-lived species. This facilitates the fine-mapping and subsequent functional characterization of genes and networks causing differences in hybrid performance and fitness.

Molecular-genetic bases of plumage coloring in chicken

The color of plumage in birds is an important feature, often determining descent to a particular species or breed. It serves as a key factor in the interaction of birds with each other due to their well-developed visual perception of the surrounding world. In poultry including chickens, the color of the plumage can be treated as a genetic marker, useful for identifying breeds, populations and breeding groups with their specific traits. The origin of diverse color plumage is the result of two interrelated physical processes, chemical and optical, due to which pigment and structural colors in the color are formed. The pigment melanin, which is presented in two forms, eumelanin and pheomelanin, is widely spread in birds. The basis for the formation of melanin is the aromatic amino acid tyrosine. The process of melano-genesis involves many loci, part of the complex expression of plumage color genes. In birds, the solid black color locus encodes the melanocortin 1 receptor (MC1R), mutations in which lead to a change in receptor activation and form different variants of the E locus. Using the GWAS analysis, possible genes affecting the formation of color in chickens were detected. The biosynthesis and types of melanin are affected by the activity of the enzyme tyrosine, and mutations in the tyrosinase gene (TYR) cause albinism in different species. The formation mechanism of brown, silver, gold, lavender and a number of other shades is determined by the influence on the work of the MC1R genes and TYR specific modifier genes. Thus, locus I currently associated with the PMEL17 gene inhibits the expression of eumelanin, and the MLPH gene affects tyrosinase function. Research on the mechanisms of formation of the secondary coloring of plumage in chickens is being actively conducted nowadays. The formation of a marble feather pattern is associated with the mutation of the endothelin B2 receptor (EDNRB2), in the coding part of the gene of which a polymorphism is found associated with the mo locus. The molecular base that causes the feather banding (locus B and autosomal recessive banding) is identified. Today, only some genes that determine the color of the plumage of chickens are studied and described. Different genes can produce similar plumage patterns, and different phenotypes can be determined by the polymorphism of a single gene. Using molecular methods, you can more accurately identify these differences. This overview shows the nature of melanin coloration in birds using the example of chickens of various breeds and also attempts to systematize knowledge about the molecular-genetic mechanisms of the appearance of various types of coloration.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

A comparison of genomic islands of differentiation across three young avian species pairs

Detailed evaluations of genomic variation between sister species often reveal distinct chromosomal regions of high relative differentiation (i.e., "islands of differentiation" in F_ST ), but there is much debate regarding the causes of this pattern. We briefly review the prominent models of genomic islands of differentiation and compare patterns of genomic differentiation in three closely related pairs of New World warblers with the goal of evaluating support for the four models. Each pair (MacGillivray's/mourning warblers; Townsend's/black-throated green warblers; and Audubon's/myrtle warblers) consists of forms that were likely separated in western and eastern North American refugia during cycles of Pleistocene glaciations and have now come into contact in western Canada, where each forms a narrow hybrid zone. We show strong differences between pairs in their patterns of genomic heterogeneity in F_ST , suggesting differing selective forces and/or differing genomic responses to similar selective forces among the three pairs. Across most of the genome, levels of within-group nucleotide diversity (π_Within ) are almost as large as levels of between-group nucleotide distance (π_Between ) within each pair, suggesting recent common ancestry and/or gene flow. In two pairs, a pattern of the F_ST peaks having low π_Between suggests that selective sweeps spread between geographically differentiated groups, followed by local differentiation. This "sweep-before-differentiation" model is consistent with signatures of gene flow within the yellow-rumped warbler species complex. These findings add to our growing understanding of speciation as a complex process that can involve phases of adaptive introgression among partially differentiated populations.

Dissection by genomic and plumage variation of a geographically complex hybrid zone between two Australian non-sister parrot species, Platycercus adscitus and Platycercus eximius

The study of hybrid zones advances understanding of the speciation process, and approaches incorporating genomic data are increasingly used to draw significant conclusions about the impact of hybridisation. Despite the progress made, the complex interplay of factors that can lead to substantially variable hybridisation outcomes are still not well understood, and many systems and/or groups remain comparatively poorly studied. Our study aims to broaden the literature on avian hybrid zones, investigating a potentially geographically and temporally complex putative hybrid zone between two native Australian non-sister parrot species, the pale-headed and eastern rosellas (Platycercus adscitus and Platycercus eximius, respectively). We analysed six plumage traits and >1400 RADseq loci and detected hybrid individuals and an unexpectedly complex geographic structure. The hybrid zone is larger than previously described due to either observer bias or its movement over recent decades. It comprises different subregions where genetic and plumage signals of admixture vary markedly in their concordance. Evidence of contemporary hybridisation (later generation and backcrossed individuals) both within and beyond the previously defined zone, when coupled with a lack of F1 hybrids and differential patterns of introgression among potentially diagnostic loci, indicates a lack of post-zygotic barriers to gene flow between species. Despite ongoing gene flow, species boundaries are likely maintained largely by strong pre-mating barriers. These findings are discussed in detail and future avenues for research into this system are proposed, which would be of benefit to the speciation and hybrid zone literature.

Admixture mapping in a hybrid zone reveals loci associated with avian feather coloration

Identifying the genetic bases for colour patterns has provided important insights into the control and expression of pigmentation and how these characteristics influence fitness. However, much more is known about the genetic bases for traits based on melanin pigments than for traits based on another major class of pigments, carotenoids. Here, we use natural admixture in a hybrid zone between Audubon's and myrtle warblers (Setophaga coronata auduboni/S. c. coronata) to identify genomic regions associated with both types of pigmentation. Warblers are known for rapid speciation and dramatic differences in plumage. For each of five plumage coloration traits, we found highly significant associations with multiple single-nucleotide polymorphisms (SNPs) across the genome and these were clustered in discrete regions. Regions near significantly associated SNPs were enriched for genes associated with keratin filaments, fibrils that make up feathers. A carotenoid-based trait that differs between the taxa-throat colour-had more than a dozen genomic regions of association. One cluster of SNPs for this trait overlaps the Scavenger Receptor Class F Member 2 (SCARF2) gene. Other scavenger receptors are presumed to be expressed at target tissues and involved in the selective movement of carotenoids into the target cells, making SCARF2 a plausible new candidate for carotenoid processing. In addition, two melanin-based plumage traits-colours of the eye line and eye spot-show very strong associations with a single genomic region mapping to chromosome 20 in the zebra finch. These findings indicate that only a subset of the genomic regions differentiated between these two warblers are associated with the plumage differences between them and demonstrate the utility of reduced-representation genomic scans in hybrid zones.