Search

The role of epistasis in shaping genetic variation and contributing to observable differences within and between populations has been the focus of much debate (1-3). In complex trait genetics, the [...], The importance of epistasis--non-additive interactions between alleles--in shaping population fitness has long been a controversial topic, hampered in part by lack of empirical evidence (1-4). Traditionally, epistasis is inferred on the basis of non-independence of genotypic values between loci for a given trait. However, epistasis for fitness should also have a genomic footprint (5-7). To capture this signal, we have developed a simple approach that relies on detecting genotype ratio distortion as a sign of epistasis, and we apply this method to a large panel of Drosophila melanogaster recombinant inbred lines (8,9). Here we confirm experimentally that instances of genotype ratio distortion represent loci with epistatic fitness effects; we conservatively estimate that any two haploid genomes in this study are expected to harbour 1.15 pairs of epistatically interacting alleles. This observation has important implications for speciation genetics, as it indicates that the raw material to drive reproductive isolation is segregating contemporaneously within species and does not necessarily require, as proposed by the Dobzhansky-Muller model, the emergence of incompatible mutations independently derived and fixed in allopatry. The relevance of our result extends beyond speciation, as it demonstrates that epistasis is widespread but that it may of ten go undetected owing to lack of statistical power or lack of genome-wide scope of the experiments. A C G A C G

Author(s): Gibran Hemani [sup.1] [sup.2] , Konstantin Shakhbazov [sup.1] [sup.2] , Harm-Jan Westra [sup.3] , Tonu Esko [sup.4] [sup.5] [sup.6] , Anjali K. Henders [sup.7] , Allan F. McRae [sup.1] [...]

The main forces directing long-term molecular evolution remain obscure. A sizable fraction of amino-acid substitutions seem to be fixed by positive selection (1-4), but it is unclear to what degree [...]

Several human genetic disorders of hemoglobin have risen in frequency because of the protection they offer against death from malaria, sickle-cell anemia being a canonical example. Here we address the issue of why this highly protective mutant, present at high frequencies in subSaharan Africa, is uncommon in Mediterranean populations that instead harbor a diverse range of thalassemic hemoglobin disorders. We demonstrate that these contrasting profiles of malaria-protective alleles can arise and be stably maintained by two well-documented phenomena: an alleviation of the clinical severity of [alpha]-and [beta]-thalassemia in compound thalassemic genotypes and a cancellation of malaria protection when [alpha]-thalassemia and the sickle-cell trait are coinherited. The complex distribution of globin mutants across Africa and the Mediterranean can therefore be explained by their specific intracellular interactions. epistasis | host genetics | malaria | thalassemia | human evolution doi/10.1073/pnas.0910840106

Evolutionary theory predicts that a population in a new environment will accumulate adaptive substitutions, but precisely how they accumulate is poorly understood. The dynamics of adaptation depend on the underlying fitness landscape. Virtually nothing is known about fitness landscapes in nature, and few methods allow us to infer the landscape from empirical data. With a view toward this inference problem, we have developed a theory that, in the weak-mutation limit, predicts how a population's mean fitness and the number of accumulated substitutions are expected to increase over time, depending on the underlying fitness landscape. We find that fitness and substitution trajectories depend not on the full distribution of fitness effects of available mutations but rather on the expected fixation probability and the expected fitness increment of mutations. We introduce a scheme that classifies landscapes in terms of the qualitative evolutionary dynamics they produce. We show that linear substitution trajectories, long considered the hallmark of neutral evolution, can arise even when mutations are strongly selected. Our results provide a basis for understanding the dynamics of adaptation and for inferring properties of an organism's fitness landscape from temporal data. Applying these methods to data from a long-term experiment, we infer the sign and strength of epistasis among beneficial mutations in the Escherichia coil genome. epistasis | fitness trajectory | substitution trajectory | weak mutation | evolution www.pnas.org/cgi/doi/10.1073/pnas.0905497106

The epistatic interactions among mutations have a large effect on the evolution of populations. In this article we provide a formalism under which epistatic interactions among pairs of mutations have a distribution whose mean can be modulated. We find that the mean epistasis is correlated to the effect of mutations or genetic robustness, which suggests that such formalism is in good agreement with most in silico models of evolution where the same pattern is observed. We further show that the evolution of epistasis is highly dependant on the intensity of drift and of how complex the organisms are, and that either positive or negative epistasis could be selected for, depending on the balance between the efficiency of selection and the intensity of drift. DOI: 10.1534/genetics.108.099127

