Premature aging in aneuploid yeast is caused in part by aneuploidy-induced defects in Ribosome Quality Control.

Premature aging is a hallmark of Down syndrome, caused by trisomy of human chromosome 21, but the reason is unclear and difficult to study in humans. We used an aneuploid model in wild yeast to show that chromosome amplification disrupts nutrient-induced cell-cycle arrest, quiescence entry, and healthy aging, across genetic backgrounds and amplified chromosomes. We discovered that these defects are due in part to aneuploidy-induced dysfunction in Ribosome Quality Control (RQC). Compared to euploids, aneuploids entering quiescence display aberrant ribosome profiles, accumulate RQC intermediates, and harbor an increased load of protein aggregates. Although they have normal proteasome capacity, aneuploids show signs of ubiquitin dysregulation, which impacts cyclin abundance to disrupt arrest. Remarkably, inducing ribosome stalling in euploids produces similar aberrations, while up-regulating limiting RQC subunits or proteins in ubiquitin metabolism alleviates many of the aneuploid defects. Our results provide implications for other aneuploidy disorders including Down syndrome.

Physiological Roles of an Acinetobacter-specific σ Factor.

The Gram-negative pathogen Acinetobacter baumannii is considered an "urgent threat" to human health due to its propensity to become antibiotic resistant. Understanding the distinct regulatory paradigms used by A. baumannii to mitigate cellular stresses may uncover new therapeutic targets. Many γ-proteobacteria use the extracytoplasmic function (ECF) σ factor, RpoE, to invoke envelope homeostasis networks in response to stress. Acinetobacter species contain the poorly characterized ECF "SigAb;" however, it is unclear if SigAb has the same physiological role as RpoE. Here, we show that SigAb is a metal stress-responsive ECF that appears unique to Acinetobacter species and distinct from RpoE. We combine promoter mutagenesis, motif scanning, and ChIP-seq to define the direct SigAb regulon, which consists of sigAb itself, the stringent response mediator, relA, and the uncharacterized small RNA, "sabS." However, RNA-seq of strains overexpressing SigAb revealed a large, indirect regulon containing hundreds of genes. Metal resistance genes are key elements of the indirect regulon, as CRISPRi knockdown of sigAb or sabS resulted in increased copper sensitivity and excess copper induced SigAb-dependent transcription. Further, we found that two uncharacterized genes in the sigAb operon, "aabA" and "aabB", have anti-SigAb activity. Finally, employing a targeted Tn-seq approach that uses CRISPR-associated transposons, we show that sigAb, aabA, and aabB are important for fitness even during optimal growth conditions. Our work reveals new physiological roles for SigAb and SabS, provides a novel approach for assessing gene fitness, and highlights the distinct regulatory architecture of A. baumannii.

Auto-MuRCiS: a streamlined software package for analysis of multiplex, randomized CRISPR interference sequencing data.

Multiplex, randomized CRISPR interference sequencing (MuRCiS) allows for the simultaneous identification of multiple gene knockouts that together influence microbial processes. Here, we report on an updated analysis tool called Auto-MuRCiS that utilizes Docker to make the analysis of these data rapid and more user-friendly.

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model.

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will only be maintained if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74-94% of the variance in aneuploid strains' growth rates is explained by the additive cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of snoRNAs and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

A Targeted Genome-scale Overexpression Platform for Proteobacteria.

Targeted, genome-scale gene perturbation screens using Clustered Regularly Interspaced Short Palindromic Repeats interference (CRISPRi) and activation (CRISPRa) have revolutionized eukaryotic genetics, advancing medical, industrial, and basic research. Although CRISPRi knockdowns have been broadly applied in bacteria, options for genome-scale overexpression face key limitations. Here, we develop a facile approach for genome-scale gene overexpression in bacteria we call, "CRISPRtOE" (CRISPR transposition and OverExpression). We create a platform for comprehensive gene targeting using CRISPR-associated transposition (CAST) and show that transposition occurs at a higher frequency in non-transcribed DNA. We then demonstrate that CRISPRtOE can upregulate gene expression in Proteobacteria with medical and industrial relevance by integrating synthetic promoters of varying strength upstream of target genes. Finally, we employ CRISPRtOE screening at the genome-scale in Escherichia coli, recovering known antibiotic targets and genes with unexplored roles in antibiotic function. We envision that CRISPRtOE will be a valuable overexpression tool for antibiotic mode of action, industrial strain optimization, and gene function discovery in bacteria.

