Parallel emergence of negative epistasis across experimental lineages

Epistatic interactions can greatly impact evolutionary phenomena, particularly the process of adaptation. Here, we leverage four parallel experimentally evolved lineages to study the emergence and trajectories of epistatic interactions in the social bacterium Myxococcus xanthus. A social gene (pilA) necessary for effective group swarming on soft agar had been deleted from the common ancestor of these lineages. During selection for competitiveness at the leading edge of growing colonies, two lineages evolved qualitatively novel mechanisms for greatly increased swarming on soft agar, whereas the other two lineages evolved relatively small increases in swarming. By reintroducing pilA into different genetic backgrounds along the four lineages, we tested whether parallel lineages showed similar patterns of epistasis. In particular, we tested whether a pattern of negative epistasis between accumulating mutations and pilA previously found in the fastest lineage would be found only in the two evolved lineages with the fastest and most striking swarming phenotypes, or rather was due to common epistatic structure across all lineages arising from the generic fixation of adaptive mutations. Our analysis reveals the emergence of negative epistasis across all four independent lineages. Further, we present results showing that the observed negative epistasis is not due exclusively to evolving populations approaching a maximum phenotypic value that inherently limits positive effects of pilA reintroduction, but rather involves direct antagonistic interactions between accumulating mutations and the reintroduced social gene.     

Experimental social evolution with Myxococcus xanthus

Genetically-based social behaviors are subject to evolutionary change in response to natural selection. Numerous microbial systems provide not only the opportunity to understand the genetic mechanisms underlying specific social interactions, but also to observe evolutionary changes in sociality over short time periods. Here we summarize experiments in which behaviors of the social bacterium Myxococcus xanthus changed extensively during evolutionary adaptation to two relatively asocial laboratory environments. M. xanthus moves cooperatively, exhibits cooperative multicellular development upon starvation and also appears to prey cooperatively on other bacteria. Replicate populations of M. xanthus were evolved in both structured (agar plate) and unstructured (liquid) environments that contained abundant resources. The importance of social cooperation for evolutionary fitness in these habitats was limited by the absence of positive selection for starvation-induced spore production or predatory efficiency. Evolved populations showed major losses in all measured categories of social proficiency- motility, predation, fruiting ability, and sporulation. Moreover, several evolved genotypes were observed to exploit the social behavior of their ancestral parent when mixed together during the developmental process. These experiments that resulted in both socially defective and socially exploitative genotypes demonstrate the power of laboratory selection experiments for studying social evolution at the microbial level. Results from additional selection experiments that place positive selection pressure on social phenotypes can be integrated with direct study of natural populations to increase our understanding of principles that underlie the evolution of microbial social behavior. This revised version was published online in August 2006 with corrections to the Cover Date.     

A SHIFT FROM MAGNITUDE TO SIGN EPISTASIS DURING ADAPTIVE EVOLUTION OF A BACTERIAL SOCIAL TRAIT

Although the importance of epistasis in evolution has long been recognized, remarkably little is known about the processes by which epistatic interactions evolve in real time in specific biological systems. Here, we have characterized how the epistatic fitness relationship between a social gene and an adapting genome changes radically over a short evolutionary time frame in the social bacterium Myxococcus xanthus. We show that a highly beneficial effect of this social gene in the ancestral genome is gradually reduced—and ultimately reversed into a deleterious effect—over the course of an experimental adaptive trajectory in which a primitive form of novel cooperation evolved. This reduction and reversal of a positive social allelic effect is driven solely by changes in the genetic context in which the gene is expressed as new mutations are sequentially fixed during adaptive evolution, and explicitly demonstrates a significant evolutionary change in the genetic architecture of an ecologically important social trait.     

Resource Level Affects Relative Performance of the Two Motility Systems of Myxococcus xanthus

