Prevalent positive epistasis in Escherichia coli and Saccharomyces cerevisiae metabolic networks

Exploring adaptation routes to cold temperatures in the Saccharomyces genus

The identification of traits that affect adaptation of microbial species to external abiotic factors, such as temperature, is key for our understanding of how biodiversity originates and can be maintained in a constantly changing environment. The Saccharomyces genus, which includes eight species with different thermotolerant profiles, represent an ideal experimental platform to study the impact of adaptive alleles in different genetic backgrounds. Previous studies identified a group of adaptive genes for maintenance of growth at lower temperatures. Here, we carried out a genus-wide assessment of the role of genes partially responsible for cold-adaptation in all eight Saccharomyces species for six candidate genes. We showed that the cold tolerance trait of S. kudriavzevii and S. eubayanus is likely to have evolved from different routes, involving genes important for the conservation of redox-balance, and for the long-chain fatty acid metabolism, respectively. For several loci, temperature- and species-dependent epistasis was detected, underscoring the plasticity and complexity of the genetic interactions. The natural isolates of S. kudriavzevii, S. jurei and S. mikatae had a significantly higher expression of the genes involved in the redox balance compared to S. cerevisiae, suggesting a role at transcriptional level. To distinguish the effects of gene expression from allelic variation, we independently replaced either the promoters or the coding sequences (CDS) of two genes in four yeast species with those derived from S. kudriavzevii. Our data consistently showed a significant fitness improvement at cold temperatures in the strains carrying the S. kudriavzevii promoter, while growth was lower upon CDS swapping. These results suggest that transcriptional strength plays a bigger role in growth maintenance at cold temperatures over the CDS and supports a model of adaptation centred on stochastic tuning of the expression network.

Pathway Evolution Through a Bottlenecking-Debottlenecking Strategy and Machine Learning-Aided Flux Balancing

The evolution of pathway enzymes enhances the biosynthesis of high-value chemicals, crucial for pharmaceutical, and agrochemical applications. However, unpredictable evolutionary landscapes of pathway genes often hinder successful evolution. Here, the presence of complex epistasis is identifued within the representative naringenin biosynthetic pathway enzymes, hampering straightforward directed evolution. Subsequently, a biofoundry-assisted strategy is developed for pathway bottlenecking and debottlenecking, enabling the parallel evolution of all pathway enzymes along a predictable evolutionary trajectory in six weeks. This study then utilizes a machine learning model, ProEnsemble, to further balance the pathway by optimizing the transcription of individual genes. The broad applicability of this strategy is demonstrated by constructing an Escherichia coli chassis with evolved and balanced pathway genes, resulting in 3.65 g L-1 naringenin. The optimized naringenin chassis also demonstrates enhanced production of other flavonoids. This approach can be readily adapted for any given number of enzymes in the specific metabolic pathway, paving the way for automated chassis construction in contemporary biofoundries.

Mitonuclear interactions modulate nutritional preference

In nature, organisms are faced with constant nutritional options which fuel key life-history traits. Studies have shown that species can actively make nutritional decisions based on internal and external cues. Metabolism itself is underpinned by complex genomic interactions involving components from both nuclear and mitochondrial genomes. Products from these two genomes must coordinate how nutrients are extracted, used and recycled. Given the complicated nature of metabolism, it is not well understood how nutritional choices are affected by mitonuclear interactions. This is under the rationale that changes in genomic interactions will affect metabolic flux and change physiological requirements. To this end we used a large Drosophila mitonuclear genetic panel, comprising nine isogenic nuclear genomes coupled to nine mitochondrial haplotypes, giving a total of 81 different mitonuclear genotypes. We use a capillary-based feeding assay to screen this panel for dietary preference between carbohydrate and protein. We find significant mitonuclear interactions modulating nutritional choices, with these epistatic interactions also being dependent on sex. Our findings support the notion that complex genomic interactions can place a constraint on metabolic flux. This work gives us deeper insights into how key metabolic interactions can have broad implications on behaviour.

Leveraging the genetic correlation between traits improves the detection of epistasis in genome-wide association studies

Epistasis, commonly defined as the interaction between genetic loci, is known to play an important role in the phenotypic variation of complex traits. As a result, many statistical methods have been developed to identify genetic variants that are involved in epistasis, and nearly all of these approaches carry out this task by focusing on analyzing one trait at a time. Previous studies have shown that jointly modeling multiple phenotypes can often dramatically increase statistical power for association mapping. In this study, we present the "multivariate MArginal ePIstasis Test" (mvMAPIT)-a multioutcome generalization of a recently proposed epistatic detection method which seeks to detect marginal epistasis or the combined pairwise interaction effects between a given variant and all other variants. By searching for marginal epistatic effects, one can identify genetic variants that are involved in epistasis without the need to identify the exact partners with which the variants interact-thus, potentially alleviating much of the statistical and computational burden associated with conventional explicit search-based methods. Our proposed mvMAPIT builds upon this strategy by taking advantage of correlation structure between traits to improve the identification of variants involved in epistasis. We formulate mvMAPIT as a multivariate linear mixed model and develop a multitrait variance component estimation algorithm for efficient parameter inference and P-value computation. Together with reasonable model approximations, our proposed approach is scalable to moderately sized genome-wide association studies. With simulations, we illustrate the benefits of mvMAPIT over univariate (or single-trait) epistatic mapping strategies. We also apply mvMAPIT framework to protein sequence data from two broadly neutralizing anti-influenza antibodies and approximately 2,000 heterogeneous stock of mice from the Wellcome Trust Centre for Human Genetics. The mvMAPIT R package can be downloaded at https://github.com/lcrawlab/mvMAPIT.

What Has Genomics Taught An Evolutionary Biologist?

Genomics, an interdisciplinary field of biology on the structure, function, and evolution of genomes, has revolutionized many subdisciplines of life sciences, including my field of evolutionary biology, by supplying huge data, bringing high-throughput technologies, and offering a new approach to biology. In this review, I describe what I have learned from genomics and highlight the fundamental knowledge and mechanistic insights gained. I focus on three broad topics that are central to evolutionary biology and beyond-variation, interaction, and selection-and use primarily my own research and study subjects as examples. In the next decade or two, I expect that the most important contributions of genomics to evolutionary biology will be to provide genome sequences of nearly all known species on Earth, facilitate high-throughput phenotyping of natural variants and systematically constructed mutants for mapping genotype-phenotype-fitness landscapes, and assist the determination of causality in evolutionary processes using experimental evolution.

Systematic evaluation of genome-wide metabolic landscapes in lactic acid bacteria reveals diet- and strain-specific probiotic idiosyncrasies

Lactic acid bacteria (LAB) are well known to elicit health benefits in humans, but their functional metabolic landscapes remain unexplored. Here, we analyze differences in growth, intestinal persistence, and postbiotic biosynthesis of six representative LAB and their interactions with 15 gut bacteria under 11 dietary regimes by combining multi-omics and in silico modeling. We confirmed predictions on short-term persistence of LAB and their interactions with commensals using cecal microbiome abundance and spent-medium experiments. Our analyses indicate that probiotic attributes are both diet and species specific and cannot be solely explained using genomics. For example, although both Lacticaseibacillus casei and Lactiplantibacillus plantarum encode similarly sized genomes with diverse capabilities, L. casei exhibits a more desirable phenotype. In addition, "high-fat/low-carb" diets more likely lead to detrimental outcomes for most LAB. Collectively, our results highlight that probiotics are not "one size fits all" health supplements and lay the foundation for personalized probiotic design.

