Genome architecture and stability in the Saccharomyces cerevisiae knockout collection

The mutagenic consequences of defective DNA repair

Multiple separate repair mechanisms safeguard the genome against various types of DNA damage, and their failure can increase the rate of spontaneous mutagenesis. The malfunction of distinct repair mechanisms leads to genomic instability through different mutagenic processes. For example, defective mismatch repair causes high base substitution rates and microsatellite instability, whereas homologous recombination deficiency is characteristically associated with deletions and chromosome instability. This review presents a comprehensive collection of all mutagenic phenotypes associated with the loss of each DNA repair mechanism, drawing on data from a variety of model organisms and mutagenesis assays, and placing greatest emphasis on systematic analyses of human cancer datasets. We describe the latest theories on the mechanism of each mutagenic process, often explained by reliance on an alternative repair pathway or the error-prone replication of unrepaired, damaged DNA. Aided by the concept of mutational signatures, the genomic phenotypes can be used in cancer diagnosis to identify defective DNA repair pathways.

Efficient biosynthesis of active hemoglobins through enhancing the import of heme in Saccharomyces cerevisiae

Hemoglobins, with heme as a cofactor, are functional proteins that have extensive applications in the fields of artificial oxygen carriers and foods. Although Saccharomyces cerevisiae is an ideal host for hemoglobin synthesis, it lacks a suitable transport system to utilize additional heme for active expression of hemoglobins, resulting in the cellular aggregation and degradation of the latter. Here, an effective heme importer, heme-responsive gene 4 (Hrg-4), was selected from six candidates through the comparison of effects on the growth rates of Δhem1 S. cerevisiae strain and the activities of various hemoglobins when supplemented with 5 mg·L-1 exogenous heme. Additionally, to counter the instability of plasmid-based expression and the metabolic burden introduced from overexpressing Hrg-4, a series of hrg-4 integrated strains were constructed and the best engineered strain with five copies of hrg-4 was chosen. We found that this engineered strain was associated with an increased binding rate of heme in monomeric leghemoglobin and multimeric human hemoglobin (76.3% and 16.5%, respectively), as well as an enhanced expression of both hemoglobins (52.8% and 17.0%, respectively). Thus, the engineered strain with improved heme uptake can be used to efficiently synthesize other heme-binding proteins and enzymes in S. cerevisiae.

Loss of mitochondrial DNA is associated with reduced DNA content variability in Saccharomyces cerevisiae

DNA content measurement by fluorescence-assisted cell sorting (FACS) provides information on cell cycle progression and DNA content variability. Saccharomyces cerevisiae mutants with DNA content variability that was reduced relative to wild-type strains had defects in mitochondrial DNA (mtDNA) maintenance and mitochondrial gene expression and were correlated with strains found to lack mtDNA ([ rho 0 ] cells) by genome sequencing and fluorescence microscopy. In contrast, mutants with increased variability had defects in cell cycle progression, which may indicate a loss of coordination between mtDNA and nuclear DNA replication. Thus, FACS measurement of DNA content variability can provide insight into cell-to-cell heterogeneity in mtDNA copy number.

Spontaneous Mutations in Saccharomyces cerevisiae mtDNA Increase Cell-to-Cell Variation in mtDNA Amount

In a eukaryotic cell, the ratio of mitochondrial DNA (mtDNA) to nuclear DNA (nDNA) is usually maintained within a specific range. This suggests the presence of a negative feedback loop mechanism preventing extensive mtDNA replication and depletion. However, the experimental data on this hypothetical mechanism are limited. In this study, we suggested that deletions in mtDNA, known to increase mtDNA abundance, can disrupt this mechanism, and thus, increase cell-to-cell variance in the mtDNA copy numbers. To test this, we generated Saccharomyces cerevisiae rho- strains with large deletions in the mtDNA and rho0 strains depleted of mtDNA. Given that mtDNA contributes to the total DNA content of exponentially growing yeast cells, we showed that it can be quantified in individual cells by flow cytometry using the DNA-intercalating fluorescent dye SYTOX green. We found that the rho- mutations increased both the levels and cell-to-cell heterogeneity in the total DNA content of G1 and G2/M yeast cells, with no association with the cell size. Furthermore, the depletion of mtDNA in both the rho+ and rho- strains significantly decreased the SYTOX green signal variance. The high cell-to-cell heterogeneity of the mtDNA amount in the rho- strains suggests that mtDNA copy number regulation relies on full-length mtDNA, whereas the rho- mtDNAs partially escape this regulation.

GenoPipe: identifying the genotype of origin within (epi)genomic datasets

Confidence in experimental results is critical for discovery. As the scale of data generation in genomics has grown exponentially, experimental error has likely kept pace despite the best efforts of many laboratories. Technical mistakes can and do occur at nearly every stage of a genomics assay (i.e. cell line contamination, reagent swapping, tube mislabelling, etc.) and are often difficult to identify post-execution. However, the DNA sequenced in genomic experiments contains certain markers (e.g. indels) encoded within and can often be ascertained forensically from experimental datasets. We developed the Genotype validation Pipeline (GenoPipe), a suite of heuristic tools that operate together directly on raw and aligned sequencing data from individual high-throughput sequencing experiments to characterize the underlying genome of the source material. We demonstrate how GenoPipe validates and rescues erroneously annotated experiments by identifying unique markers inherent to an organism's genome (i.e. epitope insertions, gene deletions and SNPs).

Role of Mitochondrial DNA in Yeast Replicative Aging

Despite the diverse manifestations of aging across different species, some common aging features and underlying mechanisms are shared. In particular, mitochondria appear to be among the most vulnerable systems in both metazoa and fungi. In this review, we discuss how mitochondrial dysfunction is related to replicative aging in the simplest eukaryotic model, the baker's yeast Saccharomyces cerevisiae. We discuss a chain of events that starts from asymmetric distribution of mitochondria between mother and daughter cells. With age, yeast mother cells start to experience a decrease in mitochondrial transmembrane potential and, consequently, a decrease in mitochondrial protein import efficiency. This induces mitochondrial protein precursors in the cytoplasm, the loss of mitochondrial DNA (mtDNA), and at the later stages - cell death. Interestingly, yeast strains without mtDNA can have either increased or decreased lifespan compared to the parental strains with mtDNA. The direction of the effect depends on their ability to activate compensatory mechanisms preventing or mitigating negative consequences of mitochondrial dysfunction. The central role of mitochondria in yeast aging and death indicates that it is one of the most complex and, therefore, deregulation-prone systems in eukaryotic cells.

Genome instability footprint under rapamycin and hydroxyurea treatments

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit

Blue mold, a postharvest disease of pome fruits, is caused by the filamentous fungus Penicillium expansum. In addition to the economic losses caused by P. expansum, food safety can be compromised, as this pathogen is mycotoxigenic. In this study, forward and reverse genetic approaches were used to identify genes involved in blue mold infection in apple fruits. For this, we generated a random T-DNA insertional mutant library. A total of 448 transformants were generated and screened for the reduced decay phenotype on apples. Of these mutants, six (T-193, T-275, T-434, T-588, T-625, and T-711) were selected for continued studies and five unique genes were identified of interest. In addition, two deletion mutants (Δt-625 and Δt-588) and a knockdown strain (t-434KD) were generated for three loci. Data show that the ∆t-588 mutant phenocopied the T-DNA insertion mutant and had virulence penalties during apple fruit decay. We hypothesize that this locus encodes a glyoxalase due to bioinformatic predictions, thus contributing to reduced colony diameter when grown in methylglyoxal (MG). This work presents novel members of signaling networks and additional genetic factors that regulate fungal virulence in the blue mold fungus during apple fruit decay.