Identification of recG genetic interactions in Escherichia coli by transposon sequencing.

IMPORTANCE: DNA damage and subsequent DNA repair processes are mutagenic in nature and an important driver of evolution in prokaryotes, including antibiotic resistance development. Genetic screening approaches, such as transposon sequencing (Tn-seq), have provided important new insights into gene function and genetic relationships. Here, we employed Tn-seq to gain insight into the function of the recG gene, which renders Escherichia coli cells moderately sensitive to a variety of DNA-damaging agents when they are absent. The reported recG genetic interactions can be used in combination with future screens to aid in a more complete reconstruction of DNA repair pathways in bacteria.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast.

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene overexpression, revealing that the fitness costs of copy-number variation can vary substantially with genetic background in a common-garden environment. But the interplay between copy-number variation tolerance and environment remains unexplored on a genomic scale. Here, we measured the tolerance to gene overexpression in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride stress. Overexpressed genes that are commonly deleterious during sodium chloride stress recapitulated those commonly deleterious under standard conditions. However, sodium chloride stress uncovered novel differences in strain responses to gene overexpression. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to sodium chloride stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to overexpression of specific genes. Although most genes were deleterious, hundreds were beneficial when overexpressed-remarkably, most of these effects were strain specific. Few beneficial genes were shared between the sodium chloride-sensitive isolates, implicating mechanistic differences behind their sodium chloride sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural copy-number variation of a sodium export pump that likely contributes to strain-specific responses to overexpression of other genes. Our results reveal extensive strain-by-environment interactions in the response to gene copy-number variation, raising important implications for the accessibility of copy-number variation-dependent evolutionary routes under times of stress.

Gene-by-environment interactions influence the fitness cost of gene copy-number variation in yeast.

Variation in gene copy number can alter gene expression and influence downstream phenotypes; thus copy-number variation (CNV) provides a route for rapid evolution if the benefits outweigh the cost. We recently showed that genetic background significantly influences how yeast cells respond to gene over-expression (OE), revealing that the fitness costs of CNV can vary substantially with genetic background in a common-garden environment. But the interplay between CNV tolerance and environment remains unexplored on a genomic scale. Here we measured the tolerance to gene OE in four genetically distinct Saccharomyces cerevisiae strains grown under sodium chloride (NaCl) stress. OE genes that are commonly deleterious during NaCl stress recapitulated those commonly deleterious under standard conditions. However, NaCl stress uncovered novel differences in strain responses to gene OE. West African strain NCYC3290 and North American oak isolate YPS128 are more sensitive to NaCl stress than vineyard BC187 and laboratory strain BY4743. Consistently, NCYC3290 and YPS128 showed the greatest sensitivities to gene OE. Although most genes were deleterious, hundreds were beneficial when overexpressed - remarkably, most of these effects were strain specific. Few beneficial genes were shared between the NaCl-sensitive isolates, implicating mechanistic differences behind their NaCl sensitivity. Transcriptomic analysis suggested underlying vulnerabilities and tolerances across strains, and pointed to natural CNV of a sodium export pump that likely contributes to strain-specific responses to OE of other genes. Our results reveal extensive strain-by-environment interaction in the response to gene CNV, raising important implications for the accessibility of CNV-dependent evolutionary routes under times of stress.

Identification of genetic interactions with priB links the PriA/PriB DNA replication restart pathway to double-strand DNA break repair in Escherichia coli.

