Evolutionary capacitance as a general feature of complex gene networks

An evolutionary capacitor buffers genotypic variation under normal conditions, thereby promoting the accumulation of hidden polymorphism. But it occasionally fails, thereby revealing this variation phenotypically. The principal example of an evolutionary capacitor is Hsp90, a molecular chaperone that targets an important set of signal transduction proteins. Experiments in Drosophila and Arabidopsis have demonstrated three key properties of Hsp90: (1) it suppresses phenotypic variation under normal conditions and releases this variation when functionally compromised; (2) its function is overwhelmed by environmental stress; and (3) it exerts pleiotropic effects on key developmental processes. But whether these properties necessarily make Hsp90 a significant and unique facilitator of adaptation is unclear. Here we use numerical simulations of complex gene networks, as well as genome-scale expression data from yeast single-gene deletion strains, to present a mechanism that extends the scope of evolutionary capacitance beyond the action of Hsp90 alone. We illustrate that most, and perhaps all, genes reveal phenotypic variation when functionally compromised, and that the availability of loss-of-function mutations accelerates adaptation to a new optimum phenotype. However, this effect does not require the mutations to be conditional on the environment. Thus, there might exist a large class of evolutionary capacitors whose effects on phenotypic variation complement the systemic, environment-induced effects of Hsp90.     

Hsp90 as a capacitor of phenotypic variation

Heat-shock protein 90 (Hsp90) chaperones the maturation of many regulatory proteins and, in the fruitfly Drosophila melanogaster, buffers genetic variation in morphogenetic pathways. Levels and patterns of genetic variation differ greatly between obligatorily outbreeding species such as fruitflies and self-fertilizing species such as the plant Arabidopsis thaliana. Also, plant development is more plastic, being coupled to environmental cues. Here we report that, in Arabidopsis accessions and recombinant inbred lines, reducing Hsp90 function produces an array of morphological phenotypes, which are dependent on underlying genetic variation. The strength and breadth of Hsp90's effects on the buffering and release of genetic variation suggests it may have an impact on evolutionary processes. We also show that Hsp90 influences morphogenetic responses to environmental cues and buffers normal development from destabilizing effects of stochastic processes. Manipulating Hsp90's buffering capacity offers a tool for harnessing cryptic genetic variation and for elucidating the interplay between genotypes, environments and stochastic events in the determination of phenotype.     

Hsp90 stress potentiates rapid cellular adaptation through induction of aneuploidy

Aneuploidy--the state of having uneven numbers of chromosomes--is a hallmark of cancer and a feature identified in yeast from diverse habitats. Recent studies have shown that aneuploidy is a form of large-effect mutation that is able to confer adaptive phenotypes under diverse stress conditions. Here we investigate whether pleiotropic stress could induce aneuploidy in budding yeast (Saccharomyces cerevisae). We show that whereas diverse stress conditions can induce an increase in chromosome instability, proteotoxic stress, caused by transient Hsp90 (also known as Hsp82 or Hsc82) inhibition or heat shock, markedly increased chromosome instability to produce a cell population with high karyotype diversity. The induced chromosome instability is linked to an evolutionarily conserved role for the Hsp90 chaperone complex in kinetochore assembly. Continued growth in the presence of an Hsp90 inhibitor resulted in the emergence of drug-resistant colonies with chromosome XV gain. This drug-resistance phenotype is a quantitative trait involving copy number increases of at least two genes located on chromosome XV. Short-term exposure to Hsp90 stress potentiated fast adaptation to unrelated cytotoxic compounds by means of different aneuploid chromosome stoichiometries. These findings demonstrate that aneuploidy is a form of stress-inducible mutation in eukaryotes, capable of fuelling rapid phenotypic evolution and drug resistance, and reveal a new role for Hsp90 in regulating the emergence of adaptive traits under stress.     

A common sex-dependent mutation in a RET enhancer underlies Hirschsprung disease risk

The identification of common variants that contribute to the genesis of human inherited disorders remains a significant challenge. Hirschsprung disease (HSCR) is a multifactorial, non-mendelian disorder in which rare high-penetrance coding sequence mutations in the receptor tyrosine kinase RET contribute to risk in combination with mutations at other genes. We have used family-based association studies to identify a disease interval, and integrated this with comparative and functional genomic analysis to prioritize conserved and functional elements within which mutations can be sought. We now show that a common non-coding RET variant within a conserved enhancer-like sequence in intron 1 is significantly associated with HSCR susceptibility and makes a 20-fold greater contribution to risk than rare alleles do. This mutation reduces in vitro enhancer activity markedly, has low penetrance, has different genetic effects in males and females, and explains several features of the complex inheritance pattern of HSCR. Thus, common low-penetrance variants, identified by association studies, can underlie both common and rare diseases.     