Identification of a common secondary mutation in the Neurospora crassa knockout collection conferring a cell fusion-defective phenotype

Gene-deletion mutants represent a powerful tool to study gene function. The filamentous fungus Neurospora crassa is a well-established model organism, and features a comprehensive gene knockout strain collection. While these mutant strains have been used in numerous studies, resulting in the functional annotation of many Neurospora genes, direct confirmation of gene-phenotype relationships is often lacking, which is particularly relevant given the possibility of background mutations, sample contamination, and/or strain mislabeling. Indeed, spontaneous mutations resulting in phenotypes resembling many cell fusion mutants have long been known to occur at relatively high frequency in N. crassa, and these secondary mutations are common in the Neurospora deletion collection. The identity of these mutations, however, is largely unknown. Here, we report that the Δada-3 strain from the N. crassa knockout collection, which exhibits a cell fusion defect, harbors a secondary mutation responsible for this phenotype. Through whole-genome sequencing and genetic analyses, we found a ~30-Kb deletion in this strain affecting a known cell fusion-related gene, so/ham-1, and show that it is the absence of this gene-and not of ada-3-that underlies its cell fusion defect. We additionally found three other knockout strains harboring the same deletion, suggesting that this mutation may be common in the collection and could have impacted previous studies. Our findings provide a cautionary note and highlight the importance of proper functional validation of strains from mutant collections. We discuss our results in the context of the spread of cell fusion-defective cheater variants in N. crassa cultures. IMPORTANCE This study emphasizes the need for careful and detailed characterization of strains from mutant collections. Specifically, we found a common deletion in various strains from the Neurospora crassa gene knockout collection that results in a cell fusion-defective phenotype. This is noteworthy because this collection is known to contain background mutations-of a largely unclear nature-that produce cell fusion-defective phenotypes. Our results describe an example of such mutations, and highlight how this common genetic defect could have impacted previous studies that have used the affected strains. Furthermore, they provide a cautionary note about the use of Neurospora strains with similar phenotypes. Lastly, these findings offer additional details relevant to our understanding of the origin and spread of cell fusion-defective cheater variants in N. crassa cultures.

Theory for High-Throughput Genetic Interaction Screening

Systematic, genome-scale genetic screens have been instrumental for elucidating genotype-phenotype relationships, but approaches for probing genetic interactions have been limited to at most ∼100 pre-selected gene combinations in mammalian cells. Here, we introduce a theory for high-throughput genetic interaction screens. The theory extends our recently developed Multiplexing using Spectral Imaging and Combinatorics (MuSIC) approach to propose ∼105 spectrally unique, genetically encoded MuSIC barcodes from 18 currently available fluorescent proteins. Simulation studies based on constraints imposed by spectral flow cytometry equipment suggest that genetic interaction screens at the human genome-scale may be possible if MuSIC barcodes can be paired to guide RNAs. While experimental testing of this theory awaits, it offers transformative potential for genetic perturbation technology and knowledge of genetic function. More broadly, the availability of a genome-scale spectral barcode library for non-destructive identification of single cells could find more widespread applications such as traditional genetic screening and high-dimensional lineage tracing.

Condition-dependent fitness effects of large synthetic chromosome amplifications

Whole-chromosome aneuploidy and large segmental amplifications can have devastating effects in multicellular organisms, from developmental disorders and miscarriage to cancer. Aneuploidy in single-celled organisms such as yeast also results in proliferative defects and reduced viability. Yet, paradoxically, CNVs are routinely observed in laboratory evolution experiments with microbes grown in stressful conditions. The defects associated with aneuploidy are often attributed to the imbalance of many differentially expressed genes on the affected chromosomes, with many genes each contributing incremental effects. An alternate hypothesis is that a small number of individual genes are large effect 'drivers' of these fitness changes when present in an altered copy number. To test these two views, we have employed a collection of strains bearing large chromosomal amplifications that we previously assayed in nutrient-limited chemostat competitions. In this study, we focus on conditions known to be poorly tolerated by aneuploid yeast-high temperature, treatment with the Hsp90 inhibitor radicicol, and growth in extended stationary phase. To identify potential genes with a large impact on fitness, we fit a piecewise constant model to fitness data across chromosome arms, filtering breakpoints in this model by magnitude to focus on regions with a large impact on fitness in each condition. While fitness generally decreased as the length of the amplification increased, we were able to identify 91 candidate regions that disproportionately impacted fitness when amplified. Consistent with our previous work with this strain collection, nearly all candidate regions were condition specific, with only five regions impacting fitness in multiple conditions.

Experimental evolution of Saccharomyces uvarum at high temperature yields elevation of maximal growth temperature and loss of the mitochondrial genome

An organism's upper thermal tolerance is a major driver of its ecology and is a complex, polygenic trait. Given the wide variance in this critical phenotype across the tree of life, it is quite striking that this trait has not proven very evolutionarily labile in experimental evolution studies of microbes. In stark contrast to recent studies, William Henry Dallinger in the 1880s reported increasing the upper thermal limit of microbes he experimentally evolved by >40°C using a very gradual temperature ramping strategy. Using a selection scheme inspired by Dallinger, we sought to increase the upper thermal limit of Saccharomyces uvarum . This species has a maximum growth temperature of 34-35°C, considerably lower than S. cerevisiae . After 136 passages on solid plates at progressively higher temperatures, we recovered a clone that can grow at 36°C, a gain of ~1.5°C. Additionally, the evolved clone lost its mitochondrial genome and cannot respire. In contrast, an induced rho 0 derivative of the ancestor shows a decrease in thermotolerance. Also, incubation of the ancestor at 34°C for 5 days increased the frequency of petite mutants drastically compared to 22°C, supporting the notion that mutation pressure rather than selection drove loss of mtDNA in the evolved clone. These results demonstrate that S. uvarum 's upper thermal limit can be elevated slightly via experimental evolution and corroborate past observations in S. cerevisiae that high temperature selection schemes can inadvertently lead to production of the potentially undesirable respiratory incompetent phenotype in yeasts.

Global analysis of the yeast knockout phenome

Genome-wide phenotypic screens in the budding yeast Saccharomyces cerevisiae, enabled by its knockout collection, have produced the largest, richest, and most systematic phenotypic description of any organism. However, integrative analyses of this rich data source have been virtually impossible because of the lack of a central data repository and consistent metadata annotations. Here, we describe the aggregation, harmonization, and analysis of ~14,500 yeast knockout screens, which we call Yeast Phenome. Using this unique dataset, we characterized two unknown genes (YHR045W and YGL117W) and showed that tryptophan starvation is a by-product of many chemical treatments. Furthermore, we uncovered an exponential relationship between phenotypic similarity and intergenic distance, which suggests that gene positions in both yeast and human genomes are optimized for function.

Overexpression profiling reveals cellular requirements in the context of genetic backgrounds and environments