Epistasis and evolutionary dependencies in human cancers

Cancer evolution is driven by the concerted action of multiple molecular alterations, which emerge and are selected during tumor progression. An alteration is selected when it provides an advantage to the tumor cell. However, the advantage provided by a specific alteration depends on the tumor lineage, cell epigenetic state, and presence of additional alterations. In this case, we say that an evolutionary dependency exists between an alteration and what influences its selection. Epistatic interactions between altered genes lead to evolutionary dependencies (EDs), by favoring or vetoing specific combinations of events. Large-scale cancer genomics studies have discovered examples of such dependencies, and showed that they influence tumor progression, disease phenotypes, and therapeutic response. In the past decade, several algorithmic approaches have been proposed to infer EDs from large-scale genomics datasets. These methods adopt diverse strategies to address common challenges and shed new light on cancer evolutionary trajectories. Here, we review these efforts starting from a simple conceptualization of the problem, presenting the tackled and still unmet needs in the field, and discussing the implications of EDs in cancer biology and precision oncology.

Past, Present, and Future of Genome Modification in Escherichia coli

Escherichia coli K-12 is one of the most well-studied species of bacteria. This species, however, is much more difficult to modify by homologous recombination (HR) than other model microorganisms. Research on HR in E. coli has led to a better understanding of the molecular mechanisms of HR, resulting in technical improvements and rapid progress in genome research, and allowing whole-genome mutagenesis and large-scale genome modifications. Developments using λ Red (exo, bet, and gam) and CRISPR-Cas have made E. coli as amenable to genome modification as other model microorganisms, such as Saccharomyces cerevisiae and Bacillus subtilis. This review describes the history of recombination research in E. coli, as well as improvements in techniques for genome modification by HR. This review also describes the results of large-scale genome modification of E. coli using these technologies, including DNA synthesis and assembly. In addition, this article reviews recent advances in genome modification, considers future directions, and describes problems associated with the creation of cells by design.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Transposon insertional mutagenesis of diverse yeast strains suggests coordinated gene essentiality polymorphisms

Due to epistasis, the same mutation can have drastically different phenotypic consequences in different individuals. This phenomenon is pertinent to precision medicine as well as antimicrobial drug development, but its general characteristics are largely unknown. We approach this question by genome-wide assessment of gene essentiality polymorphism in 16 Saccharomyces cerevisiae strains using transposon insertional mutagenesis. Essentiality polymorphism is observed for 9.8% of genes, most of which have had repeated essentiality switches in evolution. Genes exhibiting essentiality polymorphism lean toward having intermediate numbers of genetic and protein interactions. Gene essentiality changes tend to occur concordantly among components of the same protein complex or metabolic pathway and among a group of over 100 mitochondrial proteins, revealing molecular machines or functional modules as units of gene essentiality variation. Most essential genes tolerate transposon insertions consistently among strains in one or more coding segments, delineating nonessential regions within essential genes.

Emergence and propagation of epistasis in metabolic networks

Epistasis is often used to probe functional relationships between genes, and it plays an important role in evolution. However, we lack theory to understand how functional relationships at the molecular level translate into epistasis at the level of whole-organism phenotypes, such as fitness. Here, I derive two rules for how epistasis between mutations with small effects propagates from lower- to higher-level phenotypes in a hierarchical metabolic network with first-order kinetics and how such epistasis depends on topology. Most importantly, weak epistasis at a lower level may be distorted as it propagates to higher levels. Computational analyses show that epistasis in more realistic models likely follows similar, albeit more complex, patterns. These results suggest that pairwise inter-gene epistasis should be common, and it should generically depend on the genetic background and environment. Furthermore, the epistasis coefficients measured for high-level phenotypes may not be sufficient to fully infer the underlying functional relationships.

Idiosyncratic epistasis creates universals in mutational effects and evolutionary trajectories

Patterns of epistasis and shapes of fitness landscapes are of wide interest because of their bearings on a number of evolutionary theories. The common phenomena of slowing fitness increases during adaptations and diminishing returns from beneficial mutations are believed to reflect a concave fitness landscape and a preponderance of negative epistasis. Paradoxically, fitness decreases tend to decelerate and harm from deleterious mutations shrinks during the accumulation of random mutations-patterns thought to indicate a convex fitness landscape and a predominance of positive epistasis. Current theories cannot resolve this apparent contradiction. Here, we show that the phenotypic effect of a mutation varies substantially depending on the specific genetic background and that this idiosyncrasy in epistasis creates all of the above trends without requiring a biased distribution of epistasis. The idiosyncratic epistasis theory explains the universalities in mutational effects and evolutionary trajectories as emerging from randomness due to biological complexity.

Evolution of linkage and genome expansion in protocells: The origin of chromosomes

Chromosomes are likely to have assembled from unlinked genes in early evolution. Genetic linkage reduces the assortment load and intragenomic conflict in reproducing protocell models to the extent that chromosomes can go to fixation even if chromosomes suffer from a replicative disadvantage, relative to unlinked genes, proportional to their length. Here we numerically show that chromosomes spread within protocells even if recurrent deleterious mutations affecting replicating genes (as ribozymes) are considered. Dosage effect selects for optimal genomic composition within protocells that carries over to the genic composition of emerging chromosomes. Lacking an accurate segregation mechanism, protocells continue to benefit from the stochastic corrector principle (group selection of early replicators), but now at the chromosome level. A remarkable feature of this process is the appearance of multigene families (in optimal genic proportions) on chromosomes. An added benefit of chromosome formation is an increase in the selectively maintainable genome size (number of different genes), primarily due to the marked reduction of the assortment load. The establishment of chromosomes is under strong positive selection in protocells harboring unlinked genes. The error threshold of replication is raised to higher genome size by linkage due to the fact that deleterious mutations affecting protocells metabolism (hence fitness) show antagonistic (diminishing return) epistasis. This result strengthens the established benefit conferred by chromosomes on protocells allowing for the fixation of highly specific and efficient enzymes.

The emergence of chromosomes harboring several genes is a crucial ingredient of the major evolutionary transition from naked replicators to cells. Linkage of replicating genes reduces conflict between them and alleviates the problem of chance loss of genes upon stochastic protocell fission. The emerging organization of protocells maintaining several segregating chromosomes with balanced gene composition also allows for an increase in the number of gene types despite recurrent deleterious mutations. We suggest that this interim genomic organization enabled protocells to evolve specific and efficient enzymes and paved the way toward an accurate mechanism for chromosome segregation later in evolution.

Predominance of positive epistasis among drug resistance-associated mutations in HIV-1 protease

