Transcriptomics of colour patterning and coloration shifts in crows

The genetic consequences of historic climate change on the contemporary population structure of a widespread temperate North American songbird

Studies of widely distributed species can inform our understanding of how past demographic events tied to historic glaciation and ongoing population genetic processes interact to shape contemporaneous patterns of biodiversity at a continental scale. In this study, we used whole-genome resequencing to investigate the current population structure and genetic signatures of past demographic events in the widespread migratory American goldfinch (Spinus tristis). Phylogenetic relationships inferred from whole mitochondrial genomes were poorly resolved. In contrast, a genome-wide panel of > 4.5 million single nucleotide polymorphisms (SNPs) strongly supported the existence of eastern and western populations separated by western mountain ranges and additional population structuring within the western clade. Demographic modeling estimated that the eastern and western populations diverged approximately one million years ago, and both populations experienced subsequent population bottlenecks during the last glacial period. Species distribution models showed a severe contraction of suitable habitat for the American goldfinch during this period, with predicted discontinuities that are consistent with multiple, isolated glacial refugia that coincide with present-day population structure. Low overall genetic differentiation between the eastern and western populations (FST ∼ 0.01) suggests ongoing gene flow accompanied divergence, and individuals with admixed genomic signatures were sampled along a potential contact zone. Nevertheless, outlier SNPs were identified near genes associated with feather color, song, and migratory behavior and provide strong candidates for further study of the mechanisms underlying reproductive isolation and speciation in birds.

Structure, function and formation of the amniote skin pattern

From feather and hair dotted arrays to pigmented stripes and spots, the spatial distribution of skin appendages and colouration often forms visible ornaments crucial for fitness in the coat of birds and mammals. These geometrical motifs are extremely diverse in nature. Yet, phenotypic surveys evidenced common themes in variation: the orientation, appendage-specificity or pigmentation of a given region may be conserved across groups or species. Here, we review naturalist observations of natural variation in the anatomy and ecological function of the skin pattern in amniotes. We then describe several decades of genetics, mathematical modelling and experimental embryology work aiming at understanding the molecular and morphogenetic mechanisms responsible for pattern formation. We discuss how these studies provided evidence that the morphological trends and differences representative of the phenotypic landscape of skin patterns in wild amniote species is rooted in the mechanisms controlling the production of distinct compartments in the embryonic skin.

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

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

Adaptive phenotypic and genomic divergence in the common chaffinch (Fringilla coelebs) following niche expansion within a small oceanic island

According to models of ecological speciation, adaptation to adjacent, contrasting habitat types can lead to population divergence given strong enough environment-driven selection to counteract the homogenizing effect of gene flow. We tested this hypothesis in the common chaffinch (Fringilla coelebs) on the small island of La Palma, Canary Islands, where it occupies two markedly different habitats. Isotopic (δ13 C, δ15 N) analysis of feathers indicated that birds in the two habitats differed in ecosystem and/or diet, and analysis of phenotypic traits revealed significant differences in morphology and plumage colouration that are consistent with ecomorphological and ecogeographical predictions respectively. A genome-wide survey of single-nucleotide polymorphism revealed marked neutral structure that was consistent with geography and isolation by distance, suggesting low dispersal. In contrast, loci putatively under selection identified through genome-wide association and genotype-environment association analyses, revealed amarked adaptive divergence between birds in both habitats. Loci associated with phenotypic and environmental differences among habitats were distributed across the genome, as expected for polygenic traits involved in local adaptation. Our results suggest a strong role for habitat-driven local adaptation in population divergence in the chaffinches of La Palma, a process that appears to be facilitated by a strong reduction in effective dispersal distances despite the birds' high dispersal capacity.

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5 Collapse

Key Words Collapse

MESH Headings

• Animals
• Female
• Male
• Chromatin/genetics
• Crows/genetics
• Epigenesis, Genetic
• Methylation
• Dosage Compensation, Genetic
• Sex Chromosomes