Overexpression can help life adapt to stressful environments, making an examination of overexpressed genes valuable for understanding stress tolerance mechanisms. However, a systematic study of genes whose overexpression is functionally adaptive (GOFAs) under stress has yet to be conducted. We developed a new overexpression profiling method and systematically identified GOFAs in Saccharomyces cerevisiae under stress (heat, salt, and oxidative). Our results show that adaptive overexpression compensates for deficiencies and increases fitness under stress, like calcium under salt stress. We also investigated the impact of different genetic backgrounds on GOFAs, which varied among three S. cerevisiae strains reflecting differing calcium and potassium requirements for salt stress tolerance. Our study of a knockout collection also suggested that calcium prevents mitochondrial outbursts under salt stress. Mitochondria-enhancing GOFAs were only adaptive when adequate calcium was available and non-adaptive when calcium was deficient, supporting this idea. Our findings indicate that adaptive overexpression meets the cell's needs for maximizing the organism's adaptive capacity in the given environment and genetic context.

The proteomic landscape of genome-wide genetic perturbations

Functional genomic strategies have become fundamental for annotating gene function and regulatory networks. Here, we combined functional genomics with proteomics by quantifying protein abundances in a genome-scale knockout library in Saccharomyces cerevisiae, using data-independent acquisition mass spectrometry. We find that global protein expression is driven by a complex interplay of (1) general biological properties, including translation rate, protein turnover, the formation of protein complexes, growth rate, and genome architecture, followed by (2) functional properties, such as the connectivity of a protein in genetic, metabolic, and physical interaction networks. Moreover, we show that functional proteomics complements current gene annotation strategies through the assessment of proteome profile similarity, protein covariation, and reverse proteome profiling. Thus, our study reveals principles that govern protein expression and provides a genome-spanning resource for functional annotation.

Polymorphism of Saccharomyces cerevisiae Strains in DNA Metabolism Genes

Baker's yeast, S. cerevisiae, is an excellent model organism exploited for molecular genetic studies of the mechanisms of genome stability in eukaryotes. Genetic peculiarities of commonly used yeast strains impact the processes of DNA replication, repair, and recombination (RRR). We compared the genomic DNA sequence variation of the five strains that are intensively used for RRR studies. We used yeast next-generation sequencing data to detect the extent and significance of variation in 183 RRR genes. We present a detailed analysis of the differences that were found even in closely related strains. Polymorphisms of common yeast strains should be considered when interpreting the outcomes of genome stability studies, especially in cases of discrepancies between laboratories describing the same phenomena.

GenoPipe: identifying the genotype of origin within (epi)genomic datasets

Confidence in experimental results is critical for discovery. As the scale of data generation in genomics has grown exponentially, experimental error has likely kept pace despite the best efforts of many laboratories. Technical mistakes can and do occur at nearly every stage of a genomics assay (i.e., cell line contamination, reagent swapping, tube mislabelling, etc.) and are often difficult to identify post-execution. However, the DNA sequenced in genomic experiments contains certain markers (e.g., indels) encoded within and can often be ascertained forensically from experimental datasets. We developed the Genotype validation Pipeline (GenoPipe), a suite of heuristic tools that operate together directly on raw and aligned sequencing data from individual high-throughput sequencing experiments to characterize the underlying genome of the source material. We demonstrate how GenoPipe validates and rescues erroneously annotated experiments by identifying unique markers inherent to an organism’s genome (i.e., epitope insertions, gene deletions, and SNPs).

Spontaneous Mutation Rates and Spectra of Respiratory-Deficient Yeast

The yeast petite mutant was first discovered in the yeast Saccharomyces cerevisiae, which shows growth stress due to defects in genes encoding the respiratory chain. In a previous study, we described that deletion of the nuclear-encoded gene MRPL25 leads to mitochondrial genome (mtDNA) loss and the petite phenotype, which can be rescued by acquiring ATP3 mutations. The mrpl25Δ strain showed an elevated SNV (single nucleotide variant) rate, suggesting genome instability occurred during the crisis of mtDNA loss. However, the genome-wide mutation landscape and mutational signatures of mitochondrial dysfunction are unknown. In this study we profiled the mutation spectra in yeast strains with the genotype combination of MRPL25 and ATP3 in their wildtype and mutated status, along with the wildtype and cytoplasmic petite rho0 strains as controls. In addition to the previously described elevated SNV rate, we found the INDEL (insertion/deletion) rate also increased in the mrpl25Δ strain, reinforcing the occurrence of genome instability. Notably, although both are petites, the mrpl25Δ and rho0 strains exhibited different INDEL rates and transition/transversion ratios, suggesting differences in the mutational signatures underlying these two types of petites. Interestingly, the petite-related mutagenesis effect disappeared when ATP3 suppressor mutations were acquired, suggesting a cost-effective mechanism for restoring both fitness and genome stability. Taken together, we present an unbiased genome-wide characterization of the mutation rates and spectra of yeast strains with respiratory deficiency, which provides valuable insights into the impact of respiratory deficiency on genome instability.

Telomeres are shorter in wild Saccharomyces cerevisiae isolates than in domesticated ones

Telomeres are ribonucleoproteins that cap chromosome-ends and their DNA length is controlled by counteracting elongation and shortening processes. The budding yeast Saccharomyces cerevisiae has been a leading model to study telomere DNA length control and dynamics. Its telomeric DNA is maintained at a length that slightly varies between laboratory strains, but little is known about its variation at the species level. The recent publication of the genomes of over 1,000 S. cerevisiae strains enabled us to explore telomere DNA length variation at an unprecedented scale. Here, we developed a bioinformatic pipeline (YeaISTY) to estimate telomere DNA length from whole-genome sequences and applied it to the sequenced S. cerevisiae collection. Our results revealed broad natural telomere DNA length variation among the isolates. Notably, telomere DNA length is shorter in those derived from wild rather than domesticated environments. Moreover, telomere DNA length variation is associated with mitochondrial metabolism, and this association is driven by wild strains. Overall, these findings reveal broad variation in budding yeast's telomere DNA length regulation, which might be shaped by its different ecological life-styles.

Point mutations of homologs as an adaptive solution to the gene loss

Gene loss is common and influences genome evolution trajectories. Multiple adaptive strategies to compensate for gene loss have been observed, including copy number gain of paralogous genes and mutations in genes of the same pathway. By using the Ubl-specific protease 2 (ULP2) eviction model, we identify compensatory mutations in the homologous gene ULP1 by laboratory evolution and find that these mutations are capable of rescuing defects caused by the loss of ULP2. Furthermore, bioinformatics analysis of genomes of yeast gene knockout library and natural yeast isolate datasets suggests that point mutations of a homologous gene might be an additional mechanism to compensate gene loss.

OGG1 inhibition suppresses African swine fever virus replication

African swine fever virus (ASFV) is an important pathogen that causes a highly contagious and lethal disease in swine, for which neither a vaccine nor treatment is available. The DNA repair enzyme 8-oxoguanine DNA glycosylase 1 (OGG1), which excises the oxidative base lesion 8-oxo-7,8-dihydroguanine (8-oxoG), has been linked to the pathogenesis of different diseases associated with viral infections. However, the role of OGG1-base excision repair (BER) in ASFV infection has been poorly investigated. Our study aimed to characterize the alteration of host reactive oxygen species (ROS) and OGG1 and to analyse the role of OGG1 in ASFV infection. We found that ASFV infection induced high levels and dynamic changes in ROS and 8-oxoG and consistently increased the expression of OGG1. Viral yield, transcription level, and protein synthesis were reduced in ASFV-infected primary alveolar macrophages (PAMs) treated by TH5487 or SU0268 inhibiting OGG1. The expression of BER pathway associated proteins of ASFV was also suppressed in OGG1-inhibited PAMs. Furthermore, OGG1 was found to negatively regulate interferon β (IFN-β) production during ASFV infection and IFN-β could be activated by OGG1 inhibition with TH5487 and SU0268, which blocked OGG1 binding to 8-oxoG. Additionally, the interaction of OGG1 with viral MGF360-14-L protein could disturb IFN-β production to further affect ASFV replication. These results suggest that OGG1 plays the crucial role in successful viral infection and OGG1 inhibitors SU0268 or TH5487 could be used as antiviral agents for ASFV infection.