In Drosophila, the female-specific SEX-LETHAL (SXL) protein is required for oogenesis, but how Sxl interfaces with the genetic circuitry controlling oogenesis remains unknown. Here we use an allele of sans fille (snf) that specifically eliminates SXL protein in germ cells to carry out a detailed genetic and cell biological analysis of the resulting ovarian tumor phenotype. We find that tumor growth requires both Cyclin Band zero population growth, demonstrating that these mutant cells retain at least some of the essential growth-control mechanisms used by wild-type germ cells. Using a series of molecular markers, we establish that while the tumor often contains at least one apparently bona fide germline stem cell, the majority of cells exhibit an intermediate fate between a stem cell and its daughter cell fated to differentiate. In addition, snf tumors misexpress a select group of testis-enriched markers, which, remarkably, are also misexpressed in ovarian tumors that arise from the loss of bag of marbles (barn). Results of genetic epistasis experiments further reveal that barn's differentiation-promoting function depends on Sxl. Together these data demonstrate a novel role for Sxl in the lineage progression from stem cell to committed daughter cell and suggest a model in which Sxl partners with barn to facilitate this transition. DOI: 10.1534/genetics.109.100693

Biochemical and regulatory interactions central to biological networks are expected to cause extensive genetic interactions or epistasis affecting the heritability of complex traits and the distribution of genotypes in populations. However, the inference of epistasis from the observed phenotype--genotype correlation is impeded by statistical difficulties, while the theoretical understanding of the effects of epistasis remains limited, in turn limiting our ability to interpret data. Of particular interest is the biologically relevant situation of numerous interacting genetic loci with small individual contributions to fitness. Here, we present a computational model of selection dynamics involving many epistatic loci in a recombining population. We demonstrate that a large number of polymorphic interacting loci can, despite frequent recombination, exhibit cooperative behavior that locks alleles into favorable genotypes leading to a population consisting of a set of competing clones. When the recombination rate exceeds a certain critical value that depends on the strength of epistasis, this 'genotype selection' regime disappears in an abrupt transition, giving way to 'allele selection'--the regime where different loci are only weakly correlated as expected in sexually reproducing populations. We show that large populations attain highest fitness at a recombination rate just below critical. Clustering of interacting sets of genes on a chromosome leads to the emergence of an intermediate regime, where blocks of cooperating alleles lock into genetic modules. These haplotype blocks disappear in a second transition to pure allele selection. Our results demonstrate that the collective effect of many weak epistatic interactions can have dramatic effects on the population structure. gene interactions | population genetics

Epistasis has sometimes been reported in qualitative and quantitative traits of soybean [Glycine max (L.) Merr.]. Our objective was to determine the extent to which epistasis occurred in crosses between soybean lines adapted to similar latitudes in Liaoning, China, and in Ohio, USA. We crossed 'Tiefeng #27' x HS97-4534 and 'Ohio FG1' x 'Shennong #6', and developed random F4-derived lines from each cross and [BC1F.sub.3]-derived lines from the backcross to each parent. We tested the lines, along with parents and their checks, at Plain City, OH, in 2005 to 2007, and at Shenyang, Liaoning, in 2005 and 2006. Comparison of the means of parents, biparental lines, and both backcrosses revealed significant epistasis for plant height, maturity, 100-seed weight, yield, internode length, number of branches, harvest index, and content of protein and oil, but the epistasis was not expressed consistently across environments or parental combinations. Further, in some cases simple digenic epistasis did not account for the results, implying higher-order interactions among loci. Some of the results could be explained by interactions of loci controlling timing of reproductive stages and interactions between such loci and the environment. Breeders who make crosses between soybean lines adapted to northeastern China and the midwestern United States should be prepared for a high frequency of poorly adapted progeny.

Libraries of near-isogenic lines (NILs) are a powerful plant genetic resource to map quantitative trait loci (QTL). Nevertheless, QTL mapping with NILs is mostly restricted to genetic main effects. Here we propose a two-step procedure to map additive-by-additive digenic epistasis with NILs. In the first step, a generation means analysis of parents, their [F.sub.1] hybrid, and one-segment NILs and their triple testcross (TFC) progenies is used to identify in a one-dimensional scan loci exhibiting QTL-by-background interactions. In a second step, one-segment NILs with significant additive-by-additive background interactions are used to produce particular two-segment NILs to test for digenic epistatic interactions between these segments. We evaluated our approach by analyzing a random subset of a genomewide Arabidopsis thaliana NIL library for growth-related traits. The results of our experimental study illustrated the potential of the presented two-step procedure to map additive-by-additive digenic epistasis with NILs. Furthermore, our findings suggested that additive main effects as well as additive-by-additive digenic epistasis strongly influence the genetic architecture underlying growth-related traits of A. thaliana.

