Distributions of epistasis in microbes fit predictions from a fitness landscape model

How do the fitness effects of several mutations combine? Despite its simplicity, this question is central to the understanding of multilocus evolution. Epistasis (the interaction between alleles at different loci), especially epistasis for fitness traits such as reproduction and survival, influences evolutionary predictions "almost whenever multilocus genetics matters". Yet very few models have sought to predict epistasis, and none has been empirically tested. Here we show that the distribution of epistasis can be predicted from the distribution of single mutation effects, based on a simple fitness landscape model. We show that this prediction closely matches the empirical measures of epistasis that have been obtained for Escherichia coli and the RNA virus vesicular stomatitis virus. Our results suggest that a simple fitness landscape model may be sufficient to quantitatively capture the complex nature of gene interactions. This model may offer a simple and widely applicable alternative to complex metabolic network models, in particular for making evolutionary predictions.     

Pervasive joint influence of epistasis and plasticity on mutational effects in Escherichia coli

The effects of mutations on phenotype and fitness may depend on the environment (phenotypic plasticity), other mutations (genetic epistasis) or both. Here we examine the fitness effects of 18 random insertion mutations in E. coli in two resource environments and five genetic backgrounds. We tested each mutation for plasticity and epistasis by comparing its fitness effects across these ecological and genetic contexts. Some mutations had no measurable effect in any of these contexts. None of the mutations had effects on phenotypic plasticity that were independent of genetic background. However, half the mutations had epistatic interactions such that their effects differed among genetic backgrounds, usually in an environment-dependent manner. Also, the pattern of mutational effects across backgrounds indicated that epistasis had been shaped primarily by unique events in the evolutionary history of a population rather than by repeatable events associated with shared environmental history.     

Modular epistasis in yeast metabolism

Epistatic interactions, manifested in the effects of mutations on the phenotypes caused by other mutations, may help uncover the functional organization of complex biological networks. Here, we studied system-level epistatic interactions by computing growth phenotypes of all single and double knockouts of 890 metabolic genes in Saccharomyces cerevisiae, using the framework of flux balance analysis. A new scale for epistasis identified a distinctive trimodal distribution of these epistatic effects, allowing gene pairs to be classified as buffering, aggravating or noninteracting. We found that the ensuing epistatic interaction network could be organized hierarchically into function-enriched modules that interact with each other 'monochromatically' (i.e., with purely aggravating or purely buffering epistatic links). This property extends the concept of epistasis from single genes to functional units and provides a new definition of biological modularity, which emphasizes interactions between, rather than within, functional modules. Our approach can be used to infer functional gene modules from purely phenotypic epistasis measurements.     

Mechanisms and convergence of compensatory evolution in mammalian mitochondrial tRNAs

The function of protein and RNA molecules depends on complex epistatic interactions between sites. Therefore, the deleterious effect of a mutation can be suppressed by a compensatory second-site substitution. In relating a list of 86 pathogenic mutations in human tRNAs encoded by mitochondrial genes to the sequences of their mammalian orthologs, we noted that 52 pathogenic mutations were present in normal tRNAs of one or several nonhuman mammals. We found at least five mechanisms of compensation for 32 pathogenic mutations that destroyed a Watson-Crick pair in one of the four tRNA stems: restoration of the affected Watson-Crick interaction (25 cases), strengthening of another pair (4 cases), creation of a new pair (8 cases), changes of multiple interactions in the affected stem (11 cases) and changes involving the interaction between the loop and stem structures (3 cases). A pathogenic mutation and its compensating substitution are fixed in a lineage in rapid succession, and often a compensatory interaction evolves convergently in different clades. At least 10%, and perhaps as many as 50%, of all nucleotide substitutions in evolving mammalian tRNAs participate in such interactions, indicating that the evolution of tRNAs proceeds along highly epistatic fitness ridges.     

Innate partnership of HLA-B and KIR3DL1 subtypes against HIV-1

Allotypes of the natural killer (NK) cell receptor KIR3DL1 vary in both NK cell expression patterns and inhibitory capacity upon binding to their ligands, HLA-B Bw4 molecules, present on target cells. Using a sample size of over 1,500 human immunodeficiency virus (HIV)+ individuals, we show that various distinct allelic combinations of the KIR3DL1 and HLA-B loci significantly and strongly influence both AIDS progression and plasma HIV RNA abundance in a consistent manner. These genetic data correlate very well with previously defined functional differences that distinguish KIR3DL1 allotypes. The various epistatic effects observed here for common, distinct KIR3DL1 and HLA-B Bw4 combinations are unprecedented with regard to any pair of genetic loci in human disease, and indicate that NK cells may have a critical role in the natural history of HIV infection.     