Drug-resistant mutations often have deleterious impacts on replication fitness, posing a fitness cost that can only be overcome by compensatory mutations. However, the role of fitness cost in the evolution of drug resistance has often been overlooked in clinical studies or in vitro selection experiments, as these observations only capture the outcome of drug selection. In this study, we systematically profile the fitness landscape of resistance-associated sites in HIV-1 protease using deep mutational scanning. We construct a mutant library covering combinations of mutations at 11 sites in HIV-1 protease, all of which are associated with resistance to protease inhibitors in clinic. Using deep sequencing, we quantify the fitness of thousands of HIV-1 protease mutants after multiple cycles of replication in human T cells. Although the majority of resistance-associated mutations have deleterious effects on viral replication, we find that epistasis among resistance-associated mutations is predominantly positive. Furthermore, our fitness data are consistent with genetic interactions inferred directly from HIV sequence data of patients. Fitness valleys formed by strong positive epistasis reduce the likelihood of reversal of drug resistance mutations. Overall, our results support the view that strong compensatory effects are involved in the emergence of clinically observed resistance mutations and provide insights to understanding fitness barriers in the evolution and reversion of drug resistance.

Decoupling gene functions from knockout effects by evolutionary analyses

Genic functions have long been confounded by pleiotropic mutational effects. To understand such genetic effects, we examine HAP4, a well-studied transcription factor in Saccharomyces cerevisiae that functions by forming a tetramer with HAP2, HAP3 and HAP5. Deletion of HAP4 results in highly pleiotropic gene expression responses, some of which are clustered in related cellular processes (clustered effects) while most are distributed randomly across diverse cellular processes (distributed effects). Strikingly, the distributed effects that account for much of HAP4 pleiotropy tend to be non-heritable in a population, suggesting they have few evolutionary consequences. Indeed, these effects are poorly conserved in closely related yeasts. We further show substantial overlaps of clustered effects, but not distributed effects, among the four genes encoding the HAP2/3/4/5 tetramer. This pattern holds for other biochemically characterized yeast protein complexes or metabolic pathways. Examination of a set of cell morphological traits of the deletion lines yields consistent results. Hence, only some deletion effects of a gene support related biochemical understandings with the rest being often pleiotropic and evolutionarily decoupled from the gene's normal functions. This study suggests a new framework for reverse genetic analysis.

SNPeffect: identifying functional roles of SNPs using metabolic networks

Genetic sources of phenotypic variation have been a focus of plant studies aimed at improving agricultural yield and understanding adaptive processes. Genome-wide association studies identify the genetic background behind a trait by examining associations between phenotypes and single-nucleotide polymorphisms (SNPs). Although such studies are common, biological interpretation of the results remains a challenge; especially due to the confounding nature of population structure and the systematic biases thus introduced. Here, we propose a complementary analysis (SNPeffect) that offers putative genotype-to-phenotype mechanistic interpretations by integrating biochemical knowledge encoded in metabolic models. SNPeffect is used to explain differential growth rate and metabolite accumulation in A. thaliana and P. trichocarpa accessions as the outcome of SNPs in enzyme-coding genes. To this end, we also constructed a genome-scale metabolic model for Populus trichocarpa, the first for a perennial woody tree. As expected, our results indicate that growth is a complex polygenic trait governed by carbon and energy partitioning. The predicted set of functional SNPs in both species are associated with experimentally characterized growth-determining genes and also suggest putative ones. Functional SNPs were found in pathways such as amino acid metabolism, nucleotide biosynthesis, and cellulose and lignin biosynthesis, in line with breeding strategies that target pathways governing carbon and energy partition.

Recent insights into the genotype-phenotype relationship from massively parallel genetic assays

With the molecular revolution in Biology, a mechanistic understanding of the genotype-phenotype relationship became possible. Recently, advances in DNA synthesis and sequencing have enabled the development of deep mutational scanning assays, capable of scoring comprehensive libraries of genotypes for fitness and a variety of phenotypes in massively parallel fashion. The resulting empirical genotype-fitness maps pave the way to predictive models, potentially accelerating our ability to anticipate the behaviour of pathogen and cancerous cell populations from sequencing data. Besides from cellular fitness, phenotypes of direct application in industry (e.g. enzyme activity) and medicine (e.g. antibody binding) can be quantified and even selected directly by these assays. This review discusses the technological basis of and recent developments in massively parallel genetics, along with the trends it is uncovering in the genotype-phenotype relationship (distribution of mutation effects, epistasis), their possible mechanistic bases and future directions for advancing towards the goal of predictive genetics.

Synchronization of stochastic expressions drives the clustering of functionally related genes

Functionally related genes tend to be chromosomally clustered in eukaryotic genomes even after the exclusion of tandem duplicates, but the biological significance of this widespread phenomenon is unclear. We propose that stochastic expression fluctuations of neighboring genes resulting from chromatin dynamics are more or less synchronized such that their expression ratio is more stable than that for unlinked genes. Consequently, chromosomal clustering could be advantageous when the expression ratio of the clustered genes needs to stay constant, for example, because of the accumulation of toxic compounds when this ratio is altered. Evidence from manipulative experiments on the yeast GAL cluster, comprising three chromosomally adjacent genes encoding enzymes catalyzing consecutive reactions in galactose catabolism, unequivocally supports this hypothesis and elucidates how disorder in one biological phenomenon-gene expression noise-could prompt the emergence of order in another-genome organization.

Flux balance analysis with or without molecular crowding fails to predict two thirds of experimentally observed epistasis in yeast

Computational predictions of double gene knockout effects by flux balance analysis (FBA) have been used to characterized genome-wide patterns of epistasis in microorganisms. However, it is unclear how in silico predictions are related to in vivo epistasis, as FBA predicted only a minority of experimentally observed genetic interactions between non-essential metabolic genes in yeast. Here, we perform a detailed comparison of yeast experimental epistasis data to predictions generated with different constraint-based metabolic modeling algorithms. The tested methods comprise standard FBA; a variant of MOMA, which was specifically designed to predict fitness effects of non-essential gene knockouts; and two alternative implementations of FBA with macro-molecular crowding, which account approximately for enzyme kinetics. The number of interactions uniquely predicted by one method is typically larger than its overlap with any alternative method. Only 20% of negative and 10% of positive interactions jointly predicted by all methods are confirmed by the experimental data; almost all unique predictions appear to be false. More than two thirds of epistatic interactions are undetectable by any of the tested methods. The low prediction accuracies indicate that the physiology of yeast double metabolic gene knockouts is dominated by processes not captured by current constraint-based analysis methods.

Patterns and Mechanisms of Diminishing Returns from Beneficial Mutations

Diminishing returns epistasis causes the benefit of the same advantageous mutation smaller in fitter genotypes and is frequently observed in experimental evolution. However, its occurrence in other contexts, environment dependence, and mechanistic basis are unclear. Here, we address these questions using 1,005 sequenced segregants generated from a yeast cross. Under each of 47 examined environments, 66-92% of tested polymorphisms exhibit diminishing returns epistasis. Surprisingly, improving environment quality also reduces the benefits of advantageous mutations even when fitness is controlled for, indicating the necessity to revise the global epistasis hypothesis. We propose that diminishing returns originates from the modular organization of life where the contribution of each functional module to fitness is determined jointly by the genotype and environment and has an upper limit, and demonstrate that our model predictions match empirical observations. These findings broaden the concept of diminishing returns epistasis, reveal its generality and potential cause, and have important evolutionary implications.