Protein evolution is crucial for organismal adaptation and fitness. This process takes place by shaping a given 3-dimensional fold for its particular biochemical function within the metabolic requirements and constraints of the environment. The complex interplay between sequence, structure, functionality, and stability that gives rise to a particular phenotype has limited the identification of traits acquired through evolution. This is further complicated by the fact that mutations are pleiotropic, and interactions between mutations are not always understood. Antibiotic resistance mediated by [beta]-lactamases represents an evolutionary paradigm in which organismal fitness depends on the catalytic efficiency of a single enzyme. Based on this, we have dissected the structural and mechanistic features acquired by an optimized metallo-[beta]-lactamase (M[beta]L) obtained by directed evolution. We show that antibiotic resistance mediated by this enzyme is driven by 2 mutations with sign epistasis. One mutation stabilizes a catalytically relevant intermediate by fine tuning the position of 1 metal ion; whereas the other acts by augmenting the protein flexibility. We found that enzyme evolution (and the associated antibiotic resistance) occurred at the expense of the protein stability, revealing that M[beta]Ls have not exhausted their stability threshold. Our results demonstrate that flexibility is an essential trait that can be acquired during evolution on stable protein scaffolds. Directed evolution aided by a thorough characterization of the selected proteins can be successfully used to predict future evolutionary events and design inhibitors with an evolutionary perspective. antibiotic resistance | enzyme | fitness landscape | metalloproteins | epistasis

The genetic architecture of complex traits underlying physiology and disease in most organisms remains elusive. We still know little about the number of genes that underlie these traits, the magnitude of their effects, or the extent to which they interact. Chromosome substitution strains (CSSs) enable statistically powerful studies based on testing engineered inbred strains that have single, unique, and nonoverlapping genetic differences, thereby providing measures of phenotypic effects that are attributable to individual chromosomes. Here, we report a study of phenotypic effects and gene interactions for 90 blood, bone, and metabolic traits in a mouse CSS panel and 54 traits in a rat CSS panel. Two key observations emerge about the genetic architecture of these traits. First, the traits tend to be highly polygenic: across the genome, many individual chromosome substitutions each had significant phenotypic effects and, within each of the chromosomes studied, multiple distinct loci were found. Second, strong epistasis was found among the individual chromosomes. Specifically, individual chromosome substitutions often conferred surprisingly large effects (often a substantial fraction of the entire phenotypic difference between the parental strains), with the result that the sum of these individual effects often dramatically exceeded the difference between the parental strains. We suggest that strong, pervasive epistasis may reflect the presence of several phenotypically-buffered physiological states. These results have implications for identification of complex trait genes, developmental and physiological studies of phenotypic variation, and opportunities to engineer phenotypic outcomes in complex biological systems. chromosome substitution | genetic variation | quantitative trait loci

The question of what role sex plays in evolution is still open despite decades of research. It has often been assumed that sex should facilitate the increase in fitness. Hence, the fact that it may break down highly favorable genetic combinations has been seen as a problem. Here, we consider an alternative approach. We define a measure that represents the ability of alleles to perform well across different combinations and, using numerical iterations within a classical population-genetic framework, show that selection in the presence of sex favors this ability in a highly robust manner. We also show that the mechanism responsible for this effect has been out of the purview of previous theory, because it operates during the evolutionary transient, and that the breaking down of favorable genetic combinations is an integral part of it. Implications of these results and more to evolutionary theory are discussed. recombination | modularity | fitness robustness | evolvability | epistasis

Muscle histochemical characteristics are decisive determinants of meat quality. The relative percentage and diameters of the different muscular fiber types influence crucial aspects of meat such as color, tenderness, and ultimate pH. Despite its relevance, however, the information on muscle fiber genetic architecture is scant, because histochemical muscle characterization is a laborious task. Here we report a complete QTL scan of muscle fiber traits in 160 animals from a [F.sub.2] cross between Iberian and Landrace pigs using 139 markers. We identified 20 genome regions distributed along 15 porcine chromosomes (SSC1, 2, 3, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, and X) with direct and(or) epistatic effects. Epistasis was frequent and some interactions were highly significant. Chromosomes 10 and 11 seemed to behave as hubs; they harbored 2 individual QTL, but also 6 epistatic regions, Numerous individual QTL effects had cryptic alleles, with opposite effects to phenotypic pure breed differences. Many of the QTL identified here coincided with previous reports for these traits in the literature, and there was overlapping with potential candidate genes and previously reported meat quality QTL. Key words: epistasis, muscle fiber, overdominance, pig, quantitative trait locus