Systematic pathway analysis using high-resolution fitness profiling of combinatorial gene deletions

Systematic genetic interaction studies have illuminated many cellular processes. Here we quantitatively examine genetic interactions among 26 Saccharomyces cerevisiae genes conferring resistance to the DNA-damaging agent methyl methanesulfonate (MMS), as determined by chemogenomic fitness profiling of pooled deletion strains. We constructed 650 double-deletion strains, corresponding to all pairings of these 26 deletions. The fitness of single- and double-deletion strains were measured in the presence and absence of MMS. Genetic interactions were defined by combining principles from both statistical and classical genetics. The resulting network predicts that the Mph1 helicase has a role in resolving homologous recombination-derived DNA intermediates that is similar to (but distinct from) that of the Sgs1 helicase. Our results emphasize the utility of small molecules and multifactorial deletion mutants in uncovering functional relationships and pathway order.     

Functional classification of drugs by properties of their pairwise interactions

Multidrug treatments are increasingly important in medicine and for probing biological systems. Although many studies have focused on interactions between specific drugs, little is known about the system properties of a full drug interaction network. Like their genetic counterparts, two drugs may have no interaction, or they may interact synergistically or antagonistically to increase or suppress their individual effects. Here we use a sensitive bioluminescence technique to provide quantitative measurements of pairwise interactions among 21 antibiotics that affect growth rate in Escherichia coli. We find that the drug interaction network possesses a special property: it can be separated into classes of drugs such that any two classes interact either purely synergistically or purely antagonistically. These classes correspond directly to the cellular functions affected by the drugs. This network approach provides a new conceptual framework for understanding the functional mechanisms of drugs and their cellular targets and can be applied in systems intractable to mutant screening, biochemistry or microscopy.     

Epistatic interaction between KIR3DS1 and HLA-B delays the progression to AIDS

Natural killer (NK) cells provide defense in the early stages of the innate immune response against viral infections by producing cytokines and causing cytotoxicity. The killer immunoglobulin-like receptors (KIRs) on NK cells regulate the inhibition and activation of NK-cell responses through recognition of human leukocyte antigen (HLA) class I molecules on target cells KIR and HLA loci are both highly polymorphic, and some HLA class I products bind and trigger cell-surface receptors specified by KIR genes. Here we report that the activating KIR allele KIR3DS1, in combination with HLA-B alleles that encode molecules with isoleucine at position 80 (HLA-B Bw4-80Ile), is associated with delayed progression to AIDS in individuals infected with human immunodeficiency virus type 1 (HIV-1). In the absence of KIR3DS1, the HLA-B Bw4-80Ile allele was not associated with any of the AIDS outcomes measured. By contrast, in the absence of HLA-B Bw4-80Ile alleles, KIR3DS1 was significantly associated with more rapid progression to AIDS. These observations are strongly suggestive of a model involving an epistatic interaction between the two loci. The strongest synergistic effect of these loci was on progression to depletion of CD4(+) T cells, which suggests that a protective response of NK cells involving KIR3DS1 and its HLA class I ligands begins soon after HIV-1 infection.     

Epistatic buffering of fitness loss in yeast double deletion strains

Interactions between deleterious mutations have been insufficiently studied, despite the fact that their strength and direction are critical for understanding the evolution of genetic recombination and the buildup of mutational load in populations. We compiled a list of 758 yeast gene deletions causing growth defects (from the Munich Information Center for Protein Sequences database and ref. 7). Using BY4741 and BY4742 single-deletion strains, we carried out 639 random crosses and assayed growth curves of the resulting progeny. We show that the maximum growth rate averaged over strains lacking deletions and those with double deletions is higher than that of strains with single deletions, indicating a positive epistatic effect. This tendency is shared by genes belonging to a variety of functional classes. Based on our data and former theoretical work, we suggest that epistasis is likely to diminish the negative effects of mutations when the ability to produce biomass at high rates contributes significantly to fitness.     

Distribution of fitness effects among beneficial mutations before selection in experimental populations of bacteria

