Speciation by perception

Experimentally testing mate preference in an avian system with unidirectional bill color introgression

Mating behavior can play a key role in speciation by inhibiting or facilitating gene flow between closely related taxa. Hybrid zones facilitate a direct examination of mating behavior and the traits involved in establishing species barriers. The long-tailed finch (Poephila acuticauda) has two hybridizing subspecies that differ in bill color (red and yellow), and the yellow bill phenotype appears to have introgressed ~350 km eastward following secondary contact. To examine the role of mate choice on bill color introgression, we performed behavioral assays using natural and manipulated bill colors. We found an assortative female mating preference for males of their own subspecies when bill color was not manipulated. However, we did not find this assortative preference in trials based on artificially manipulated bill color. This could suggest that assortative preference is not fixed entirely on bill color and instead may be based on a different trait (e.g., song) or a combination of traits, or alternatively may be due to lower statistical power alongside the bill manipulations being unconvincing to the female choosers. Intriguingly, we find a bias in the inheritance of bill color in captive bred F1 hybrid females. Previous modeling suggests that assortative mate preference and this kind of partial dominance in the underlying genes may together contribute to introgression, making the genetic architecture of bill color in this system a priority for future research.

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.

Extra-pair paternity as a strategy to reduce the costs of heterospecific reproduction? Insights from the crow hybrid zone

Within hybrid zones of socially monogamous species, the number of mating opportunities with a conspecific can be limited. As a consequence, individuals may mate with a heterospecific (social) partner despite possible fitness costs to their hybrid offspring. Extra-pair copulations with a conspecific may thus arise as a possible post hoc strategy to reduce the costs of hybridization. We here assessed the rate of extra-pair paternity in the hybrid zone between all-black carrion crows (Corvus (corone) corone) and grey hooded crows (C. (c.) cornix) and tested whether extra-pair paternity (EPP) was more likely in broods where parents differed in plumage colour. The proportion of broods with at least one extra-pair offspring and the proportion of extra-pair offspring were low overall (6.98% and 2.90%, respectively) with no evidence of hybrid broods having higher EPP rates than purebred nests.

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.

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.

An extensive candidate gene approach to speciation: diversity, divergence and linkage disequilibrium in candidate pigmentation genes across the European crow hybrid zone

Colouration patterns have an important role in adaptation and speciation. The European crow system, in which all-black carrion crows and grey-coated hooded crows meet in a narrow hybrid zone, is a prominent example. The marked phenotypic difference is maintained by assortative mating in the absence of neutral genetic divergence, suggesting the presence of few pigmentation genes of major effect. We made use of the rich phenotypic and genetic resources in mammals and identified a comprehensive panel of 95 candidate pigmentation genes for birds. Based on functional annotation, we chose a subset of the most promising 37 candidates, for which we developed a marker system that demonstrably works across the avian phylogeny. In total, we sequenced 107 amplicons (∼3 loci per gene, totalling 60 kb) in population samples of crows (n=23 for each taxon). Tajima's D, Fu's FS, DHEW and HKA (Hudson-Kreitman-Aguade) statistics revealed several amplicons that deviated from neutrality; however, none of these showed significantly elevated differentiation between the two taxa. Hence, colour divergence in this system may be mediated by uncharacterized pigmentation genes or regulatory regions outside genes. Alternatively, the observed high population recombination rate (4Ner∼0.03), with overall linkage disequilibrium dropping rapidly within the order of few 100 bp, may compromise the power to detect causal loci with nearby markers. Our results add to the debate as to the utility of candidate gene approaches in relation to genomic features and the genetic architecture of the phenotypic trait in question.

Nucleotide divergence vs. gene expression differentiation: comparative transcriptome sequencing in natural isolates from the carrion crow and its hybrid zone with the hooded crow

Recent advances in sequencing technology promise to provide new strategies for studying population differentiation and speciation phenomena in their earliest phases. We focus here on the black carrion crow (Corvus [corone] corone), which forms a zone of hybridization and overlap with the grey coated hooded crow (Corvus [corone] cornix). However, although these semispecies are taxonomically distinct, previous analyses based on several types of genetic markers did not reveal significant molecular differentiation between them. We here corroborate this result with sequence data obtained from a set of 25 nuclear intronic loci. Thus, the system represents a case of a very early phase of species divergence that requires new molecular approaches for its description. We have therefore generated RNAseq expression profiles using barcoded massively parallel pyrosequencing of brain mRNA from six individuals of the carrion crow and five individuals from a hybrid zone with the hooded crow. We obtained 856 675 reads from two runs, with average read length of 270 nt and coverage of 8.44. Reads were assembled de novo into 19 552 contigs, 70% of which could be assigned to annotated genes in chicken and zebra finch. This resulted in a total of 7637 orthologous genes and a core set of 1301 genes that could be compared across all individuals. We find a clear clustering of expression profiles for the pure carrion crow animals and disperse profiles for the animals from the hybrid zone. These results suggest that gene expression differences may indeed be a sensitive indicator of initial species divergence.

Geographic and temporal variation in moth chemical communication

In moth pheromone communication signals, both quantitative and qualitative intraspecific differences have been found across geographic regions. Such variation has generally been hypothesized to be due to selection, but evidence of genetic control of these differences is largely lacking. To explore the patterns of variation in pheromone signals, we quantified variation in the female sex pheromone blend and male responses of two closely related noctuid moth species in five different geographic regions for 2-3 consecutive years. We found significant variation in the ratios of sex pheromone blend components as well as in male response, not only between geographic regions but also within a region between consecutive years. The temporal variation was of a similar magnitude as the geographic variation. As far as we know, this is the first study reporting such temporal variation in moth chemical communication systems. The geographic variation seems to at least partly be controlled by genetic factors, and to be correlated with the quality of the local chemical environment. However, the pattern of temporal variation within populations suggests that optimization of the pheromonal signal also may be driven by within-generation physiological adjustments by the moths in response to their experience of the local chemical environment.