Collapse

Grants

• HORIZON EUROPE European Research Council
• Swedish Research Council
• Knut och Alice Wallenbergs Stiftelse

Collapse

Affiliation(s)

Ana Catalán Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Justin Merondun Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Ulrich Knief Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany Evolutionary Biology & Ecology,Faculty of Biology, University of Freiburg, Freiburg, Germany
Jochen B. W. Wolf Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany

Collapse

A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility

Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.

A conserved molecular template underlies color pattern diversity in estrildid finches

The color patterns that adorn animals' coats not only exhibit extensive diversity linked to various ecological functions but also display recurrences in geometry, orientation, or body location. How processes of pattern formation shape such phenotypic trends remains a mystery. Here, we surveyed plumage color patterns in passerine finches displaying extreme apparent variation and identified a conserved set of color domains. We linked these domains to putative embryonic skin regions instructed by early developmental tissues and outlined by the combinatory expression of few genetic markers. We found that this embryonic prepattern is largely conserved in birds displaying drastic color differences in the adult, interspecies variation resulting from the masking or display of each domain depending on their coloration. This work showed that a simple molecular landscape serves as common spatial template to extensive color pattern variation in finches, revealing that early conserved landmarks and molecular pathways are a major cause of phenotypic trends.

The evolutionary history and mechanistic basis of female ornamentation in a tropical songbird

Ornamentation, such as the showy plumage of birds, is widespread among female vertebrates, yet the evolutionary pressures shaping female ornamentation remain uncertain. In part this is due to a poor understanding of the mechanistic route to ornamentation in females. To address this issue, we evaluated the evolutionary history of ornament expression in a tropical passerine bird, the White-shouldered Fairywren, whose females, but not males, strongly vary between populations in occurrence of ornamented black-and-white plumage. We first use phylogenomic analysis to demonstrate that female ornamentation is derived and that female ornamentation evolves independently of changes in male plumage. We then use exogenous testosterone in a field experiment to induce partial ornamentation in naturally unornamented females. By sequencing the transcriptome of experimentally induced ornamented and natural feathers, we identify genes expressed during ornament production and evaluate the degree to which female ornamentation in this system is associated with elevated testosterone, as is common in males. We reveal that some ornamentation in females is linked to testosterone and that sexes differ in ornament-linked gene expression. Lastly, using genomic outlier analysis we identify a candidate melanogenesis gene that lies in a region of high genomic divergence among populations that is also differentially expressed in feather follicles of different female plumages. Taken together, these findings are consistent with sex-specific selection favoring the evolution of female ornaments and demonstrate a key role for testosterone in generating population divergence in female ornamentation through gene regulation. More broadly, our work highlights similarities and differences in how ornamentation evolves in the sexes.

Ineffective silencing of transposable elements on an avian W Chromosome

One of the defining features of transposable elements (TEs) is their ability to move to new locations in the host genome. To minimise the potentially deleterious effects of de novo TE insertions, hosts have evolved several mechanisms to control TE activity, including recombination-mediated removal and epigenetic silencing; however, increasing evidence suggests that silencing of TEs is often incomplete. The crow family experienced a recent radiation of LTR retrotransposons (LTRs), offering an opportunity to gain insight into the regulatory control of young, potentially still active TEs. We quantified the abundance of TE-derived transcripts across several tissues in 15 Eurasian crows (Corvus (corone) spp.) raised under common garden conditions and find evidence for ineffective TE suppression on the female-specific W Chromosome. Using RNA-seq data, we show that ~ 9.5% of all transcribed TEs had considerably greater (average: 16-fold) transcript abundance in female crows, and that more than 85% of these female-biased TEs originated on the W Chromosome. After accounting for differences in TE density among chromosomal classes, W-linked TEs were significantly more highly expressed than TEs residing on other chromosomes, consistent with ineffective silencing on the former. Together, our results suggest that the crow W Chromosome acts as a source of transcriptionally active TEs, with possible negative fitness consequences for female birds analogous to Drosophila (an X/Y system), where overexpression of Y-linked TEs is associated with male-specific aging and fitness loss ('toxic Y').