The extent to which a population diverges from its ancestor through adaptive evolution depends on variation supplied by novel beneficial mutations. Extending earlier work, recent theory makes two predictions that seem to be robust to biological details: the distribution of fitness effects among beneficial mutations before selection should be (i) exponential and (ii) invariant, meaning it is always exponential regardless of the fitness rank of the wild-type allele. Here we test these predictions by assaying the fitness of 665 independently derived single-step mutations in the bacterium Pseudomonas fluorescens across a range of environments. We show that the distribution of fitness effects among beneficial mutations is indistinguishable from an exponential despite marked variation in the fitness rank of the wild type across environments. These results suggest that the initial step in adaptive evolution--the production of novel beneficial mutants from which selection sorts--is very general, being characterized by an approximately exponential distribution with many mutations of small effect and few of large effect. We also document substantial variation in the pleiotropic costs of antibiotic resistance, a result that may have implications for strategies aimed at eliminating resistant pathogens in animal and human populations.     

Perspectives en biologie moléculaire appliquée au suivi thérapeutique des personnes infectées par le VIH

The determination of plasma HIV-1 RNA concentration in infected individuals is now a key part of the assessment of the stage of the disease and the effectiveness of treatment. CD4 cell counts provide independent information, but many studies have demonstrated the prognostic value of baseline HIV-1 RNA levels and reductions in HIV-1 RNA in response to antiretroviral therapy. Moreover, several lines of evidence suggest that a long-lasting response to treatment is associated with HIV-RNA levels being reduced to levels undetectable by very sensitive assays within 6 months of initiating treatment. Recent progress in monitoring plasma RNA levels includes (i) the use of new primers for detecting and quantifying non-clade B HIV-1 subtypes with analytical performances similar to those of the clade B subtype, (ii) the development of more sensitive procedures and (iii) the automating of amplification and detection steps. Other molecular forms of HIV, such as cell-associated RNA and proviral DNA, can be determined in patients with plasma RNA levels below the detection limit who have been given active antiretroviral therapy. Finally, several PCR-based assays are now available for detecting viral mutations conferring drug resistance. Nevertheless the clinical value of these techniques for guiding clinical decisions requires further evaluation.     

Intérêt de la biologie moléculaire dans le diagnostic des infections vlrales du système nerveux central

The use of PCR and hybridization for the detection of viral nucleic acids in cerebrospinal fluid (CSF) currently represents the most reliable method for the meningoencephalitis diagnosis due to herpes viruses, enteroviruses, JC polyomavirus and HIV-1 primary infection. The molecular tools are less suitable for the postinfectious encephalitis associated with a systemic viral infection such as neurologic complications of measles, rubella, influenza and varicella, The speed and high sensitivity allow earlier diagnosis and specific treatment. The monitoring of the antiviral therapy is assessed by the development of quantitative amplification assays in CSF or by a negative PCR result. The pattern mutations conferring resistance and the genetic variability of viruses (HIV-1, Dengue virus) could be determined from different locations in the genome. Theses methods performed over the past several years as research tools have now a large scale application with the commercialization of some tests. But before becoming the first line diagnostic test, PCR of CSF could be conducted with control quality and evaluated protocols for avoiding false positive and negative results.     

Pleiotropic fitness effects of the Tre1-Gr5a region in Drosophila melanogaster

The abundance of transposable elements and DNA repeat sequences in mammalian genomes raises the question of whether such insertions represent passive evolutionary baggage or may influence the expression of complex traits. We addressed this question in Drosophila melanogaster, in which the effects of single transposable elements on complex traits can be assessed in genetically identical individuals reared in controlled environments. Here we demonstrate that single P-element insertions in the intergenic region between the gustatory receptor 5a (Gr5a, also known as Tre) and trapped in endoderm 1 (Tre1), which encodes an orphan receptor, exert complex pleiotropic effects on fitness traits, including selective nutrient intake, life span, and resistance to starvation and heat stress. Mutations in this region interact epistatically with downstream components of the insulin signaling pathway. Transposon-induced sex-specific and sex-antagonistic effects further accentuate the complex influences that intergenic transposable elements can contribute to quantitative trait phenotypes.     

Life extension through neurofibromin mitochondrial regulation and antioxidant therapy for neurofibromatosis-1 in Drosophila melanogaster