Biocontainment strategies for in vivo applications of Saccharomyces boulardii

The human gastrointestinal tract is a complex and dynamic environment, playing a crucial role in human health. Microorganisms engineered to express a therapeutic activity have emerged as a novel modality to manage numerous diseases. Such advanced microbiome therapeutics (AMTs) must be contained within the treated individual. Hence safe and robust biocontainment strategies are required to prevent the proliferation of microbes outside the treated individual. Here we present the first biocontainment strategy for a probiotic yeast, demonstrating a multi-layered strategy combining an auxotrophic and environmental-sensitive strategy. We knocked out the genes THI6 and BTS1, causing thiamine auxotrophy and increased sensitivity to cold, respectively. The biocontained Saccharomyces boulardii showed restricted growth in the absence of thiamine above 1 ng/ml and exhibited a severe growth defect at temperatures below 20°C. The biocontained strain was well tolerated and viable in mice and demonstrated equal efficiency in peptide production as the ancestral non-biocontained strain. In combination, the data support that thi6∆ and bts1∆ enable biocontainment of S. boulardii, which could be a relevant chassis for future yeast-based AMTs.

Fungal Strains with Identical Genomes Were Found at a Distance of 2000 Kilometers after 40 Years

Heredity and variation are inherent characteristics of species and are mainly reflected in the stability and variation of the genome; the former is relative, while the latter is continuous. However, whether life has both stable genomes and extremely diverse genomes at the same time is unknown. In this study, we isolated Sclerotinia sclerotiorum strains from sclerotium samples in Quincy, Washington State, USA, and found that four single-sclerotium-isolation strains (PB4, PB273, PB615, and PB623) had almost identical genomes to the reference strain 1980 isolated in the west of Nebraska 40 years ago. The genome of strain PB4 sequenced by the next-generation sequencing (NGS) and Pacific Biosciences (PacBio) sequencing carried only 135 single nucleotide polymorphisms (SNPs) and 18 structural variations (SVs) compared with the genome of strain 1980 and 48 SNPs were distributed on Contig_20. Based on data generated by NGS, three other strains, PB273, PB615, and PB623, had 256, 275, and 262 SNPs, respectively, against strain 1980, which were much less than in strain PB4 (532 SNPs) and none of them occurred on Contig_20, suggesting much closer genomes to strain 1980 than to strain PB4. All other strains from America and China are rich in SNPs with a range of 34,391-77,618 when compared with strain 1980. We also found that there were 39-79 SNPs between strain PB4 and its sexual offspring, 53.1% of which also occurred on Contig_20. Our discoveries show that there are two types of genomes in S. sclerotiorum, one is very stable and the other tends to change constantly. Investigating the mechanism of such genome stability will enhance our understanding of heredity and variation.

A comparative analysis of telomere length maintenance circuits in fission and budding yeast

The natural ends of the linear eukaryotic chromosomes are protected by telomeres, which also play an important role in aging and cancer development. Telomere length varies between species, but it is strictly controlled in all organisms. The process of Telomere Length Maintenance (TLM) involves many pathways, protein complexes and interactions that were first discovered in budding and fission yeast model organisms (Saccharomyces cerevisiae, Schizosaccharomyces pombe). In particular, large-scale systematic genetic screens in budding yeast uncovered a network of ≈500 genes that, when mutated, cause telomeres to lengthen or to shorten. In contrast, the TLM network in fission yeast remains largely unknown and systematic data is still lacking. In this work we try to close this gap and develop a unified interpretable machine learning framework for TLM gene discovery and phenotype prediction in both species. We demonstrate the utility of our framework in pinpointing the pathways by which TLM homeostasis is maintained and predicting novel TLM genes in fission yeast. The results of this study could be used for better understanding of telomere biology and serve as a step towards the adaptation of computational methods based on telomeric data for human prognosis.

Genetically controlled mtDNA deletions prevent ROS damage by arresting oxidative phosphorylation

Deletion of mitochondrial DNA in eukaryotes is currently attributed to rare accidental events associated with mitochondrial replication or repair of double-strand breaks. We report the discovery that yeast cells arrest harmful intramitochondrial superoxide production by shutting down respiration through genetically controlled deletion of mitochondrial oxidative phosphorylation genes. We show that this process critically involves the antioxidant enzyme superoxide dismutase 2 and two-way mitochondrial-nuclear communication through Rtg2 and Rtg3. While mitochondrial DNA homeostasis is rapidly restored after cessation of a short-term superoxide stress, long-term stress causes maladaptive persistence of the deletion process, leading to complete annihilation of the cellular pool of intact mitochondrial genomes and irrevocable loss of respiratory ability. This shows that oxidative stress-induced mitochondrial impairment may be under strict regulatory control. If the results extend to human cells, the results may prove to be of etiological as well as therapeutic importance with regard to age-related mitochondrial impairment and disease.

Changes to the mtDNA copy number during yeast culture growth

We show that the mitochondrial DNA (mtDNA) copy number in growing cultures of the yeast Saccharomyces cerevisiae increases by a factor of up to 4, being lowest (approx. 10 per haploid genome) and stable during rapid fermentative growth, and highest at the end of the respiratory phase. When yeast are grown on glucose, the onset of the mtDNA copy number increase coincides with the early stages of the diauxic shift, and the increase continues through respiration. A lesser yet still substantial copy number increase occurs when yeast are grown on a nonfermentable carbon source, i.e. when there is no diauxic shift. The mtDNA copy number increase during and for some time after the diauxic shift is not driven by an increase in cell size. The copy number increase occurs in both haploid and diploid strains but is markedly attenuated in a diploid wild isolate that is a ready sporulator. Strain-to-strain differences in mtDNA copy number are least apparent in fermentation and most apparent in late respiration or stationary phase. While changes in mitochondrial morphology and function were previously known to accompany changes in physiological state, it had not been previously shown that the mtDNA copy number changes substantially over time in a clonal growing culture. The mtDNA copy number in yeast is therefore a highly dynamic phenotype.

Directed yeast genome evolution by controlled introduction of trans-chromosomic structural variations

Naturally occurring structural variations (SVs) are a considerable source of genomic variation that can reshape the 3D architecture of chromosomes. Controllable methods aimed at introducing the complex SVs and their related molecular mechanisms have remained farfetched. In this study, an SV-prone yeast strain was developed using Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) technology with two synthetic chromosomes, namely synV and synX. The biosynthesis of astaxanthin is used as a readout and a proof of concept for the application of SVs in industries. Our findings showed that complex SVs, including a pericentric inversion and a trans-chromosome translocation between synV and synX, resulted in two neo-chromosomes and a 2.7-fold yield of astaxanthin. Also, genetic targets were mapped, which resulted in a higher astaxanthin yield, thus demonstrating the SVs' ability to reorganize genetic information along the chromosomes. The rational design of trans-chromosome translocation and pericentric inversion enabled precise induction of these phenomena. Collectively, this study provides an effective tool to not only accelerate the directed genome evolution but also to reveal the mechanistic insight of complex SVs for altering phenotypes.

rDNA array length is a major determinant of replicative lifespan in budding yeast

