Coadapted genomes and selection on hybrids: Fisher's geometric model explains a variety of empirical patterns

Natural selection plays a variety of roles in hybridization, speciation, and admixture. Most research has focused on two extreme cases: crosses between closely related inbred lines, where hybrids are fitter than their parents, or crosses between effectively isolated species, where hybrids suffer severe breakdown. But many natural populations must fall into intermediate regimes, with multiple types of gene interaction, and these are more difficult to study. Here, we develop a simple fitness landscape model, and show that it naturally interpolates between previous modeling approaches, which were designed for the extreme cases, and invoke either mildly deleterious recessives, or discrete hybrid incompatibilities. Our model yields several new predictions, which we test with genomic data from Mytilus mussels, and published data from plants (Zea, Populus, and Senecio) and animals (Mus, Teleogryllus, and Drosophila). The predictions are generally supported, and the model explains a number of surprising empirical patterns. Our approach enables novel and complementary uses of genome-wide datasets, which do not depend on identifying outlier loci, or "speciation genes" with anomalous effects. Given its simplicity and flexibility, and its predictive successes with a wide range of data, the approach should be readily extendable to other outstanding questions in the study of hybridization.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Heterogeneous Patterns of Genetic Diversity and Differentiation in European and Siberian Chiffchaff (Phylloscopus collybita abietinus/P. tristis)

Identification of candidate genes for trait variation in diverging lineages and characterization of mechanistic underpinnings of genome differentiation are key steps toward understanding the processes underlying the formation of new species. Hybrid zones provide a valuable resource for such investigations, since they allow us to study how genomes evolve as species exchange genetic material and to associate particular genetic regions with phenotypic traits of interest. Here, we use whole-genome resequencing of both allopatric and hybridizing populations of the European (Phylloscopus collybita abietinus) and the Siberian chiffchaff (P. tristis)—two recently diverged species which differ in morphology, plumage, song, habitat, and migration—to quantify the regional variation in genome-wide genetic diversity and differentiation, and to identify candidate regions for trait variation. We find that the levels of diversity, differentiation, and divergence are highly heterogeneous, with significantly reduced global differentiation, and more pronounced differentiation peaks in sympatry than in allopatry. This pattern is consistent with regional differences in effective population size and recurrent background selection or selective sweeps reducing the genetic diversity in specific regions prior to lineage divergence, but the data also suggest that postdivergence selection has resulted in increased differentiation and fixed differences in specific regions. We find that hybridization and backcrossing is common in sympatry, and that phenotype is a poor predictor of the genomic composition of sympatric birds. The combination of a differentiation scan approach with identification of fixed differences pinpoint a handful of candidate regions that might be important for trait variation between the two species.