Abnormal pattern detected in fragile-X patients by pulsed-field gel electrophoresis.

The fragile-X syndrome is the most frequent inherited form of mental retardation, with an incidence of 1 in 1,500 males. It is characterized by the presence of a fragile site at Xq27.3 induced in vitro by folate deprivation or by inhibitors of deoxynucleotide synthesis. Its mode of inheritance is unusual for an X-linked trait, with incomplete penetrance in both males and females. Some phenotypically normal males transmit the mutation to all their daughters who rarely express any symptoms, but penetrance is high in sons and daughters of these carrier women. Genetic and physical mapping of the Xq27-q28 region has confirmed that the disease locus is located at or very near the fragile site. Hypotheses proposed to account for the abnormalities in the inheritance of the disease include sequence rearrangements by meiotic recombination or a mutation that affects reactivation of an inactive X chromosome during differentiation of female germ cells. To detect such rearrangements, or methylation changes that may reflect a locally inactive X chromosome, we used pulsed-field gel analysis of DNA from fragile-X patients with probes close to the fragile-X locus. The probe Do33 (DXS465) detected abnormal patterns in fragile-X patients, but not in normal controls or in non-expressing male transmitters.     

Familial predisposition to Wilms' tumour does not map to the short arm of chromosome 11

Wilms' tumour of the kidney usually occurs sporadically, but can also segregate as an autosomal dominant trait with incomplete penetrance. Patients with the WAGR syndrome of aniridia, genitourinary anomalies, mental retardation and high risk of Wilms' tumour have overlapping deletions of chromosome 11p13 which has suggested a possible location for a Wilms' tumour locus. Moreover, many sporadic tumours have lost a portion of chromosome 11p. A second locus at 11p15 is implicated by association of the tumour with the Wiedemann-Beckwith syndrome and by tumour-specific losses of chromosome 11 confined to 11p15. Here we report a multipoint linkage analysis of a family segregating for Wilms' tumour, using polymorphic DNA markers mapped to chromosome 11p. The results exclude the predisposing mutation from both locations. In a second family, the 11p15 alleles lost in the tumour were derived from the affected parent, thus precluding this region as the location of the inherited mutation. These findings imply an aetiological heterogeneity for Wilms' tumour and raise questions concerning the general applicability of the carcinogenesis model that has been useful in the understanding of retinoblastoma.     

A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig

Most traits and disorders have a multifactorial background indicating that they are controlled by environmental factors as well as an unknown number of quantitative trait loci (QTLs). The identification of mutations underlying QTLs is a challenge because each locus explains only a fraction of the phenotypic variation. A paternally expressed QTL affecting muscle growth, fat deposition and size of the heart in pigs maps to the IGF2 (insulin-like growth factor 2) region. Here we show that this QTL is caused by a nucleotide substitution in intron 3 of IGF2. The mutation occurs in an evolutionarily conserved CpG island that is hypomethylated in skeletal muscle. The mutation abrogates in vitro interaction with a nuclear factor, probably a repressor, and pigs inheriting the mutation from their sire have a threefold increase in IGF2 messenger RNA expression in postnatal muscle. Our study establishes a causal relationship between a single-base-pair substitution in a non-coding region and a QTL effect. The result supports the long-held view that regulatory mutations are important for controlling phenotypic variation.     

Non-adaptive origins of interactome complexity

The boundaries between prokaryotes, unicellular eukaryotes and multicellular eukaryotes are accompanied by orders-of-magnitude reductions in effective population size, with concurrent amplifications of the effects of random genetic drift and mutation. The resultant decline in the efficiency of selection seems to be sufficient to influence a wide range of attributes at the genomic level in a non-adaptive manner. A key remaining question concerns the extent to which variation in the power of random genetic drift is capable of influencing phylogenetic diversity at the subcellular and cellular levels. Should this be the case, population size would have to be considered as a potential determinant of the mechanistic pathways underlying long-term phenotypic evolution. Here we demonstrate a phylogenetically broad inverse relation between the power of drift and the structural integrity of protein subunits. This leads to the hypothesis that the accumulation of mildly deleterious mutations in populations of small size induces secondary selection for protein-protein interactions that stabilize key gene functions. By this means, the complex protein architectures and interactions essential to the genesis of phenotypic diversity may initially emerge by non-adaptive mechanisms.     

Direct estimation of per nucleotide and genomic deleterious mutation rates in Drosophila

