The importance, extent, and mode of interspecific gene flow for the evolution of species has long been debated. Characterization of genomic differentiation in a classic example of hybridization between all-black carrion crows and gray-coated hooded crows identified genome-wide introgression extending far beyond the morphological hybrid zone. Gene expression divergence was concentrated in pigmentation genes expressed in gray versus black feather follicles. Only a small number of narrow genomic islands exhibited resistance to gene flow. One prominent genomic region (<2 megabases) harbored 81 of all 82 fixed differences (of 8.4 million single-nucleotide polymorphisms in total) linking genes involved in pigmentation and in visual perception-a genomic signal reflecting color-mediated prezygotic isolation. Thus, localized genomic selection can cause marked heterogeneity in introgression landscapes while maintaining phenotypic divergence.

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.

Transcriptome Shotgun Sequencing (RNA-seq) has been readily embraced by geneticists and molecular ecologists alike. As with all high-throughput technologies, it is critical to understand which analytic strategies are best suited and which parameters may bias the interpretation of the data. Here we use a comprehensive simulation approach to explore how various features of the transcriptome (complexity, degree of polymorphism π, alternative splicing), technological processing (sequencing error ε, library normalization) and bioinformatic workflow (de novo vs. mapping assembly, reference genome quality) impact transcriptome quality and inference of differential gene expression (DE). We find that transcriptome assembly and gene expression profiling (EdgeR vs. BaySeq software) works well even in the absence of a reference genome and is robust across a broad range of parameters. We advise against library normalization and in most situations advocate mapping assemblies to an annotated genome of a divergent sister clade, which generally outperformed de novo assembly (Trans-Abyss, Trinity, Soapdenovo-Trans). Transcriptome complexity (size, paralogs, alternative splicing isoforms) negatively affected the assembly and DE profiling, whereas the effects of sequencing error and polymorphism were almost negligible. Finally, we highlight the challenge of gene name assignment for de novo assemblies, the importance of mapping strategies and raise awareness of challenges associated with the quality of reference genomes. Overall, our results have significant practical and methodological implications and can provide guidance in the design and analysis of RNA-seq experiments, particularly for organisms where genomic background information is lacking.

Animal coloration is one of the most conspicuous phenotypic traits in natural populations and has important implications for adaptation and speciation. Changes in coloration can occur over surprisingly short evolutionary timescales, while recurrence of similar colour patterns across large phylogenetic distances is also common. Even though the genetic basis of pigment production is well understood, little is known about the mechanisms regulating colour patterning. In this study, we shed light on the molecular elements regulating regional pigment production in two genetically near-identical crow taxa with striking differences in a eumelanin-based phenotype: black carrion and grey-coated hooded crows. We produced a high-quality genome annotation and analysed transcriptome data from a 2 × 2 design of active melanogenic feather follicles from head (black in both taxa) and torso (black in carrion and grey in hooded crow). Extensive, parallel expression differences between body regions in both taxa, enriched for melanogenesis genes (e.g. ASIP, CORIN, and ALDH6), indicated the presence of cryptic prepatterning also in all-black carrion crows. Meanwhile, colour-specific expression (grey vs. black) was limited to a small number of melanogenesis genes in close association with the central transcription factor MITF (most notably HPGDS, NDP and RASGRF1). We conclude that colour pattern differences between the taxa likely result from an interaction between divergence in upstream elements of the melanogenesis pathway and genes that provide an underlying prepattern across the body through positional information. A model of evolutionary stable prepatterns that can be exposed and masked through simple regulatory changes may explain the phylogenetically independent recurrence of colour patterns that is observed across corvids and many other vertebrate groups.

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.

To investigate the origins and stages of vertebrate adaptive radiation, we reconstructed the spatial and temporal histories of adaptive alleles underlying major phenotypic axes of diversification from the genomes of 202 Caribbean pupfishes. On a single Bahamian island, ancient standing variation from disjunct geographic sources was reassembled into new combinations under strong directional selection for adaptation to the novel trophic niches of scale-eating and molluscivory. We found evidence for two longstanding hypotheses of adaptive radiation: hybrid swarm origins and temporal stages of adaptation. Using a combination of population genomics, transcriptomics, and genome-wide association mapping, we demonstrate that this microendemic adaptive radiation of novel trophic specialists on San Salvador Island, Bahamas experienced twice as much adaptive introgression as generalist populations on neighboring islands and that adaptive divergence occurred in stages. First, standing regulatory variation in genes associated with feeding behavior (prlh, cfap20, and rmi1) were swept to fixation by selection, then standing regulatory variation in genes associated with craniofacial and muscular development (itga5, ext1, cyp26b1, and galr2) and finally the only de novo nonsynonymous substitution in an osteogenic transcription factor and oncogene (twist1) swept to fixation most recently. Our results demonstrate how ancient alleles maintained in distinct environmental refugia can be assembled into new adaptive combinations and provide a framework for reconstructing the spatiotemporal landscape of adaptation and speciation.

As is true of virtually every realm of the biological sciences, our understanding of speciation is increasingly informed by the genomic revolution of the past decade. Investigators can ask detailed questions relating to both the extrinsic (e.g. inter-and intra-population and ecological interactions) and intrinsic (e.g. genome content and architecture) forces that drive speciation. Technologies ranging from restriction-site associated DNA sequencing (RADseq), to whole genome sequencing and assembly, to transcriptomics, to CRISPR are revolutionizing the means by which investigators can both frame and test hypotheses of lineage diversification. Our review aims to examine both extrinsic and intrinsic aspects of speciation. Genome-scale data have already served to fundamentally clarify the role of gene flow during (and after) speciation, although we predict that the differential propensity for speciation among phylogenetic lineages will be one of the most exciting frontiers for future genomic investigation. We propose that a unified theory of speciation will take into account the idiosyncratic features of genomic architecture examined in the light of each organism's biology and ecology drawn from across the full breadth of the Tree of Life.

Song divergence between closely related taxa may play a critical role in the evolutionary processes of speciation and hybridization. We explored song variation between two Ecuadorian subspecies of the gray-breasted wood-wren (Henicorhina leucophrys) and tested the impact of song divergence on response behaviors. Songs were significantly different between the two subspecies, even between two parapatric populations 10 km apart. Playback experiments revealed an asymmetric response pattern to these divergent subspecies specific songs; one subspecies responded more to songs of its own subspecies than to the other subspecies' songs, whereas the second responded equally strongly to songs of both subspecies. While song parameters revealed a mixed pattern of divergence between allopatric and parapatric populations, the majority of spectral characteristics showed increased divergence in parapatry, suggestive of character displacement. This increased song divergence in parapatry appeared to affect behavioral responses to playback as discriminating responses were most prominent in parapatry and against parapatric songs. The clear behavioral impact of subspecies-specific song differences supports a potential role for song as an acoustic barrier to gene flow. The asymmetric nature of the responses suggests that song divergence could affect the direction of gene flow and the position of the subspecies-specific transition.