The substrate specificity of eukaryotic cytosolic chaperonin CCT

The cytosolic chaperonin CCT (chaperonin containing TCP-1) is an ATP-dependent double-ring protein machine mediating the folding of members of the eukaryotic cytoskeletal protein families. The actins and tubulins are obligate substrates of CCT because they are completely dependent on CCT activity to reach their native states. Genetic and proteomic analysis of the CCT interactome in the yeast Saccharomyces cerevisiae revealed a CCT network of approximately 300 genes and proteins involved in many fundamental biological processes. We classified network members into sets such as substrates, CCT cofactors and CCT-mediated assembly processes. Many members of the 7-bladed propeller family of proteins are commonly found tightly bound to CCT isolated from human and plant cells and yeasts. The anaphase promoting complex (APC/C) cofactor propellers, Cdh1p and Cdc20p, are also obligate substrates since they both require CCT for folding and functional activation. In vitro translation analysis in prokaryotic and eukaryotic cell extracts of a set of yeast propellers demonstrates their highly differential interactions with CCT and GroEL (another chaperonin). Individual propeller proteins have idiosyncratic interaction modes with CCT because they emerged independently with neo-functions many times throughout eukaryotic evolution. We present a toy model in which cytoskeletal protein biogenesis and folding flux through CCT couples cell growth and size control to time dependent cell cycle mechanisms.This article is part of a discussion meeting issue 'Allostery and molecular machines'.

Metabolic Determinants of Enzyme Evolution in a Genome-Scale Bacterial Metabolic Network

Different genes and proteins evolve at very different rates. To identify the factors that explain these differences is an important aspect of research in molecular evolution. One such factor is the role a protein plays in a large molecular network. Here, we analyze the evolutionary rates of enzyme-coding genes in the genome-scale metabolic network of Escherichia coli to find the evolutionary constraints imposed by the structure and function of this complex metabolic system. Central and highly connected enzymes appear to evolve more slowly than less connected enzymes, but we find that they do so as a by-product of their high abundance, and not because of their position in the metabolic network. In contrast, enzymes catalyzing reactions with high metabolic flux-high substrate to product conversion rates-evolve slowly even after we account for their abundance. Moreover, enzymes catalyzing reactions that are difficult to by-pass through alternative pathways, such that they are essential in many different genetic backgrounds, also evolve more slowly. Our analyses show that an enzyme's role in the function of a metabolic network affects its evolution more than its place in the network's structure. They highlight the value of a system-level perspective for studies of molecular evolution.

Multiplexed deactivated CRISPR-Cas9 gene expression perturbations deter bacterial adaptation by inducing negative epistasis

The ever-increasing threat of multi-drug resistant bacteria, a shrinking antibiotic pipeline, and the innate ability of microorganisms to adapt necessitates long-term strategies to slow the evolution of antibiotic resistance. Here we develop an approach, dubbed Controlled Hindrance of Adaptation of OrganismS or CHAOS, involving induction of epistasis between gene perturbations to deter adaption. We construct a combinatorial library of multiplexed, deactivated CRISPR-Cas9 devices to systematically perturb gene expression in Escherichia coli. While individual perturbations improved fitness during antibiotic exposure, multiplexed perturbations caused large fitness loss in a significant epistatic fashion. Strains exhibiting epistasis adapted significantly more slowly over three to fourteen days, and loss in adaptive potential was shown to be sustainable. Finally, we show that multiplexed peptide nucleic acids increase the antibiotic susceptibility of clinically isolated Carbapenem-resistant E. coli in an epistatic fashion. Together, these results suggest a new therapeutic strategy for restricting the evolution of antibiotic resistance.

Peter Otoupal et al. present CHAOS, an approach for preventing the development of antibiotic resistance in E. coli through CRISPR-Cas9-based perturbation of gene expression. They show that multiplexed perturbations decrease fitness of clinically-isolated Carbapenem-resistant E. coli upon antibiotic exposure.

An Ultrahigh-Dimensional Mapping Model of High-order Epistatic Networks for Complex Traits

BACKGROUND

Genetic interactions involving more than two loci have been thought to affect quantitatively inherited traits and diseases more pervasively than previously appreciated. However, the detection of such high-order interactions to chart a complete portrait of genetic architecture has not been well explored.

METHODS

We present an ultrahigh-dimensional model to systematically characterize genetic main effects and interaction effects of various orders among all possible markers in a genetic mapping or association study. The model was built on the extension of a variable selection procedure, called iFORM, derived from forward selection. The model shows its unique power to estimate the magnitudes and signs of high-order epistatic effects, in addition to those of main effects and pairwise epistatic effects.

RESULTS

The statistical properties of the model were tested and validated through simulation studies. By analyzing a real data for shoot growth in a mapping population of woody plant, mei (Prunus mume), we demonstrated the usefulness and utility of the model in practical genetic studies. The model has identified important high-order interactions that contribute to shoot growth for mei.

CONCLUSION

The model provides a tool to precisely construct genotype-phenotype maps for quantitative traits by identifying any possible high-order epistasis which is often ignored in the current genetic literature.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

Sex and genotype effects on nutrient-dependent fitness landscapes in Drosophila melanogaster

The sexes perform different reproductive roles and have evolved sometimes strikingly different phenotypes. One focal point of adaptive divergence occurs in the context of diet and metabolism, and males and females of a range of species have been shown to require different nutrients to maximize their fitness. Biochemical analyses in Drosophila melanogaster have confirmed that dimorphism in dietary requirements is associated with molecular sex differences in metabolite titres. In addition, they also showed significant within-sex genetic variation in the metabolome. To date however, it is unknown whether this metabolic variation translates into differences in reproductive fitness. The answer to this question is crucial to establish whether genetic variation is selectively neutral or indicative of constraints on sex-specific physiological adaptation and optimization. Here we assay genetic variation in consumption and metabolic fitness effects by screening male and female fitness of thirty D. melanogaster genotypes across four protein-to-carbohydrate ratios. In addition to confirming sexual dimorphism in consumption and fitness, we find significant genetic variation in male and female dietary requirements. Importantly, these differences are not explained by feeding responses and probably reflect metabolic variation that, in turn, suggests the presence of genetic constraints on metabolic dimorphism.

Dietary choices are influenced by genotype, mating status, and sex in Drosophila melanogaster

Mating causes many changes in physiology, behavior, and gene expression in a wide range of organisms. These changes are predicted to be sex specific, influenced by the divergent reproductive roles of the sexes. In female insects, mating is associated with an increase in egg production which requires high levels of nutritional input with direct consequences for the physiological needs of individual females. Consequently, females alter their nutritional acquisition in line with the physiological demands imposed by mating. Although much is known about the female mating-induced nutritional response, far less is known about changes in males. In addition, it is unknown whether variation between genotypes translates into variation in dietary behavioral responses. Here we examine mating-induced shifts in male and female dietary preferences across genotypes of Drosophila melanogaster. We find sex- and genotype-specific effects on both the quantity and quality of the chosen diet. These results contribute to our understanding of sex-specific metabolism and reveal genotypic variation that influences responses to physiological demands.