Spontaneous mutations are the source of genetic variation required for evolutionary change, and are therefore important for many aspects of evolutionary biology. For example, the divergence between taxa at neutrally evolving sites in the genome is proportional to the per nucleotide mutation rate, u (ref. 1), and this can be used to date speciation events by assuming a molecular clock. The overall rate of occurrence of deleterious mutations in the genome each generation (U) appears in theories of nucleotide divergence and polymorphism, the evolution of sex and recombination, and the evolutionary consequences of inbreeding. However, estimates of U based on changes in allozymes or DNA sequences and fitness traits are discordant. Here we directly estimate u in Drosophila melanogaster by scanning 20 million bases of DNA from three sets of mutation accumulation lines by using denaturing high-performance liquid chromatography. From 37 mutation events that we detected, we obtained a mean estimate for u of 8.4 x 10(-9) per generation. Moreover, we detected significant heterogeneity in u among the three mutation-accumulation-line genotypes. By multiplying u by an estimate of the fraction of mutations that are deleterious in natural populations of Drosophila, we estimate that U is 1.2 per diploid genome. This high rate suggests that selection against deleterious mutations may have a key role in explaining patterns of genetic variation in the genome, and help to maintain recombination and sexual reproduction.     

Climate-driven population divergence in sex-determining systems

Sex determination is a fundamental biological process, yet its mechanisms are remarkably diverse. In vertebrates, sex can be determined by inherited genetic factors or by the temperature experienced during embryonic development. However, the evolutionary causes of this diversity remain unknown. Here we show that live-bearing lizards at different climatic extremes of the species' distribution differ in their sex-determining mechanisms, with temperature-dependent sex determination in lowlands and genotypic sex determination in highlands. A theoretical model parameterized with field data accurately predicts this divergence in sex-determining systems and the consequence thereof for variation in cohort sex ratios among years. Furthermore, we show that divergent natural selection on sex determination across altitudes is caused by climatic effects on lizard life history and variation in the magnitude of between-year temperature fluctuations. Our results establish an adaptive explanation for intra-specific divergence in sex-determining systems driven by phenotypic plasticity and ecological selection, thereby providing a unifying framework for integrating the developmental, ecological and evolutionary basis for variation in vertebrate sex determination.     

Cryptic genetic variation promotes rapid evolutionary adaptation in an RNA enzyme

Cryptic variation is caused by the robustness of phenotypes to mutations. Cryptic variation has no effect on phenotypes in a given genetic or environmental background, but it can have effects after mutations or environmental change. Because evolutionary adaptation by natural selection requires phenotypic variation, phenotypically revealed cryptic genetic variation may facilitate evolutionary adaptation. This is possible if the cryptic variation happens to be pre-adapted, or "exapted", to a new environment, and is thus advantageous once revealed. However, this facilitating role for cryptic variation has not been proven, partly because most pertinent work focuses on complex phenotypes of whole organisms whose genetic basis is incompletely understood. Here we show that populations of RNA enzymes with accumulated cryptic variation adapt more rapidly to a new substrate than a population without cryptic variation. A detailed analysis of our evolving RNA populations in genotype space shows that cryptic variation allows a population to explore new genotypes that become adaptive only in a new environment. Our observations show that cryptic variation contains new genotypes pre-adapted to a changed environment. Our results highlight the positive role that robustness and epistasis can have in adaptive evolution.     

Genetic linkage of Werner's syndrome to five markers on chromosome 8.

Werner's syndrome (WS) is a rare autosomal recessive disease in which the affected individuals display symptoms of premature ageing. The substantial phenotypic overlap between WS and normal ageing indicates that these two conditions may have pathogenetic mechanisms in common. The WS mutation has pleiotropic effects, and patients and their cells show many differences compared with normals. Despite extensive study of the clinical and biochemical features of this disorder, the primary genetic defect remains unknown. We have undertaken a genetic linkage study in an effort to identify the locus of the primary defect. Here we report close genetic linkage of the WS mutation to a group of markers on chromosome 8.     

A map of human genome variation from population-scale sequencing

The 1000 Genomes Project aims to provide a deep characterization of human genome sequence variation as a foundation for investigating the relationship between genotype and phenotype. Here we present results of the pilot phase of the project, designed to develop and compare different strategies for genome-wide sequencing with high-throughput platforms. We undertook three projects: low-coverage whole-genome sequencing of 179 individuals from four populations; high-coverage sequencing of two mother-father-child trios; and exon-targeted sequencing of 697 individuals from seven populations. We describe the location, allele frequency and local haplotype structure of approximately 15 million single nucleotide polymorphisms, 1 million short insertions and deletions, and 20,000 structural variants, most of which were previously undescribed. We show that, because we have catalogued the vast majority of common variation, over 95% of the currently accessible variants found in any individual are present in this data set. On average, each person is found to carry approximately 250 to 300 loss-of-function variants in annotated genes and 50 to 100 variants previously implicated in inherited disorders. We demonstrate how these results can be used to inform association and functional studies. From the two trios, we directly estimate the rate of de novo germline base substitution mutations to be approximately 10(-8) per base pair per generation. We explore the data with regard to signatures of natural selection, and identify a marked reduction of genetic variation in the neighbourhood of genes, due to selection at linked sites. These methods and public data will support the next phase of human genetic research.     