The adventurous (A) and social (S) motility systems of the microbial predator Myxococcus xanthus show differential swarming performance on distinct surface types. Under standard laboratory conditions, A-motility performs well on hard agar but poorly on soft agar, whereas the inverse pattern is shown by S-motility. These properties may allow M. xanthus to swarm effectively across a greater diversity of natural surfaces than would be possible with one motility system alone. Nonetheless, the range of ecological conditions under which dual motility enhances effective swarming across distinct surfaces and how ecological parameters affect the complementarity of A-motility and S-motility remain unclear. Here we have examined the role of nutrient concentration in determining swarming patterns driven by dual motility on distinct agar surfaces, as well as the relative contributions of A-motility and S-motility to these patterns. Swarm expansion rates of dually motile (A⁺S⁺), solely A-motile (A⁺S⁻), and solely S-motile (A⁻S⁺) strains were compared on hard and soft agar across a wide range of casitone concentrations. At low casitone concentrations (0–0.1%), swarming on soft agar driven by S-motility is very poor, and is significantly slower than swarming on hard agar driven by A-motility. This reverses at high casitone concentration (1–3.2%) such that swarming on soft agar is much faster than swarming on hard agar. This pattern greatly constrained the ability of M. xanthus to encounter patches of prey bacteria on a soft agar surface when nutrient levels between the patches were low. The swarming patterns of a strain that is unable to produce extracellular fibrils indicate that these appendages are responsible for the elevated swarming of S-motility at high resource levels. Together, these data suggest that large contributions by S-motility to predatory swarming in natural soils may be limited to soft, wet, high-nutrient conditions that may be uncommon. Several likely benefits of S-motility to the M. xanthus life cycle are discussed, including synergistic interactions with A-motility across a wide variety of conditions.     

Epistatic pleiotropy and the genetic architecture of covariation within early and late-developing skull trait complexes in mice

The role of epistasis as a source of trait variation is well established, but its role as a source of covariation among traits (i.e., as a source of "epistatic pleiotropy") is rarely considered. In this study we examine the relative importance of epistatic pleiotropy in producing covariation within early and late-developing skull trait complexes in a population of mice derived from an intercross of the Large and Small inbred strains. Significant epistasis was found for several pairwise combinations of the 21 quantitative trait loci (QTL) affecting early developing traits and among the 20 QTL affecting late-developing traits. The majority of the epistatic effects were restricted to single traits but epistatic pleiotropy still contributed significantly to covariances. Because of their proportionally larger effects on variances than on covariances, epistatic effects tended to reduce within-group correlations of traits and reduce their overall degree of integration. The expected contributions of single-locus and two-locus epistatic pleiotropic QTL effects to the genetic covariance between traits were analyzed using a two-locus population genetic model. The model demonstrates that, for single-locus or epistatic pleiotropy to contribute to trait covariances in the study population, both traits must show the same pattern of single-locus or epistatic effects. As a result, a large number of the cases where loci show pleiotropic effects do not contribute to the covariance between traits in this population because the loci show a different pattern of effect on the different traits. In general, covariance patterns produced by single-locus and epistatic pleiotropy predicted by the model agreed well with actual values calculated from the QTL analysis. Nearly all single-locus and epistatic pleiotropic effects contributed positive components to covariances between traits, suggesting that genetic integration in the skull is achieved by a complex combination of pleiotropic effects.     

Genetic variation in pleiotropy: differential epistasis as a source of variation in the allometric relationship between long bone lengths and body weight

Pleiotropy is an aspect of genetic architecture underlying the phenotypic covariance structure. The presence of genetic variation in pleiotropy is necessary for natural selection to shape patterns of covariation between traits. We examined the contribution of differential epistasis to variation in the intertrait relationship and the nature of this variation. Genetic variation in pleiotropy was revealed by mapping quantitative trait loci (QTLs) affecting the allometry of mouse limb and tail length relative to body weight in the mouse-inbred strain LG/J by SM/J intercross. These relationship QTLs (rQTLs) modify relationships between the traits affected by a common pleiotropic locus. We detected 11 rQTLs, mostly affecting allometry of multiple bones. We further identified epistatic interactions responsible for the observed allometric variation. Forty loci that interact epistatically with the detected rQTLs were identified. We demonstrate how these epistatic interactions differentially affect the body size variance and the covariance of traits with body size. We conclude that epistasis, by differentially affecting both the canalization and mean values of the traits of a pleiotropic domain, causes variation in the covariance structure. Variation in pleiotropy maintains evolvability of the genetic architecture, in particular the evolvability of its modular organization.     

The contribution of epistatic pleiotropy to the genetic architecture of covariation among polygenic traits in mice