We consider the dynamics of a nonrecombining haploid population of finite size that accumulates deleterious mutations irreversibly. This ratchet-like process occurs at a finite speed in the absence of epistasis, but it has been suggested that synergistic epistasis can halt the ratchet. Using a diffusion theory, we find explicit analytical expressions for the typical time between successive clicks of the ratchet for both nonepistatic and epistatic fitness functions. Our calculations show that the interclick time is of a scaling form that in the absence of epistasis gives a speed that is determined by size of the least-loaded class and the selection coefficient. With synergistic interactions, the ratchet speed is found to approach zero rapidly for arbitrary epistasis. Our analytical results are in good agreement with the numerical simulations.

A new measure of directional linkage disequilibrium is developed for detecting epistatic selection on interacting genes. Simulations show that by orienting the direction of linkage disequilibrium on the basis of the ancestral-derived status of alleles, the new measure indeed improves the power to detect a positive fitness interaction between two new mutations.

In this contribution, we study the genetic mechanisms leading to differences in the observed growth patterns of domesticated White Leghorn chickens and their wild ancestor the red jungle fowl. An epistatic QTL analysis for several body-weight measures from hatch to adulthood confirms earlier findings that polymorphisms at >15 loci contribute to body-weight determination in an F2 intercross between these populations and that many loci are involved in complex genetic interactions. Here, we use a new genetic model to decompose the genetic effects of this multilocus epistatic genetic network. The results show how the functional modeling of genetic effects provides new insights into how genetic interactions in a large set of loci jointly contribute to phenotypic expression. By exploring the functional effects of QTL alleles, we show that some alleles can display temporal shifts in the expression of genetic effects due to their dependencies on the genetic background. Our results demonstrate that the effects of many genes are dependent on genetic interactions with other loci and how their involvement in the domestication process relies on these interactions.

Describing at a genomic scale how mutations in different genes influence one another is essential to the understanding of how genotype correlates with phenotype and remains a major challenge in biology. Previous studies pointed out the need for accurate measurements of not only synthetic but also buffering interactions in the characterization of genetic networks and functional modules. We developed a sensitive and efficient method that allows such measurements at a genomic scale in yeast. In a pilot experiment (41 genome-wide screens), we quantified the fitness of 140,000 double deletion strains relative to the corresponding single mutants and identified many genetic interactions. In addition to synthetic growth defects (validated experimentally with factors newly identified as genetically interfering with mRNA degradation), most of the identified genetic interactions measured weak epistatic effects. These weak effects, rarely meaningful when considered individually, were crucial to defining specific signatures for many gene deletions and had a major contribution in defining clusters of functionally related genes. epistasis | functional genomics | genetic screen | mRNA decapping | Saccharomyces cerevisiae

The CpxAR two-component signal transduction system in Escherichia coli and other pathogens senses diverse envelope stresses and promotes the transcription of a variety of genes that remedy these stresses. An important member of the CpxAR regulon is cpxP. The CpxA-dependent transcription of cpxP has been linked to stresses such as misfolded proteins and alkaline pH. It also has been proposed that acetyl phosphate, the intermediate of the phosphotransacetylase (Pta)-acetate kinase (AckA) pathway, can activate the transcription of cpxP in a CpxA-independent manner by donating its phosphoryl group to CpxR. We tested this hypothesis by measuring the transcription of cpxP using mutants with mutations in the CpxAR pathway, mutants with mutations in the Pta-AckA pathway, and mutants with a combination of both types of mutations. From this epistasis analysis, we learned that CpxR integrates diverse stimuli. The stimuli that originate in the envelope depend on CpxA, while those associated with growth and central metabolism depend on the Pta-AckA pathway. While CpxR could receive a phosphoryi group from acetyl phosphate, this global signal was not the primary trigger for CpxR activation associated with the Pta-AckA pathway. On the strength of these results, we contend that the interactions between central metabolism and signal transduction can be quite complex and that successful investigations of such interactions must include a complete epistatic analysis.

Sometimes mutations in two genes produce a phenotype that is surprising in light of each mutation's individual effects. This phenomenon, which defines genetic interaction, can reveal functional relationships between genes and pathways. For example, double mutants with surprisingly slow growth define synergistic interactions that can identify compensatory pathways or protein complexes. Recent studies have used four mathematically distinct definitions of genetic interaction (here termed Product, Additive, Log, and Min). Whether this choice holds practical consequences has not been clear, because the definitions yield identical results under some conditions. Here, we show that the choice among alternative definitions can have profound consequences. Although 52% of known synergistic genetic interactions in Saccharomyces cerevisiae were inferred according to the Min definition, we find that both Product and Log definitions (shown here to be practically equivalent) are better than Min for identifying functional relationships. Additionally, we show that the Additive and Log definitions, each commonly used in population genetics, lead to differing conclusions related to the selective advantages of sexual reproduction. epistasis | fitness | gene function