Genes encoding ribosomal RNA (rDNA) are organized into a repetitive array in eukaryotic genomes. The copy number of these genes often varies and is responsive to genetics and environment. Here, we show that variation in copy number at the rDNA locus is capable of altering replicative lifespan in yeast. These results indicate that considering rDNA copy number, and conditions that could potentially change this dynamic chromosome locus, is critical for evaluating lifespan. We propose that this rDNA locus represents the kind of repeated element in eukaryotic genomes that escapes easy detection, yet has phenotypic consequences, in this case lifespan.

The complex processes and interactions that regulate aging and determine lifespan are not fully defined for any organism. Here, taking advantage of recent technological advances in studying aging in budding yeast, we discovered a previously unappreciated relationship between the number of copies of the ribosomal RNA gene present in its chromosomal array and replicative lifespan (RLS). Specifically, the chromosomal ribosomal DNA (rDNA) copy number (rDNA CN) positively correlated with RLS and this interaction explained over 70% of variability in RLS among a series of wild-type strains. In strains with low rDNA CN, SIR2 expression was attenuated and extrachromosomal rDNA circle (ERC) accumulation was increased, leading to shorter lifespan. Suppressing ERC formation by deletion of FOB1 eliminated the relationship between rDNA CN and RLS. These data suggest that previously identified rDNA CN regulatory mechanisms limit lifespan. Importantly, the RLSs of reported lifespan-enhancing mutations were significantly impacted by rDNA CN, suggesting that changes in rDNA CN might explain the magnitude of some of those reported effects. We propose that because rDNA CN is modulated by environmental, genetic, and stochastic factors, considering rDNA CN is a prerequisite for accurate interpretation of lifespan data.

Modulating DNA Repair Pathways to Diversify Genomic Alterations in Saccharomyces cerevisiae

Nuclease based genome editing systems have emerged as powerful tools to drive genomic alterations and enhance genome evolution via precise engineering in the various human and microbial cells. However, error-prone DNA repair has not been well studied previously to generate diverse genomic alterations and novel phenotypes. Here, we systematically investigated the potential interplay between DNA double strand break (DSB) repair and genome editing tools, and found that modulating the DSB end resection proteins could significantly improve mutational efficiency and diversity without exogenous DNA template in yeast. Deleting SAE2, EXO1, or FUN30, or overexpressing MRE11-H125N (nuclease-dead allele of MRE11), for DSB end resection markedly increased the efficiency of CRISPR/SpCas9 (more than 22-fold) and CRISPR/AsCpf1 (more than 30-fold)-induced mutagenesis. Deleting SAE2 or overexpressing MRE11-H125N substantially diversified CRISPR/SpCas9 or AsCpf1-induced mutation 2–3-fold at URA3 locus, and 3–5-fold at ADE2 locus. Thus, the error-prone DNA repair protein was employed to develop a novel mutagenic genome editing (mGE) strategy, which can increase the mutation numbers and effectively improve the ethanol/glycerol ratio of Saccharomyces cerevisiae through modulating the expression of FPS1 and GPD1. This study highlighted the feasibility of potentially reshaping the capability of genome editing by regulating the different DSB repair proteins and can thus expand the application of genome editing in diversifying gene expression and enhancing genome evolution.

IMPORTANCE Most of the published papers about nuclease-assisted genome editing focused on precision engineering in human cells. However, the topic of inducing mutagenesis via error-prone repair has often been ignored in yeast. In this study, we reported that perturbing DNA repair, especially modifications of the various DSB end resection-related proteins, could greatly improve the mutational efficiency and diversity, and thus functionally reshape the capability of the different genome editing tools without requiring an exogenous DNA template in yeast. Specifically, mutagenic genome editing (mGE) was developed based on CRISPR/AsCpf1 and MRE11-H125N overexpression, and used to generate promoters of different strengths more efficiently. Thus, this work provides a novel method to diversify gene expression and enhance genome evolution.

Adaptive laboratory evolution in S. cerevisiae highlights role of transcription factors in fungal xenobiotic resistance

In vitro evolution and whole genome analysis were used to comprehensively identify the genetic determinants of chemical resistance in Saccharomyces cerevisiae. Sequence analysis identified many genes contributing to the resistance phenotype as well as numerous amino acids in potential targets that may play a role in compound binding. Our work shows that compound-target pairs can be conserved across multiple species. The set of 25 most frequently mutated genes was enriched for transcription factors, and for almost 25 percent of the compounds, resistance was mediated by one of 100 independently derived, gain-of-function SNVs found in a 170 amino acid domain in the two Zn₂C₆ transcription factors YRR1 and YRM1 (p < 1 × 10⁻¹⁰⁰). This remarkable enrichment for transcription factors as drug resistance genes highlights their important role in the evolution of antifungal xenobiotic resistance and underscores the challenge to develop antifungal treatments that maintain potency.

Ottilie et al. employ an experimental evolution approach to investigate the role of transcription factors in yeast chemical resistance. Most emergent mutations in resistant strains were enriched in transcription factor coding genes, highlighting their importance in drug resistance.

Effects of aneuploidy on cell behaviour and function

Aneuploidy, a genomic alternation characterized by deviations in the copy number of chromosomes, affects organisms from early development through to aging. Although it is a main cause of human pregnancy loss and a hallmark of cancer, how aneuploidy affects cellular function has been elusive. The last two decades have seen rapid advances in the understanding of the causes and consequences of aneuploidy at the molecular and cellular levels. These studies have uncovered effects of aneuploidy that can be beneficial or detrimental to cells and organisms in an environmental context-dependent and karyotype-dependent manner. Aneuploidy also imposes general stress on cells that stems from an imbalanced genome and, consequently, also an imbalanced proteome. These insights provide the fundamental framework for understanding the impact of aneuploidy in genome evolution, human pathogenesis and drug resistance.

VID22 counteracts G-quadruplex-induced genome instability

Genome instability is a condition characterized by the accumulation of genetic alterations and is a hallmark of cancer cells. To uncover new genes and cellular pathways affecting endogenous DNA damage and genome integrity, we exploited a Synthetic Genetic Array (SGA)-based screen in yeast. Among the positive genes, we identified VID22, reported to be involved in DNA double-strand break repair. vid22Δ cells exhibit increased levels of endogenous DNA damage, chronic DNA damage response activation and accumulate DNA aberrations in sequences displaying high probabilities of forming G-quadruplexes (G4-DNA). If not resolved, these DNA secondary structures can block the progression of both DNA and RNA polymerases and correlate with chromosome fragile sites. Vid22 binds to and protects DNA at G4-containing regions both in vitro and in vivo. Loss of VID22 causes an increase in gross chromosomal rearrangement (GCR) events dependent on G-quadruplex forming sequences. Moreover, the absence of Vid22 causes defects in the correct maintenance of G4-DNA rich elements, such as telomeres and mtDNA, and hypersensitivity to the G4-stabilizing ligand TMPyP4. We thus propose that Vid22 is directly involved in genome integrity maintenance as a novel regulator of G4 metabolism.

DNA Damage Responses during the Cell Cycle: Insights from Model Organisms and Beyond

Genome damage is a threat to all organisms. To respond to such damage, DNA damage responses (DDRs) lead to cell cycle arrest, DNA repair, and cell death. Many DDR components are highly conserved, whereas others have adapted to specific organismal needs. Immense progress in this field has been driven by model genetic organism research. This review has two main purposes. First, we provide a survey of model organism-based efforts to study DDRs. Second, we highlight how model organism study has contributed to understanding how specific DDRs are influenced by cell cycle stage. We also look forward, with a discussion of how future study can be expanded beyond typical model genetic organisms to further illuminate how the genome is protected.

