Spontaneous single-nucleotide substitutions and microsatellite mutations have distinct distributions of fitness effects.

The fitness effects of new mutations determine key properties of evolutionary processes. Beneficial mutations drive evolution, yet selection is also shaped by the frequency of small-effect deleterious mutations, whose combined effect can burden otherwise adaptive lineages and alter evolutionary trajectories and outcomes in clonally evolving organisms such as viruses, microbes, and tumors. The small effect sizes of these important mutations have made accurate measurements of their rates difficult. In microbes, assessing the effect of mutations on growth can be especially instructive, as this complex phenotype is closely linked to fitness in clonally evolving organisms. Here, we perform high-throughput time-lapse microscopy on cells from mutation-accumulation strains to precisely infer the distribution of mutational effects on growth rate in the budding yeast, Saccharomyces cerevisiae. We show that mutational effects on growth rate are overwhelmingly negative, highly skewed towards very small effect sizes, and frequent enough to suggest that deleterious hitchhikers may impose a significant burden on evolving lineages. By using lines that accumulated mutations in either wild-type or slippage repair-defective backgrounds, we further disentangle the effects of 2 common types of mutations, single-nucleotide substitutions and simple sequence repeat indels, and show that they have distinct effects on yeast growth rate. Although the average effect of a simple sequence repeat mutation is very small (approximately 0.3%), many do alter growth rate, implying that this class of frequent mutations has an important evolutionary impact.

Extent and context dependence of pleiotropy revealed by high-throughput single-cell phenotyping.

Pleiotropy-when a single mutation affects multiple traits-is a controversial topic with far-reaching implications. Pleiotropy plays a central role in debates about how complex traits evolve and whether biological systems are modular or are organized such that every gene has the potential to affect many traits. Pleiotropy is also critical to initiatives in evolutionary medicine that seek to trap infectious microbes or tumors by selecting for mutations that encourage growth in some conditions at the expense of others. Research in these fields, and others, would benefit from understanding the extent to which pleiotropy reflects inherent relationships among phenotypes that correlate no matter the perturbation (vertical pleiotropy). Alternatively, pleiotropy may result from genetic changes that impose correlations between otherwise independent traits (horizontal pleiotropy). We distinguish these possibilities by using clonal populations of yeast cells to quantify the inherent relationships between single-cell morphological features. Then, we demonstrate how often these relationships underlie vertical pleiotropy and how often these relationships are modified by genetic variants (quantitative trait loci [QTL]) acting via horizontal pleiotropy. Our comprehensive screen measures thousands of pairwise trait correlations across hundreds of thousands of yeast cells and reveals ample evidence of both vertical and horizontal pleiotropy. Additionally, we observe that the correlations between traits can change with the environment, genetic background, and cell-cycle position. These changing dependencies suggest a nuanced view of pleiotropy: biological systems demonstrate limited pleiotropy in any given context, but across contexts (e.g., across diverse environments and genetic backgrounds) each genetic change has the potential to influence a larger number of traits. Our method suggests that exploiting pleiotropy for applications in evolutionary medicine would benefit from focusing on traits with correlations that are less dependent on context.

Decanalizing thinking on genetic canalization.