Non-invasive prenatal measurement of the fetal genome

The vast majority of prenatal genetic testing requires invasive sampling. However, this poses a risk to the fetus, so one must make a decision that weighs the desire for genetic information against the risk of an adverse outcome due to hazards of the testing process. These issues are not required to be coupled, and it would be desirable to discover genetic information about the fetus without incurring a health risk. Here we demonstrate that it is possible to non-invasively sequence the entire prenatal genome. Our results show that molecular counting of parental haplotypes in maternal plasma by shotgun sequencing of maternal plasma DNA allows the inherited fetal genome to be deciphered non-invasively. We also applied the counting principle directly to each allele in the fetal exome by performing exome capture on maternal plasma DNA before shotgun sequencing. This approach enables non-invasive exome screening of clinically relevant and deleterious alleles that were paternally inherited or had arisen as de novo germline mutations, and complements the haplotype counting approach to provide a comprehensive view of the fetal genome. Non-invasive determination of the fetal genome may ultimately facilitate the diagnosis of all inherited and de novo genetic disease.     

Germ-line transmission of a mutated p53 gene in a cancer-prone family with Li-Fraumeni syndrome

Tumour suppressor genes, whose usual function seems to be controlling normal cell proliferation, have been implicated in many inherited and sporadic forms of malignancies Much evidence supports the concept of tumour formation by loss-of-function mutations in suppressor genes, as predicted by the two-hit model of Knudson and DeMars. The suppressor gene, p53, is affected in such a manner by numerous mutations, which occur in a variety of human tumours. These mutations usually represent the loss of one allele and the substitution of a single base in the other. We have now analysed the p53 gene in a family affected by Li-Fraumeni syndrome, a rare autosomal dominant syndrome characterized by the occurrence of diverse mesenchymal and epithelial neoplasms at multiple sites. In some instances the neoplasms seem to be related to exposure to carcinogens, including ionizing radiation. The Li-Fraumeni family that we studied had noncancerous skin fibroblasts (NSF) with an unusual radiation-resistant phenotype. DNA derived from the NSF cells of four family members, spanning two generations, had the same point mutation in codon 245 (GGC----GAC) of the p53 gene. This mutation leads to substitution of aspartic acid for glycine in one of the regions identified as a frequent target of point mutations in p53. The NSF cell lines with the mutation also retained the normal p53 allele. This inherited p53 mutation may predispose the members of this family to increased susceptibility to cancer.     

Mutation and sex in a competitive world

How do deleterious mutations interact to affect fitness? The answer to this question has substantial implications for a variety of important problems in population biology, including the evolution of sex, the rate of adaptation and the conservation of small populations. Here we analyse a mathematical model of competition for food in which deleterious mutations affect competitive ability. We show that, if individuals usually compete in small groups, then competition can easily lead to a type of genetic interaction known as synergistic epistasis. This means that a deleterious mutation is most damaging in a genome that already has many other deleterious mutations. We also show that competition in small groups can produce a large advantage for sexual populations, both in mean fitness and in ability to resist invasion by asexual lineages. One implication of our findings is that experimental efforts to demonstrate synergistic epistasis may not succeed unless the experiments are redesigned to make them much more naturalistic.     

Antigenic variation in trypanosomes

In its mammalian host, Trypanosoma brucei is able to change the antigenic character of its glycoprotein surface coat and so evade the host's immune response. This phenotypic change seems to occur spontaneously in 1 in 10,000 individuals but is not due to genetic mutation: host antibody is not necessary for its induction but plays a selective part in bringing about the gross changes in parasite numbers and antigenic character observed in the bloodstream by destroying the main component of what is actually a heterogeneous population. The infecting trypanosome population injected into the mammalian host by the tsetse fly vector may also be heterogeneous. Such heterogeneity complicates plans to vaccinate cattle and people against the African trypanosomes based on the premise that the metacyclic trypanosomes of a clone bear the same surface antigen.     

A polymorphic DNA marker linked to cystic fibrosis is located on chromosome 7

Although cystic fibrosis (CF) is among the most common inherited diseases in Caucasian populations, the basic biochemical defect is not yet known. CF is inherited as an autosomal recessive trait apparently due to mutations in a single gene, whence the efforts made to identify the genetic locus responsible by linkage studies. Two markers have recently been identified that are genetically linked to CF: one is a genetic variation in serum level of activity of the enzyme paraoxonase, and the other is a restriction fragment length polymorphism (RFLP) identified with a randomly isolated DNA probe. We report here that the genetic locus DOCRI-917 defined by the cloned DNA probe is located on chromosome 7.