Divergent selection for oil and protein concentration in the corn (Zea mays L.) kernel was initiated at the University of Illinois in 1896 by C.G. Hopkins. In 2005, 106 generations of selection had been completed for high oil and 105 for high protein Limits to selection for low oil and low protein were reached but not for high oil or high protein. Over the more than 100 yr of the existence of the program a number of attempts have been made to analyze the experiment using quantitative genetic tools. The purpose of this paper is to trace the use of Quantitative genetic techniques to analyze the results of divergent long-term selection for oil and protein, to relate results to the question of the need for divergent parents for quantitative trait locus (QTL) analysis, and to provide a new look at the reasons for long-continued progress from selection. Key findings include (i) progress from selection was much greater than could nave been predicted; (ii) based on both classical quantitative genetic analysis and QTL studies, a large number of QTL are involved in control of the three traits; (iii) the number of QTL identified in a given study cannot be predicted by the magnitude of genetic variance or the divergence of the parents but is a function of the number of markers used and the number of lines evaluated; and (iv) epistasis may be an important factor in explaining long-term response to selection.

Epistasis has been proposed as a possible explanation for the long continued progress from selection in the Illinois long-term corn (Zea mays L.) selection strains. From the crosses of Illinois High Oil (IHO) x Illinois Low Oil (ILO) and of Illinois High Protein (IHP) x Illinois Low Protein (ILP), 500 [S.sub.2] lines were developed. Each IHOxILO line was genotyped for 479 single nucleotide polymorphism (SNP) markers, and each IHPxILP line was genotyped for 499 SNP markers. Per se and testcross progenies of each [S.sub.2] line were evaluated for oil, protein, and starch. A large number of QTL were found to control these traits. The objective of this paper is to report the results of analysis of these crosses for two-way epistatic interactions, More epistatic interactions were significant than expected by chance. The proportion of significant interactions that involved a marker significant in a single marker analysis (SMA) was no greater than expected by chance. The number of markers associated only with significant epistatic effects ranged from 46.3 to 72.2% of the total number of markers significant for either an interaction effect or from SMA, The number of QTL controlling a trait is much greater than will be found by analyzing for significant QTL main effects. Thus, epistasis could contribute to the long continued response to selection in the Illinois long-term selection strains and also may help explain the continued success of commercial corn breeding.

Inwardly rectifying potassium (Kir) channels form gates in the cell membrane that regulate the flow of [K.sup.+] ions into and out of the cell, thereby influencing the membrane potential and electrical signaling of many cell types, including neurons and cardiomyocytes. Kir-channel function depends on other cellular proteins that aid in the folding of channel subunits, assembly into tetrameric complexes, trafficking of quality-controlled channels to the plasma membrane, and regulation of channel activity at the cell surface. We used the yeast Saccharomyces cerevisiae as a model system to identify proteins necessary for the functional expression of a mammalian Kir channel at the cell surface. A screen of 376 yeast strains, each lacking one nonessential protein localized to the early secretory pathway, identified seven deletion strains in which functional expression of the Kir channel at the plasma membrane was impaired. Six deletions were of genes with known functions in trafficking and lipid biosynthesis (sur4[DELTA], csg2[DELTA], erv14[DELTA], emp24[DELTA], erv25[DELTA], and bst1[DELTA]), and one deletion was of an uncharacterized gene (yil39w[DELTA]). We provide genetic and functional evidence that Yil039wp, a conserved, phosphoesterase domain-containing protein, which we named 'trafficking of Emp24p/ Erv25p-dependent cargo disrupted 1' (Ted1p), acts together with Emp24p/Erv25p in cargo exit from the endoplasmic reticulum (ER). The seven yeast proteins identified in our screen likely impact Kir-channel functional expression at the level of vesicle budding from the ER and/or the local lipid environment at the plasma membrane. yeast deletion library | YIL039W/TED1 | GPI anchor | GIRK2 | epistasis mini array profile

