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.     

Whole-genome sequencing of multiple Arabidopsis thaliana populations

The plant Arabidopsis thaliana occurs naturally in many different habitats throughout Eurasia. As a foundation for identifying genetic variation contributing to adaptation to diverse environments, a 1001 Genomes Project to sequence geographically diverse A. thaliana strains has been initiated. Here we present the first phase of this project, based on population-scale sequencing of 80 strains drawn from eight regions throughout the species' native range. We describe the majority of common small-scale polymorphisms as well as many larger insertions and deletions in the A. thaliana pan-genome, their effects on gene function, and the patterns of local and global linkage among these variants. The action of processes other than spontaneous mutation is identified by comparing the spectrum of mutations that have accumulated since A. thaliana diverged from its closest relative 10 million years ago with the spectrum observed in the laboratory. Recent species-wide selective sweeps are rare, and potentially deleterious mutations are more common in marginal populations.     

Germline mutations in DIS3L2 cause the Perlman syndrome of overgrowth and Wilms tumor susceptibility

Perlman syndrome is a congenital overgrowth syndrome inherited in an autosomal recessive manner that is associated with Wilms tumor susceptibility. We mapped a previously unknown susceptibility locus to 2q37.1 and identified germline mutations in DIS3L2, a homolog of the Schizosaccharomyces pombe dis3 gene, in individuals with Perlman syndrome. Yeast dis3 mutant strains have mitotic abnormalities. Yeast Dis3 and its human homologs, DIS3 and DIS3L1, have exoribonuclease activity and bind to the core RNA exosome complex. DIS3L2 has a different intracellular localization and lacks the PIN domain found in DIS3 and DIS3L1; nevertheless, we show that DIS3L2 has exonuclease activity. DIS3L2 inactivation was associated with mitotic abnormalities and altered expression of mitotic checkpoint proteins. DIS3L2 overexpression suppressed the growth of human cancer cell lines, and knockdown enhanced the growth of these cells. We also detected evidence of DIS3L2 mutations in sporadic Wilms tumor. These observations suggest that DIS3L2 has a critical role in RNA metabolism and is essential for the regulation of cell growth and division.     

Metabolic gene-deletion strains of Escherichia coli evolve to computationally predicted growth phenotypes

Genome-scale metabolic models have a promising ability to describe cellular phenotypes accurately. Here we show that strains of Escherichia coli carrying a deletion of a single metabolic gene increase their growth rates (by 87% on average) during adaptive evolution and that the endpoint growth rates can be predicted computationally in 39 of 50 (78%) strains tested. These results show that computational models can be used to predict the eventual effects of genetic modifications.     

Inherited variants of MYH associated with somatic G:C-->T:A mutations in colorectal tumors

Inherited defects of base excision repair have not been associated with any human genetic disorder, although mutations of the genes mutM and mutY, which function in Escherichia coli base excision repair, lead to increased transversions of G:C to T:A. We have studied family N, which is affected with multiple colorectal adenomas and carcinoma but lacks an inherited mutation of the adenomatous polyposis coli gene (APC) that is associated with familial adenomatous polyposis. Here we show that 11 tumors from 3 affected siblings contain 18 somatic inactivating mutations of APC and that 15 of these mutations are G:C-->A transversions--a significantly greater proportion than is found in sporadic tumors or in tumors associated with familial adenomatous polyposis. Analysis of the human homolog of mutY, MYH, showed that the siblings were compound heterozygotes for the nonconservative missense variants Tyr165Cys and Gly382Asp. These mutations affect residues that are conserved in mutY of E. coli (Tyr82 and Gly253). Tyrosine 82 is located in the pseudo-helix-hairpin-helix (HhH) motif and is predicted to function in mismatch specificity. Assays of adenine glycosylase activity of the Tyr82Cys and Gly253Asp mutant proteins with 8-oxoG:A and G:A substrates show that their activity is reduced significantly. Our findings link the inherited variants in MYH to the pattern of somatic APC mutation in family N and implicate defective base excision repair in predisposition to tumors in humans.     

Evidence that positive selection drives Y-chromosome degeneration in Drosophila miranda