The metabolic growth limitations of petite cells lacking the mitochondrial genome

Eukaryotic cells can survive the loss of their mitochondrial genome, but consequently suffer from severe growth defects. 'Petite yeasts', characterized by mitochondrial genome loss, are instrumental for studying mitochondrial function and physiology. However, the molecular cause of their reduced growth rate remains an open question. Here we show that petite cells suffer from an insufficient capacity to synthesize glutamate, glutamine, leucine and arginine, negatively impacting their growth. Using a combination of molecular genetics and omics approaches, we demonstrate the evolution of fast growth overcomes these amino acid deficiencies, by alleviating a perturbation in mitochondrial iron metabolism and by restoring a defect in the mitochondrial tricarboxylic acid cycle, caused by aconitase inhibition. Our results hence explain the slow growth of mitochondrial genome-deficient cells with a partial auxotrophy in four amino acids that results from distorted iron metabolism and an inhibited tricarboxylic acid cycle.

Ploidy and recombination proficiency shape the evolutionary adaptation to constitutive DNA replication stress

In haploid budding yeast, evolutionary adaptation to constitutive DNA replication stress alters three genome maintenance modules: DNA replication, the DNA damage checkpoint, and sister chromatid cohesion. We asked how these trajectories depend on genomic features by comparing the adaptation in three strains: haploids, diploids, and recombination deficient haploids. In all three, adaptation happens within 1000 generations at rates that are correlated with the initial fitness defect of the ancestors. Mutations in individual genes are selected at different frequencies in populations with different genomic features, but the benefits these mutations confer are similar in the three strains, and combinations of these mutations reproduce the fitness gains of evolved populations. Despite the differences in the selected mutations, adaptation targets the same three functional modules in strains with different genomic features, revealing a common evolutionary response to constitutive DNA replication stress.

Recombination-mediated genome rearrangements

Homologous recombination (HR) is a universal DNA double-strand break (DSB) repair pathway that uses an intact DNA molecule as a template. Signature HR reactions are homology search and DNA strand invasion catalyzed by the prototypical RecA-ssDNA filament (Rad51 and Dmc1 in eukaryotes), which produces heteroduplex DNA-containing joint molecules (JMs). These reactions uniquely infringe on the DNA strands association established at replication, on the basis of substantial sequence similarity. For that reason, and despite the high fidelity of its templated nature, DSB repair by HR authorizes the alteration of genome structure, guided by repetitive DNA elements. The resulting structural variations (SVs) can involve vast genomic regions, potentially affecting multiple coding sequences and regulatory elements at once, with possible pathological consequences. Here, we discuss recent advances in our understanding of genetic and molecular vulnerabilities of HR leading to SVs, and of the various fidelity-enforcing factors acting across scales on the balancing act of this complex pathway. An emphasis is put on extra-chomosomal DNAs, both product of, and substrate for HR-mediated chromosomal rearrangements.

Methylation of histone H3 at lysine 37 by Set1 and Set2 prevents spurious DNA replication

DNA replication initiates at genomic locations known as origins of replication, which, in S. cerevisiae, share a common DNA consensus motif. Despite being virtually nucleosome-free, origins of replication are greatly influenced by the surrounding chromatin state. Here, we show that histone H3 lysine 37 mono-methylation (H3K37me1) is catalyzed by Set1p and Set2p and that it regulates replication origin licensing. H3K37me1 is uniformly distributed throughout most of the genome, but it is scarce at replication origins, where it increases according to the timing of their firing. We find that H3K37me1 hinders Mcm2 interaction with chromatin, maintaining low levels of MCM outside of conventional replication origins. Lack of H3K37me1 results in defective DNA replication from canonical origins while promoting replication events at inefficient and non-canonical sites. Collectively, our results indicate that H3K37me1 ensures correct execution of the DNA replication program by protecting the genome from inappropriate origin licensing and spurious DNA replication.

A case of convergent-gene interference in the budding yeast knockout library causing chromosome instability

To maintain genome stability, organisms depend on faithful chromosome segregation, a process affected by diverse genetic pathways, some of which are not directly linked to mitosis. In this study, we set out to explore one such pathway represented by an undercharacterized gene, SNO1, identified previously in screens of the yeast knockout (YKO) library for mitotic fidelity genes. We found that the causative factor increasing mitotic error rate in the sno1Δ mutant is not loss of the Sno1 protein, but rather perturbation to the mRNA of the neighboring convergent gene, CTF13, encoding an essential component for forming the yeast kinetochore. This is caused by a combination of the Kanamycin resistance gene and the transcriptional terminator used in the YKO library affecting the CTF13 mRNA level and quality . We further provide a list of gene pairs potentially subjected to this artifact, which may be useful for accurate phenotypic interpretation of YKO mutants.

Protection of the C. elegans germ cell genome depends on diverse DNA repair pathways during normal proliferation

Maintaining genome integrity is particularly important in germ cells to ensure faithful transmission of genetic information across generations. Here we systematically describe germ cell mutagenesis in wild-type and 61 DNA repair mutants cultivated over multiple generations. ~44% of the DNA repair mutants analysed showed a >2-fold increased mutagenesis with a broad spectrum of mutational outcomes. Nucleotide excision repair deficiency led to higher base substitution rates, whereas polh-1(Polη) and rev-3(Polζ) translesion synthesis polymerase mutants resulted in 50-400 bp deletions. Signatures associated with defective homologous recombination fall into two classes: 1) brc-1/BRCA1 and rad-51/RAD51 paralog mutants showed increased mutations across all mutation classes, 2) mus-81/MUS81 and slx-1/SLX1 nuclease, and him-6/BLM, helq-1/HELQ or rtel-1/RTEL1 helicase mutants primarily accumulated structural variants. Repetitive and G-quadruplex sequence-containing loci were more frequently mutated in specific DNA repair backgrounds. Tandem duplications embedded in inverted repeats were observed in helq-1 helicase mutants, and a unique pattern of 'translocations' involving homeologous sequences occurred in rip-1 recombination mutants. atm-1/ATM checkpoint mutants harboured structural variants specifically enriched in subtelomeric regions. Interestingly, locally clustered mutagenesis was only observed for combined brc-1 and cep-1/p53 deficiency. Our study provides a global view of how different DNA repair pathways contribute to prevent germ cell mutagenesis.

Mutagenic mechanisms of cancer-associated DNA polymerase ϵ alleles

A single amino acid residue change in the exonuclease domain of human DNA polymerase ϵ, P286R, is associated with the development of colorectal cancers, and has been shown to impart a mutator phenotype. The corresponding Pol ϵ allele in the yeast Saccharomyces cerevisiae (pol2-P301R), was found to drive greater mutagenesis than an entirely exonuclease-deficient Pol ϵ (pol2-4), an unexpected phenotype of ultra-mutagenesis. By studying the impact on mutation frequency, type, replication-strand bias, and sequence context, we show that ultra-mutagenesis is commonly observed in yeast cells carrying a range of cancer-associated Pol ϵ exonuclease domain alleles. Similarities between mutations generated by these alleles and those generated in pol2-4 cells indicate a shared mechanism of mutagenesis that yields a mutation pattern similar to cancer Signature 14. Comparison of POL2 ultra-mutator with pol2-M644G, a mutant in the polymerase domain decreasing Pol ϵ fidelity, revealed unexpected analogies in the sequence context and strand bias of mutations. Analysis of mutational patterns unique to exonuclease domain mutant cells suggests that backtracking of the polymerase, when the mismatched primer end cannot be accommodated in the proofreading domain, results in the observed insertions and T>A mutations in specific sequence contexts.