Weak Epistasis May Drive Adaptation in Recombining Bacteria

The impact of epistasis on the evolution of multi-locus traits depends on recombination. While sexually reproducing eukaryotes recombine so frequently that epistasis between polymorphisms is not considered to play a large role in short-term adaptation, many bacteria also recombine, some to the degree that their populations are described as "panmictic" or "freely recombining." However, whether this recombination is sufficient to limit the ability of selection to act on epistatic contributions to fitness is unknown. We quantify homologous recombination in five bacterial pathogens and use these parameter estimates in a multilocus model of bacterial evolution with additive and epistatic effects. We find that even for highly recombining species (e.g., Streptococcus pneumoniae or Helicobacter pylori), selection on weak interactions between distant mutations is nearly as efficient as for an asexual species, likely because homologous recombination typically transfers only short segments. However, for strong epistasis, bacterial recombination accelerates selection, with the dynamics dependent on the amount of recombination and the number of loci. Epistasis may thus play an important role in both the short- and long-term adaptive evolution of bacteria, and, unlike in eukaryotes, is not limited to strong effect sizes, closely linked loci, or other conditions that limit the impact of recombination.

Evolutionary adaptations to new environments generally reverse plastic phenotypic changes

Organismal adaptation to a new environment may start with plastic phenotypic changes followed by genetic changes, but whether the plastic changes are stepping stones to genetic adaptation is debated. Here we address this question by investigating gene expression and metabolic flux changes in the two-phase adaptation process using transcriptomic data from multiple experimental evolution studies and computational metabolic network analysis, respectively. We discover that genetic changes more frequently reverse than reinforce plastic phenotypic changes in virtually every adaptation. Metabolic network analysis reveals that, even in the presence of plasticity, organismal fitness drops after environmental shifts, but largely recovers through subsequent evolution. Such fitness trajectories explain why plastic phenotypic changes are genetically compensated rather than strengthened. In conclusion, although phenotypic plasticity may serve as an emergency response to a new environment that is necessary for survival, it does not generally facilitate genetic adaptation by bringing the organismal phenotype closer to the new optimum.

Phenotypic plasticity has been suggested to facilitate survival in new environments and subsequent adaptation. Here, the authors reanalyze transcriptomic data from experimental evolution studies in combination with computational metabolic network analysis and show that genetic adaptation tends to reverse plastic changes in order to recover fitness.

Negative Epistasis in Experimental RNA Fitness Landscapes

Mutations and their effects on fitness are a fundamental component of evolution. The effects of some mutations change in the presence of other mutations, and this is referred to as epistasis. Epistasis can occur between mutations in different genes or within the same gene. A systematic study of epistasis requires the analysis of numerous mutations and their combinations, which has recently become feasible with advancements in DNA synthesis and sequencing. Here we review the mutational effects and epistatic interactions within RNA molecules revealed by several recent high-throughput mutational studies involving two ribozymes studied in vitro, as well as a tRNA and a snoRNA studied in yeast. The data allow an analysis of the distribution of fitness effects of individual mutations as well as combinations of two or more mutations. Two different approaches to measuring epistasis in the data both reveal a predominance of negative epistasis, such that higher combinations of two or more mutations are typically lower in fitness than expected from the effect of each individual mutation. These data are in contrast to past studies of epistasis that used computationally predicted secondary structures of RNA that revealed a predominance of positive epistasis. The RNA data reviewed here are more similar to that found from mutational experiments on individual protein enzymes, suggesting that a common thermodynamic framework may explain negative epistasis between mutations within macromolecules.

Enhanced identification of synergistic and antagonistic emergent interactions among three or more drugs

Interactions among drugs play a critical role in the killing efficacy of multi-drug treatments. Recent advances in theory and experiment for three-drug interactions enable the search for emergent interactions-ones not predictable from pairwise interactions. Previous work has shown it is easier to detect synergies and antagonisms among pairwise interactions when a rescaling method is applied to the interaction metric. However, no study has carefully examined whether new types of normalization might be needed for emergence. Here, we propose several rescaling methods for enhancing the classification of the higher order drug interactions based on our conceptual framework. To choose the rescaling that best separates synergism, antagonism and additivity, we conducted bacterial growth experiments in the presence of single, pairwise and triple-drug combinations among 14 antibiotics. We found one of our rescaling methods is far better at distinguishing synergistic and antagonistic emergent interactions than any of the other methods. Using our new method, we find around 50% of emergent interactions are additive, much less than previous reports of greater than 90% additivity. We conclude that higher order emergent interactions are much more common than previously believed, and we argue these findings for drugs suggest that appropriate rescaling is crucial to infer higher order interactions.

Deciphering the Genic Basis of Yeast Fitness Variation by Simultaneous Forward and Reverse Genetics

The budding yeast Saccharomyces cerevisiae is the best studied eukaryote in molecular and cell biology, but its utility for understanding the genetic basis of phenotypic variation in natural populations is limited by inefficient association mapping due to strong and complex population structure. To overcome this challenge, we generated genome sequences for 85 strains and performed a comprehensive population genomic survey of a total of 190 diverse strains. We identified considerable variation in population structure among chromosomes and identified 181 genes that are absent from the reference genome. Many of these nonreference genes are expressed and we functionally confirmed that two of these genes confer increased resistance to antifungals. Next, we simultaneously measured the growth rates of over 4,500 laboratory strains, each of which lacks a nonessential gene, and 81 natural strains across multiple environments using unique DNA barcode present in each strain. By combining the genome-wide reverse genetic information gained from the gene deletion strains with a genome-wide association analysis from the natural strains, we identified genomic regions associated with fitness variation in natural populations. To experimentally validate a subset of these associations, we used reciprocal hemizygosity tests, finding that while the combined forward and reverse genetic approaches can identify a single causal gene, the phenotypic consequences of natural genetic variation often follow a complicated pattern. The resources and approach provided outline an efficient and reliable route to association mapping in yeast and significantly enhance its value as a model for understanding the genetic mechanisms underlying phenotypic variation and evolution in natural populations.

Influence of epistasis on response to genomic selection using complete sequence data

BACKGROUND

The effect of epistasis on response to selection is a highly debated topic. Here, we investigated the impact of epistasis on response to sequence-based selection via genomic best linear prediction (GBLUP) in a regime of strong non-symmetrical epistasis under divergent selection, using real Drosophila sequence data. We also explored the possible advantage of including epistasis in the evaluation model and/or of knowing the causal mutations.

RESULTS

Response to selection was almost exclusively due to changes in allele frequency at a few loci with a large effect. Response was highly asymmetric (about four phenotypic standard deviations higher for upward than downward selection) due to the highly skewed site frequency spectrum. Epistasis accentuated this asymmetry and affected response to selection by modulating the additive genetic variance, which was sustained for longer under upward selection whereas it eroded rapidly under downward selection. Response to selection was quite insensitive to the evaluation model, especially under an additive scenario. Nevertheless, including epistasis in the model when there was none eventually led to lower accuracies as selection proceeded. Accounting for epistasis in the model, if it existed, was beneficial but only in the medium term. There was not much gain in response if causal mutations were known, compared to using sequence data, which is likely due to strong linkage disequilibrium, high heritability and availability of phenotypes on candidates.