The concept of genetic canalization has had an abiding influence on views of complex-trait evolution. A genetically canalized system has evolved to become less sensitive to the effects of mutation. When a gene product that supports canalization is compromised, the phenotypic impacts of a mutation should be more pronounced. This expected increase in mutational effects not only has important consequences for evolution, but has also motivated strategies to treat disease. However, recent studies demonstrate that, when putative agents of genetic canalization are impaired, systems do not behave as expected. Here, we review the evidence that is used to infer whether particular gene products are agents of genetic canalization. Then we explain how such inferences often succumb to a converse error. We go on to show that several candidate agents of genetic canalization increase the phenotypic impacts of some mutations while decreasing the phenotypic impacts of others. These observations suggest that whether a gene product acts as a 'buffer' (lessening mutational effects) or a 'potentiator' (increasing mutational effects) is not a fixed property of the gene product but instead differs for the different mutations with which it interacts. To investigate features of genetic interactions that might predispose them toward buffering versus potentiation, we explore simulated gene-regulatory networks. Similarly to putative agents of genetic canalization, the gene products in simulated networks also modify the phenotypic effects of mutations in other genes without a strong overall tendency towards lessening or increasing these effects. In sum, these observations call into question whether complex traits have evolved to become less sensitive (i.e., are canalized) to genetic change, and the degree to which trends exist that predict how one genetic change might alter another's impact. We conclude by discussing approaches to address these and other open questions that are brought into focus by re-thinking genetic canalization.

Control of nongenetic heterogeneity in growth rate and stress tolerance of Saccharomyces cerevisiae by cyclic AMP-regulated transcription factors.

Genetically identical cells exhibit extensive phenotypic variation even under constant and benign conditions. This so-called nongenetic heterogeneity has important clinical implications: within tumors and microbial infections, cells show nongenetic heterogeneity in growth rate and in susceptibility to drugs or stress. The budding yeast, Saccharomyces cerevisiae, shows a similar form of nongenetic heterogeneity in which growth rate correlates positively with susceptibility to acute heat stress at the single-cell level. Using genetic and chemical perturbations, combined with high-throughput single-cell assays of yeast growth and gene expression, we show here that heterogeneity in intracellular cyclic AMP (cAMP) levels acting through the conserved Ras/cAMP/protein kinase A (PKA) pathway and its target transcription factors, Msn2 and Msn4, underlies this nongenetic heterogeneity. Lower levels of cAMP correspond to slower growth, as shown by direct comparison of cAMP concentration in subpopulations enriched for slower vs. faster growing cells. Concordantly, an endogenous reporter of this pathway's activity correlates with growth in individual cells. The paralogs Msn2 and Msn4 differ in their roles in nongenetic heterogeneity in a way that demonstrates slow growth and stress tolerance are not inevitably linked. Heterogeneity in growth rate requires each, whereas only Msn2 is required for heterogeneity in expression of Tsl1, a subunit of trehalose synthase that contributes to acute-stress tolerance. Perturbing nongenetic heterogeneity by mutating genes in this pathway, or by culturing wild-type cells with the cell-permeable cAMP analog 8-bromo-cAMP or the PKA inhibitor H89, significantly impacts survival of acute heat stress. Perturbations that increase intracellular cAMP levels reduce the slower-growing subpopulation and increase susceptibility to acute heat stress, whereas PKA inhibition slows growth and decreases susceptibility to acute heat stress. Loss of Msn2 reduces, but does not completely eliminate, the correlation in individual cells between growth rate and acute-stress survival, suggesting a major role for the Msn2 pathway in nongenetic heterogeneity but also a residual benefit of slow growth. Our results shed light on the genetic control of nongenetic heterogeneity and suggest a possible means of defeating bet-hedging pathogens or tumor cells by making them more uniformly susceptible to treatment.

Deep sequencing of natural and experimental populations of Drosophila melanogaster reveals biases in the spectrum of new mutations.

Mutations provide the raw material of evolution, and thus our ability to study evolution depends fundamentally on having precise measurements of mutational rates and patterns. We generate a data set for this purpose using (1) de novo mutations from mutation accumulation experiments and (2) extremely rare polymorphisms from natural populations. The first, mutation accumulation (MA) lines are the product of maintaining flies in tiny populations for many generations, therefore rendering natural selection ineffective and allowing new mutations to accrue in the genome. The second, rare genetic variation from natural populations allows the study of mutation because extremely rare polymorphisms are relatively unaffected by the filter of natural selection. We use both methods in Drosophila melanogaster, first generating our own novel data set of sequenced MA lines and performing a meta-analysis of all published MA mutations (∼2000 events) and then identifying a high quality set of ∼70,000 extremely rare (≤0.1%) polymorphisms that are fully validated with resequencing. We use these data sets to precisely measure mutational rates and patterns. Highlights of our results include: a high rate of multinucleotide mutation events at both short (∼5 bp) and long (∼1 kb) genomic distances, showing that mutation drives GC content lower in already GC-poor regions, and using our precise context-dependent mutation rates to predict long-term evolutionary patterns at synonymous sites. We also show that de novo mutations from independent MA experiments display similar patterns of single nucleotide mutation and well match the patterns of mutation found in natural populations.

