Heterogeneity and concordance in locus-specific differentiation and introgression between species of towhees.

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well-differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping-by-sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome-wide expectation using a Bayesian genomic cline model. We also estimated locus-specific FST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus-specific introgression in the two independent hybrid zones (6-20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.

Phenotypic differentiation is associated with divergent sexual selection among closely related barn swallow populations.

Sexual selection plays a key role in the diversification of numerous animal clades and may accelerate trait divergence during speciation. However, much of our understanding of this process comes from phylogenetic comparative studies, which rely on surrogate measures such as dimorphism that may not represent selection in wild populations. In this study, we assess sexual selection pressures for multiple male visual signals across four barn swallow (Hirundo rustica) populations. Our sample encompassed 2400 linear km and two described subspecies: European H. r. rustica (in the Czech Republic and Romania) and eastern Mediterranean H. r. transitiva (in Israel), as well as a potential area of contact (in Turkey). We demonstrate significant phenotypic differentiation in four sexual signalling axes, despite very low-level genomic divergence and no comparable divergence in an ecological trait. Moreover, the direction of phenotypic divergence is consistent with differences in sexual selection pressures among subspecies. Thus, H. r. transitiva, which have the darkest ventral plumage of any population, experience directional selection for darker plumage. Similarly, H. r. rustica, which have the longest tail feathers of any population, experience directional selection for elongated tail feathers and disruptive selection for ventral plumage saturation. These results suggest that sexual selection is the primary driver of phenotypic differentiation in this species. Our findings add to growing evidence of phenotypic divergence with gene flow. However, to our knowledge, this is the first study to relate direct measures of the strength and targets of sexual selection to phenotypic divergence among closely related wild populations.

Genome-wide differentiation in closely related populations: the roles of selection and geographic isolation.

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as 'isolation by adaptation' (IBA). IBA can be complementary to the well-known pattern of 'isolation by distance' (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome-wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole-genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome-wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population-level differentiation in this group of closely related populations.

When directional selection reduces geographic variation in traits mediating species interactions.

Although we often focus on the causes of geographic variation, understanding processes that act to reduce geographic variation is also important. Here, we consider a process whereby adaptive foraging across the landscape and directional selection exerted by a conifer seed predator, the common crossbill (Loxia curvirostra), potentially act to homogenize geographic variation in the defensive traits of its prey. We measured seed predation and phenotypic selection exerted by crossbills on black pine (Pinus nigra) at two sites in the Pindos Mountains, Greece. Seed predation by crossbills was over an order of magnitude higher at the site where cone scale thickness was significantly thinner, which was also the cone trait that was the target of selection at the high predation site. Additional comparisons of selection differentials demonstrate that crossbills exert selection on black pine that is consistent in form across space and time, and increases in strength with increasing seed predation. If predators distribute themselves in relation to the defensive traits of their prey and the strength of selection predators exert is proportional to the amount of predation, then predators may act to homogenize trait variation among populations of their prey in a process analogous to coevolutionary alternation with escalation.

Development of genetic diversity, differentiation and structure over 500 years in four ponderosa pine populations.

Population history plays an important role in shaping contemporary levels of genetic variation and geographic structure. This is especially true in small, isolated range-margin populations, where effects of inbreeding, genetic drift and gene flow may be more pronounced than in large continuous populations. Effects of landscape fragmentation and isolation distance may have implications for persistence of range-margin populations if they are demographic sinks. We studied four small, disjunct populations of ponderosa pine over a 500-year period. We coupled demographic data obtained through dendroecological methods with microsatellite data to discern how and when contemporary levels of allelic diversity, among and within-population levels of differentiation, and geographic structure, arose. Alleles accumulated rapidly following initial colonization, demonstrating proportionally high levels of gene flow into the populations. At population sizes of approximately 100 individuals, allele accumulation saturated. Levels of genetic differentiation among populations (F(ST) and Jost's D(est)) and diversity within populations (F(IS)) remained stable through time. There was no evidence of geographic genetic structure at any time in the populations' history. Proportionally, high gene flow in the early stages of population growth resulted in rapid accumulation of alleles and quickly created relatively homogenous genetic patterns among populations. Our study demonstrates that contemporary levels of genetic diversity were formed quickly and early in population development. How contemporary genetic diversity accumulates over time is a key facet of understanding population growth and development. This is especially relevant given the extent and speed at which species ranges are predicted to shift in the coming century.

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins.

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation (FST) and introgression (genomic cline center [α] and rate [ß]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of FST, α, and ß were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated FST, α, and ß were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.

The n = 1 constraint in population genomics.

A key objective of population genomics is to identify portions of the genome that have been shaped by natural selection rather than by neutral divergence. A previously recognized but underappreciated challenge to this objective is that observations of allele frequencies across genomes in natural populations often correspond to a single, unreplicated instance of the outcome of evolution. This is because the composition of each individual genomic region and population is expected to be the outcome of a unique array of evolutionary processes. Given a single observation, inference of the evolutionary processes that led to the observed state of a locus is associated with considerable uncertainty. This constraint on inference can be ameliorated by utilizing multi-allelic (e.g. DNA haplotypes) rather than bi-allelic markers, by analysing two or more populations with certain models and by utilizing studies of replicated experimental evolution. Future progress in population genomics will follow from research that recognizes the 'n = 1 constraint' and that utilizes appropriate and explicit evolutionary models for analysis.

Coevolution between crossbills and black pine: the importance of competitors, forest area and resource stability.

Studies of predator-prey interactions have found that geographically structured coevolution has played an important role in the adaptive diversification of crossbills (Loxia spp.). We extend those studies by considering common crossbills (L. curvirostra) in the Mediterranean where they rely on seeds in the cones of black pine (Pinus nigra). On the continent, where tree squirrels (Sciurus vulgaris) are present, enhanced defenses against crossbills were most evident in larger areas of pine forest. On islands in the absence of tree squirrels, crossbills and black pine have coevolved in a predator-prey arms race on Cyprus but not Corsica. In contrast to other conifers that island endemic crossbills rely upon, black pine does not hold seeds in its cones year round. Consequently, key to the strong crossbill-pine interaction on Cyprus is likely the presence of an alternative conifer that provides seeds during early summer when black pine seeds are scarce.