Unraveling patterns of disrupted gene expression across a complex tissue

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

Rapid radiation in a highly diverse marine environment

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

A ROR2 coding variant is associated with craniofacial variation in domestic pigeons

Vertebrate craniofacial morphogenesis is a highly orchestrated process that is directed by evolutionarily conserved developmental pathways.1,2 Within species, canalized development typically produces modest morphological variation. However, as a result of millennia of artificial selection, the domestic pigeon displays radical craniofacial variation within a single species. One of the most striking cases of pigeon craniofacial variation is the short-beak phenotype, which has been selected in numerous breeds. Classical genetic experiments suggest that pigeon beak length is regulated by a small number of genetic factors, one of which is sex linked (Ku2 locus).3-5 However, the genetic underpinnings of pigeon craniofacial variation remain unknown. Using geometric morphometrics and quantitative trait locus (QTL) mapping on an F2 intercross between a short-beaked Old German Owl (OGO) and a medium-beaked Racing Homer (RH), we identified a single Z chromosome locus that explains a majority of the variation in beak morphology in the F2 population. Complementary comparative genomic analyses revealed that the same locus is strongly differentiated between breeds with short and medium beaks. Within the Ku2 locus, we identified an amino acid substitution in the non-canonical Wnt receptor ROR2 as a putative regulator of pigeon beak length. The non-canonical Wnt pathway serves critical roles in vertebrate neural crest cell migration and craniofacial morphogenesis.6,7 In humans, ROR2 mutations cause Robinow syndrome, a congenital disorder characterized by skeletal abnormalities, including a widened and shortened facial skeleton.8,9 Our results illustrate how the extraordinary craniofacial variation among pigeons can reveal genetic regulators of vertebrate craniofacial diversity.

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

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

Beyond the Chicken: Alternative Avian Models for Developmental Physiological Research

Biomedical research focusing on physiological, morphological, behavioral, and other aspects of development has long depended upon the chicken (Gallus gallus domesticus) as a key animal model that is presumed to be typical of birds and generally applicable to mammals. Yet, the modern chicken in its many forms is the result of artificial selection more intense than almost any other domesticated animal. A consequence of great variation in genotype and phenotype is that some breeds have inherent aberrant physiological and morphological traits that may show up relatively early in development (e.g., hypertension, hyperglycemia, and limb defects in the broiler chickens). While such traits can be useful as models of specific diseases, this high degree of specialization can color general experimental results and affect their translational value. Against this background, in this review we first consider the characteristics that make an animal model attractive for developmental research (e.g., accessibility, ease of rearing, size, fecundity, development rates, genetic variation, etc.). We then explore opportunities presented by the embryo to adult continuum of alternative bird models, including quail, ratites, songbirds, birds of prey, and corvids. We conclude by indicating that expanding developmental studies beyond the chicken model to include additional avian groups will both validate the chicken model as well as potentially identify even more suitable avian models for answering questions applicable to both basic biology and the human condition.

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.

Pooled Sequencing Analysis of Geese (Anser cygnoides) Reveals Genomic Variations Associated With Feather Color

During the domestication of the goose a change in its feather color took place, however, the molecular mechanisms responsible for this change are not completely understood. Here, we performed whole-genome resequencing on three pooled samples of geese (feral and domestic geese), with two distinct feather colors, to identify genes that might regulate feather color. We identified around 8 million SNPs within each of the three pools and validated allele frequencies for a subset of these SNPs using PCR and Sanger sequencing. Several genomic regions with signatures of differential selection were found when we compared the gray and white feather color populations using the F_ST and Hp approaches. When we combined previous functional studies with our genomic analyses we identified 26 genes (KITLG, MITF, TYRO3, KIT, AP3B1, SMARCA2, ROR2, CSNK1G3, CCDC112, VAMP7, SLC16A2, LOC106047519, RLIM, KIAA2022, ST8SIA4, LOC106044163, TRPM6, TICAM2, LOC106038556, LOC106038575, LOC106038574, LOC106038594, LOC106038573, LOC106038604, LOC106047489, and LOC106047492) that potentially regulate feather color in geese. These results substantially expand the catalog of potential feather color regulators in geese and provide a basis for further studies on domestication and avian feather coloration.

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.