Author Correction: Testing the key assumption of heritability estimates based on genome-wide genetic relatedness.

In the original paper, we used the variable "URBRUR08," from the 2008 survey wave as a measure of childhood urbanicity. Upon further investigation we realized that this variable actually measured Beale urban-rural code during the respondent's adulthood.  Thus, we reran our analysis of the pseudo-heritability of childhood urbanicity using the variable. The original results hold such that even with the first 20 principal components held constant, childhood urban-rural status appears to be ~20% "heritable" in GREML models-a figure that is actually higher than the original estimate reported in the paper (14% controlling for 25 PCs, 15% controlling for 10 PCs, and 29% controlling for two PCs). Meanwhile, the heritabilities of the other phenotypes-height, BMI and education-still do not change when they are residualized on childhood urbanicity. In other words, the original results of the paper do not change.

Selection Transforms the Landscape of Genetic Variation Interacting with Hsp90.

The protein-folding chaperone Hsp90 has been proposed to buffer the phenotypic effects of mutations. The potential for Hsp90 and other putative buffers to increase robustness to mutation has had major impact on disease models, quantitative genetics, and evolutionary theory. But Hsp90 sometimes contradicts expectations for a buffer by potentiating rapid phenotypic changes that would otherwise not occur. Here, we quantify Hsp90's ability to buffer or potentiate (i.e., diminish or enhance) the effects of genetic variation on single-cell morphological features in budding yeast. We corroborate reports that Hsp90 tends to buffer the effects of standing genetic variation in natural populations. However, we demonstrate that Hsp90 tends to have the opposite effect on genetic variation that has experienced reduced selection pressure. Specifically, Hsp90 tends to enhance, rather than diminish, the effects of spontaneous mutations and recombinations. This result implies that Hsp90 does not make phenotypes more robust to the effects of genetic perturbation. Instead, natural selection preferentially allows buffered alleles to persist and thereby creates the false impression that Hsp90 confers greater robustness.

Shifting sugars and shifting paradigms.

No organism lives in a constant environment. Based on classical studies in molecular biology, many have viewed microbes as following strict rules for shifting their metabolic activities when prevailing conditions change. For example, students learn that the bacterium Escherichia coli makes proteins for digesting lactose only when lactose is available and glucose, a better sugar, is not. However, recent studies, including three PLOS Biology papers examining sugar utilization in the budding yeast Saccharomyces cerevisiae, show that considerable heterogeneity in response to complex environments exists within and between populations. These results join similar recent results in other organisms that suggest that microbial populations anticipate predictable environmental changes and hedge their bets against unpredictable ones. The classical view therefore represents but one special case in a range of evolutionary adaptations to environmental changes that all organisms face.

Precise estimates of mutation rate and spectrum in yeast.

Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.

Essential gene disruptions reveal complex relationships between phenotypic robustness, pleiotropy, and fitness.

