A gene's eye view of epistasis,selection and speciation

In this mini‐review, I discuss the effects of gene interaction or epistasis from a `gene's eye view.' By a `gene's eye view' of epistasis, I mean that I will consider a single, bi‐allelic locus, A , whose effects on fitness result only from its interactions with alleles of another, unknown locus, X . I will show how changes in the frequencies of alleles at the background locus affect the relationship of alleles at the A ‐locus to fitness. Changing the genetic background changes the fundamental characteristics of the A ‐locus, such as the magnitude and sign of allelic effects on fitness, and, consequently, it changes the strength and pattern of selection. I consider each of the four kinds of pair–wise interactions between two loci and show that some kinds of epistasis are more sensitive than others to population genetic subdivision. Lastly, I show that some kinds of epistasis are more likely than others to affect the process of speciation and contribute to or be responsible for general genetic features of interspecific hybrids, such as Haldane's rule.     

Detecting the genomic signal of polygenic adaptation and the role of epistasis in evolution

Over the last decade, the genomic revolution has offered the possibility to generate tremendous amounts of data that contain valuable information on the genetic basis of phenotypic traits, such as those linked to human diseases or those that allow for species to adapt to a changing environment. Most ecologically relevant traits are controlled by a large number of genes with small individual effects on trait variation, but that are connected with one another through complex developmental, metabolic and biochemical networks. As a result, it has recently been suggested that most adaptation events in natural populations are reached via correlated changes at multiple genes at a time, for which the name polygenic adaptation has been coined. The current challenge is to develop methods to extract the relevant information from genomic data to detect the signature of polygenic evolutionary change. The symposium entitled “Detecting the Genomic Signal of Polygenic Adaptation and the Role of Epistasis in Evolution” held in 2017 at the University of Zürich aimed at reviewing our current state of knowledge. In this review, we use the talks of the invited speakers to summarize some of the most recent developments in this field.     

Effects of genetic drift on variance components under a general model of epistasis

