Honors & Awards


  • Life Science Alliance Exchange Grant Recipient, Life Science Alliance/EMBL (2023)
  • Hanna Gray Fellow, Howard Hughes Medical Institute (2022-ongoing)
  • Postdoctoral NRSA Fellowship awardee, F32 NIAID Ruth L. Kirschstein Postdoctoral Individual National Research Service Award (2021-2024)
  • Postdoctoral Fellowship awardee, School of Medicine Dean's Postdoctoral Fellowship (2021)
  • Postdoctoral Fellow, Stanford Center for Computational, Evolutionary and Human Genomics (CEHG) (2020-2021)
  • Postdoctoral Training Grant awardee, Stanford Genomics Training Program NHGRI training grant (2020-2021)
  • Graduate Research Fellowship Program (GRFP) awardee, National Science Foundation (2015-2019)
  • Predoctoral Training Grant awardee, University of Washington NHGRI Genome Training Grant (2014-2015)

Stanford Advisors


Lab Affiliations


All Publications


  • The Viral K1 Killer Yeast System: Toxicity, Immunity, and Resistance. Yeast (Chichester, England) Chan, A., Hays, M., Sherlock, G. 2025

    Abstract

    Killer yeasts, such as the K1 killer strain of S. Cerevisiae, express a secreted anti-competitive toxin whose production and propagation require the presence of two vertically-transmitted dsRNA viruses. In sensitive cells lacking killer virus infection, toxin binding to the cell wall results in ion pore formation, disruption of osmotic homeostasis, and cell death. However, the exact mechanism(s) of K1 toxin killing activity, how killer yeasts are immune to their own toxin, and which factors could influence adaptation and resistance to K1 toxin within formerly sensitive populations are still unknown. Here, we describe the state of knowledge about K1 killer toxin, including current models of toxin processing and killing activity, and a summary of known modifiers of K1 toxin immunity and resistance. In addition, we discuss two key signaling pathways, HOG (high osmolarity glycerol) and CWI (cell wall integrity), whose involvement in an adaptive response to K1 killer toxin in sensitive cells has been previously documented but requires further study. As both host-virus and sensitive-killer competition have been documented in killer systems like K1, further characterization of K1 killer yeasts may provide a useful model system for study of both intracellular genetic conflict and counter-adaptation between competing sensitive and killer populations.

    View details for DOI 10.1002/yea.3987

    View details for PubMedID 39853823

  • Genetic conflicts in budding yeast: The 2mu plasmid as a model selfish element. Seminars in cell & developmental biology Hays, M. 2024; 161-162: 31-41

    Abstract

    Antagonistic coevolution, arising from genetic conflict, can drive rapid evolution and biological innovation. Conflict can arise both between organisms and within genomes. This review focuses on budding yeasts as a model system for exploring intra- and inter-genomic genetic conflict, highlighting in particular the 2-micron (2mu) plasmid as a model selfish element. The 2mu is found widely in laboratory strains and industrial isolates of Saccharomyces cerevisiae and has long been known to cause host fitness defects. Nevertheless, the plasmid is frequently ignored in the context of genetic, fitness, and evolution studies. Here, I make a case for further exploring the evolutionary impact of the 2mu plasmid as well as other selfish elements of budding yeasts, discuss recent advances, and, finally, future directions for the field.

    View details for DOI 10.1016/j.semcdb.2024.04.002

    View details for PubMedID 38598944

  • Paths to adaptation under fluctuating nitrogen starvation: The spectrum of adaptive mutations in Saccharomyces cerevisiae is shaped by retrotransposons and microhomology-mediated recombination. PLoS genetics Hays, M., Schwartz, K., Schmidtke, D. T., Aggeli, D., Sherlock, G. 2023; 19 (5): e1010747

    Abstract

    There are many mechanisms that give rise to genomic change: while point mutations are often emphasized in genomic analyses, evolution acts upon many other types of genetic changes that can result in less subtle perturbations. Changes in chromosome structure, DNA copy number, and novel transposon insertions all create large genomic changes, which can have correspondingly large impacts on phenotypes and fitness. In this study we investigate the spectrum of adaptive mutations that arise in a population under consistently fluctuating nitrogen conditions. We specifically contrast these adaptive alleles and the mutational mechanisms that create them, with mechanisms of adaptation under batch glucose limitation and constant selection in low, non-fluctuating nitrogen conditions to address if and how selection dynamics influence the molecular mechanisms of evolutionary adaptation. We observe that retrotransposon activity accounts for a substantial number of adaptive events, along with microhomology-mediated mechanisms of insertion, deletion, and gene conversion. In addition to loss of function alleles, which are often exploited in genetic screens, we identify putative gain of function alleles and alleles acting through as-of-yet unclear mechanisms. Taken together, our findings emphasize that how selection (fluctuating vs. non-fluctuating) is applied also shapes adaptation, just as the selective pressure (nitrogen vs. glucose) does itself. Fluctuating environments can activate different mutational mechanisms, shaping adaptive events accordingly. Experimental evolution, which allows a wider array of adaptive events to be assessed, is thus a complementary approach to both classical genetic screens and natural variation studies to characterize the genotype-to-phenotype-to-fitness map.

    View details for DOI 10.1371/journal.pgen.1010747

    View details for PubMedID 37192196

  • A natural variant of the essential host gene MMS21 restricts the parasitic 2-micron plasmid in Saccharomyces cerevisiae. eLife Hays, M. n., Young, J. M., Levan, P. F., Malik, H. S. 2020; 9

    Abstract

    Antagonistic coevolution with selfish genetic elements (SGEs) can drive evolution of host resistance. Here, we investigated host suppression of 2-micron (2m) plasmids, multicopy nuclear parasites that have co-evolved with budding yeasts. We developed SCAMPR (Single-Cell Assay for Measuring Plasmid Retention) to measure copy number heterogeneity and 2m plasmid loss in live cells. We identified three S. cerevisiae strains that lack endogenous 2m plasmids and reproducibly inhibit mitotic plasmid stability. Focusing on the Y9 ragi strain, we determined that plasmid restriction is heritable and dominant. Using bulk segregant analysis, we identified a high-confidence Quantitative Trait Locus (QTL) with a single variant of MMS21 associated with increased 2m instability. MMS21 encodes a SUMO E3 ligase and an essential component of the Smc5/6 complex, involved in sister chromatid cohesion, chromosome segregation, and DNA repair. Our analyses leverage natural variation to uncover a novel means by which budding yeasts can overcome highly successful genetic parasites.

    View details for DOI 10.7554/eLife.62337

    View details for PubMedID 33063663

  • Independent Origins of Yeast Associated with Coffee and Cacao Fermentation. Current biology : CB Ludlow, C. L., Cromie, G. A., Garmendia-Torres, C., Sirr, A., Hays, M., Field, C., Jeffery, E. W., Fay, J. C., Dudley, A. M. 2016; 26 (7): 965-71

    Abstract

    Modern transportation networks have facilitated the migration and mingling of previously isolated populations of plants, animals, and insects. Human activities can also influence the global distribution of microorganisms. The best-understood example is yeasts associated with winemaking. Humans began making wine in the Middle East over 9,000 years ago [1, 2]. Selecting favorable fermentation products created specialized strains of Saccharomyces cerevisiae [3, 4] that were transported along with grapevines. Today, S. cerevisiae strains residing in vineyards around the world are genetically similar, and their population structure suggests a common origin that followed the path of human migration [3-7]. Like wine, coffee and cacao depend on microbial fermentation [8, 9] and have been globally dispersed by humans. Theobroma cacao originated in the Amazon and Orinoco basins of Colombia and Venezuela [10], was cultivated in Central America by Mesoamerican peoples, and was introduced to Europeans by Hernán Cortés in 1530 [11]. Coffea, native to Ethiopia, was disseminated by Arab traders throughout the Middle East and North Africa in the 6(th) century and was introduced to European consumers in the 17(th) century [12]. Here, we tested whether the yeasts associated with coffee and cacao are genetically similar, crop-specific populations or genetically diverse, geography-specific populations. Our results uncovered populations that, while defined by niche and geography, also bear signatures of admixture between major populations in events independent of the transport of the plants. Thus, human-associated fermentation and migration may have affected the distribution of yeast involved in the production of coffee and chocolate.

    View details for DOI 10.1016/j.cub.2016.02.012

    View details for PubMedID 27020745

    View details for PubMedCentralID PMC4821677

  • Aneuploidy underlies a multicellular phenotypic switch. Proceedings of the National Academy of Sciences of the United States of America Tan, Z., Hays, M., Cromie, G. A., Jeffery, E. W., Scott, A. C., Ahyong, V., Sirr, A., Skupin, A., Dudley, A. M. 2013; 110 (30): 12367-72

    Abstract

    Although microorganisms are traditionally used to investigate unicellular processes, the yeast Saccharomyces cerevisiae has the ability to form colonies with highly complex, multicellular structures. Colonies with the "fluffy" morphology have properties reminiscent of bacterial biofilms and are easily distinguished from the "smooth" colonies typically formed by laboratory strains. We have identified strains that are able to reversibly toggle between the fluffy and smooth colony-forming states. Using a combination of flow cytometry and high-throughput restriction-site associated DNA tag sequencing, we show that this switch is correlated with a change in chromosomal copy number. Furthermore, the gain of a single chromosome is sufficient to switch a strain from the fluffy to the smooth state, and its subsequent loss to revert the strain back to the fluffy state. Because copy number imbalance of six of the 16 S. cerevisiae chromosomes and even a single gene can modulate the switch, our results support the hypothesis that the state switch is produced by dosage-sensitive genes, rather than a general response to altered DNA content. These findings add a complex, multicellular phenotype to the list of molecular and cellular traits known to be altered by aneuploidy and suggest that chromosome missegregation can provide a quick, heritable, and reversible mechanism by which organisms can toggle between phenotypes.

    View details for DOI 10.1073/pnas.1301047110

    View details for PubMedID 23812752

    View details for PubMedCentralID PMC3725063

  • Use of pleiotropy to model genetic interactions in a population. PLoS genetics Carter, G. W., Hays, M., Sherman, A., Galitski, T. 2012; 8 (10): e1003010

    Abstract

    Systems-level genetic studies in humans and model systems increasingly involve both high-resolution genotyping and multi-dimensional quantitative phenotyping. We present a novel method to infer and interpret genetic interactions that exploits the complementary information in multiple phenotypes. We applied this approach to a population of yeast strains with randomly assorted perturbations of five genes involved in mating. We quantified pheromone response at the molecular level and overall mating efficiency. These phenotypes were jointly analyzed to derive a network of genetic interactions that mapped mating-pathway relationships. To determine the distinct biological processes driving the phenotypic complementarity, we analyzed patterns of gene expression to find that the pheromone response phenotype is specific to cellular fusion, whereas mating efficiency was a combined measure of cellular fusion, cell cycle arrest, and modifications in cellular metabolism. We applied our novel method to global gene expression patterns to derive an expression-specific interaction network and demonstrate applicability to global transcript data. Our approach provides a basis for interpretation of genetic interactions and the generation of specific hypotheses from populations assayed for multiple phenotypes.

    View details for DOI 10.1371/journal.pgen.1003010

    View details for PubMedID 23071457

    View details for PubMedCentralID PMC3469415

  • Spatiotemporal patterns of gene expression during development of a complex colony morphology. PloS one Cromie, G. A., Tan, Z., Hays, M., Sirr, A., Dudley, A. M. 2024; 19 (12): e0311061

    Abstract

    Clonal communities of single celled organisms, such as bacterial or fungal colonies and biofilms, are spatially structured, with subdomains of cells experiencing differing environmental conditions. In the development of such communities, cell specialization is not only important to respond and adapt to the local environment but has the potential to increase the fitness of the clonal community through division of labor. Here, we examine colony development in a yeast strain (F13) that produces colonies with a highly structured "ruffled" phenotype in the colony periphery and an unstructured "smooth" phenotype in the colony center. We demonstrate that in the F13 genetic background deletions of transcription factors can either increase (dig1D, sfl1D) or decrease (tec1D) the degree of colony structure. To investigate the development of colony structure, we carried out gene expression analysis on F13 and the three deletion strains using RNA-seq. Samples were taken early in colony growth (day2), which precedes ruffled phenotype development in F13, and from the peripheral and central regions of colonies later in development (day5), at which time these regions are structured and unstructured (respectively) in F13. We identify genes responding additively and non-additively to the genotype and spatiotemporal factors and cluster these genes into a number of different expression patterns. We identify clusters whose expression correlates closely with the degree of colony structure in each sample and include genes with known roles in the development of colony structure. Individual deletion of 26 genes sampled from different clusters identified 5 with strong effects on colony morphology (BUD8, CIS3, FLO11, MSB2 and SFG1), all of which eliminated or greatly reduced the structure of the F13 outer region.

    View details for DOI 10.1371/journal.pone.0311061

    View details for PubMedID 39637084

  • Data-driven multiscale modeling reveals the role of metabolic coupling for the spatio-temporal growth dynamics of yeast colonies. BMC molecular and cell biology Intosalmi, J., Scott, A. C., Hays, M., Flann, N., Yli-Harja, O., Lähdesmäki, H., Dudley, A. M., Skupin, A. 2019; 20 (1): 59

    Abstract

    Multicellular entities like mammalian tissues or microbial biofilms typically exhibit complex spatial arrangements that are adapted to their specific functions or environments. These structures result from intercellular signaling as well as from the interaction with the environment that allow cells of the same genotype to differentiate into well-organized communities of diversified cells. Despite its importance, our understanding how this cell-cell and metabolic coupling lead to functionally optimized structures is still limited.Here, we present a data-driven spatial framework to computationally investigate the development of yeast colonies as such a multicellular structure in dependence on metabolic capacity. For this purpose, we first developed and parameterized a dynamic cell state and growth model for yeast based on on experimental data from homogeneous liquid media conditions. The inferred model is subsequently used in a spatially coarse-grained model for colony development to investigate the effect of metabolic coupling by calibrating spatial parameters from experimental time-course data of colony growth using state-of-the-art statistical techniques for model uncertainty and parameter estimations. The model is finally validated by independent experimental data of an alternative yeast strain with distinct metabolic characteristics and illustrates the impact of metabolic coupling for structure formation.We introduce a novel model for yeast colony formation, present a statistical methodology for model calibration in a data-driven manner, and demonstrate how the established model can be used to generate predictions across scales by validation against independent measurements of genetically distinct yeast strains.

    View details for DOI 10.1186/s12860-019-0234-z

    View details for PubMedID 31856706

    View details for PubMedCentralID PMC6923950

  • Transcriptional Profiling of Biofilm Regulators Identified by an Overexpression Screen in Saccharomyces cerevisiae G3-GENES GENOMES GENETICS Cromie, G. A., Tan, Z., Hays, M., Sirr, A., Jeffery, E. W., Dudley, A. M. 2017; 7 (8): 2845–54

    Abstract

    Biofilm formation by microorganisms is a major cause of recurring infections and removal of biofilms has proven to be extremely difficult given their inherent drug resistance . Understanding the biological processes that underlie biofilm formation is thus extremely important and could lead to the development of more effective drug therapies, resulting in better infection outcomes. Using the yeast Saccharomyces cerevisiae as a biofilm model, overexpression screens identified DIG1, SFL1, HEK2, TOS8, SAN1, and ROF1/YHR177W as regulators of biofilm formation. Subsequent RNA-seq analysis of biofilm and nonbiofilm-forming strains revealed that all of the overexpression strains, other than DIG1 and TOS8, were adopting a single differential expression profile, although induced to varying degrees. TOS8 adopted a separate profile, while the expression profile of DIG1 reflected the common pattern seen in most of the strains, plus substantial DIG1-specific expression changes. We interpret the existence of the common transcriptional pattern seen across multiple, unrelated overexpression strains as reflecting a transcriptional state, that the yeast cell can access through regulatory signaling mechanisms, allowing an adaptive morphological change between biofilm-forming and nonbiofilm states.

    View details for DOI 10.1534/g3.117.042440

    View details for Web of Science ID 000407314000038

    View details for PubMedID 28673928

    View details for PubMedCentralID PMC5555487

  • Dissecting Gene Expression Changes Accompanying a Ploidy-Based Phenotypic Switch. G3 (Bethesda, Md.) Cromie, G. A., Tan, Z., Hays, M., Jeffery, E. W., Dudley, A. M. 2017; 7 (1): 233-246

    Abstract

    Aneuploidy, a state in which the chromosome number deviates from a multiple of the haploid count, significantly impacts human health. The phenotypic consequences of aneuploidy are believed to arise from gene expression changes associated with the altered copy number of genes on the aneuploid chromosomes. To dissect the mechanisms underlying altered gene expression in aneuploids, we used RNA-seq to measure transcript abundance in colonies of the haploid Saccharomyces cerevisiae strain F45 and two aneuploid derivatives harboring disomies of chromosomes XV and XVI. F45 colonies display complex "fluffy" morphologies, while the disomic colonies are smooth, resembling laboratory strains. Our two disomes displayed similar transcriptional profiles, a phenomenon not driven by their shared smooth colony morphology nor simply by their karyotype. Surprisingly, the environmental stress response (ESR) was induced in F45, relative to the two disomes. We also identified genes whose expression reflected a nonlinear interaction between the copy number of a transcriptional regulatory gene on chromosome XVI, DIG1, and the copy number of other chromosome XVI genes. DIG1 and the remaining chromosome XVI genes also demonstrated distinct contributions to the effect of the chromosome XVI disome on ESR gene expression. Expression changes in aneuploids appear to reflect a mixture of effects shared between different aneuploidies and effects unique to perturbing the copy number of particular chromosomes, including nonlinear copy number interactions between genes. The balance between these two phenomena is likely to be genotype- and environment-specific.

    View details for DOI 10.1534/g3.116.036160

    View details for PubMedID 27836908

    View details for PubMedCentralID PMC5217112

  • Identification and characterization of a drug-sensitive strain enables puromycin-based translational assays in Saccharomyces cerevisiae. Yeast (Chichester, England) Cary, G. A., Yoon, S. H., Torres, C. G., Wang, K., Hays, M., Ludlow, C., Goodlett, D. R., Dudley, A. M. 2014; 31 (5): 167-78

    Abstract

    Puromycin is an aminonucleoside antibiotic with structural similarity to aminoacyl tRNA. This structure allows the drug to bind the ribosomal A site and incorporate into nascent polypeptides, causing chain termination, ribosomal subunit dissociation and widespread translational arrest at high concentrations. In contrast, at sufficiently low concentrations, puromycin incorporates primarily at the C-terminus of proteins. While a number of techniques utilize puromycin incorporation as a tool for probing translational activity in vivo, these methods cannot be applied in yeasts that are insensitive to puromycin. Here, we describe a mutant strain of the yeast Saccharomyces cerevisiae that is sensitive to puromycin and characterize the cellular response to the drug. Puromycin inhibits the growth of yeast cells mutant for erg6∆, pdr1∆ and pdr3∆ (EPP) on both solid and liquid media. Puromycin also induces the aggregation of the cytoplasmic processing body component Edc3 in the mutant strain. We establish that puromycin is rapidly incorporated into yeast proteins and test the effects of puromycin on translation in vivo. This study establishes the EPP strain as a valuable tool for implementing puromycin-based assays in yeast, which will enable new avenues of inquiry into protein production and maturation.

    View details for DOI 10.1002/yea.3007

    View details for PubMedID 24610064

    View details for PubMedCentralID PMC4013229

  • High-throughput tetrad analysis. Nature methods Ludlow, C. L., Scott, A. C., Cromie, G. A., Jeffery, E. W., Sirr, A., May, P., Lin, J., Gilbert, T. L., Hays, M., Dudley, A. M. 2013; 10 (7): 671-5

    Abstract

    Tetrad analysis has been a gold-standard genetic technique for several decades. Unfortunately, the need to manually isolate, disrupt and space tetrads has relegated its application to small-scale studies and limited its integration with high-throughput DNA sequencing technologies. We have developed a rapid, high-throughput method, called barcode-enabled sequencing of tetrads (BEST), that uses (i) a meiosis-specific GFP fusion protein to isolate tetrads by FACS and (ii) molecular barcodes that are read during genotyping to identify spores derived from the same tetrad. Maintaining tetrad information allows accurate inference of missing genetic markers and full genotypes of missing (and presumably nonviable) individuals. An individual researcher was able to isolate over 3,000 yeast tetrads in 3 h, an output equivalent to that of almost 1 month of manual dissection. BEST is transferable to other microorganisms for which meiotic mapping is significantly more laborious.

    View details for DOI 10.1038/nmeth.2479

    View details for PubMedID 23666411

    View details for PubMedCentralID PMC3696418

  • Predicting the effects of copy-number variation in double and triple mutant combinations. Pacific Symposium on Biocomputing. Pacific Symposium on Biocomputing Carter, G. W., Hays, M., Li, S., Galitski, T. 2012: 19-30

    Abstract

    The study of genetic interactions is a powerful tool in inferring structure and function of biological networks. To date, genetic interaction studies have been dominated by pair-wise gene deletion screens. However, classical genetic analysis and natural genetic variation involve diverse gene forms ranging from null alleles to copy number variants. Moreover, genetic variation is typically multifactorial. Addressing multiple combinatorial genetic variations ranging in gene activity is therefore of critical value. We approach this problem using genetic network modeling that quantitatively encodes how genes influence the activity of one another and phenotype outcomes. A network model was initially inferred from linear decomposition of gene expression data. We used this network to predict the effects of combining multi-copy and deletion mutations of specific gene pairs and a gene triplet. Predicted expression patterns across hundreds of genes were experimentally validated. Prediction success was critically dependent on how a multi-copy gene interacted with other genes in the network model. This strategy provides a template for the inference, prediction, and testing of genetically complex hypotheses involving diverse genetic variation.

    View details for PubMedID 22174259

    View details for PubMedCentralID PMC3334851