The concept of robustness in biology has gained much attention recently, but a mechanistic understanding of how genetic networks regulate phenotypic variation has remained elusive. One approach to understand the genetic architecture of variability has been to analyze dispensable gene deletions in model organisms; however, the most important genes cannot be deleted. Here, we have utilized two systems in yeast whereby essential genes have been altered to reduce expression. Using high-throughput microscopy and image analysis, we have characterized a large number of morphological phenotypes, and their associated variation, for the majority of essential genes in yeast. Our results indicate that phenotypic robustness is more highly dependent upon the expression of essential genes than on the presence of dispensable genes. Morphological robustness appears to be a general property of a genotype that is closely related to pleiotropy. While the fitness profile across a range of expression levels is idiosyncratic to each gene, the global pattern indicates that there is a window in which phenotypic variation can be released before fitness effects are observable.

Histone variant HTZ1 shows extensive epistasis with, but does not increase robustness to, new mutations.

Biological systems produce phenotypes that appear to be robust to perturbation by mutations and environmental variation. Prior studies identified genes that, when impaired, reveal previously cryptic genetic variation. This result is typically interpreted as evidence that the disrupted gene normally increases robustness to mutations, as such robustness would allow cryptic variants to accumulate. However, revelation of cryptic genetic variation is not necessarily evidence that a mutationally robust state has been made less robust. Demonstrating a difference in robustness requires comparing the ability of each state (with the gene perturbed or intact) to suppress the effects of new mutations. Previous studies used strains in which the existing genetic variation had been filtered by selection. Here, we use mutation accumulation (MA) lines that have experienced minimal selection, to test the ability of histone H2A.Z (HTZ1) to increase robustness to mutations in the yeast Saccharomyces cerevisiae. HTZ1, a regulator of chromatin structure and gene expression, represents a class of genes implicated in mutational robustness. It had previously been shown to increase robustness of yeast cell morphology to fluctuations in the external or internal microenvironment. We measured morphological variation within and among 79 MA lines with and without HTZ1. Analysis of within-line variation confirms that HTZ1 increases microenvironmental robustness. Analysis of between-line variation shows the morphological effects of eliminating HTZ1 to be highly dependent on the line, which implies that HTZ1 interacts with mutations that have accumulated in the lines. However, lines without HTZ1 are, as a group, not more phenotypically diverse than lines with HTZ1 present. The presence of HTZ1, therefore, does not confer greater robustness to mutations than its absence. Our results provide experimental evidence that revelation of cryptic genetic variation cannot be assumed to be caused by loss of robustness, and therefore force reevaluation of prior claims based on that assumption.

Polygenic cis-regulatory adaptation in the evolution of yeast pathogenicity.

The acquisition of new genes, via horizontal transfer or gene duplication/diversification, has been the dominant mechanism thus far implicated in the evolution of microbial pathogenicity. In contrast, the role of many other modes of evolution--such as changes in gene expression regulation-remains unknown. A transition to a pathogenic lifestyle has recently taken place in some lineages of the budding yeast Saccharomyces cerevisiae. Here we identify a module of physically interacting proteins involved in endocytosis that has experienced selective sweeps for multiple cis-regulatory mutations that down-regulate gene expression levels in a pathogenic yeast. To test if these adaptations affect virulence, we created a panel of single-allele knockout strains whose hemizygous state mimics the genes' adaptive down-regulations, and measured their virulence in a mammalian host. Despite having no growth advantage in standard laboratory conditions, nearly all of the strains were more virulent than their wild-type progenitor, suggesting that these adaptations likely played a role in the evolution of pathogenicity. Furthermore, genetic variants at these loci were associated with clinical origin across 88 diverse yeast strains, suggesting the adaptations may have contributed to the virulence of a wide range of clinical isolates. We also detected pleiotropic effects of these adaptations on a wide range of morphological traits, which appear to have been mitigated by compensatory mutations at other loci. These results suggest that cis-regulatory adaptation can occur at the level of physically interacting modules and that one such polygenic adaptation led to increased virulence during the evolution of a pathogenic yeast.

On the Nature and Evolutionary Impact of Phenotypic Robustness Mechanisms.