Evolutionary biologists have identified several factors that could explain the widespread phenomena of sex and recombination. One hypothesis is that host-parasite interactions favor sex and recombination because they favor the production of rare genotypes. A problem with many of the early models of this socalled Red Queen hypothesis is that several factors are acting together: directional selection, fluctuating epistasis, and drift. It is thus difficult to identify what exactly is selecting for sex in these models. Is one factor more important than the others or is it the synergistic action of these different factors that really matters? Here we focus on the analysis of a simple model with a single mechanism that might select for sex: fluctuating epistasis. We first analyze the evolution of sex and recombination when the temporal fluctuations are driven by the abiotic environment. We then analyze the evolution of sex and recombination in a two-species coevolutionary model, where directional selection is absent (allele frequencies remain fixed) and temporal variation in epistasis is induced by coevolution with the antagonist species. In both cases we contrast situations with weak and strong selection and derive the evolutionarily stable (ES) recombination rate. The ES recombination rate is most sensitive to the period of the cycles, which in turn depends on the strength of epistasis. In particular, more virulent parasites cause more rapid cycles and consequently increase the ES recombination rate of the host. Although the ES strategy is maximized at an intermediate period, some recombination is favored even when fluctuations are very slow. By contrast, the amplitude of the cycles has no effect on the ES level of sex and recombination, unless sex and recombination are costly, in which case higher-amplitude cycles allow the evolution of higher rates of sex and recombination. In the coevolutionary model, the amount of recombination in the interacting species also has a large effect on the ES, with evolution favoring higher rates of sex and recombination than in the interacting species. In general, the ES recombination rate is less than or equal to the recombination rate that would maximize mean fitness. We also discuss the effect of migration when sex and recombination evolve in a metapopulation. We find that intermediate parasite migration rates maximize the degree of local adaptation of the parasite and lead to a higher ES recombination rate in the host.

Quantitative trait loci (QTL) method was used to study the genetic architecture of traits associated with immune defense and reproduction in male Bombus terrestris bumblebees. The study revealed phenotypic variance based on epistatic interactions.

Why are most genes dispensable? The impact of gene deletions may depend on the environment (plasticity), the presence of compensatory mechanisms (mutational robustness), or both. Here, we analyze the interaction between these two forces by exploring the condition-dependence of synthetic genetic interactions that define redundant functions and alternative pathways. We performed systems-level flux balance analysis of the yeast (Saccharomyces cerevisiae) metabolic network to identify genetic interactions and then tested the model's predictions with in vivo gene-deletion studies. We found that the majority of synthetic genetic interactions are restricted to certain environmental conditions, partly because of the lack of compensation under some (but not all) nutrient conditions. Moreover, the phylogenetic cooccurrence of synthetically interacting pairs is not significantly different from random expectation. These findings suggest that these gene pairs have at least partially independent functions, and, hence, compensation is only a byproduct of their evolutionary history. Experimental analyses that used multiple gene deletion strains not only confirmed predictions of the model but also showed that investigation of false predictions may both improve functional annotation within the model and also lead to the discovery of higher-order genetic interactions. Our work supports the view that functional redundancy may be more apparent than real, and it offers a unified framework for the evolution of environmental adaptation and mutational robustness. epistasis | genetic robustness | Saccharomyces cerevisiae | environmental dependence I flux balance analysis

Vibrio cholerae undergoes phenotypic variation that generates two morphologically different variants, termed smooth and rugose. The transcriptional profiles of the two variants differ greatly, and many of the differentially regulated genes are controlled by a complex regulatory circuitry that includes the transcriptional regulators VpsR, VpsT, and HapR. In this study, we identified the VpsT regulon and compared the VpsT and VpsR regulons to elucidate the contribution of each positive regulator to the rugose variant transcriptional profile and associated phenotypes. We have found that although the VpsT and VpsR regulons are very similar, the magnitude of the gene regulation accomplished by each regulator is different. We also determined that cdgA, which encodes a GGDEF domain protein, is partially responsible for the altered vps gene expression between the vpsT and vpsR mutants. Analysis of epistatic relationships among hapR, vpsT, and vpsR with respect to a whole-genome expression profile, colony morphology, and biofiim formation revealed that vpsR is epistatic to hapR and vpsT. Expression of virulence genes was increased in a vpsR hapR double mutant relative to a hapR mutant, suggesting that VpsR negatively regulates virulence gene expression in the hapR mutant. These results show that a complex regulatory interplay among VpsT, VpsR, HapR, and GGDEF/EAL family proteins controls transcription of the genes required for Vibrio polysaccharide and virulence factor production in V. cholerae.

The influence of temperature on cyto-nuclear effects on egg-to-adult development time in seed beetles is assessed. The findings show that environmental factors such as temperature do exert selection on cytoplasmic genes by favoring specific cyto-nuclear genetic combinations, and are consistent with the suggestion that complex genotype-by-environment interactions might promote the maintenance of polymorphism in mitochondrial genes.

Functional dependencies between genes are a defining characteristic of gene networks underlying quantitative traits. However, recent studies show that the proportion of the genetic variation that can be attributed to statistical epistasis varies from almost zero to very high. It is thus of fundamental as well as instrumental importance to better understand whether different functional dependency patterns among polymorphic genes give rise to distinct statistical interaction patterns or not. Here we address this issue by combining a quantitative genetic model approach with genotype-phenotype models capable of translating allelic variation and regulatory principles into phenotypic variation at the level of gene expression. We show that gene regulatory networks with and without feedback motifs can exhibit a wide range of possible statistical genetic architectures with regard to both type of effect explaining phenotypic variance and number of apparent loci underlying the observed phenotypic effect. Although all motifs are capable of harboring significant interactions, positive feedback gives rise to higher amounts and more types of statistical epistasis. The results also suggest that the inclusion of statistical interaction terms in genetic models will increase the chance to detect additional QTL as well as functional dependencies between genetic loci over a broad range of regulatory regimes. This article illustrates how statistical genetic methods can fruitfully be combined with nonlinear systems dynamics to elucidate biological issues beyond reach of each methodology in isolation.