The contribution that pleiotropic effects of individual loci make to covariation among traits is well understood theoretically and is becoming well documented empirically. However, little is known about the role of epistasis in determining patterns of covariation among traits. To address this problem we combine a quantitative trait locus (QTL) analysis with a two-locus model to assess the contribution of epistasis to the genetic architecture of variation and covariation of organ weights and limb bone lengths in a backcross population of mice created from the M16i and CAST/Ei strains. Significant epistasis was exhibited by 14 pairwise combinations of QTL for organ weights and 10 combinations of QTL for limb bone lengths, which contributed, on average, about 5% of the variation in organ weights and 8% in limb bone lengths beyond that of single-locus QTL effects. Epistatic pleiotropy was much more common in the limb bones (seven of 10 epistatic combinations affecting limb bone lengths were pleiotropic) than the organs (three of the 14 epistatic combinations affecting organ weights were pleiotropic). In both cases, epistatic pleiotropy was less common than single-locus pleiotropy. Epistatic pleiotropy accounted for an average of 6% of covariation among organ weights and 21% of covariation among limb bone lengths, which represented an average of one-fifth (for organ weights) and one-third (for limb bone lengths) of the total genetic covariance between traits. Thus, although epistatic pleiotropy made a smaller contribution than single-locus pleiotropy, it clearly made a significant contribution to the genetic architecture of variation/covariation.     

Stalk size and altruism investment within and among populations of the social amoeba

Reproductive division of labour is common in many societies, including those of eusocial insects, cooperatively breeding vertebrates, and most forms of multicellularity. However, conflict over what is best for the individual vs. the group can prevent an optimal division of labour from being achieved. In the social amoeba Dictyostelium discoideum, cells aggregate to become multicellular and a fraction behaves altruistically, forming a dead stalk that supports the rest. Theory suggests that intra‐organismal conflict over spore–stalk cell fate can drive rapid evolutionary change in allocation traits, leading to polymorphisms within populations or rapid divergence between them. Here, we assess several proxies for stalk size and spore–stalk allocation as metrics of altruism investment among strains and across geographic regions. We observe geographic divergence in stalk height that can be partly explained by differences in multicellular size, as well as variation among strains in clonal spore–stalk allocation, suggesting within‐population variation in altruism investment. Analyses of chimeras comprised of strains from the same vs. different populations indicated genotype‐by‐genotype epistasis, where the morphology of the chimeras deviated significantly from the average morphology of the strains developed clonally. The significantly negative epistasis observed for allopatric pairings suggests that populations are diverging in their spore–stalk allocation behaviours, generating incompatibilities when they encounter one another. Our results demonstrate divergence in microbial social traits across geographically separated populations and demonstrate how quantification of genotype‐by‐genotype interactions can elucidate the trajectory of social trait evolution in nature.     

Nonadditive genetic components in resistance of the red flour beetle Tribolium castanaeum against parasite infection

Genetically coupled antagonistic coevolution between host and parasites can select for the maintenance of recombination in the host. Mechanistically, maintenance of recombination relies on epistatic interactions between resistance genes creating linkage disequilibria (LD). The role of epistasis in host resistance traits is however only partly understood. Therefore, we applied the joint scaling principle to assess epistasis and other nonadditive genetic components of two resistance traits, survival, and parasite spore load, in population crosses of the red flour beetle Tribolium castanaeum under infections with the microsporidian Nosema whitei. We found nonadditive components only in infected populations but not in control populations. The genetic architecture underlying survival under parasite infection was more complex than that of spore load. Accordingly, the observed negative correlation between survival and spore load was mainly based on a correlation between shared additive components. Breakdown of resistance was especially strong in F2 crosses between resistant lines indicating that multiple epistatic routes can lead to the same adaptation. In general, the wide range of nonoverlapping genetic components between crosses indicated that parasite resistance in T. castanaeum can be understood as a multi peaked fitness landscape with epistasis contributing substantially to phenotypic differentiation in resistance.     

Pleiotropic effects on mandibular morphology II: differential epistasis and genetic variation in morphological integration

The evolution of morphological modularity through the sequestration of pleiotropy to sets of functionally and developmentally related traits requires genetic variation in the relationships between traits. Genetic variation in relationships between traits can result from differential epistasis, where epistatic relationships for pairs of loci are different for different traits. This study maps relationship quantitative trait loci (QTLs), specifically QTLs that affect the relationship between individual mandibular traits and mandible length, across the genome in an F2 intercross of the LG/J and SM/J inbred mouse strains (N = 1045). We discovered 23 relationship QTLs scattered throughout the genome. All mandibular traits were involved in one or more relationship QTL. When multiple traits were affected at a relationship QTL, the traits tended to come from a developmentally restricted region of the mandible, either the muscular processes or the alveolus. About one-third of the relationship QTLs correspond to previously located trait QTLs affecting the same traits. These results comprise examples of genetic variation necessary for an evolutionary response to selection on the range of pleiotropic effects.     