Untangling the structural and molecular mechanisms underlying colour and rapid colour change in a lizard, Agama atra

With functions as diverse as communication, protection and thermoregulation, coloration is one of the most important traits in lizards. The ability to change colour as a function of varying social and environmental conditions is thus an important innovation. While colour change is present in animals ranging from squids, to fish and reptiles, not much is known about the mechanisms behind it. Traditionally, colour change was attributed to migration of pigments, in particular melanin. More recent work has shown that the changes in nanostructural configuration inside iridophores are able to produce a wide palette of colours. However, the genetic mechanisms underlying colour, and colour change in particular, remain unstudied. Here we use a combination of transcriptomic and microscopic data to show that melanin, iridophores and pteridines are the main colour-producing mechanisms in Agama atra, and provide molecular and structural data suggesting that rapid colour change is achieved via melanin dispersal in combination with iridophore organization. This work demonstrates the power of combining genotypic (gene expression) and phenotypic (microscopy) information for addressing physiological questions, providing a basis for future studies of colour change.

Rapid speciation via the evolution of pre-mating isolation in the Iberá Seedeater

Behavioral isolation can catalyze speciation and permit the slow accumulation of additional reproductive barriers between co-occurring organisms. We illustrate how this process occurs by examining the genomic and behavioral bases of pre-mating isolation between two bird species (Sporophila hypoxantha and the recently discovered S. iberaensis) that belong to the southern capuchino seedeaters, a recent, rapid radiation characterized by variation in male plumage coloration and song. Although these two species co-occur without obvious ecological barriers to reproduction, we document behaviors indicating species recognition by song and plumage traits and strong assortative mating associated with genomic regions underlying male plumage patterning. Plumage differentiation likely originated through the reassembly of standing genetic variation, indicating how novel sexual signals may quickly arise and maintain species boundaries.

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.

Discovery and population genomics of structural variation in a songbird genus

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping.

Differential divergence in autosomes and sex chromosomes is associated with intra-island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.

Epistatic mutations under divergent selection govern phenotypic variation in the crow hybrid zone

The evolution of genetic barriers opposing inter-specific gene flow is key to the origin of new species. Drawing from information of over 400 admixed genomes sourced from replicate transects across the European hybrid zone between all-black carrion crows and grey-coated hooded crows, we decipher the interplay between phenotypic divergence and selection at the molecular level. Over 68% of plumage variation was explained by epistasis between the gene NDP and a ~2.8 Mb region on chromosome 18 with suppressed recombination. Both pigmentation loci showed evidence for divergent selection resisting introgression. This study reveals how few, large-effect loci can govern prezygotic isolation and shield phenotypic divergence from gene flow.

In situ quantification of individual mRNA transcripts in melanocytes discloses gene regulation of relevance to speciation

Functional validation of candidate genes involved in adaptation and speciation remains challenging. Here, we exemplify the utility of a method quantifying individual mRNA transcripts in revealing the molecular basis of divergence in feather pigment synthesis during early-stage speciation in crows. Using a padlock probe assay combined with rolling circle amplification, we quantified cell-type-specific gene expression in the histological context of growing feather follicles. Expression of Tyrosinase Related Protein 1 (TYRP1), Solute Carrier Family 45 member 2 (SLC45A2) and Hematopoietic Prostaglandin D Synthase (HPGDS) was melanocyte-limited and significantly reduced in follicles from hooded crow, explaining the substantially lower eumelanin content in grey versus black feathers. The central upstream Melanocyte Inducing Transcription Factor (MITF) only showed differential expression specific to melanocytes - a feature not captured by bulk RNA-seq. Overall, this study provides insight into the molecular basis of an evolutionary young transition in pigment synthesis, and demonstrates the power of histologically explicit, statistically substantiated single-cell gene expression quantification for functional genetic inference in natural populations.

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