Biologists have long observed that physiological and developmental processes are insensitive, or robust, to many genetic and environmental perturbations. A complete understanding of the evolutionary causes and consequences of this robustness is lacking. Recent progress has been made in uncovering the regulatory mechanisms that underlie environmental robustness in particular. Less is known about robustness to the effects of mutations, and indeed the evolution of mutational robustness remains a controversial topic. The controversy has spread to related topics, in particular the evolutionary relevance of cryptic genetic variation. This review aims to synthesize current understanding of robustness mechanisms and to cut through the controversy by shedding light on what is and is not known about mutational robustness. Some studies have confused mutational robustness with non-additive interactions between mutations (epistasis). We conclude that a profitable way forward is to focus investigations (and rhetoric) less on mutational robustness and more on epistasis.

Bet hedging in yeast by heterogeneous, age-correlated expression of a stress protectant.

Genetically identical cells grown in the same culture display striking cell-to-cell heterogeneity in gene expression and other traits. A crucial challenge is to understand how much of this heterogeneity reflects the noise tolerance of a robust system and how much serves a biological function. In bacteria, stochastic gene expression results in cell-to-cell heterogeneity that might serve as a bet-hedging mechanism, allowing a few cells to survive through an antimicrobial treatment while others perish. Despite its clinical importance, the molecular mechanisms underlying bet hedging remain unclear. Here, we investigate the mechanisms of bet hedging in Saccharomyces cerevisiae using a new high-throughput microscopy assay that monitors variable protein expression, morphology, growth rate, and survival outcomes of tens of thousands of yeast microcolonies simultaneously. We find that clonal populations display broad distributions of growth rates and that slow growth predicts resistance to heat killing in a probabalistic manner. We identify several gene products that are likely to play a role in bet hedging and confirm that Tsl1, a trehalose-synthesis regulator, is an important component of this resistance. Tsl1 abundance correlates with growth rate and replicative age and predicts survival. Our results suggest that yeast bet hedging results from multiple epigenetic growth states determined by a combination of stochastic and deterministic factors.

Sibling genes as environment: Sibling dopamine genotypes and adolescent health support frequency dependent selection.

While research consistently suggests siblings matter for individual outcomes, it remains unclear why. At the same time, studies of genetic effects on health typically correlate variants of a gene with the average level of behavioral or health measures, ignoring more complicated genetic dynamics. Using National Longitudinal Study of Adolescent Health data, we investigate whether sibling genes moderate individual genetic expression. We compare twin variation in health-related absences and self-rated health by genetic differences at three locations related to dopamine regulation and transport to test sibship-level cross-person gene-gene interactions. Results suggest effects of variation at these genetic locations are moderated by sibling genes. Although the mechanism remains unclear, this evidence is consistent with frequency dependent selection and suggests much genetic research may violate the stable unit treatment value assumption.

Genetic and nongenetic determinants of cell growth variation assessed by high-throughput microscopy.

In microbial populations, growth initiation and proliferation rates are major components of fitness and therefore likely targets of selection. We used a high-throughput microscopy assay, which enables simultaneous analysis of tens of thousands of microcolonies, to determine the sources and extent of growth rate variation in the budding yeast (Saccharomyces cerevisiae) in different glucose environments. We find that cell growth rates are regulated by the extracellular concentration of glucose as proposed by Monod (1949), but that significant heterogeneity in growth rates is observed among genetically identical individuals within an environment. Yeast strains isolated from different geographic locations and habitats differ in their growth rate responses to different glucose concentrations. Inheritance patterns suggest that the genetic determinants of growth rates in different glucose concentrations are distinct. In addition, we identified genotypes that differ in the extent of variation in growth rate within an environment despite nearly identical mean growth rates, providing evidence that alleles controlling phenotypic variability segregate in yeast populations. We find that the time to reinitiation of growth (lag) is negatively correlated with growth rate, yet this relationship is strain-dependent. Between environments, the respirative activity of individual cells negatively correlates with glucose abundance and growth rate, but within an environment respirative activity and growth rate show a positive correlation, which we propose reflects differences in protein expression capacity. Our study quantifies the sources of genetic and nongenetic variation in cell growth rates in different glucose environments with unprecedented precision, facilitating their molecular genetic dissection.