Ecological Variables Affecting Predatory Success in <Emphasis Type="BoldItalic">Myxococcus xanthus</Emphasis>

The feeding efficiency of microbial predators depends on both the availability of various prey species and abiotic variables. Myxococcus xanthus is a bacterial predator that searches for microbial prey by gliding motility, and then kills and lyses its prey with secreted compounds. We manipulated three ecological variables to examine their effects on the predatory performance of M. xanthus to better understand its behavior and how it affects prey populations. Experiments were designed to determine how surface solidity (hard vs soft agar), density of prey patches (1 vs 2 cm grids), and type of prey (Gram-positive Micrococcus luteus vs Gram-negative Escherichia coli) affect predatory swarming and prey killing by M. xanthus. The prey were dispersed in patches on a buffered agar surface. M. xanthus swarms attacked a greater proportion of prey patches when patches were densely arranged on a hard-agar surface, compared with either soft-agar surfaces or low-patch-density arrangements. These ecological variables did not significantly influence the rate of killing of individual prey within a patch, although a few surviving prey were more likely to be recovered on soft agar than on hard agar. These results indicate that M. xanthus quickly kills most nearby E. coli or M. luteus regardless of the surface. However, the ability of M. xanthus to search out patches of these prey is affected by surface hardness, the density of prey patches, and the prey species.     

Long-term adaptation of epistatic genetic networks

Gene networks are likely to govern most traits in nature. Mutations at these genes often show functional epistatic interactions that lead to complex genetic architectures and variable fitness effects in different genetic backgrounds. Understanding how epistatic genetic systems evolve in nature remains one of the great challenges in evolutionary biology. Here we combine an analytical framework with individual-based simulations to generate novel predictions about long-term adaptation of epistatic networks. We find that relative to traits governed by independently evolving genes, adaptation with epistatic gene networks is often characterized by longer waiting times to selective sweeps, lower standing genetic variation, and larger fitness effects of adaptive mutations. This may cause epistatic networks to either adapt more slowly or more quickly relative to a nonepistatic system. Interestingly, epistatic networks may adapt faster even when epistatic effects of mutations are on average deleterious. Further, we study the evolution of epistatic properties of adaptive mutations in gene networks. Our results show that adaptive mutations with small fitness effects typically evolve positive synergistic interactions, whereas adaptive mutations with large fitness effects evolve positive synergistic and negative antagonistic interactions at approximately equal frequencies. These results provide testable predictions for adaptation of traits governed by epistatic networks and the evolution of epistasis within networks.     

WHY ONTOGENY MATTERS DURING ADAPTATION: DEVELOPMENTAL NICHE CONSTRUCTION AND PLEIOTORPY ACROSS THE LIFE CYCLE IN ARABIDOPSIS THALIANA

This case study of adaptation in Arabidopsis thaliana shows that natural selection on early life stages can be intense and can influence the evolution of subsequent traits. Two mechanisms contribute to this influence: pleiotropy across developmental stages and developmental niche construction. Examples are given of pleiotropy of environmentally cued development across life stages, and potential ways that pleiotropy can be relieved are discussed. In addition, this case study demonstrates how the timing of prior developmental transitions determines the seasonal environment experienced subsequently, and that such developmental niche construction alters phenotypic expression of subsequent traits, the expression of genetic variation of those traits, and natural selection on those traits and alleles associated with them. As such, developmental niche construction modifies pleiotropic relationships across the life cycle in ways that influence the dynamics of adaptation. Understanding the genetic basis of life‐cycle variation therefore requires consideration of environmental effects on pleiotropy.     

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.     

The genetic architecture of domestication in the chicken: effects of pleiotropy and linkage

The extent of pleiotropy and epistasis in quantitative traits remains equivocal. In the case of pleiotropy, multiple quantitative trait loci are often taken to be pleiotropic if their confidence intervals overlap, without formal statistical tests being used to ascertain if these overlapping loci are statistically significantly pleiotropic. Additionally, the degree to which the genetic correlations between phenotypic traits are reflected in these pleiotropic quantitative trait loci is often variable, especially in the case of antagonistic pleiotropy. Similarly, the extent of epistasis in various morphological, behavioural and life-history traits is also debated, with a general problem being the sample sizes required to detect such effects. Domestication involves a large number of trade-offs, which are reflected in numerous behavioural, morphological and life-history traits which have evolved as a consequence of adaptation to selective pressures exerted by humans and captivity. The comparison between wild and domestic animals allows the genetic analysis of the traits that differ between these population types, as well as being a general model of evolution. Using a large F(2) intercross between wild and domesticated chickens, in combination with a dense SNP and microsatellite marker map, both pleiotropy and epistasis were analysed. The majority of traits were found to segregate in 11 tight 'blocks' and reflected the trade-offs associated with domestication. These blocks were shown to have a pleiotropic 'core' surrounded by more loosely linked loci. In contrast, epistatic interactions were almost entirely absent, with only six pairs identified over all traits analysed. These results give insights both into the extent of such blocks in evolution and the development of domestication itself.     