Recombination between loci underlying mate choice and ecological traits is a major evolutionary force acting against speciation with gene flow. The evolution of linkage disequilibrium between such loci is therefore a fundamental step in the origin of species. Here, we show that this process can take place in the absence of physical linkage in hamlets-a group of closely related reef fishes from the wider Caribbean that differ essentially in colour pattern and are reproductively isolated through strong visually-based assortative mating. Using full-genome analysis, we identify four narrow genomic intervals that are consistently differentiated among sympatric species in a backdrop of extremely low genomic divergence. These four intervals include genes involved in pigmentation (sox10), axial patterning (hoxc13a), photoreceptor development (casz1) and visual sensitivity (SWS and LWS opsins) that develop islands of long-distance and inter-chromosomal linkage disequilibrium as species diverge. The relatively simple genomic architecture of species differences facilitates the evolution of linkage disequilibrium in the presence of gene flow.

Introgression of regulatory alleles and a missense coding mutation drive plumage pattern diversity in the rock pigeon

Birds and other vertebrates display stunning variation in pigmentation patterning, yet the genes controlling this diversity remain largely unknown. Rock pigeons (Columba livia) are fundamentally one of four color pattern phenotypes, in decreasing order of melanism: T-check, checker, bar (ancestral), or barless. Using whole-genome scans, we identified NDP as a candidate gene for this variation. Allele-specific expression differences in NDP indicate cis-regulatory divergence between ancestral and melanistic alleles. Sequence comparisons suggest that derived alleles originated in the speckled pigeon (Columba guinea), providing a striking example of introgression. In contrast, barless rock pigeons have an increased incidence of vision defects and, like human families with hereditary blindness, carry start-codon mutations in NDP. In summary, we find that both coding and regulatory variation in the same gene drives wing pattern diversity, and post-domestication introgression supplied potentially advantageous melanistic alleles to feral populations of this ubiquitous urban bird.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.

Dynamic regulation of mRNA and miRNA associated with the developmental stages of skin pigmentation in Japanese ornamental carp

The Japanese ornamental carp (Cyprinus carpio var. Koi) is famous for multifarious colors and patterns, making it commonly culture and trade across the world. Although functional genes and inheritance of color traits have been commonly studied, seldom attentions were focused on the genetic regulation during the developmental process of pigmentation. To better understand the mechanism of skin color development, we observed the morphogenesis of pigment cells during the post-embryonic stages and analysed the temporal expression pattern of mRNAs/miRNAs profiles in four distinct developmental stages. 59 and 103 differentially expressed genes/miRNAs (DEGs/DEMs) associated with pigmentation and skin were identified, including pax7, mitf, tyr, tyrp1, etc., and the highest DEGs were detected at 11 days post hatching (dph). In addition, the functional characteristics of mRNAs/miRNAs associated with pteridine and carotenoid pathway were also examined. Furthermore, 65 miRNA-mRNA interaction pairs related to pigmentation, pteridines and carotenoids metabolism were detected between different stages. Interestingly, the largest pairs appeared in the transition from 11 dph to 48 dph, which had the similar trend with DEGs further manifesting the importance of 11 dph. This study produced a comprehensive programme of DEGs/DEMs during color development, which will provide resources to understand the regulation mechanism in color formation. The understanding of genetic basis in color formation might promote the production and breeding of the Koi carp.

Genomics of coloration in natural animal populations