CONCLUSIONS

Epistatic interactions affect the response to genomic selection by modulating the additive genetic variance used for selection. Epistasis releases additive variance that may increase response to selection compared to a pure additive genetic action. Furthermore, genomic evaluation models and, in particular, GBLUP are robust, i.e. adding complexity to the model did not modify substantially the response (for a given architecture).

Systems level analysis of the Chlamydomonas reinhardtii metabolic network reveals variability in evolutionary co-conservation

Metabolic networks, which are mathematical representations of organismal metabolism, are reconstructed to provide computational platforms to guide metabolic engineering experiments and explore fundamental questions on metabolism. Systems level analyses, such as interrogation of phylogenetic relationships within the network, can provide further guidance on the modification of metabolic circuitries. Chlamydomonas reinhardtii, a biofuel relevant green alga that has retained key genes with plant, animal, and protist affinities, serves as an ideal model organism to investigate the interplay between gene function and phylogenetic affinities at multiple organizational levels. Here, using detailed topological and functional analyses, coupled with transcriptomics studies on a metabolic network that we have reconstructed for C. reinhardtii, we show that network connectivity has a significant concordance with the co-conservation of genes; however, a distinction between topological and functional relationships is observable within the network. Dynamic and static modes of co-conservation were defined and observed in a subset of gene-pairs across the network topologically. In contrast, genes with predicted synthetic interactions, or genes involved in coupled reactions, show significant enrichment for both shorter and longer phylogenetic distances. Based on our results, we propose that the metabolic network of C. reinhardtii is assembled with an architecture to minimize phylogenetic profile distances topologically, while it includes an expansion of such distances for functionally interacting genes. This arrangement may increase the robustness of C. reinhardtii's network in dealing with varied environmental challenges that the species may face. The defined evolutionary constraints within the network, which identify important pairings of genes in metabolism, may offer guidance on synthetic biology approaches to optimize the production of desirable metabolites.

Detecting epistasis with the marginal epistasis test in genetic mapping studies of quantitative traits

Epistasis, commonly defined as the interaction between multiple genes, is an important genetic component underlying phenotypic variation. Many statistical methods have been developed to model and identify epistatic interactions between genetic variants. However, because of the large combinatorial search space of interactions, most epistasis mapping methods face enormous computational challenges and often suffer from low statistical power due to multiple test correction. Here, we present a novel, alternative strategy for mapping epistasis: instead of directly identifying individual pairwise or higher-order interactions, we focus on mapping variants that have non-zero marginal epistatic effects-the combined pairwise interaction effects between a given variant and all other variants. By testing marginal epistatic effects, we can identify candidate variants that are involved in epistasis without the need to identify the exact partners with which the variants interact, thus potentially alleviating much of the statistical and computational burden associated with standard epistatic mapping procedures. Our method is based on a variance component model, and relies on a recently developed variance component estimation method for efficient parameter inference and p-value computation. We refer to our method as the "MArginal ePIstasis Test", or MAPIT. With simulations, we show how MAPIT can be used to estimate and test marginal epistatic effects, produce calibrated test statistics under the null, and facilitate the detection of pairwise epistatic interactions. We further illustrate the benefits of MAPIT in a QTL mapping study by analyzing the gene expression data of over 400 individuals from the GEUVADIS consortium.

Pervasive antagonistic interactions among hybrid incompatibility loci

Species barriers, expressed as hybrid inviability and sterility, are often due to epistatic interactions between divergent loci from two lineages. Theoretical models indicate that the strength, direction, and complexity of these genetic interactions can strongly affect the expression of interspecific reproductive isolation and the rates at which new species evolve. Nonetheless, empirical analyses have not quantified the frequency with which loci are involved in interactions affecting hybrid fitness, and whether these loci predominantly interact synergistically or antagonistically, or preferentially involve loci that have strong individual effects on hybrid fitness. We systematically examined the prevalence of interactions between pairs of short chromosomal regions from one species (Solanum habrochaites) co-introgressed into a heterospecific genetic background (Solanum lycopersicum), using lines containing pairwise combinations of 15 chromosomal segments from S. habrochaites in the background of S. lycopersicum (i.e., 95 double introgression lines). We compared the strength of hybrid incompatibility (either pollen sterility or seed sterility) expressed in each double introgression line to the expected additive effect of its two component single introgressions. We found that epistasis was common among co-introgressed regions. Interactions for hybrid dysfunction were substantially more prevalent in pollen fertility compared to seed fertility phenotypes, and were overwhelmingly antagonistic (i.e., double hybrids were less unfit than expected from additive single introgression effects). This pervasive antagonism is expected to attenuate the rate at which hybrid infertility accumulates among lineages over time (i.e., giving diminishing returns as more reproductive isolation loci accumulate), as well as decouple patterns of accumulation of sterility loci and hybrid incompatibility phenotypes. This decoupling effect might explain observed differences between pollen and seed fertility in their fit to theoretical predictions of the accumulation of isolation loci, including the ‘snowball’ effect.

A characteristic feature of new species is their inability to produce fertile or viable hybrids with other lineages. This post-zygotic reproductive isolation is caused by dysfunctional interactions between genes that have newly evolved changes in the diverging lineages. Whether these interactions occur between pairs of divergent alleles, or involve more complex networks of genes, can have strong effects on how rapidly reproductive isolation—and therefore new species—evolve. The complexity of these interactions, however, is poorly understood in empirical systems. We examined the fertility of hybrids that carried one or two chromosomal regions from a close relative, finding that hybrids with two of these heterospecific regions were frequently less sterile than would be expected from the joint fitness of hybrids that have the same regions singly. This ‘less-than-additive’ effect on hybrid sterility was widespread (observed in 20% of pairwise combinations), and especially pronounced for male sterility. We infer that genes contributing to male sterility form a more tightly connected network than previously thought, implying that reproductive isolation is evolving by incremental dysfunction of complex interactions rather than by independent pairwise incompatibilities. We use simulations to illustrate these expected patterns of accumulation of reproductive isolation when it involves highly interconnected gene networks.

HiSeeker: Detecting High-Order SNP Interactions Based on Pairwise SNP Combinations

Detecting single nucleotide polymorphisms’ (SNPs) interaction is one of the most popular approaches for explaining the missing heritability of common complex diseases in genome-wide association studies. Many methods have been proposed for SNP interaction detection, but most of them only focus on pairwise interactions and ignore high-order ones, which may also contribute to complex traits. Existing methods for high-order interaction detection can hardly handle genome-wide data and suffer from low detection power, due to the exponential growth of search space. In this paper, we proposed a flexible two-stage approach (called HiSeeker) to detect high-order interactions. In the screening stage, HiSeeker employs the chi-squared test and logistic regression model to efficiently obtain candidate pairwise combinations, which have intermediate or significant associations with the phenotype for interaction detection. In the search stage, two different strategies (exhaustive search and ant colony optimization-based search) are utilized to detect high-order interactions from candidate combinations. The experimental results on simulated datasets demonstrate that HiSeeker can more efficiently and effectively detect high-order interactions than related representative algorithms. On two real case-control datasets, HiSeeker also detects several significant high-order interactions, whose individual SNPs and pairwise interactions have no strong main effects or pairwise interaction effects, and these high-order interactions can hardly be identified by related algorithms.