Epistatic interaction maps relative to multiple metabolic phenotypes

An epistatic interaction between two genes occurs when the phenotypic impact of one gene depends on another gene, often exposing a functional association between them. Due to experimental scalability and to evolutionary significance, abundant work has been focused on studying how epistasis affects cellular growth rate, most notably in yeast. However, epistasis likely influences many different phenotypes, affecting our capacity to understand cellular functions, biochemical networks adaptation, and genetic diseases. Despite its broad significance, the extent and nature of epistasis relative to different phenotypes remain fundamentally unexplored. Here we use genome-scale metabolic network modeling to investigate the extent and properties of epistatic interactions relative to multiple phenotypes. Specifically, using an experimentally refined stoichiometric model for Saccharomyces cerevisiae, we computed a three-dimensional matrix of epistatic interactions between any two enzyme gene deletions, with respect to all metabolic flux phenotypes. We found that the total number of epistatic interactions between enzymes increases rapidly as phenotypes are added, plateauing at approximately 80 phenotypes, to an overall connectivity that is roughly 8-fold larger than the one observed relative to growth alone. Looking at interactions across all phenotypes, we found that gene pairs interact incoherently relative to different phenotypes, i.e. antagonistically relative to some phenotypes and synergistically relative to others. Specific deletion-deletion-phenotype triplets can be explained metabolically, suggesting a highly informative role of multi-phenotype epistasis in mapping cellular functions. Finally, we found that genes involved in many interactions across multiple phenotypes are more highly expressed, evolve slower, and tend to be associated with diseases, indicating that the importance of genes is hidden in their total phenotypic impact. Our predictions indicate a pervasiveness of nonlinear effects in how genetic perturbations affect multiple metabolic phenotypes. The approaches and results reported could influence future efforts in understanding metabolic diseases and the role of biochemical regulation in the cell.     

Epistasis and the evolutionary dynamics of measured genotypic values during simulated serial bottlenecks

The evolutionary effects of epistasis have been primarily explored analytically and most empirical studies have utilized yeast, viral and bacterial populations. Empirical analyses in multi‐cellular organisms are rare because of experimental constraints. Here, we report the results of a genome‐wide scan for two‐way epistasis in 16 traits related to body size and composition in F₂ mice from the LG/J by SM/J intercross. We analyze two‐locus genotypic values at quantitative trait loci (QTL), which provides an especially detailed view of epistatic architectures, to evaluate their predicted evolutionary consequences via Monte Carlo simulations. Epistatic profiles vary, but all traits show complicated genetic architectures which are largely hidden in single locus QTL scans. On average, detected epistatic effects are comparable in size to marginal effects. Simulations demonstrate an expected preservation, and often inflation, of heritable variance across several generations of small effective population size for many identified epistatic pairs over a range of starting allele frequencies.     

Complex genotype interactions influence social fitness during the developmental phase of the social amoeba Dictyostelium discoideum

When individuals interact, phenotypic variation can be partitioned into direct genetic effects (DGEs) of the individuals’ own genotypes, indirect genetic effects (IGEs) of their social partners’ genotypes and epistatic interactions between the genotypes of interacting individuals (‘genotype‐by‐genotype (G×G) epistasis’). These components can all play important roles in evolutionary processes, but few empirical studies have examined their importance. The social amoeba Dictyostelium discoideum provides an ideal system to measure these effects during social interactions and development. When starved, free‐living amoebae aggregate and differentiate into a multicellular fruiting body with a dead stalk that holds aloft viable spores. By measuring interactions among a set of natural strains, we quantify DGEs, IGEs and G×G epistasis affecting spore formation. We find that DGEs explain most of the phenotypic variance (57.6%) whereas IGEs explain a smaller (13.3%) but highly significant component. Interestingly, G×G epistasis explains nearly a quarter of the variance (23.0%), highlighting the complex nature of genotype interactions. These results demonstrate the large impact that social interactions can have on development and suggest that social effects should play an important role in developmental evolution in this system.