Collisions between DNA replication complexes (replisomes) and impediments such as damaged DNA or proteins tightly bound to the chromosome lead to premature dissociation of replisomes at least once per cell cycle in Escherichia coli. Left unrepaired, these events produce incompletely replicated chromosomes that cannot be properly partitioned into daughter cells. DNA replication restart, the process that reloads replisomes at prematurely terminated sites, is therefore essential in E. coli and other bacteria. Three replication restart pathways have been identified in E. coli: PriA/PriB, PriA/PriC, and PriC/Rep. A limited number of genetic interactions between replication restart and other genome maintenance pathways have been defined, but a systematic study placing replication restart reactions in a broader cellular context has not been performed. We have utilized transposon-insertion sequencing to identify new genetic interactions between DNA replication restart pathways and other cellular systems. Known genetic interactors with the priB replication restart gene (uniquely involved in the PriA/PriB pathway) were confirmed and several novel priB interactions were discovered. Targeted genetic and imaging-based experiments with priB and its genetic partners revealed significant double-strand DNA break accumulation in strains with mutations in dam, rep, rdgC, lexA, or polA. Modulating the activity of the RecA recombinase partially suppressed the detrimental effects of rdgC or lexA mutations in ΔpriB cells. Taken together, our results highlight roles for several genes in double-strand DNA break homeostasis and define a genetic network that facilitates DNA repair/processing upstream of PriA/PriB-mediated DNA replication restart in E. coli.

Transcriptomic Data Sets for Zymomonas mobilis 2032 during Fermentation of Ammonia Fiber Expansion (AFEX)-Pretreated Corn Stover and Switchgrass Hydrolysates.

The transcriptomes of Zymomonas mobilis 2032 were captured during the fermentation of ammonia fiber expansion (AFEX)-pretreated corn stover and switchgrass hydrolysates containing different concentrations of glucose and xylose. RNA samples were collected when Z. mobilis was fermenting glucose or xylose. Here, we present transcriptome sequencing (RNA-Seq) data obtained during separate phases of glucose or xylose consumption.

Comparative chemical genomic profiling across plant-based hydrolysate toxins reveals widespread antagonism in fitness contributions.

The budding yeast Saccharomyces cerevisiae has been used extensively in fermentative industrial processes, including biofuel production from sustainable plant-based hydrolysates. Myriad toxins and stressors found in hydrolysates inhibit microbial metabolism and product formation. Overcoming these stresses requires mitigation strategies that include strain engineering. To identify shared and divergent mechanisms of toxicity and to implicate gene targets for genetic engineering, we used a chemical genomic approach to study fitness effects across a library of S. cerevisiae deletion mutants cultured anaerobically in dozens of individual compounds found in different types of hydrolysates. Relationships in chemical genomic profiles identified classes of toxins that provoked similar cellular responses, spanning inhibitor relationships that were not expected from chemical classification. Our results also revealed widespread antagonistic effects across inhibitors, such that the same gene deletions were beneficial for surviving some toxins but detrimental for others. This work presents a rich dataset relating gene function to chemical compounds, which both expands our understanding of plant-based hydrolysates and provides a useful resource to identify engineering targets.

Crabtree/Warburg-like aerobic xylose fermentation by engineered Saccharomyces cerevisiae.

Bottlenecks in the efficient conversion of xylose into cost-effective biofuels have limited the widespread use of plant lignocellulose as a renewable feedstock. The yeast Saccharomyces cerevisiae ferments glucose into ethanol with such high metabolic flux that it ferments high concentrations of glucose aerobically, a trait called the Crabtree/Warburg Effect. In contrast to glucose, most engineered S. cerevisiae strains do not ferment xylose at economically viable rates and yields, and they require respiration to achieve sufficient xylose metabolic flux and energy return for growth aerobically. Here, we evolved respiration-deficient S. cerevisiae strains that can grow on and ferment xylose to ethanol aerobically, a trait analogous to the Crabtree/Warburg Effect for glucose. Through genome sequence comparisons and directed engineering, we determined that duplications of genes encoding engineered xylose metabolism enzymes, as well as TKL1, a gene encoding a transketolase in the pentose phosphate pathway, were the causative genetic changes for the evolved phenotype. Reengineered duplications of these enzymes, in combination with deletion mutations in HOG1, ISU1, GRE3, and IRA2, increased the rates of aerobic and anaerobic xylose fermentation. Importantly, we found that these genetic modifications function in another genetic background and increase the rate and yield of xylose-to-ethanol conversion in industrially relevant switchgrass hydrolysate, indicating that these specific genetic modifications may enable the sustainable production of industrial biofuels from yeast. We propose a model for how key regulatory mutations prime yeast for aerobic xylose fermentation by lowering the threshold for overflow metabolism, allowing mutations to increase xylose flux and to redirect it into fermentation products.