Indole-3-acetic acid is a physiological inhibitor of TORC1 in yeast

Indole-3-acetic acid (IAA) is the most common, naturally occurring phytohormone that regulates cell division, differentiation, and senescence in plants. The capacity to synthesize IAA is also widespread among plant-associated bacterial and fungal species, which may use IAA as an effector molecule to define their relationships with plants or to coordinate their physiological behavior through cell-cell communication. Fungi, including many species that do not entertain a plant-associated life style, are also able to synthesize IAA, but the physiological role of IAA in these fungi has largely remained enigmatic. Interestingly, in this context, growth of the budding yeast Saccharomyces cerevisiae is sensitive to extracellular IAA. Here, we use a combination of various genetic approaches including chemical-genetic profiling, SAturated Transposon Analysis in Yeast (SATAY), and genetic epistasis analyses to identify the mode-of-action by which IAA inhibits growth in yeast. Surprisingly, these analyses pinpointed the target of rapamycin complex 1 (TORC1), a central regulator of eukaryotic cell growth, as the major growth-limiting target of IAA. Our biochemical analyses further demonstrate that IAA inhibits TORC1 both in vivo and in vitro. Intriguingly, we also show that yeast cells are able to synthesize IAA and specifically accumulate IAA upon entry into stationary phase. Our data therefore suggest that IAA contributes to proper entry of yeast cells into a quiescent state by acting as a metabolic inhibitor of TORC1.

Auxins are a major group of plant phytohormones that are critical for growth and development. Amongst the auxins, indole-3-acetic acid (IAA) is the most common, naturally occurring phytohormone that regulates cell division, differentiation, and senescence in plants. Interestingly, the capacity to synthesize and secrete IAA is also widespread among fungi, including the budding yeast Saccharomyces cerevisiae, but the role of IAA in fungi has largely remained unknown. Here, we confirm an earlier observation that IAA inhibits growth of budding yeast and show by diverse genetic and biochemical means that IAA restrains budding yeast growth by inhibiting the target of rapamycin complex 1 (TORC1), a highly conserved eukaryotic regulator of growth. Intriguingly, budding yeast cells accumulate IAA specifically when limited for nutrients, which suggests that IAA plays a hitherto unknown physiological role in contributing to the establishment of cellular quiescence by acting as a metabolic inhibitor of TORC1.

Auxin-mediated protein depletion for metabolic engineering in terpene-producing yeast

In metabolic engineering, loss-of-function experiments are used to understand and optimise metabolism. A conditional gene inactivation tool is required when gene deletion is lethal or detrimental to growth. Here, we exploit auxin-inducible protein degradation as a metabolic engineering approach in yeast. We demonstrate its effectiveness using terpenoid production. First, we target an essential prenyl-pyrophosphate metabolism protein, farnesyl pyrophosphate synthase (Erg20p). Degradation successfully redirects metabolic flux toward monoterpene (C10) production. Second, depleting hexokinase-2, a key protein in glucose signalling transduction, lifts glucose repression and boosts production of sesquiterpene (C15) nerolidol to 3.5 g L-1 in flask cultivation. Third, depleting acetyl-CoA carboxylase (Acc1p), another essential protein, delivers growth arrest without diminishing production capacity in nerolidol-producing yeast, providing a strategy to decouple growth and production. These studies demonstrate auxin-mediated protein degradation as an advanced tool for metabolic engineering. It also has potential for broader metabolic perturbation studies to better understand metabolism.

A High-Throughput Method for Screening for Genes Controlling Bacterial Conjugation of Antibiotic Resistance

The rapid transmission of antibiotic resistance genes on conjugative plasmids between bacterial host cells is a major cause of the accelerating antibiotic resistance crisis. There are currently no experimental platforms for fast and cost-efficient screening of genetic effects on antibiotic resistance transmission by conjugation, which prevents understanding and targeting conjugation.

The rapid horizontal transmission of antibiotic resistance genes on conjugative plasmids between bacterial host cells is a major cause of the accelerating antibiotic resistance crisis. There are currently no experimental platforms for fast and cost-efficient screening of genetic effects on antibiotic resistance transmission by conjugation, which prevents understanding and targeting conjugation. We introduce a novel experimental framework to screen for conjugation-based horizontal transmission of antibiotic resistance between >60,000 pairs of cell populations in parallel. Plasmid-carrying donor strains are constructed in high-throughput. We then mix the resistance plasmid-carrying donors with recipients in a design where only transconjugants can reproduce, measure growth in dense intervals, and extract transmission times as the growth lag. As proof-of-principle, we exhaustively explore chromosomal genes controlling F-plasmid donation within Escherichia coli populations, by screening the Keio deletion collection in high replication. We recover all seven known chromosomal gene mutants affecting conjugation as donors and identify many novel mutants, all of which diminish antibiotic resistance transmission. We validate nine of the novel genes’ effects in liquid mating assays and complement one of the novel genes’ effect on conjugation (rseA). The new framework holds great potential for exhaustive disclosing of candidate targets for helper drugs that delay resistance development in patients and societies and improve the longevity of current and future antibiotics. Further, the platform can easily be adapted to explore interspecies conjugation, plasmid-borne factors, and experimental evolution and be used for rapid construction of strains.

IMPORTANCE The rapid transmission of antibiotic resistance genes on conjugative plasmids between bacterial host cells is a major cause of the accelerating antibiotic resistance crisis. There are currently no experimental platforms for fast and cost-efficient screening of genetic effects on antibiotic resistance transmission by conjugation, which prevents understanding and targeting conjugation. We introduce a novel experimental framework to screen for conjugation-based horizontal transmission of antibiotic resistance between >60,000 pairs of cell populations in parallel. As proof-of-principle, we exhaustively explore chromosomal genes controlling F-plasmid donation within E. coli populations. We recover all previously known and many novel chromosomal gene mutants that affect conjugation efficiency. The new framework holds great potential for rapid screening of compounds that decrease transmission. Further, the platform can easily be adapted to explore interspecies conjugation, plasmid-borne factors, and experimental evolution and be used for rapid construction of strains.

Trajectory and uniqueness of mutational signatures in yeast mutators

The acquisition of mutations plays critical roles in adaptation, evolution, senescence, and tumorigenesis. Massive genome sequencing has allowed extraction of specific features of many mutational landscapes but it remains difficult to retrospectively determine the mechanistic origin(s), selective forces, and trajectories of transient or persistent mutations and genome rearrangements. Here, we conducted a prospective reciprocal approach to inactivate 13 single or multiple evolutionary conserved genes involved in distinct genome maintenance processes and characterize de novo mutations in 274 diploid Saccharomyces cerevisiae mutation accumulation lines. This approach revealed the diversity, complexity, and ultimate uniqueness of mutational landscapes, differently composed of base substitutions, small insertions/deletions (InDels), structural variants, and/or ploidy variations. Several landscapes parallel the repertoire of mutational signatures in human cancers while others are either novel or composites of subsignatures resulting from distinct DNA damage lesions. Notably, the increase of base substitutions in the homologous recombination-deficient Rad51 mutant, specifically dependent on the Polζ translesion polymerase, yields COSMIC signature 3 observed in BRCA1/BRCA2-mutant breast cancer tumors. Furthermore, "mutome" analyses in highly polymorphic diploids and single-cell bottleneck lineages revealed a diverse spectrum of loss-of-heterozygosity (LOH) signatures characterized by interstitial and terminal chromosomal events resulting from interhomolog mitotic cross-overs. Following the appearance of heterozygous mutations, the strong stimulation of LOHs in the rad27/FEN1 and tsa1/PRDX1 backgrounds leads to fixation of homozygous mutations or their loss along the lineage. Overall, these mutomes and their trajectories provide a mechanistic framework to understand the origin and dynamics of genome variations that accumulate during clonal evolution.

