Pronounced differentiation on the Z chromosome and parts of the autosomes in crowned sparrows contrasts with mitochondrial paraphyly: implications for speciation.

When a single species evolves into multiple descendent species, some parts of the genome can play a key role in the evolution of reproductive isolation while other parts flow between the evolving species via interbreeding. Genomic evolution during the speciation process is particularly interesting when major components of the genome-for instance, sex chromosomes vs. autosomes vs. mitochondrial DNA-show widely differing patterns of relationships between three diverging populations. The golden-crowned sparrow (Zonotrichia atricapilla) and the white-crowned sparrow (Zonotrichia leucophrys) are phenotypically differentiated sister species that are largely reproductively isolated despite possessing similar mitochondrial genomes, likely due to recent introgression. We assessed variation in more than 45,000 single nucleotide polymorphisms to determine the structure of nuclear genomic differentiation between these species and between two hybridizing subspecies of Z. leucophrys. The two Z. leucophrys subspecies show moderate levels of relative differentiation and patterns consistent with a history of recurrent selection in both ancestral and daughter populations, with much of the sex chromosome Z and a large region on the autosome 1A showing increased differentiation compared to the rest of the genome. The two species Z. leucophrys and Z. atricapilla show high relative differentiation and strong heterogeneity in the level of differentiation among various chromosomal regions, with a large portion of the sex chromosome (Z) showing highly divergent haplotypes between these species. Studies of speciation often emphasize mitochondrial DNA differentiation, but speciation between Z. atricapilla and Z. leucophrys appears primarily associated with Z chromosome divergence and more moderately associated with autosomal differentiation, whereas mitochondria are highly similar due apparently to recent introgression. These results add to the growing body of evidence for highly heterogeneous patterns of genomic differentiation during speciation, with some genomic regions showing a lack of gene flow between populations many hundreds of thousands of years before other genomic regions.

Geographic variability of hybridization between red-breasted and red-naped sapsuckers.

Hybrid zones reveal the strength of reproductive isolation between populations undergoing speciation and are a key tool in evolutionary biology research. Multiple replicate transects across the same hybrid zone offer insight into the dynamics of hybridization in different environments, clarifying the role of extrinsic forces on the speciation process. Red-breasted and red-naped sapsuckers (Sphyrapicus ruber and Sphyrapicus nuchalis) have a long zone of contact over approximately 1,600 km from central British Columbia, Canada to central California, USA. We used Genotyping-by-Sequencing data from three independent sapsucker hybrid zone transects to compare hybridization dynamics between these species under variable geoclimatic conditions. We generated geographic clines of the genomic data to compare hybrid zone widths and used random forests models and linear regression to assess the relationship between climate and sapsucker ancestry along each transect. Our results show variation in the directionality of backcrossing, often indicative of moving hybrid zones. We note variable cline widths among transects, indicating differences in selection maintaining hybrid zone dynamics. Furthermore, random forests models identified different variables in close association with sapsucker ancestry across each transect. These results indicate a lack of repeatability across replicate transects and a strong influence of the local environment on hybrid zone dynamics.

Population genomics of an emergent tri-species hybrid zone.

Isolating barriers that drive speciation are commonly studied in the context of two-species hybrid zones. There is, however, evidence that more complex introgressive relationships are common in nature. Here, we use field observations and genomic analysis, including the sequencing and assembly of a novel reference genome, to study an emergent hybrid zone involving two colliding hybrid zones of three woodpecker species: red-breasted, red-naped, and yellow-bellied sapsuckers (Sphyrapicus ruber, S. nuchalis, and S. varius). Surveys of the area surrounding Prince George, British Columbia, Canada, show that all three species are sympatric, and Genotyping-by-Sequencing identifies hybrids from each species pair and birds with ancestry from all three species. Observations of phenotypes and genotypes of mated pairs provide evidence for assortative mating, though there is some heterospecific pairing. Hybridization is more extensive in this tri-species hybrid zone than in two di-species hybrid zones. However, there is no evidence of a hybrid swarm and admixture is constrained to contact zones, so we classify this region as a tension zone and invoke selection against hybrids as a likely mechanism maintaining species boundaries. Analysis of sapsucker age classes does not show disadvantages in hybrid survival to adulthood, so we speculate the selection upholding the tension zone may involve hybrid fecundity. Gene flow among all sapsuckers in di-species hybrid zones suggests introgression probably occurred before the formation of this tri-species hybrid zone, and might result from bridge hybridization, vagrancies, or other three-species interactions.

Population Genetics and Speciation of Yellow-Bellied, Red-Naped, and Red-Breasted Sapsuckers (Sphyrapicus varius, S. nuchalis, and S. ruber).

The root of understanding speciation lies in determining the forces which drive it. In many closely-related species, including Sphyrapicus varius, S. nuchalis, and S. ruber, it is assumed that speciation occurred due to isolation in multiple Pleistocene refugia. We used genetic data from 457 samples at the control region (CR), cytochrome oxidase I (COI), and chromo-helicase DNA binding protein (CHD1Z) to examine range-wide population genetic structure and differentiation amongst these 3 species across each species' breeding range. In addition, we modeled these species' ecological niches for the Holocene (~6000 years ago), Last Glacial Maximum (~22000 years ago), and Last Interglacial (~120000-140000 years ago) to determine if Pleistocene glaciations could have contributed to allopatric distributions, therefore allowing these groups to differentiate. Population genetic data show a potential Pleistocene refugium in Haida Gwaii, an east-west split among S. varius, and low genetic differentiation within each species. Our CR data show some polyphyly, while COI and CHD1Z data show differentiation among species using composite genotypes. Ecological Niche Modeling shows a large amount of niche overlap at each time period suggesting that S. varius, S. nuchalis, and S. ruber may not have been completely allopatric, and these species likely had repeated intermittent contact. Our data support the growing body of research that suggests differentiation despite gene flow.

Population genetic isolation and limited connectivity in the purple finch (Haemorhous purpureus).

Using a combination of mitochondrial and z-linked sequences, microsatellite data, and spatio-geographic modeling, we examined historical and contemporary factors influencing the population genetic structure of the purple finch (Haemorhous purpureus). Mitochondrial DNA data show the presence of two distinct groups corresponding to the two subspecies, H. p. purpureus and H. p. californicus. The two subspecies likely survived in separate refugia during the last glacial maximum, one on the Pacific Coast and one east of the Rocky Mountains, and now remain distinct lineages with little evidence of gene flow between them. Southwestern British Columbia is a notable exception, as subspecies mixing between central British Columbia and Vancouver Island populations suggests a possible contact zone in this region. Z-linked data support two mitochondrial groups; however, Coastal Oregon and central British Columbia sites show evidence of mixing. Contemporary population structure based on microsatellite data identified at least six genetic clusters: three H. p. purpureus clusters, two H. p. californicus clusters, and one mixed cluster, which likely resulted from high site fidelity and isolation by distance, combined with sexual selection on morphological characters reinforcing subspecies differences.