To systematically explore complex genetic interactions, we constructed ~200,000 yeast triple mutants and scored negative trigenic interactions. We selected double mutant query genes across a broad spectrum of biological processes, spanning a range of quantitative features of the global digenic interaction network and tested for a genetic interaction with a third mutation. Trigenic interactions often occurred among functionally related genes and essential genes were hubs on the trigenic network. Despite their functional enrichment, trigenic interactions tended to link genes in distant bioprocesses and display a weaker magnitude than digenic interactions. We estimate that the global trigenic interaction network is ~100-fold larger than the global digenic network, highlighting the potential for complex genetic interactions to impact the biology of inheritance, including the genotype to phenotype relationship.

Recombination often differs markedly between males and females. Here we present the first analysis of sex-specific recombination in Gasterosteus sticklebacks. Using whole-genome sequencing of 15 crosses between G. aculeatus and G. nipponicus, we localized 698 crossovers with a median resolution of 2.3 kb. We also used a bioinformatic approach to infer historical sex-averaged recombination patterns for both species. Recombination is greater in females than males on all chromosomes, and overall map length is 1.64 times longer in females. The locations of crossovers differ strikingly between sexes. Crossovers cluster toward chromosome ends in males, but are distributed more evenly across chromosomes in females. Suppression of recombination near the centromeres in males causes crossovers to cluster at the ends of long arms in acrocentric chromosomes, and greatly reduces crossing over on short arms. The effect of centromeres on recombination is much weaker in females. Genomic differentiation between G. aculeatus and G. nipponicus is strongly correlated with recombination rate, and patterns of differentiation along chromosomes are strongly influenced by male-specific telomere and centromere effects. We found no evidence for fine-scale correlations between recombination and local gene content in either sex. We discuss hypotheses for the origin of sexual dimorphism in recombination and its consequences for sexually antagonistic selection and sex chromosome evolution.

The evolution of postzygotic reproductive isolation is an important component of speciation. But before isolation is complete there is sometimes a phase of heterosis in which hybrid fitness exceeds that of the two parental species. The genetics and evolution of heterosis and postzygotic isolation have typically been studied in isolation, precluding the development of a unified theory of speciation. Here, we develop a model that incorporates both positive and negative gene interactions, and accounts for the evolution of both heterosis and postzygotic isolation. We parameterize the model with recent data on the fitness effects of 10,000 mutations in yeast, singly and in pairwise epistatic combinations. The model makes novel predictions about the types of interactions that contribute to declining hybrid fitness. We reproduce patterns familiar from earlier models of speciation (e.g. Haldane’s Rule and Darwin’s Corollary) and identify new mechanisms that may underlie these patterns. Our approach provides a general framework for integrating experimental data from gene interaction networks into speciation theory and makes new predictions about the genetic mechanisms of speciation.

Speciation is characterized by the development of reproductive isolating barriers between diverging groups. Intrinsic post-zygotic barriers of the type envisioned by Bateson, Dobzhansky, and Muller are deleterious epistatic interactions among loci that reduce hybrid fitness, leading to reproductive isolation. The first formal population genetic model of the development of these barriers was published by Orr in 1995, and here we develop a more general model of this process by incorporating finite protein–protein interaction networks, which reduce the probability of deleterious interactions in vivo. Our model shows that the development of deleterious interactions is limited by the density of the protein–protein interaction network. We have confirmed our analytical predictions of the number of possible interactions given the number of allele substitutions by using simulations on the Saccharomyces cerevisiae protein–protein interaction network. These results allow us to define the rate at which deleterious interactions are expected to form, and hence the speciation rate, for any protein–protein interaction network.

Chromosomal inversions are frequently implicated in isolating species. Models have shown how inversions can evolve in the context of postmating isolation. Inversions are also frequently associated with mating preferences, a topic that has not been studied theoretically. Here, we show how inversions can spread by capturing a mating preference locus and one or more loci involved with epistatic incompatibilities. Inversions can be established under broad conditions ranging from near panmixis to nearly complete speciation. These results provide a hypothesis to explain the growing number of examples of inversions associated with premating isolating mechanisms.

Understanding the factors that produce and maintain genetic variation is a central goal of evolutionary biology. Despite a century of genetic analysis, the evolutionary history underlying patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic and phenotypic variation is partitioned across global D. melanogaster populations, and specifically in its putative ancestral range in Subtropical Africa, remains unresolved. Here, we integrate genomic and behavioral analyses to assess patterns of population genetic structure, admixture, mate preference, and genetic incompatibility throughout the range of this model organism. Our analysis includes 174 new accessions from novel and under-sampled regions within Subtropical Africa. We find that while almost all Out of Africa genomes correspond to a single genetic ancestry, different geographic regions within Africa contain multiple distinct ancestries, including the presence of substantial cryptic diversity within Subtropical Africa. We find evidence for significant admixture- and variation in admixture rates-between geographic regions within Africa, as well as between African and non-African lineages. By combining behavioral analysis with population genomics, we demonstrate that female mate choice is highly polymorphic, behavioral types are not monophyletic, and that genomic differences between behavioral types correspond to many regions across the genome. These include regions associated with neurological development, behavior, olfactory perception, and learning. Finally, we discovered that many individual pairs of putative incompatibility loci likely evolved during or after the expansion of D. melanogaster out of Africa. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the distribution of reproductive barriers that are polymorphic within species.

With the rise of affordable next-generation sequencing technology, introgression-or the exchange of genetic materials between taxa-has become widely perceived to be a ubiquitous phenomenon in nature. Although this claim is supported by several keystone studies, no thorough assessment of the frequency of introgression across eukaryotes in nature has been performed to date. In this manuscript, we aim to address this knowledge gap by examining patterns of introgression across eukaryotes. We collated a single statistic, Patterson's D, which can be used as a test for introgression across 123 studies to further assess how taxonomic group, divergence time, and sequencing technology influence reports of introgression. Overall, introgression has mostly been measured in plants and vertebrates, with less attention given to the rest of the Eukaryotes. We find that the most frequently used metrics to detect introgression are difficult to compare across studies and even more so across biological systems due to differences in study effort, reporting standards, and methodology. Nonetheless, our analyses reveal several intriguing patterns, including the observation that differences in sequencing technologies may bias values of Patterson's D and that introgression may differ throughout the course of the speciation process. Together, these results suggest the need for a unified approach to quantifying introgression in natural communities and highlight important areas of future research that can be better assessed once this unified approach is met.

Intralocus sexually antagonistic selection occurs when an allele is beneficial to one sex but detrimental to the other. This form of selection is thought to be key to the evolution of sex chromosomes but is hard to detect. Here we perform an analysis of phased young sex chromosomes to look for signals of sexually antagonistic selection in the Japan Sea stickleback ( Gasterosteus nipponicus ). Phasing allows us to date the suppression of recombination on the sex chromosome and provides unprecedented resolution to identify sexually antagonistic selection in the recombining region of the chromosome. We identify four windows with elevated divergence between the X and Y in the recombining region, all in or very near genes associated with phenotypes potentially under sexually antagonistic selection in humans. We are unable, however, to rule out the alternative hypothesis that the peaks of divergence result from demographic effects. Thus, although sexually antagonistic selection is a key hypothesis for the formation of supergenes on sex chromosomes, it remains challenging to detect. This article is part of the theme issue ‘Genomic architecture of supergenes: causes and evolutionary consequences’.