Why does the Y chromosome harbor so few functional loci? Evolutionary theory predicts that Y chromosomes degenerate because they lack genetic recombination. Both positive and negative selection models have been invoked to explain this degeneration, as both can result in the recurrent fixation of linked deleterious mutations on a nonrecombining Y chromosome. To distinguish between these models, I investigated patterns of nucleotide variability along 37 kb of the recently formed neo-Y chromosome in Drosophila miranda. Levels of nucleotide variability on this chromosome are 30 times lower than in highly recombining portions of the genome. Both positive and negative selection models can result in reduced variability levels, but their effects on the frequency spectrum of mutations differ. Using coalescent simulations, I show that the patterns of nucleotide variability on the neo-Y chromosome are unlikely under deleterious mutation models (including background selection and Muller's ratchet) but are expected under recent positive selection. These results implicate positive selection as an important force driving the degeneration of Y chromosomes; adaptation at a few loci, possibly increasing male fitness, occurs at the cost of most other genes on this chromosome.     

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.     

AXIN1 mutations in hepatocellular carcinomas, and growth suppression in cancer cells by virus-mediated transfer of AXIN1

The Wnt signaling pathway is essential for development and organogenesis. Wnt signaling stabilizes beta-catenin, which accumulates in the cytoplasm, binds to 1-cell factor (TCF; also known as lymphocyte enhancer-binding factor, LEF) and then upregulates downstream genes. Mutations in CTNNB1 (encoding beta-catenin) or APC (adenomatous polyposis coli) have been reported in human neoplasms including colon cancers and hepatocellular carcinomas (HCCs). Because HCC5 tend to show accumulation of beta-catenin more often than mutations in CTNNB1, we looked for mutations in AXIN1, encoding a key factor for Wnt signaling, in 6 HCC cell lines and 100 primary HCC5. Among the 4 cell lines and 87 HCC5 in which we did not detect CTNNB1 mutations, we identified AXIN1 mutations in 3 cell lines and 6 mutations in 5 of the primary HCCs. In cell lines containing mutations in either gene, we observed increased DNA binding of TCF associated with beta-catenin in nuclei. Adenovirus mediated gene transfer of wild-type AXINI induced apoptosis in hepatocellular and colorectal cancer cells that had accumulated beta-catenin as a consequence of either APC, CTNNB1 or AXIN1 mutation, suggesting that axin may be an effective therapeutic molecule for suppressing growth of hepatocellular and colorectal cancers.     

Predicting phenotypic variation in yeast from individual genome sequences

A central challenge in genetics is to predict phenotypic variation from individual genome sequences. Here we construct and evaluate phenotypic predictions for 19 strains of Saccharomyces cerevisiae. We use conservation-based methods to predict the impact of protein-coding variation within genes on protein function. We then rank strains using a prediction score that measures the total sum of function-altering changes in different sets of genes reported to influence over 100 phenotypes in genome-wide loss-of-function screens. We evaluate our predictions by comparing them with the observed growth rate and efficiency of 15 strains tested across 20 conditions in quantitative experiments. The median predictive performance, as measured by ROC AUC, was 0.76, and predictions were more accurate when the genes reported to influence a trait were highly connected in a functional gene network.     

Use of whole genome sequencing to estimate the mutation rate of Mycobacterium tuberculosis during latent infection

Tuberculosis poses a global health emergency, which has been compounded by the emergence of drug-resistant Mycobacterium tuberculosis (Mtb) strains. We used whole-genome sequencing to compare the accumulation of mutations in Mtb isolated from cynomolgus macaques with active, latent or reactivated disease. We sequenced 33 Mtb isolates from nine macaques with an average genome coverage of 93% and an average read depth of 117×. Based on the distribution of SNPs observed, we calculated the mutation rates for these disease states. We found a similar mutation rate during latency as during active disease or in a logarithmically growing culture over the same period of time. The pattern of polymorphisms suggests that the mutational burden in vivo is because of oxidative DNA damage. We show that Mtb continues to acquire mutations during disease latency, which may explain why isoniazid monotherapy for latent tuberculosis is a risk factor for the emergence of isoniazid resistance.     

Natural selection has driven population differentiation in modern humans

The considerable range of observed phenotypic variation in human populations may reflect, in part, distinctive processes of natural selection and adaptation to variable environmental conditions. Although recent genome-wide studies have identified candidate regions under selection, it is not yet clear how natural selection has shaped population differentiation. Here, we have analyzed the degree of population differentiation at 2.8 million Phase II HapMap single-nucleotide polymorphisms. We find that negative selection has globally reduced population differentiation at amino acid-altering mutations, particularly in disease-related genes. Conversely, positive selection has ensured the regional adaptation of human populations by increasing population differentiation in gene regions, primarily at nonsynonymous and 5'-UTR variants. Our analyses identify a fraction of loci that have contributed, and probably still contribute, to the morphological and disease-related phenotypic diversity of current human populations.     

Mutations in ORC1, encoding the largest subunit of the origin recognition complex, cause microcephalic primordial dwarfism resembling Meier-Gorlin syndrome

Studies into disorders of extreme growth failure (for example, Seckel syndrome and Majewski osteodysplastic primordial dwarfism type II) have implicated fundamental cellular processes of DNA damage response signaling and centrosome function in the regulation of human growth. Here we report that mutations in ORC1, encoding a subunit of the origin recognition complex, cause microcephalic primordial dwarfism resembling Meier-Gorlin syndrome. We establish that these mutations disrupt known ORC1 functions including pre-replicative complex formation and origin activation. ORC1 deficiency perturbs S-phase entry and S-phase progression. Additionally, we show that Orc1 depletion in zebrafish is sufficient to markedly reduce body size during rapid embryonic growth. Our data suggest a model in which ORC1 mutations impair replication licensing, slowing cell cycle progression and consequently impeding growth during development, particularly at times of rapid proliferation. These findings establish a novel mechanism for the pathogenesis of microcephalic dwarfism and show a surprising but important developmental impact of impaired origin licensing.     

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.     

Apc modulates embryonic stem-cell differentiation by controlling the dosage of beta-catenin signaling

The Wnt signal-transduction pathway induces the nuclear translocation of membrane-bound beta-catenin (Catnb) and has a key role in cell-fate determination. Tight somatic regulation of this signal is essential, as uncontrolled nuclear accumulation of beta-catenin can cause developmental defects and tumorigenesis in the adult organism. The adenomatous polyposis coli gene (APC) is a major controller of the Wnt pathway and is essential to prevent tumorigenesis in a variety of tissues and organs. Here, we have investigated the effect of different mutations in Apc on the differentiation potential of mouse embryonic stem (ES) cells. We provide genetic and molecular evidence that the ability and sensitivity of ES cells to differentiate into the three germ layers is inhibited by increased doses of beta-catenin by specific Apc mutations. These range from a severe differentiation blockade in Apc alleles completely deficient in beta-catenin regulation to more specific neuroectodermal, dorsal mesodermal and endodermal defects in more hypomorphic alleles. Accordingly, a targeted oncogenic mutation in Catnb also affects the differentiation potential of ES cells. Expression profiling of wildtype and Apc-mutated teratomas supports the differentiation defects at the molecular level and pinpoints a large number of downstream structural and regulating genes. Chimeric experiments showed that this effect is cell-autonomous. Our results imply that constitutive activation of the Apc/beta-catenin signaling pathway results in differentiation defects in tissue homeostasis, and possibly underlies tumorigenesis in the colon and other self-renewing tissues.     

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.     

X-linked Cornelia de Lange syndrome owing to SMC1L1 mutations

Cornelia de Lange syndrome is a multisystem developmental disorder characterized by facial dysmorphisms, upper limb abnormalities, growth delay and cognitive retardation. Mutations in the NIPBL gene, a component of the cohesin complex, account for approximately half of the affected individuals. We report here that mutations in SMC1L1 (also known as SMC1), which encodes a different subunit of the cohesin complex, are responsible for CdLS in three male members of an affected family and in one sporadic case.     

Mutations in the CCN gene family member WISP3 cause progressive pseudorheumatoid dysplasia.

Members of the CCN (for CTGF, cyr61/cef10, nov) gene family encode cysteine-rich secreted proteins with roles in cell growth and differentiation. Cell-specific and tissue-specific differences in the expression and function of different CCN family members suggest they have non-redundant roles. Using a positional-candidate approach, we found that mutations in the CCN family member WISP3 are associated with the autosomal recessive skeletal disorder progressive pseudorheumatoid dysplasia (PPD; MIM 208230). PPD is an autosomal recessive disorder that may be initially misdiagnosed as juvenile rheumatoid arthritis. Its population incidence has been estimated at 1 per million in the United Kingdom, but it is likely to be higher in the Middle East and Gulf States. Affected individuals are asymptomatic in early childhood. Signs and symptoms of disease typically develop between three and eight years of age. Clinically and radiographically, patients experience continued cartilage loss and destructive bone changes as they age, in several instances necessitating joint replacement surgery by the third decade of life. Extraskeletal manifestations have not been reported in PPD. Cartilage appears to be the primary affected tissue, and in one patient, a biopsy of the iliac crest revealed abnormal nests of chondrocytes and loss of normal cell columnar organization in growth zones. We have identified nine different WISP3 mutations in unrelated, affected individuals, indicating that the gene is essential for normal post-natal skeletal growth and cartilage homeostasis.     

Indian Hedgehog is an antagonist of Wnt signaling in colonic epithelial cell differentiation

Wnt signaling defines the colonic epithelial progenitor cell phenotype, and mutations in the gene adenomatous polyposis coli (APC) that activate the Wnt pathway cause the familial adenomatous polyposis coli (FAP) syndrome and most sporadic colon cancers. The mechanisms that regulate the transition of epithelial precursor cells into their differentiated derivatives are poorly characterized. We report that Indian hedgehog (Ihh) is expressed by mature colonocytes and regulates their differentiation in vitro and in vivo. Hedgehog (Hh) signaling restricts the expression of Wnt targets to the base of the colonic crypt in vivo, and transfection of Ihh into colon cancer cells leads to a downregulation of both components of the nuclear TCF4-beta-catenin complex and abrogates endogenous Wnt signaling in vitro. In turn, expression of Ihh is downregulated in polyps of individuals with FAP and expression of doxycycline-inducible dominant negative TCF4 (dnTCF4) restores Ihh expression in APC mutant DLD-1 colon cancer cells. These data identify a new Wnt-Hh axis in colonic epithelial renewal.     

MYCN haploinsufficiency is associated with reduced brain size and intestinal atresias in Feingold syndrome

Feingold syndrome is characterized by variable combinations of esophageal and duodenal atresias, microcephaly, learning disability, syndactyly and cardiac defect. We show here that heterozygous mutations in the gene MYCN are present in Feingold syndrome. All mutations are predicted to disrupt both the full-length protein and a new shortened MYCN isoform, suggesting that multiple aspects of early embryogenesis and postnatal brain growth in humans are tightly regulated by MYCN dosage.