When more is less: Emergent suppressive interactions in three-drug combinations

BACKGROUND

In drug-drug interactions, there are surprising cases in which the growth inhibition of bacteria by a single antibiotic decreases when a second antibiotic is added. These interactions are termed suppressive and have been argued to have the potential to limit the evolution of resistance. Nevertheless, little attention has been given to suppressive interactions because clinical studies typically search for increases in killing efficiency and because suppressive interactions are believed to be rare based on pairwise studies.

RESULTS

Here, we quantify the effects of single-, double-, and triple-drug combinations from a set of 14 antibiotics and 3 bacteria strains, totaling 364 unique three-drug combinations per bacteria strain. We find that increasing the number of drugs can increase the prevalence of suppressive interactions: 17% of three-drug combinations are suppressive compared to 5% of two-drug combinations in this study. Most cases of suppression we find (97%) are "hidden" cases for which the triple-drug bacterial growth is less than the single-drug treatments but exceeds that of a pairwise combination.

CONCLUSIONS

We find a surprising number of suppressive interactions in higher-order drug combinations. Without examining lower-order (pairwise) bacterial growth, emergent suppressive effects would be missed, potentially affecting our understanding of evolution of resistance and treatment strategies for resistant pathogens. These findings suggest that careful examination of the full factorial of drug combinations is needed to uncover suppressive interactions in higher-order combinations.

Functional characterisation of long intergenic non-coding RNAs through genetic interaction profiling in Saccharomyces cerevisiae

BACKGROUND

Transcriptome studies have revealed that many eukaryotic genomes are pervasively transcribed producing numerous long non-coding RNAs (lncRNAs). However, only a few lncRNAs have been ascribed a cellular role thus far, with most regulating the expression of adjacent genes. Even less lncRNAs have been annotated as essential hence implying that the majority may be functionally redundant. Therefore, the function of lncRNAs could be illuminated through systematic analysis of their synthetic genetic interactions (GIs).

RESULTS

Here, we employ synthetic genetic array (SGA) in Saccharomyces cerevisiae to identify GIs between long intergenic non-coding RNAs (lincRNAs) and protein-coding genes. We first validate this approach by demonstrating that the telomerase RNA TLC1 displays a GI network that corresponds to its well-described function in telomere length maintenance. We subsequently performed SGA screens on a set of uncharacterised lincRNAs and uncover their connection to diverse cellular processes. One of these lincRNAs, SUT457, exhibits a GI profile associating it to telomere organisation and we consistently demonstrate that SUT457 is required for telomeric overhang homeostasis through an Exo1-dependent pathway. Furthermore, the GI profile of SUT457 is distinct from that of its neighbouring genes suggesting a function independent to its genomic location. Accordingly, we show that ectopic expression of this lincRNA suppresses telomeric overhang accumulation in sut457Δ cells assigning a trans-acting role for SUT457 in telomere biology.

CONCLUSIONS

Overall, our work proposes that systematic application of this genetic approach could determine the functional significance of individual lncRNAs in yeast and other complex organisms.

Genetic variants in Alzheimer disease - molecular and brain network approaches

Genetic studies in late-onset Alzheimer disease (LOAD) are aimed at identifying core disease mechanisms and providing potential biomarkers and drug candidates to improve clinical care of AD. However, owing to the complexity of LOAD, including pathological heterogeneity and disease polygenicity, extraction of actionable guidance from LOAD genetics has been challenging. Past attempts to summarize the effects of LOAD-associated genetic variants have used pathway analysis and collections of small-scale experiments to hypothesize functional convergence across several variants. In this Review, we discuss how the study of molecular, cellular and brain networks provides additional information on the effects of LOAD-associated genetic variants. We then discuss emerging combinations of these omic data sets into multiscale models, which provide a more comprehensive representation of the effects of LOAD-associated genetic variants at multiple biophysical scales. Furthermore, we highlight the clinical potential of mechanistically coupling genetic variants and disease phenotypes with multiscale brain models.

Network of epistatic interactions within a yeast snoRNA

Epistatic interactions play a fundamental role in molecular evolution, but little is known about the spatial distribution of these interactions within genes. To systematically survey a model landscape of intragenic epistasis, we quantified the fitness of ~60,000 Saccharomyces cerevisiae strains expressing randomly mutated variants of the 333-nucleotide-long U3 small nucleolar RNA (snoRNA). The fitness effects of individual mutations were correlated with evolutionary conservation and structural stability. Many mutations had small individual effects but had large effects in the context of additional mutations, which indicated negative epistasis. Clusters of negative interactions were explained by local thermodynamic threshold effects, whereas positive interactions were enriched among large-effect sites and between base-paired nucleotides. We conclude that high-throughput mapping of intragenic epistasis can identify key structural and functional features of macromolecules.

Genome scale models of yeast: towards standardized evaluation and consistent omic integration

Genome scale models (GEMs) have enabled remarkable advances in systems biology, acting as functional databases of metabolism, and as scaffolds for the contextualization of high-throughput data. In the case of Saccharomyces cerevisiae (budding yeast), several GEMs have been published and are currently used for metabolic engineering and elucidating biological interactions. Here we review the history of yeast's GEMs, focusing on recent developments. We study how these models are typically evaluated, using both descriptive and predictive metrics. Additionally, we analyze the different ways in which all levels of omics data (from gene expression to flux) have been integrated in yeast GEMs. Relevant conclusions and current challenges for both GEM evaluation and omic integration are highlighted.

Adaptive Genetic Robustness of Escherichia coli Metabolic Fluxes

Genetic robustness refers to phenotypic invariance in the face of mutation and is a common characteristic of life, but its evolutionary origin is highly controversial. Genetic robustness could be an intrinsic property of biological systems, a result of direct natural selection, or a byproduct of selection for environmental robustness. To differentiate among these hypotheses, we analyze the metabolic network of Escherichia coli and comparable functional random networks. Treating the flux of each reaction as a trait and computationally predicting trait values upon mutations or environmental shifts, we discover that 1) genetic robustness is greater for the actual network than the random networks, 2) the genetic robustness of a trait increases with trait importance and this correlation is stronger in the actual network than in the random networks, and 3) the above result holds even after the control of environmental robustness. These findings demonstrate an adaptive origin of genetic robustness, consistent with the theoretical prediction that, under certain conditions, direct selection is sufficiently powerful to promote genetic robustness in cellular organisms.

The fitness landscape of a tRNA gene

