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.     

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.     

Epistasis analysis with global transcriptional phenotypes

Classical epistasis analysis can determine the order of function of genes in pathways using morphological, biochemical and other phenotypes. It requires knowledge of the pathway's phenotypic output and a variety of experimental expertise and so is unsuitable for genome-scale analysis. Here we used microarray profiles of mutants as phenotypes for epistasis analysis. Considering genes that regulate activity of protein kinase A in Dictyostelium, we identified known and unknown epistatic relationships and reconstructed a genetic network with microarray phenotypes alone. This work shows that microarray data can provide a uniform, quantitative tool for large-scale genetic network analysis.     

Biochemical networking contributes more to genetic buffering in human and mouse metabolic pathways than does gene duplication

During evolution different genes evolve at unequal rates, reflecting the varying functional constraints on phenotype. An important contributor to this variation is genetic buffering, which reduces the potential detrimental effects of mutations. We studied whether gene duplication and redundant metabolic networks affect genetic buffering by comparing the evolutionary rate of 242 human and mouse orthologous genes involved in metabolic pathways. A gene with a redundant network is defined as one for which the structural layout of metabolic pathways provides an alternative metabolic route that can, in principle, compensate for the loss of a protein function encoded by the gene. We found that genes with redundant networks evolve at similar rates as did genes without redundant networks, [corrected] but no significant difference was detected between single-copy genes and gene families. This implies that redundancy in metabolic networks provides significantly more genetic buffering than do gene families. We also found that genes encoding proteins involved in glycolysis and gluconeogenesis showed as a group a distinct pattern of variation, in contrast to genes involved in other pathways. These results suggest that redundant networks provide genetic buffering and contribute to the functional diversification of metabolic pathways.     

Multiple knockout analysis of genetic robustness in the yeast metabolic network

Genetic robustness characterizes the constancy of the phenotype in face of heritable perturbations. Previous investigations have used comprehensive single and double gene knockouts to study gene essentiality and pairwise gene interactions in the yeast Saccharomyces cerevisiae. Here we conduct an in silico multiple knockout investigation of a flux balance analysis model of the yeast's metabolic network. Cataloging gene sets that provide mutual functional backup, we identify sets of up to eight interacting genes and characterize the 'k robustness' (the depth of backup interactions) of each gene. We find that 74% (360) of the metabolic genes participate in processes that are essential to growth in a standard laboratory environment, compared with only 13% previously found to be essential using single knockouts. The genes' k robustness is shown to be a solid indicator of their biological buffering capacity and is correlated with both the genes' environmental specificity and their evolutionary retention.     

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.     

Identification of acquired somatic mutations in the gene encoding chromatin-remodeling factor ATRX in the alpha-thalassemia myelodysplasia syndrome (ATMDS)

Inherited mutations of specific genes have elucidated the normal roles of the proteins they encode by relating specific mutations to particular phenotypes. But many potentially informative mutations in such genes are lethal early in development. Consequently, inherited mutations may not reflect all the functional roles of such proteins. Acquired, somatic defects should reflect a wider spectrum of mutations because they are not prone to negative selection in development. It has been difficult to identify such mutations so far, but microarray analysis provides a new opportunity to do so. Using this approach, we have shown that in individuals with myelodysplasia associated with alpha-thalassemia (ATMDS), somatic mutations of the gene encoding the chromatin remodeling factor ATRX cause an unexpectedly severe hematological phenotype compared with the wide spectrum of inherited mutations affecting this gene. These findings cast new light on this pleiotropic cofactor, which appears to be an essential component rather than a mere facilitator of globin gene expression.     

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.     

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.     

Robustness against mutations in genetic networks of yeast

There are two principal mechanisms that are responsible for the ability of an organism's physiological and developmental processes to compensate for mutations. In the first, genes have overlapping functions, and loss-of-function mutations in one gene will have little phenotypic effect if there are one or more additional genes with similar functions. The second mechanism has its origin in interactions between genes with unrelated functions, and has been documented in metabolic and regulatory gene networks. Here I analyse, on a genome-wide scale, which of these mechanisms of robustness against mutations is more prevalent. I used functional genomics data from the yeast Saccharomyces cerevisiae to test hypotheses related to the following: if gene duplications are mostly responsible for robustness, then a correlation is expected between the similarity of two duplicated genes and the effect of mutations in one of these genes. My results demonstrate that interactions among unrelated genes are the major cause of robustness against mutations. This type of robustness is probably an evolved response of genetic networks to stabilizing selection.     

Complex interactions between genes controlling trafficking in primary cilia

Cilia-associated human genetic disorders are striking in the diversity of their abnormalities and their complex inheritance. Inactivation of the retrograde ciliary motor by mutations in DYNC2H1 causes skeletal dysplasias that have strongly variable expressivity. Here we define previously unknown genetic relationships between Dync2h1 and other genes required for ciliary trafficking. Mutations in mouse Dync2h1 disrupt cilia structure, block Sonic hedgehog signaling and cause midgestation lethality. Heterozygosity for Ift172, a gene required for anterograde ciliary trafficking, suppresses cilia phenotypes, Sonic hedgehog signaling defects and early lethality of Dync2h1 homozygotes. Ift122, like Dync2h1, is required for retrograde ciliary trafficking, but reduction of Ift122 gene dosage also suppresses the Dync2h1 phenotype. These genetic interactions illustrate the cell biology underlying ciliopathies and argue that mutations in intraflagellar transport genes cause their phenotypes because of their roles in cilia architecture rather than direct roles in signaling.     