We analyze the changes in the mean and variance components of a quantitative trait caused by changes in allele frequencies, concentrating on the effects of genetic drift. We use a general representation of epistasis and dominance that allows an arbitrary relation between genotype and phenotype for any number of diallelic loci. We assume initial and final Hardy-Weinberg and linkage equilibrium in our analyses of drift-induced changes. Random drift generates transient linkage disequilibria that cause correlations between allele frequency fluctuations at different loci. However, we show that these have negligible effects, at least for interactions among small numbers of loci. Our analyses are based on diffusion approximations that summarize the effects of drift in terms of F, the inbreeding coefficient, interpreted as the expected proportional decrease in heterozygosity at each locus. For haploids, the variance of the trait mean after a population bottleneck is var(delta(z)) = sigma(n)k=1 FkV(A(k)), where n is the number of loci contributing to the trait variance, V(A(1)) = V(A) is the additive genetic variance, and V(A(k)) is the kth-order additive epistatic variance. The expected additive genetic variance after the bottleneck, denoted (V*(A)), is closely related to var(delta(z)); (V*(A)) = (1 - F) sigma(n)k=1 kFk-1V(A(k)). Thus, epistasis inflates the expected additive variance above V(A)(1 - F), the expectation under additivity. For haploids (and diploids without dominance), the expected value of every variance component is inflated by the existence of higher order interactions (e.g., third-order epistasis inflates (V*(AA. This is not true in general with diploidy, because dominance alone can reduce (V*(A)) below V(A)(1 - F) (e.g., when dominant alleles are rare). Without dominance, diploidy produces simple expressions: var(delta(z)) = sigma(n)k=1 (2F)kV(A(k)) and (V(A)) = (1 - F) sigma(n)k=1 k(2F)k-1V(A(k)). With dominance (and even without epistasis), var(delta(z)) and (V*(A)) no longer depend solely on the variance components in the base population. For small F, the expected additive variance simplifies to (V*(A)) approximately equal to (1 - F)V(A) + 4FV(AA) + 2FV(D) + 2FC(AD), where C(AD) is a sum of two terms describing covariances between additive effects and dominance and additive X dominance interactions. Whether population bottlenecks lead to expected increases in additive variance depends primarily on the ratio of nonadditive to additive genetic variance in the base population, but dominance precludes simple predictions based solely on variance components. We illustrate these results using a model in which genotypic values are drawn at random, allowing extreme and erratic epistatic interactions. Although our analyses clarify the conditions under which drift is expected to increase V(A), we question the evolutionary importance of such increases.     

Nematode–bacteria mutualism: Selection within the mutualism supersedes selection outside of the mutualism

The coevolution of interacting species can lead to codependent mutualists. Little is known about the effect of selection on partners within verses apart from the association. Here, we determined the effect of selection on bacteria (Xenorhabdus nematophila) both within and apart from its mutualistic partner (a nematode, Steinernema carpocapsae). In nature, the two species cooperatively infect and kill arthropods. We passaged the bacteria either together with (M+), or isolated from (M?), nematodes under two different selection regimes: random selection (S?) and selection for increased virulence against arthropod hosts (S+). We found that the isolated bacteria evolved greater virulence under selection for greater virulence (M?S+) than under random selection (M?S?). In addition, the response to selection in the isolated bacteria (M?S+) caused a breakdown of the mutualism following reintroduction to the nematode. Finally, selection for greater virulence did not alter the evolutionary trajectories of bacteria passaged within the mutualism (M+S+ = M+S?), indicating that selection for the maintenance of the mutualism was stronger than selection for increased virulence. The results show that selection on isolated mutualists can rapidly breakdown beneficial interactions between species, but that selection within a mutualism can supersede external selection, potentially generating codependence over time.     

Integrating the ‘genomic mosaic’ view of species into studies of biotic interactions: a comment on Bernhardsson et al. ()

The co‐occurrence of geographical structure in herbivore communities, metabolomes and defence genes in forest trees has been analysed in the context of ‘geographical mosaics’ of coevolution. A deeper understanding of these important issues will require full integration of a ‘genomic mosaic’ view of species into community ecology.     

Signature of sexual conflict is actually conflict resolved

There has been substantial interest of late in using population genetic methods to study sexual conflict, where an allele increases the fitness of one sex at some cost to the other (Mank, 2017). Population genomic scans for sexual conflict offer an important advance given the difficulties of identifying antagonistic alleles from more traditional methods, and could greatly increase our understanding of the extent and loci of sexual conflict. This is particularly true for studies in natural populations, for which obtaining accurate fitness measurements for each sex can be challenging. In this issue of Molecular Ecology, Bissegger, Laurentino, Roesti, and Berner (2019) present a cautionary tale about how to interpret these population genomic data.     

THE INTENSITY OF SELECTION ACTING ON THE COUCH POTATO GENE—SPATIAL–TEMPORAL VARIATION IN A DIAPAUSE CLINE

Cosmopolitan populations of Drosophila melanogaster have co‐opted a form of reproductive diapause to overwinter in northern populations. Polymorphism in the couch potato gene has been implicated in genetic variation for this diapause trait. Using a collection of 20 populations from Florida to Canada and 11 collections from 3 years in a Pennsylvania orchard, we estimated the allele frequencies for 15 single nucleotide polymorphisms (SNPs) in the couch potato gene. These include the specific polymorphism associated with diapause inducability. We find that the SNP polymorphism, 48034(A/T), is correlated with latitude and its frequencies are predicted by the incidence of diapause trait. We find that the clinal patterns for cpo SNPs sampled in 1997 are similar to the same SNPs sampled in 2009–2010. SNPs that show apparent associations with cpo expression are also clinal with the low‐expression allele increasing in frequency, as would be predicted from functional knockout studies of cpo. Finally, we see a significant pattern where the frequency of the diapause‐causing allele drops in frequency during the summer season, consistent with the drop in the incidence of the diapause trait. The selection required to drive this response is large, roughly 24% to 59% per generation depending on the degree of dominance.     

Predicting epistasis: an experimental test of metabolic control theory with bacterial transcription and translation

Epistatic interactions between mutations are thought to play a crucial role in a number of evolutionary processes, including adaptation and sex. Evidence for epistasis is abundant, but tests of general theoretical models that can predict epistasis are lacking. In this study, I test the ability of metabolic control theory to predict epistasis using a novel experimental approach that combines phenotypic and genetic perturbations of enzymes involved in gene expression and protein synthesis in the bacterium Pseudomonas aeruginosa. These experiments provide experimental support for two key predictions of metabolic control theory: (i) epistasis between genes involved in the same pathway is antagonistic; (ii) epistasis becomes increasingly antagonistic as mutational severity increases. Metabolic control theory is a general theory that applies to any set of genes that are involved in the same linear processing chain, not just metabolic pathways, and I argue that this theory is likely to have important implications for predicting epistasis between functionally coupled genes, such as those involved in antibiotic resistance. Finally, this study highlights the fact that phenotypic manipulations of gene activity provide a powerful method for studying epistasis that complements existing genetic methods.     

How can distinct egg polymorphism be maintained in the rufescent prinia (Prinia rufescens)–plaintive cuckoo (Cacomantis merulinus) interaction—a modeling approach

In avian brood parasitism, both the host and the parasite are expected to develop various conflicting adaptations; hosts develop a defense against parasitism, such as an ability to recognize and reject parasitic eggs that look unlike their own, while parasites evolve egg mimicry to counter this host defense. Hosts may further evolve to generate various egg phenotypes that are not mimicked by parasites. Difference in egg phenotype critically affects the successful reproduction of hosts and parasites. Recent studies have shown that clear polymorphism in egg phenotype is observed in several host–parasite interactions, which suggests that egg polymorphism may be a more universal phenomenon than previously thought. We examined the mechanism for maintaining egg polymorphism in the rufescent prinia (Prinia rufescens) that is parasitized by the plaintive cuckoo (Cacomantis merulinus) from a theoretical viewpoint based on a mathematical model. The prinia has four distinct egg phenotypes: immaculate white, immaculate blue, white with spots, and blue with spots. Only two egg phenotypes, white with spots and blue with spots, are found in the cuckoo population. We show that the observed prinia and cuckoo phenotypes cannot be at an equilibrium and that egg polymorphism can be maintained either at stationary equilibrium or with dynamic, frequency oscillations, depending on the mutation rates of the background color and spottiness. Long‐term monitoring of the prinia–cuckoo interaction over a wide geographic range is needed to test the results of the model analyses.     

Environmental robustness and the adaptability of populations

Recent work has shown that genetic robustness can either facilitate or impede adaptation. But the impact of environmental robustness on adaptation remains unclear. Environmental robustness helps ensure that organisms consistently develop the same phenotype in the face of "environmental noise" during development. Under purifying selection, those genotypes that express the optimal phenotype most reliably will be selectively favored. The resulting reduction in genetic variation tends to slow adaptation when the population is faced with a novel target phenotype. However, environmental noise sometimes induces the expression of an alternative advantageous phenotype, which may speed adaptation by genetic assimilation. Here, we use a population-genetic model to explore how these two opposing effects of environmental noise influence the capacity of a population to adapt. We analyze how the rate of adaptation depends on the frequency of environmental noise, the degree of environmental robustness in the population, the distribution of environmental robustness across genotypes, the population size, and the strength of selection for a newly adaptive phenotype. Over a broad regime, we find that environmental noise can either facilitate or impede adaptation. Our analysis uncovers several surprising insights about the relationship between environmental noise and adaptation, and it provides a general framework for interpreting empirical studies of both genetic and environmental robustness.     

The evolutionary interplay of intergroup conflict and altruism in humans: a review of parochial altruism theory and prospects for its extension

Drawing on an idea proposed by Darwin, it has recently been hypothesized that violent intergroup conflict might have played a substantial role in the evolution of human cooperativeness and altruism. The central notion of this argument, dubbed ‘parochial altruism’, is that the two genetic or cultural traits, aggressiveness against the out-groups and cooperativeness towards the in-group, including self-sacrificial altruistic behaviour, might have coevolved in humans. This review assesses the explanatory power of current theories of ‘parochial altruism’. After a brief synopsis of the existing literature, two pitfalls in the interpretation of the most widely used models are discussed: potential direct benefits and high relatedness between group members implicitly induced by assumptions about conflict structure and frequency. Then, a number of simplifying assumptions made in the construction of these models are pointed out which currently limit their explanatory power. Next, relevant empirical evidence from several disciplines which could guide future theoretical extensions is reviewed. Finally, selected alternative accounts of evolutionary links between intergroup conflict and intragroup cooperation are briefly discussed which could be integrated with parochial altruism in the future.     

INTROGRESSION IN THE DROSOPHILA SUBOBSCURA—D. MADEIRENSIS SISTER SPECIES: EVIDENCE OF GENE FLOW IN NUCLEAR GENES DESPITE MITOCHONDRIAL DIFFERENTIATION

Species hybridization, and thus the potential for gene flow, was once viewed as reproductive mistake. However, recent analysis based on large datasets and newly developed models suggest that gene exchange is not as rare as originally suspected. To investigate the history and speciation of the closely related species Drosophila subobscura, D. madeirensis, and D. guanche, we obtained polymorphism and divergence data for 26 regions throughout the genome, including the Y chromosome and mitochondrial DNA. We found that the D. subobscura X/autosome ratio of silent nucleotide diversity is significantly smaller than the 0.75 expected under neutrality. This pattern, if held genomewide, may reflect a faster accumulation of beneficial mutations on the X chromosome than on autosomes. We also detected evidence of gene flow in autosomal regions, while sex chromosomes remain distinct. This is consistent with the large X effect on hybrid male sterility seen in this system and the presence of two X chromosome inversions fixed between species. Overall, our data conform to chromosomal speciation models in which rearrangements are proposed to serve as gene flow barriers. Contrary to other observations in Drosophila, the mitochondrial genome appears resilient to gene flow in the presence of nuclear exchange.     

Sexual antagonism leads to a mosaic of X-autosome conflict

Males and females have different optimal values for some traits, such as body size. When the same genes control these traits in both sexes, selection pushes in opposite directions in males and females. Alleles at autosomal loci spend equal amounts of time in males and females, suggesting that the sexually antagonistic selective forces may approximately balance between the opposing optima. Frank and Crespi noted that alleles on the X chromosome spend twice as much time in diploid females as in haploid males. That distinction between the sexes may tend to favor X-linked genes that push more strongly toward the female optimum than the male optimum. The female bias of X-linked genes opposes the intermediate optimum of autosomal genes, potentially creating a difference between the direction of selection on traits favored by X chromosomes and autosomes. Patten has recently argued that explicit genetic assumptions about dominance and the relative magnitude of allelic effects may lead X-linked genes to favor the male rather than the female optimum, contradicting Frank and Crespi. This article combines the insights of those prior analyses into a new, more general theory. We find some parameter combinations for X-linked loci that favor a female bias and other parameter combinations that favor a male bias. We conclude that the X likely contains a mosaic pattern of loci that differ with autosomes over sexually antagonistic traits. The overall tendency for a female or male bias on the X depends on prior assumptions about the distribution of key parameters across X-linked loci. Those parameters include the dominance coefficient and the way in which ploidy influences the magnitude of allelic effects.     

HETEROGENEOUS GENOMIC DIFFERENTIATION IN MARINE THREESPINE STICKLEBACKS: ADAPTATION ALONG AN ENVIRONMENTAL GRADIENT

Evolutionary divergence among populations occupying ecologically distinct environments can occur even in the face of on‐going gene flow. However, the genetic underpinnings, as well as the scale and magnitude at which this differentiation occurs in marine habitats are not well understood. We investigated the patterns and degree of genomic heterogeneity in threespine sticklebacks (Gasterosteus aculeatus) by assessing genetic variability in 20 nongenic and 20 genic (associated with genes important for freshwater adaptation) microsatellite loci in samples collected from 38 locations spanning the entire Baltic Sea coast to the North Sea boundary. Population divergence (F_ST ≈ 0.026) and structuring (five genetic clusters) was significantly more pronounced in the genic as compared to nongenic markers (F_ST ≈ 0.008; no genetic clusters). Patterns of divergence in the genic markers—45% of which were identified as outliers—correlated with local differences in salinity. Yet, a strong positive correlation between divergence in genic and nongenic markers, and their association with environmental factors suggests that adaptive divergence is reducing gene flow across the genome. Apart from providing a clear demonstration of heterogeneous genomic patterns of differentiation in a marine species, the results are indicative of adaptive population structuring across the relatively young Baltic Sea in spite of ample opportunities for gene flow.     

Extensive gene flow along the urban–rural gradient in a migratory colonial bird

Urban colonization by wildlife involves a combination of several different mechanisms, including phenotype or genotype sorting, phenotypic plasticity and microevolutionary adaptation. Combination of these processes can produce a rapid phenotypic, but also genetic divergence of urban versus rural populations. Here, we examined the pattern of genetic differentiation between urban and rural populations of a colonial migratory bird, the black‐headed gull Chroicocephalus ridibundus. To this end, we sampled ca 170 individuals from six (two urban and four rural) colonies in northern Poland, and genotyped them at ten microsatellite loci. Our analysis provided evidence for negligible genetic divergence of urban and rural colonies, as assessed with fixation index F_ST and Nei's unbiased genetic distance D (mean pairwise urban‐rural comparisons: F_ST = 0.003 ± 0.001 [SE] and D = 0.012 ± 0.006 [SE]). Bayesian clustering methods provided support for homogeneous genetic structure across all urban and rural populations. Also, we found no support for reduced allelic diversity in urban versus rural colonies. These results stand in a stark contrast to the previous findings on the genetic consequences of urbanization in birds. We hypothesize that this pattern could possibly be attributed to the important life‐history characters of the black‐headed gull, including coloniality, migratoriness, and high dispersal propensity. Our study provides a novel insight into the urban landscape genetics, underlining large variation in the mechanisms of urban colonization and its genetic consequences in wild animal populations.     

Domestication and fitness in the wild: A multivariate view

Domesticated species continually escaping and interbreeding with wild relatives impose a migration load on wild populations. As domesticated stocks become increasingly different as a result of artificial and natural selection in captivity, fitness of escapees in the wild is expected to decline, reducing the effective rate of migration into wild populations. Recent theory suggest that this may alleviate and eventually eliminate the resulting migration load. I develop a multivariate model of trait and wild fitness evolution resulting from the joint effects of artificial and natural selection in the captive environment. Initially, the evolutionary trajectory is dominated by the effects of artificial selection causing a fast initial decline in fitness of escapees in the wild. In later phases, through the counteracting effects of correlational multivariate natural selection in captivity, the mean phenotype is pushed in directions of weak stabilizing selection, allowing a sustained response in the trait subject to artificial selection. Provided that there is some alignment between the adaptive landscapes in the wild and in captivity, these phases are associated with slower rates of decline in wild fitness of the domesticated stock, suggesting that detrimental effects on wild populations are likely to remain a concern in the foreseeable future.