Comparing the utility of in vivo transposon mutagenesis approaches in yeast species to infer gene essentiality

In vivo transposon mutagenesis, coupled with deep sequencing, enables large-scale genome-wide mutant screens for genes essential in different growth conditions. We analyzed six large-scale studies performed on haploid strains of three yeast species (Saccharomyces cerevisiae, Schizosaccaromyces pombe, and Candida albicans), each mutagenized with two of three different heterologous transposons (AcDs, Hermes, and PiggyBac). Using a machine-learning approach, we evaluated the ability of the data to predict gene essentiality. Important data features included sufficient numbers and distribution of independent insertion events. All transposons showed some bias in insertion site preference because of jackpot events, and preferences for specific insertion sequences and short-distance vs long-distance insertions. For PiggyBac, a stringent target sequence limited the ability to predict essentiality in genes with few or no target sequences. The machine learning approach also robustly predicted gene function in less well-studied species by leveraging cross-species orthologs. Finally, comparisons of isogenic diploid versus haploid S. cerevisiae isolates identified several genes that are haplo-insufficient, while most essential genes, as expected, were recessive. We provide recommendations for the choice of transposons and the inference of gene essentiality in genome-wide studies of eukaryotic haploid microbes such as yeasts, including species that have been less amenable to classical genetic studies.

Systematic analysis of nuclear gene function in respiratory growth and expression of the mitochondrial genome in S. cerevisiae

The production of metabolic energy in form of ATP by oxidative phosphorylation depends on the coordinated action of hundreds of nuclear-encoded mitochondrial proteins and a handful of proteins encoded by the mitochondrial genome (mtDNA). We used the yeast Saccharomyces cerevisiae as a model system to systematically identify the genes contributing to this process. Integration of genome-wide high-throughput growth assays with previously published large data sets allowed us to define with high confidence a set of 254 nuclear genes that are indispensable for respiratory growth. Next, we induced loss of mtDNA in the yeast deletion collection by growth on ethidium bromide-containing medium and identified twelve genes that are essential for viability in the absence of mtDNA (i.e. petite-negative). Replenishment of mtDNA by cytoduction showed that respiratory-deficient phenotypes are highly variable in many yeast mutants. Using a mitochondrial genome carrying a selectable marker, ARG8^m, we screened for mutants that are specifically defective in maintenance of mtDNA and mitochondrial protein synthesis. We found that up to 176 nuclear genes are required for expression of mitochondria-encoded proteins during fermentative growth. Taken together, our data provide a comprehensive picture of the molecular processes that are required for respiratory metabolism in a simple eukaryotic cell.

Candida albicans Genetic Background Influences Mean and Heterogeneity of Drug Responses and Genome Stability during Evolution in Fluconazole

Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

The importance of within-species diversity in determining the evolutionary potential of a population to evolve drug resistance or tolerance is not well understood, including in eukaryotic pathogens. To examine the influence of genetic background, we evolved replicates of 20 different clinical isolates of Candida albicans, a human fungal pathogen, in fluconazole, the commonly used antifungal drug. The isolates hailed from the major C. albicans clades and had different initial levels of drug resistance and tolerance to the drug. The majority of replicates rapidly increased in fitness in the evolutionary environment, with the degree of improvement inversely correlated with parental strain fitness in the drug. Improvement was largely restricted to up to the evolutionary level of drug: only 4% of the evolved replicates increased resistance (MIC) above the evolutionary level of drug. Prevalent changes were altered levels of drug tolerance (slow growth of a subpopulation of cells at drug concentrations above the MIC) and increased diversity of genome size. The prevalence and predominant direction of these changes differed in a strain-specific manner, but neither correlated directly with parental fitness or improvement in fitness. Rather, low parental strain fitness was correlated with high levels of heterogeneity in fitness, tolerance, and genome size among evolved replicates. Thus, parental strain background is an important determinant in mean improvement to the evolutionary environment as well as the diversity of evolved phenotypes, and the range of possible responses of a pathogen to an antimicrobial drug cannot be captured by in-depth study of a single strain background.

IMPORTANCE Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

Role of the 40S beak ribosomal protein eS12 in ribosome biogenesis and function in Saccharomyces cerevisiae

In eukaryotes, the beak structure of 40S subunits is formed by the protrusion of the 18S rRNA helix 33 and three ribosomal proteins: eS10, eS12 and eS31. The exact role of these proteins in ribosome biogenesis is not well understood. While eS10 is an essential protein encoded by two paralogous genes in Saccharomyces cerevisiae, eS12 and eS31 are not essential proteins encoded by the single-copy genes RPS12 and UBI3, respectively. Here, we have analysed the contribution of yeast eS12 to ribosome biogenesis and compared it with that of eS31. Polysome analysis reveals that deletion of either RPS12 or UBI3 results in equivalent 40S deficits. Analysis of pre-rRNA processing indicates that eS12, akin to eS31, is required for efficient processing of 20S pre-rRNA to mature 18S rRNA. Moreover, we show that the 20S pre-rRNA accumulates within cytoplasmic pre-40S particles, as deduced from FISH experiments and the lack of nuclear retention of 40S subunit reporter proteins, in rps12∆ and ubi3∆ cells. However, these particles containing 20S pre-rRNA are not efficiently incorporated into polyribosomes. We also provide evidence for a genetic interaction between eS12 or eS31 and the late-acting 40S assembly factors Enp1 and Ltv1, which appears not to be linked to the dynamics of their association with or release from pre-40S particles in the absence of either eS12 or eS31. Finally, we show that eS12- and eS31-deficient ribosomes exhibit increased levels of translational misreading. Altogether, our data highlight distinct important roles of the beak region during ribosome assembly and function.

Discordant evolution of mitochondrial and nuclear yeast genomes at population level

Background

Mitochondria are essential organelles partially regulated by their own genomes. The mitochondrial genome maintenance and inheritance differ from the nuclear genome, potentially uncoupling their evolutionary trajectories. Here, we analysed mitochondrial sequences obtained from the 1011 Saccharomyces cerevisiae strain collection and identified pronounced differences with their nuclear genome counterparts.

Results

In contrast with pre-whole genome duplication fungal species, S. cerevisiae mitochondrial genomes show higher genetic diversity compared to the nuclear genomes. Strikingly, mitochondrial genomes appear to be highly admixed, resulting in a complex interconnected phylogeny with a weak grouping of isolates, whereas interspecies introgressions are very rare. Complete genome assemblies revealed that structural rearrangements are nearly absent with rare inversions detected. We tracked intron variation in COX1 and COB to infer gain and loss events throughout the species evolutionary history. Mitochondrial genome copy number is connected with the nuclear genome and linearly scale up with ploidy. We observed rare cases of naturally occurring mitochondrial DNA loss, petite, with a subset of them that do not suffer the expected growth defect in fermentable rich media.

Conclusions

Overall, our results illustrate how differences in the biology of two genomes coexisting in the same cells can lead to discordant evolutionary histories.