Fitness landscapes describe the genotype-fitness relationship and represent major determinants of evolutionary trajectories. However, the vast genotype space, coupled with the difficulty of measuring fitness, has hindered the empirical determination of fitness landscapes. Combining precise gene replacement and next-generation sequencing, we quantified Darwinian fitness under a high-temperature challenge for more than 65,000 yeast strains, each carrying a unique variant of the single-copy tRNA(CCU)(Arg) gene at its native genomic location. Approximately 1% of single point mutations in the gene were beneficial and 42% were deleterious. Almost half of all mutation pairs exhibited statistically significant epistasis, which had a strong negative bias, except when the mutations occurred at Watson-Crick paired sites. Fitness was broadly correlated with the predicted fraction of correctly folded transfer RNA (tRNA) molecules, thereby revealing a biophysical basis of the fitness landscape.

Thoughts Toward a Theory of Natural Selection: The Importance of Microbial Experimental Evolution

Natural selection should no longer be thought of simply as a primitive (magical) concept that can be used to support all kinds of evolutionary theorizing. We need to develop causal theories of natural selection; how it arises. Because the factors contributing to the creation of natural selection are expected to be complex and intertwined, theories explaining the causes of natural selection can only be developed through the experimental method. Microbial experimental evolution provides many benefits that using other organisms does not. Microorganisms are small, so millions can be housed in a test tube; they have short generation times, so evolution over hundreds of generations can be easily studied; they can grow in chemically defined media, so the environment can be precisely defined; and they can be frozen, so the fitness of strains or populations can be directly compared across time. Microbial evolution experiments can be divided into two types. The first is to measure the selection coefficient of two known strains over the first 50 or so generations, before advantageous mutations rise to high frequency. This type of experiment can be used to directly test hypotheses. The second is to allow microbial cultures to evolve over many hundreds or thousands of generations and follow the genetic changes, to infer what phenotypes are selected. In the last section of this article, I propose that selection coefficients are not constant, but change as the population becomes fitter, introducing the idea of the selection space.

A Bayesian model for detection of high-order interactions among genetic variants in genome-wide association studies

BACKGROUND

A central question for disease studies and crop improvements is how genetics variants drive phenotypes. Genome Wide Association Study (GWAS) provides a powerful tool for characterizing the genotype-phenotype relationships in complex traits and diseases. Epistasis (gene-gene interaction), including high-order interaction among more than two genes, often plays important roles in complex traits and diseases, but current GWAS analysis usually just focuses on additive effects of single nucleotide polymorphisms (SNPs). The lack of effective computational modelling of high-order functional interactions often leads to significant under-utilization of GWAS data.

RESULTS

We have developed a novel Bayesian computational method with a Markov Chain Monte Carlo (MCMC) search, and implemented the method as a Bayesian High-order Interaction Toolkit (BHIT) for detecting epistatic interactions among SNPs. BHIT first builds a Bayesian model on both continuous data and discrete data, which is capable of detecting high-order interactions in SNPs related to case--control or quantitative phenotypes. We also developed a pipeline that enables users to apply BHIT on different species in different use cases.

CONCLUSIONS

Using both simulation data and soybean nutritional seed composition studies on oil content and protein content, BHIT effectively detected some high-order interactions associated with phenotypes, and it outperformed a number of other available tools. BHIT is freely available for academic users at http://digbio.missouri.edu/BHIT/.

Puzzling role of genetic risk factors in human longevity: "risk alleles" as pro-longevity variants

Complex diseases are major contributors to human mortality in old age. Paradoxically, many genetic variants that have been associated with increased risks of such diseases are found in genomes of long-lived people, and do not seem to compromise longevity. Here we argue that trade-off-like and conditional effects of genes can play central role in this phenomenon and in determining longevity. Such effects may occur as result of: (i) antagonistic influence of gene on the development of different health disorders; (ii) change in the effect of gene on vulnerability to death with age (especially, from “bad” to “good”); (iii) gene–gene interaction; and (iv) gene–environment interaction, among other factors. A review of current knowledge provides many examples of genetic factors that may increase the risk of one disease but reduce chances of developing another serious health condition, or improve survival from it. Factors that may increase risk of a major disease but attenuate manifestation of physical senescence are also discussed. Overall, available evidence suggests that the influence of a genetic variant on longevity may be negative, neutral or positive, depending on a delicate balance of the detrimental and beneficial effects of such variant on multiple health and aging related traits. This balance may change with age, internal and external environments, and depend on genetic surrounding. We conclude that trade-off-like and conditional genetic effects are very common and may result in situations when a disease “risk allele” can also be a pro-longevity variant, depending on context. We emphasize importance of considering such effects in both aging research and disease prevention.

Epistatic interactions among metabolic genes depend upon environmental conditions

When the effect of the state of one gene is dependent on the state of another gene in more than an additive or a neutral way, the phenomenon is termed epistasis. In particular, positive epistasis signifies that the impact of the double deletion is less severe than the neutral combination, while negative epistasis signifies that the double deletion is more severe. Epistatic interactions between genes affect the fitness landscape of an organism in its environment and are believed to be important for the evolution of sex and the evolution of recombination. Here we use large-scale computational metabolic models of microorganisms to study epistasis computationally using Flux Balance Analysis (FBA). We study what the effects of the environment are on epistatic interactions between metabolic genes in three different microorganisms: the model bacterium E. coli, the cyanobacteria Synechocystis PCC6803 and the model green algae, C. reinhardtii. Prior studies have shown that under standard laboratory conditions epistatic interactions between metabolic genes are dominated by positive epistasis. We show here that epistatic interactions depend strongly upon environmental conditions, i.e. the source of carbon, the carbon/oxygen ratio, and for photosynthetic organisms, the intensity of light. By a comparative analysis of flux distributions under different conditions, we show that whether epistatic interactions are positive or negative depends upon the topology of the carbon flow between the reactions affected by the pair of genes being considered. Thus complex metabolic networks can show epistasis even without explicit interactions between genes, and the direction and the scale of epistasis are dependent on network flows. Our results suggest that the path of evolutionary adaptation in fluctuating environments is likely to be very history dependent because of the strong effect of the environment on epistasis.

Dynamic epistasis under varying environmental perturbations

Epistasis describes the phenomenon that mutations at different loci do not have independent effects with regard to certain phenotypes. Understanding the global epistatic landscape is vital for many genetic and evolutionary theories. Current knowledge for epistatic dynamics under multiple conditions is limited by the technological difficulties in experimentally screening epistatic relations among genes. We explored this issue by applying flux balance analysis to simulate epistatic landscapes under various environmental perturbations. Specifically, we looked at gene-gene epistatic interactions, where the mutations were assumed to occur in different genes. We predicted that epistasis tends to become more positive from glucose-abundant to nutrient-limiting conditions, indicating that selection might be less effective in removing deleterious mutations in the latter. We also observed a stable core of epistatic interactions in all tested conditions, as well as many epistatic interactions unique to each condition. Interestingly, genes in the stable epistatic interaction network are directly linked to most other genes whereas genes with condition-specific epistasis form a scale-free network. Furthermore, genes with stable epistasis tend to have similar evolutionary rates, whereas this co-evolving relationship does not hold for genes with condition-specific epistasis. Our findings provide a novel genome-wide picture about epistatic dynamics under environmental perturbations.