Common variation in three genes, including a noncoding variant in CFH, strongly influences risk of age-related macular degeneration

Age-related macular degeneration (AMD) is a common, late-onset disease with seemingly typical complexity: recurrence ratios for siblings of an affected individual are three- to sixfold higher than in the general population, and family-based analysis has resulted in only modestly significant evidence for linkage. In a case-control study drawn from a US-based population of European descent, we have identified a previously unrecognized common, noncoding variant in CFH, the gene encoding complement factor H, that substantially increases the influence of this locus on AMD, and we have strongly replicated the associations of four other previously reported common alleles in three genes (P values ranging from 10(-6) to 10(-70)). Despite excellent power to detect epistasis, we observed purely additive accumulation of risk from alleles at these genes. We found no differences in association of these loci with major phenotypic categories of advanced AMD. Genotypes at these five common SNPs define a broad spectrum of interindividual disease risk and explain about half of the classical sibling risk of AMD in our study population.     

Genome-wide association study and mouse model identify interaction between RET and EDNRB pathways in Hirschsprung disease

Genetic studies of Hirschsprung disease, a common congenital malformation, have identified eight genes with mutations that can be associated with this condition. Mutations at individual loci are, however, neither necessary nor sufficient to cause clinical disease. We conducted a genome-wide association study in 43 Mennonite family trios using 2,083 microsatellites and single-nucleotide polymorphisms and a new multipoint linkage disequilibrium method that searches for association arising from common ancestry. We identified susceptibility loci at 10q11, 13q22 and 16q23; the gene at 13q22 is EDNRB, encoding a G protein-coupled receptor (GPCR) and the gene at 10q11 is RET, encoding a receptor tyrosine kinase (RTK). Statistically significant joint transmission of RET and EDNRB alleles in affected individuals and non-complementation of aganglionosis in mouse intercrosses between Ret null and the Ednrb hypomorphic piebald allele are suggestive of epistasis between EDNRB and RET. Thus, genetic interaction between mutations in RET and EDNRB is an underlying mechanism for this complex disorder.     

Exposing the fitness contribution of duplicated genes

Duplicate genes from the whole-genome duplication (WGD) in yeast are often dispensable--removing one copy has little or no phenotypic consequence. It is unknown, however, whether such dispensability reflects insignificance of the ancestral function or compensation from paralogs. Here, using precise competition-based measurements of the fitness cost of single and double deletions, we estimate the exposed fitness contribution of WGD duplicate genes in metabolism and bound the importance of their ancestral pre-duplication function. We find that the functional overlap between paralogs sufficiently explains the apparent dispensability of individual WGD genes. Furthermore, the lower bound on the fitness value of the ancestral function, which is estimated by the degree of synergistic epistasis, is at least as large as the average fitness cost of deleting single non-WGD genes. These results suggest that most metabolic functions encoded by WGD genes are important today and were also important at the time of duplication.     

Mutations in origin recognition complex gene ORC4 cause Meier-Gorlin syndrome

Meier-Gorlin syndrome is a rare autosomal recessive genetic condition whose primary clinical hallmarks include small stature, small external ears and small or absent patellae. Using marker-assisted mapping in multiple families from a founder population and traditional coding exon sequencing of positional candidate genes, we identified three different mutations in the gene encoding ORC4, a component of the eukaryotic origin recognition complex, in five individuals with Meier-Gorlin syndrome. In two such individuals that were negative for mutations in ORC4, we found potential mutations in ORC1 and CDT1, two other genes involved in origin recognition. ORC4 is well conserved in eukaryotes, and the yeast equivalent of the human ORC4 missense mutation was shown to be pathogenic in functional assays of cell growth. This is the first report, to our knowledge, of a germline mutation in any gene of the origin recognition complex in a vertebrate organism.     

Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.     

Disruption of Bardet-Biedl syndrome ciliary proteins perturbs planar cell polarity in vertebrates

The evolutionarily conserved planar cell polarity (PCP) pathway (or noncanonical Wnt pathway) drives several important cellular processes, including epithelial cell polarization, cell migration and mitotic spindle orientation. In vertebrates, PCP genes have a vital role in polarized convergent extension movements during gastrulation and neurulation. Here we show that mice with mutations in genes involved in Bardet-Biedl syndrome (BBS), a disorder associated with ciliary dysfunction, share phenotypes with PCP mutants including open eyelids, neural tube defects and disrupted cochlear stereociliary bundles. Furthermore, we identify genetic interactions between BBS genes and a PCP gene in both mouse (Ltap, also called Vangl2) and zebrafish (vangl2). In zebrafish, the augmented phenotype results from enhanced defective convergent extension movements. We also show that Vangl2 localizes to the basal body and axoneme of ciliated cells, a pattern reminiscent of that of the BBS proteins. These data suggest that cilia are intrinsically involved in PCP processes.     

Restriction enzyme-generated siRNA (REGS) vectors and libraries

Small interfering RNA (siRNA) technology facilitates the study of loss of gene function in mammalian cells and animal models, but generating multiple siRNA vectors using oligonucleotides is slow, inefficient and costly. Here we describe a new, enzyme-mediated method for generating numerous functional siRNA constructs from any gene of interest or pool of genes. To test our restriction enzyme-generated siRNA (REGS) system, we silenced a transgene and two endogenous genes and obtained the predicted phenotypes. REGS generated on average 34 unique siRNAs per kilobase of sequence. REGS enabled us to create enzymatically a complex siRNA library (>4 x 10(5) clones) from double-stranded cDNA encompassing known and unknown genes with 96% of the clones containing inserts of the appropriate size.