We investigated the pathophysiology of neurofibromatosis-1 (NF1) in Drosophila melanogaster by inactivation or overexpression of the NF1 gene. NF1 gene mutants had shortened life spans and increased vulnerability to heat and oxidative stress in association with reduced mitochondrial respiration and elevated reactive oxygen species (ROS) production. Flies overexpressing NF1 had increased life spans, improved reproductive fitness, increased resistance to oxidative and heat stress in association with increased mitochondrial respiration and a 60% reduction in ROS production. These phenotypic effects proved to be modulated by the adenylyl cyclase/cyclic AMP (cAMP)/protein kinase A pathway, not the Ras/Raf pathway. Treatment of wild-type D. melanogaster with cAMP analogs increased their life span, and treatment of NF1 mutants with metalloporphyrin catalytic antioxidant compounds restored their life span. Thus, neurofibromin regulates longevity and stress resistance through cAMP regulation of mitochondrial respiration and ROS production, and NF1 may be treatable using catalytic antioxidants.     

The 'evolvability' of promiscuous protein functions

How proteins with new functions (e.g., drug or antibiotic resistance or degradation of man-made chemicals) evolve in a matter of months or years is still unclear. This ability is dependent on the induction of new phenotypic traits by a small number of mutations (plasticity). But mutations often have deleterious effects on functions that are essential for survival. How are these seemingly conflicting demands met at the single-protein level? Results from directed laboratory evolution experiments indicate that the evolution of a new function is driven by mutations that have little effect on the native function but large effects on the promiscuous functions that serve as starting point. Thus, an evolving protein can initially acquire increased fitness for a new function without losing its original function. Gene duplication and the divergence of a completely new protein may then follow.     

Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes

Epidemics of drug-resistant bacteria emerge worldwide, even as resistant strains frequently have reduced fitness compared to their drug-susceptible counterparts. Data from model systems suggest that the fitness cost of antimicrobial resistance can be reduced by compensatory mutations; however, there is limited evidence that compensatory evolution has any significant role in the success of drug-resistant bacteria in human populations. Here we describe a set of compensatory mutations in the RNA polymerase genes of rifampicin-resistant M. tuberculosis, the etiologic agent of human tuberculosis (TB). M. tuberculosis strains harboring these compensatory mutations showed a high competitive fitness in vitro. Moreover, these mutations were associated with high fitness in vivo, as determined by examining their relative clinical frequency across patient populations. Of note, in countries with the world's highest incidence of multidrug-resistant (MDR) TB, more than 30% of MDR clinical isolates had this form of mutation. Our findings support a role for compensatory evolution in the global epidemics of MDR TB.     

Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is the most common primary liver malignancy. Here, we performed high-resolution copy-number analysis on 125 HCC tumors and whole-exome sequencing on 24 of these tumors. We identified 135 homozygous deletions and 994 somatic mutations of genes with predicted functional consequences. We found new recurrent alterations in four genes (ARID1A, RPS6KA3, NFE2L2 and IRF2) not previously described in HCC. Functional analyses showed tumor suppressor properties for IRF2, whose inactivation, exclusively found in hepatitis B virus (HBV)-related tumors, led to impaired TP53 function. In contrast, inactivation of chromatin remodelers was frequent and predominant in alcohol-related tumors. Moreover, association of mutations in specific genes (RPS6KA3-AXIN1 and NFE2L2-CTNNB1) suggested that Wnt/β-catenin signaling might cooperate in liver carcinogenesis with both oxidative stress metabolism and Ras/mitogen-activated protein kinase (MAPK) pathways. This study provides insight into the somatic mutational landscape in HCC and identifies interactions between mutations in oncogene and tumor suppressor gene mutations related to specific risk factors.     

Trans-dominant inhibition of RNA viral replication can slow growth of drug-resistant viruses

The high error rates of viral RNA-dependent RNA polymerases create heterogeneous viral populations whose disparate RNA genomes affect each other's survival. We systematically screened the poliovirus genome and identified four sets of dominant mutations. Mutated alleles in capsid- and polymerase-coding regions resulted in dominant negative phenotypes, probably due to the proteins' oligomeric properties. We also identified dominant mutations in an RNA element required for priming RNA synthesis (CRE) and in the protein primer (VPg), suggesting that nonproductive priming intermediates are inhibitory. Mutations that inhibit the activity of viral proteinase 2A were dominant, arguing that inhibition of its known intramolecular activity creates a toxic product. Viral products that, when defective, dominantly interfere with growth of nondefective viruses will probably be excellent drug targets because drug-sensitive viruses should be dominant over drug-resistant variants. Accordingly, a virus sensitive to anticapsid compound WIN51711 dominantly inhibited the intracellular growth of a drug-resistant virus. Therefore, dominant inhibitor screening should validate or predict targets for antiviral therapy with reduced risk for drug resistance.