The female-specific doublesex isoform regulates pleiotropic transcription factors to pattern genital development in Drosophila.

Regulatory networks driving morphogenesis of animal genitalia must integrate sexual identity and positional information. Although the genetic hierarchy that controls somatic sexual identity in the fly Drosophila melanogaster is well understood, there are very few cases in which the mechanism by which it controls tissue-specific gene activity is known. In flies, the sex-determination hierarchy terminates in the doublesex (dsx) gene, which produces sex-specific transcription factors via alternative splicing of its transcripts. To identify sex-specifically expressed genes downstream of dsx that drive the sexually dimorphic development of the genitalia, we performed genome-wide transcriptional profiling of dissected genital imaginal discs of each sex at three time points during early morphogenesis. Using a stringent statistical threshold, we identified 23 genes that have sex-differential transcript levels at all three time points, of which 13 encode transcription factors, a significant enrichment. We focus here on three sex-specifically expressed transcription factors encoded by lozenge (lz), Drop (Dr) and AP-2. We show that, in female genital discs, Dsx activates lz and represses Dr and AP-2. We further show that the regulation of Dr by Dsx mediates the previously identified expression of the fibroblast growth factor Branchless in male genital discs. The phenotypes we observe upon loss of lz or Dr function in genital discs explain the presence or absence of particular structures in dsx mutant flies and thereby clarify previously puzzling observations. Our time course of expression data also lays the foundation for elucidating the regulatory networks downstream of the sex-specifically deployed transcription factors.

Testing the key assumption of heritability estimates based on genome-wide genetic relatedness.

Comparing genetic and phenotypic similarity among unrelated individuals seems a promising way to quantify the genetic component of traits while avoiding the problematic assumptions plaguing twin- and other kin-based estimates of heritability. One approach uses a Genetic Relatedness Estimation through Maximum Likelihood (GREML) model for individuals who are related at less than 0.025 to predict their phenotypic similarity by their genetic similarity. Here we test the key underlying assumption of this approach: that genetic relatedness is orthogonal to environmental similarity. Using data from the Health and Retirement Study (and two other surveys), we show two unrelated individuals may be more likely to have been reared in a similar environment (urban versus nonurban setting) if they are genetically similar. This effect is not eliminated by controls for population structure. However, when we include this environmental confound in GREML models, heritabilities do not change substantially and thus potential bias in estimates of most biological phenotypes is probably minimal.

Sperm-storage defects and live birth in Drosophila females lacking spermathecal secretory cells.

Male Drosophila flies secrete seminal-fluid proteins that mediate proper sperm storage and fertilization, and that induce changes in female behavior. Females also produce reproductive-tract secretions, yet their contributions to postmating physiology are poorly understood. Large secretory cells line the female's spermathecae, a pair of sperm-storage organs. We identified the regulatory regions controlling transcription of two genes exclusively expressed in these spermathecal secretory cells (SSC): Spermathecal endopeptidase 1 (Send1), which is expressed in both unmated and mated females, and Spermathecal endopeptidase 2 (Send2), which is induced by mating. We used these regulatory sequences to perform precise genetic ablations of the SSC at distinct time points relative to mating. We show that the SSC are required for recruiting sperm to the spermathecae, but not for retaining sperm there. The SSC also act at a distance in the reproductive tract, in that their ablation: (1) reduces sperm motility in the female's other sperm-storage organ, the seminal receptacle; and (2) causes ovoviviparity--the retention and internal development of fertilized eggs. These results establish the reproductive functions of the SSC, shed light on the evolution of live birth, and open new avenues for studying and manipulating female fertility in insects.