Natural variation in the consequences of gene overexpression and its implications for evolutionary trajectories.

Copy number variation through gene or chromosome amplification provides a route for rapid phenotypic variation and supports the long-term evolution of gene functions. Although the evolutionary importance of copy-number variation is known, little is understood about how genetic background influences its tolerance. Here, we measured fitness costs of over 4000 overexpressed genes in 15 Saccharomyces cerevisiae strains representing different lineages, to explore natural variation in tolerating gene overexpression (OE). Strain-specific effects dominated the fitness costs of gene OE. We report global differences in the consequences of gene OE, independent of the amplified gene, as well as gene-specific effects that were dependent on the genetic background. Natural variation in the response to gene OE could be explained by several models, including strain-specific physiological differences, resource limitations, and regulatory sensitivities. This work provides new insight on how genetic background influences tolerance to gene amplification and the evolutionary trajectories accessible to different backgrounds.

Delila-PY, a Pipeline for Utilizing the Delila Suite of Software to Identify Potential DNA Binding Motifs.

Predicting potential DNA binding motifs is a critical part of understanding gene expression across all domains of life. Here, we report the development of Delila-PY, an easy-to-use pipeline for utilizing the Delila suite of software to identify DNA binding motifs.

CRISpy-Pop: A Web Tool for Designing CRISPR/Cas9-Driven Genetic Modifications in Diverse Populations.