The symbiotic relationship between the mitochondrial and nuclear genomes coordinates metabolic energy production and is fundamental to life among eukaryotes. Consequently, there is potential for strong selection to shape interactions between these two genomes. Substantial research attention has focused on the possibility that within-population sequence polymorphism in mitochondrial DNA (mtDNA) is maintained by mitonuclear fitness interactions. Early theory predicted that selection will often eliminate mitochondrial polymorphisms. However, recent models demonstrate that intergenomic interactions can promote the maintenance of polymorphism, especially if the nuclear genes involved are linked to the X chromosome. Most empirical studies to date that have assessed cytonuclear fitness interactions have studied variation across populations and it is still unclear how general and strong such interactions are within populations. We experimentally tested for cytonuclear interactions within a laboratory population of Drosophila melanogaster using 25 randomly sampled cytoplasmic genomes, expressed in three different haploid nuclear genetic backgrounds, while eliminating confounding effects of intracellular bacteria (e.g., Wolbachia). We found sizable cytonuclear fitness interactions within this population and present limited evidence suggesting that these effects were sex specific. Moreover, the relative fitness of cytonuclear genotypes was environment specific. Sequencing of mtDNA (2752 bp) revealed polymorphism within the population, suggesting that the observed cytoplasmic genetic effects may be mitochondrial in origin.

Genetic interactions in fitness are studied by using modifier theory. The effects on fitness of two linked genes are perturbed by alleles at a third linked locus that controls the extent of epistasis in fitness between the first two. This epistasis is determined by a symmetric interaction matrix, and it is shown that a modifier allele that increases epistasis will invade when the linkage between the other two genes is sufficiently tight and these genes are in linkage disequilibrium. With linkage equilibrium among the major loci, increased or decreased epistasis may evolve depending on the allele frequencies at these loci. modifier theory | linkage disequilibrium | symmetricviabilities | interaction matrix | external stability

The LIN-3-LET-23-mediated inductive signaling pathway plays a major role during vulval development in C. elegans. Studies on the components of this pathway have revealed positive as well as negative regulators that function to modulate the strength and specificity of the signal transduction cascade. We have carried out genetic screens to identify new regulators of this pathway by screening for suppressors of lin-3 vulvaless phenotype. The screens recovered three loci including alleles of gap-1 and a new gene represented by sli-3. Our genetic epistasis experiments suggest that sli-3 functions either downstream or in parallel to nuclear factors lin-1 and sur-2. sli-3 synergistically interacts with the previously identified negative regulators of the let-23 signaling pathway and causes excessive cell proliferation. However, in the absence of any other mutation sli-3 mutant animals display wild-type vulval induction and morphology. We propose that sli-3 functions as a negative regulator of vulval induction and defines a branch of the inductive signaling pathway. We provide evidence that sli-3 interacts with the EGF signaling pathway components during vulval induction but not during viability and ovulation processes. Thus, sli-3 helps define specificity of the EGF signaling to induce the vulva.

Additive and nonadditive genetic effects on preweaning weight gain (PWG) of a commercial crossbred population were estimated using different genetic models and estimation methods. The data set consisted of 103,445 records on purebred and crossbred Nelore-Hereford calves raised under pasture conditions on farms located in south, southeast, and middle west Brazilian regions. In addition to breed additive and dominance effects, the models including different epistasis covariables were tested. Models considering joint additive and environment (latitude) by genetic effects interactions were also applied. In a first step, analyses were carried out under animal models. In a second step, preadjusted records were analyzed using ordinary least squares (OLS) and ridge regression (RR). The results reinforced evidence that breed additive and dominance effects are not sufficient to explain the observed variability in preweaning traits of Bos taurus x Bos indicus calves, and that genotype x environment interaction plays an important role in the evaluation of crossbred calves. Data were ill-conditioned to estimate the effects of genotype x environment interactions. Models including these effects presented multicolinearity problems. In this case, RR seemed to be a powerful tool for obtaining more plausible and stable estimates. Estimated prediction error variances and variance inflation factors were drastically reduced, and many effects that were not significant under ordinary least squares became significant under RR. Predictions of PWG based on RR estimates were more acceptable from a biological perspective. In temperate and subtropical regions, calves with intermediate genetic compositions (close to 1/2 Nelore) exhibited greater predicted PWG. In the tropics, predicted PWG increased linearly as genotype got closer to Nelore. Key words: crossbreeding, epistasis, genotype x environment interaction, heterosis, multicolinearity, ridge regression