Animal coloration has traditionally been the target of genetic and evolutionary studies. However, until very recently, the study of the genetic basis of animal coloration has been mainly restricted to model species, whereas research on non-model species has been either neglected or mainly based on candidate approaches, and thereby limited by the knowledge obtained in model species. Recent high-throughput sequencing technologies allow us to overcome previous limitations, and open new avenues to study the genetic basis of animal coloration in a broader number of species and colour traits, and to address the general relevance of different genetic structures and their implications for the evolution of colour. In this review, we highlight aspects where genome-wide studies could be of major utility to fill in the gaps in our understanding of the biology and evolution of animal coloration. The new genomic approaches have been promptly adopted to study animal coloration although substantial work is still needed to consider a larger range of species and colour traits, such as those exhibiting continuous variation or based on reflective structures. We argue that a robust advancement in the study of animal coloration will also require large efforts to validate the functional role of the genes and variants discovered using genome-wide tools.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.

Challenges and advances for transcriptome assembly in non-model species

Analyses of high-throughput transcriptome sequences of non-model organisms are based on two main approaches: de novo assembly and genome-guided assembly using mapping to assign reads prior to assembly. Given the limits of mapping reads to a reference when it is highly divergent, as is frequently the case for non-model species, we evaluate whether using blastn would outperform mapping methods for read assignment in such situations (>15% divergence). We demonstrate its high performance by using simulated reads of lengths corresponding to those generated by the most common sequencing platforms, and over a realistic range of genetic divergence (0% to 30% divergence). Here we focus on gene identification and not on resolving the whole set of transcripts (i.e. the complete transcriptome). For simulated datasets, the transcriptome-guided assembly based on blastn recovers 94.8% of genes irrespective of read length at 0% divergence; however, assignment rate of reads is negatively correlated with both increasing divergence level and reducing read lengths. Nevertheless, we still observe 92.6% of recovered genes at 30% divergence irrespective of read length. This analysis also produces a categorization of genes relative to their assignment, and suggests guidelines for data processing prior to analyses of comparative transcriptomics and gene expression to minimize potential inferential bias associated with incorrect transcript assignment. We also compare the performances of de novo assembly alone vs in combination with a transcriptome-guided assembly based on blastn both via simulation and empirically, using data from a cyprinid fish species and from an oak species. For any simulated scenario, the transcriptome-guided assembly using blastn outperforms the de novo approach alone, including when the divergence level is beyond the reach of traditional mapping methods. Combining de novo assembly and a related reference transcriptome for read assignment also addresses the bias/error in contigs caused by the dependence on a related reference alone. Empirical data corroborate these findings when assembling transcriptomes from the two non-model organisms: Parachondrostoma toxostoma (fish) and Quercus pubescens (plant). For the fish species, out of the 31,944 genes known from D. rerio, the guided and de novo assemblies recover respectively 20,605 and 20,032 genes but the performance of the guided assembly approach is much higher for both the contiguity and completeness metrics. For the oak, out of the 29,971 genes known from Vitis vinifera, the transcriptome-guided and de novo assemblies display similar performance, but the new guided approach detects 16,326 genes where the de novo assembly only detects 9,385 genes.

Genomewide patterns of variation in genetic diversity are shared among populations, species and higher-order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST , PBS, d_xy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.

The genomic mosaicism of hybrid speciation

Hybridization is widespread in nature and, in some instances, can result in the formation of a new hybrid species. We investigate the genetic foundation of this poorly understood process through whole-genome analysis of the hybrid Italian sparrow and its progenitors. We find overall balanced yet heterogeneous levels of contribution from each parent species throughout the hybrid genome and identify areas of novel divergence in the hybrid species exhibiting signals consistent with balancing selection. High-divergence areas are disproportionately located on the Z chromosome and overrepresented in gene networks relating to key traits separating the focal species, which are likely involved in reproductive barriers and/or species-specific adaptations. Of special interest are genes and functional groups known to affect body patterning, beak morphology, and the immune system, which are important features of diversification and fitness. We show that a combination of mosaic parental inheritance and novel divergence within the hybrid lineage has facilitated the origin and maintenance of an avian hybrid species.