CRISPR/Cas9 is a powerful tool for editing genomes, but design decisions are generally made with respect to a single reference genome. With population genomic data becoming available for an increasing number of model organisms, researchers are interested in manipulating multiple strains and lines. CRISpy-pop is a web application that generates and filters guide RNA sequences for CRISPR/Cas9 genome editing for diverse yeast and bacterial strains. The current implementation designs and predicts the activity of guide RNAs against more than 1000 Saccharomyces cerevisiae genomes, including 167 strains frequently used in bioenergy research. Zymomonas mobilis, an increasingly popular bacterial bioenergy research model, is also supported. CRISpy-pop is available as a web application (https://CRISpy-pop.glbrc.org/) with an intuitive graphical user interface. CRISpy-pop also cross-references the human genome to allow users to avoid the selection of guide RNAs with potential biosafety concerns. Additionally, CRISpy-pop predicts the strain coverage of each guide RNA within the supported strain sets, which aids in functional population genetic studies. Finally, we validate how CRISpy-pop can accurately predict the activity of guide RNAs across strains using population genomic data.

Phosphoproteome Response to Dithiothreitol Reveals Unique Versus Shared Features of Saccharomyces cerevisiae Stress Responses.

To cope with sudden changes in the external environment, the budding yeast Saccharomyces cerevisiae orchestrates a multifaceted response that spans many levels of physiology. Several studies have interrogated the transcriptome response to endoplasmic reticulum (ER) stress and the role of regulators such as the Ire1 kinase and Hac1 transcription factors. However, less is known about responses to ER stress at other levels of physiology. Here, we used quantitative phosphoproteomics and computational network inference to uncover the yeast phosphoproteome response to the reducing agent dithiothreitol (DTT) and the upstream signaling network that controls it. We profiled wild-type cells and mutants lacking IRE1 or MAPK kinases MKK1 and MKK2, before and at various times after DTT treatment. In addition to revealing downstream targets of these kinases, our inference approach predicted new regulators in the DTT response, including cell-cycle regulator Cdc28 and osmotic-response kinase Rck2, which we validated computationally. Our results also revealed similarities and surprising differences in responses to different stress conditions, especially in the response of protein kinase A targets. These results have implications for the breadth of signaling programs that can give rise to common stress response signatures.

COnTORT: COmprehensive Transcriptomic ORganizational Tool for Simultaneously Retrieving and Organizing Numerous Gene Expression Data Sets from the NCBI Gene Expression Omnibus Database.

We introduce COnTORT (COmprehensive Transcriptomic ORganizational Tool), a publicly available program that retrieves all available gene expression data and associated metadata for an organism from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. The data are compiled into text files that can be used for downstream bioinformatic applications.

Network inference reveals novel connections in pathways regulating growth and defense in the yeast salt response.

Cells respond to stressful conditions by coordinating a complex, multi-faceted response that spans many levels of physiology. Much of the response is coordinated by changes in protein phosphorylation. Although the regulators of transcriptome changes during stress are well characterized in Saccharomyces cerevisiae, the upstream regulatory network controlling protein phosphorylation is less well dissected. Here, we developed a computational approach to infer the signaling network that regulates phosphorylation changes in response to salt stress. We developed an approach to link predicted regulators to groups of likely co-regulated phospho-peptides responding to stress, thereby creating new edges in a background protein interaction network. We then use integer linear programming (ILP) to integrate wild type and mutant phospho-proteomic data and predict the network controlling stress-activated phospho-proteomic changes. The network we inferred predicted new regulatory connections between stress-activated and growth-regulating pathways and suggested mechanisms coordinating metabolism, cell-cycle progression, and growth during stress. We confirmed several network predictions with co-immunoprecipitations coupled with mass-spectrometry protein identification and mutant phospho-proteomic analysis. Results show that the cAMP-phosphodiesterase Pde2 physically interacts with many stress-regulated transcription factors targeted by PKA, and that reduced phosphorylation of those factors during stress requires the Rck2 kinase that we show physically interacts with Pde2. Together, our work shows how a high-quality computational network model can facilitate discovery of new pathway interactions during osmotic stress.

Genome-wide association across Saccharomyces cerevisiae strains reveals substantial variation in underlying gene requirements for toxin tolerance.

Cellulosic plant biomass is a promising sustainable resource for generating alternative biofuels and biochemicals with microbial factories. But a remaining bottleneck is engineering microbes that are tolerant of toxins generated during biomass processing, because mechanisms of toxin defense are only beginning to emerge. Here, we exploited natural diversity in 165 Saccharomyces cerevisiae strains isolated from diverse geographical and ecological niches, to identify mechanisms of hydrolysate-toxin tolerance. We performed genome-wide association (GWA) analysis to identify genetic variants underlying toxin tolerance, and gene knockouts and allele-swap experiments to validate the involvement of implicated genes. In the process of this work, we uncovered a surprising difference in genetic architecture depending on strain background: in all but one case, knockout of implicated genes had a significant effect on toxin tolerance in one strain, but no significant effect in another strain. In fact, whether or not the gene was involved in tolerance in each strain background had a bigger contribution to strain-specific variation than allelic differences. Our results suggest a major difference in the underlying network of causal genes in different strains, suggesting that mechanisms of hydrolysate tolerance are very dependent on the genetic background. These results could have significant implications for interpreting GWA results and raise important considerations for engineering strategies for industrial strain improvement.

Electrostatic confinement and manipulation of DNA molecules for genome analysis.

Very large DNA molecules enable comprehensive analysis of complex genomes, such as human, cancer, and plants because they span across sequence repeats and complex somatic events. When physically manipulated, or analyzed as single molecules, long polyelectrolytes are problematic because of mechanical considerations that include shear-mediated breakage, dealing with the massive size of these coils, or the length of stretched DNAs using common experimental techniques and fluidic devices. Accordingly, we harness analyte "issues" as exploitable advantages by our invention and characterization of the "molecular gate," which controls and synchronizes formation of stretched DNA molecules as DNA dumbbells within nanoslit geometries. Molecular gate geometries comprise micro- and nanoscale features designed to synergize very low ionic strength conditions in ways we show effectively create an "electrostatic bottle." This effect greatly enhances molecular confinement within large slit geometries and supports facile, synchronized electrokinetic loading of nanoslits, even without dumbbell formation. Device geometries were considered at the molecular and continuum scales through computer simulations, which also guided our efforts to optimize design and functionalities. In addition, we show that the molecular gate may govern DNA separations because DNA molecules can be electrokinetically triggered, by varying applied voltage, to enter slits in a size-dependent manner. Lastly, mapping the Mesoplasmaflorum genome, via synchronized dumbbell formation, validates our nascent approach as a viable starting point for advanced development that will build an integrated system capable of large-scale genome analysis.