The effect of inbreeding depression on survival in two populations of the rare species Brassica insularis is investigated using both controlled crosses and a marker-based approach. Results show that although combining the approaches is not necessary to obtain simple point estimates of inbreeding depression, using molecular markers might give insight into the genetic basis of inbreeding depression, such as the occurrence of epistatic interactions among deleterious alleles or purging.

Whether systematic genetic interactions (epistasis) occur at the genomic scale remains a challenging topic in evolutionary biology. Epistasis should make a significant contribution to variation in complex traits and influence the evolution of genetic systems as sex, diploidy, dominance, or the contamination of genomes with deleterious mutations. We have collected data from widely different organisms and quantified epistasis in a common, per-generation scale. Simpler genomes, such as those of RNA viruses, display antagonistic epistasis (mutations have smaller effects together than expected); bacterial microorganisms do not apparently deviate from independent effects, whereas in multicellular eukaryotes, a transition toward synergistic epistasis occurs (mutations have larger effects together than expected). We propose that antagonistic epistasis might be a property of compact genomes with few nonpleiotropic biological functions, whereas in complex genomes, synergism might emerge from mutational robustness. mutation | robustness | antagonism | evolution | synergism

We examine ionizing radiation (IR) sensitivity and epistasis relationships of several Saccharomyces mutants affecting post-translational modifications of histones H2B and H3. Mutants bre1[DELTA], lge1[DELTA], and rtf1[DELTA], defective in histone H2B lysine 123 ubiquitination, show IR sensitivity equivalent to that of the dot1[DELTA] mutant that we reported on earlier, consistent with published findings that Dot1p requires H2B K123 ubiquitination to fully methylate histone H3 K79. This implicates progressive K79 methylation rather than mono-methylation in IR resistance. The set2[DELTA] mutant, defective in H3 K36 methylation, shows mild IR sensitivity whereas mutants that abolish H3 K4 methylation resemble wild type. The dot1[DELTA], bre1[DELTA], and lge1[DELTA] mutants show epistasis for IR sensitivity. The paf1[DELTA] mutant, also reportedly defective in H2B K123 ubiquitination, confers no sensitivity. The rad6[DELTA], rad51null rad50[DELTA], and rad9A mutations are epistatic to bre1[DELTA] and dot1[DELTA], but rad18[DELTA] and rad5[DELTA] show additivity with bre1[DELTA], dot1[DELTA], and each other. The bre1[DELTA] rad18[DELTA] double mutant resembles rad6[DELTA] in sensitivity; thus the role of Rad6p in ubiquitinating H2B accounts for its extra sensitivity compared to rad18[DELTA]. We conclude that IR resistance conferred by BRE1 and DOT1 is mediated through homologous recombinational repair, not postreplication repair, and confirm findings of a Ga checkpoint role for the RAD6/BRE1/DOT1 pathway.

Motile bacteria follow gradients of attractant and repellent chemicals with high sensitivity. Their chemoreceptors are physically clustered, which may enable them to function as a cooperative array. Although native chemoreceptor molecules are typically transmembrane homodimers, they appear to associate through their cytoplasmic tips to form trimers of dimers, which may be an important architectural element in the assembly and operation of receptor clusters. The five receptors of Escherichia coli that mediate most of its chemotactic and aerotactic behaviors have identical trimer contact residues and have been shown by in vivo crosslinking methods to form mixed trimers of dimers. Mutations at the trimer contact sites of Tsr, the serine chemoreceptor, invariably abrogate Tsr function, but some of those lesions (designated Tsr*) are epistatic and block the function of heterologous chemoreceptors. We isolated and characterized mutations (designated Tar^) in the aspartate chemoreceptor that restored function to Tsr* receptors. The suppressors arose at or near the Tar trimer contact sites and acted in an allele-specific fashion on Tsr* partners. Alone, many Tar^ receptors were unable to mediate chemotactic responses to aspartate, but all formed clusters with varying efficiencies. Most of those Tar^ receptors were epistatic to WT Tsr, but some regained Tar function in combination with a suppressible Tsr* partner. Tar^-Tsr* suppression most likely occurs through compensatory changes in the conformation or dynamics of a mixed receptor signaling complex, presumably based on trimer-of-dimer interactions. These collaborative teams may be responsible for the high-gain signaling properties of bacterial chemoreceptors. chemotaxis | epistasis | receptor clustering | signaling teams | trimers of dimers