Combination of short-read, long-read, and optical mapping assemblies reveals large-scale tandem repeat arrays with population genetic implications

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

Covariation in levels of nucleotide diversity in homologous regions of the avian genome long after completion of lineage sorting

Closely related species may show similar levels of genetic diversity in homologous regions of the genome owing to shared ancestral variation still segregating in the extant species. However, after completion of lineage sorting, such covariation is not necessarily expected. On the other hand, if the processes that govern genetic diversity are conserved, diversity may potentially covary even among distantly related species. We mapped regions of conserved synteny between the genomes of two divergent bird species-collared flycatcher and hooded crow-and identified more than 600 Mb of homologous regions (66% of the genome). From analyses of whole-genome resequencing data in large population samples of both species we found nucleotide diversity in 200 kb windows to be well correlated (Spearman's ρ = 0.407). The correlation remained highly similar after excluding coding sequences. To explain this covariation, we suggest that a stable avian karyotype and a conserved landscape of recombination rate variation render the diversity-reducing effects of linked selection similar in divergent bird lineages. Principal component regression analysis of several potential explanatory variables driving heterogeneity in flycatcher diversity levels revealed the strongest effects from recombination rate variation and density of coding sequence targets for selection, consistent with linked selection. It is also possible that a stable karyotype is associated with a conserved genomic mutation environment contributing to covariation in diversity levels between lineages. Our observations imply that genetic diversity is to some extent predictable.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a 'brown' ancestral population, the dominant 'grey' allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.

Identification of Submergence-Responsive MicroRNAs and Their Targets Reveals Complex MiRNA-Mediated Regulatory Networks in Lotus (Nelumbo nucifera Gaertn)

MicroRNAs (miRNAs) are endogenous non-coding RNAs with important regulatory functions in plant development and stress responses. However, their population abundance in lotus (Nelumbo nucifera Gaertn) has so far been poorly described, particularly in response to stresses. In this work, submergence-related miRNAs and their target genes were systematically identified, compared, and validated at the transcriptome-wide level using high-throughput sequencing data of small RNA, Mrna, and the degradome. A total of 128 known and 20 novel miRNAs were differentially expressed upon submergence. We identified 629 target transcripts for these submergence-responsive miRNAs. Based on the miRNA expression profiles and GO and KEGG annotation of miRNA target genes, we suggest possible molecular responses and physiological changes of lotus in response to submergence. Several metabolic, physiological and morphological adaptations-related miRNAs, i.e., NNU_far-miR159, NNU_gma-miR393h, and NNU_aly-miR319c-3p, were found to play important regulatory roles in lotus response to submergence. This work will contribute to a better understanding of miRNA-regulated adaption responses of lotus to submergence stress.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

Evolution of heterogeneous genome differentiation across multiple contact zones in a crow species complex

Uncovering the genetic basis of species diversification is a central goal in evolutionary biology. Yet, the link between the accumulation of genomic changes during population divergence and the evolutionary forces promoting reproductive isolation is poorly understood. Here, we analysed 124 genomes of crow populations with various degrees of genome-wide differentiation, with parallelism of a sexually selected plumage phenotype, and ongoing hybridization. Overall, heterogeneity in genetic differentiation along the genome was best explained by linked selection exposed on a shared genome architecture. Superimposed on this common background, we identified genomic regions with signatures of selection specific to independent phenotypic contact zones. Candidate pigmentation genes with evidence for divergent selection were only partly shared, suggesting context-dependent selection on a multigenic trait architecture and parallelism by pathway rather than by repeated single-gene effects. This study provides insight into how various forms of selection shape genome-wide patterns of genomic differentiation as populations diverge.