Quantitative phenotyping via deep barcode sequencing

Improving the Accuracy of Bulk Fitness Assays by Correcting Barcode Processing Biases

Measuring the fitnesses of genetic variants is a fundamental objective in evolutionary biology. A standard approach for measuring microbial fitnesses in bulk involves labeling a library of genetic variants with unique sequence barcodes, competing the labeled strains in batch culture, and using deep sequencing to track changes in the barcode abundances over time. However, idiosyncratic properties of barcodes can induce nonuniform amplification or uneven sequencing coverage that causes some barcodes to be over- or under-represented in samples. This systematic bias can result in erroneous read count trajectories and misestimates of fitness. Here, we develop a computational method, named REBAR (Removing the Effects of Bias through Analysis of Residuals), for inferring the effects of barcode processing bias by leveraging the structure of systematic deviations in the data. We illustrate this approach by applying it to two independent data sets, and demonstrate that this method estimates and corrects for bias more accurately than standard proxies, such as GC-based corrections. REBAR mitigates bias and improves fitness estimates in high-throughput assays without introducing additional complexity to the experimental protocols, with potential applications in a range of experimental evolution and mutation screening contexts.

Functional analysis of the Aspergillus fumigatus kinome identifies a druggable DYRK kinase that regulates septal plugging

More than 10 million people suffer from lung diseases caused by the pathogenic fungus Aspergillus fumigatus. Azole antifungals represent first-line therapeutics for most of these infections but resistance is rising, therefore the identification of antifungal targets whose inhibition synergises with the azoles could improve therapeutic outcomes. Here, we generate a library of 111 genetically barcoded null mutants of Aspergillus fumigatus in genes encoding protein kinases, and show that loss of function of kinase YakA results in hypersensitivity to the azoles and reduced pathogenicity. YakA is an orthologue of Candida albicans Yak1, a TOR signalling pathway kinase involved in modulation of stress responsive transcriptional regulators. We show that YakA has been repurposed in A. fumigatus to regulate blocking of the septal pore upon exposure to stress. Loss of YakA function reduces the ability of A. fumigatus to penetrate solid media and to grow in mouse lung tissue. We also show that 1-ethoxycarbonyl-beta-carboline (1-ECBC), a compound previously shown to inhibit C. albicans Yak1, prevents stress-mediated septal spore blocking and synergises with the azoles to inhibit A. fumigatus growth.

Thiolutin has complex effects in vivo but is a direct inhibitor of RNA polymerase II in vitro

Thiolutin is a natural product transcription inhibitor with an unresolved mode of action. Thiolutin and the related dithiolopyrrolone holomycin chelate Zn2+ and previous studies have concluded that RNA Polymerase II (Pol II) inhibition in vivo is indirect. Here, we present chemicogenetic and biochemical approaches to investigate thiolutin's mode of action in Saccharomyces cerevisiae. We identify mutants that alter sensitivity to thiolutin. We provide genetic evidence that thiolutin causes oxidation of thioredoxins in vivo and that thiolutin both induces oxidative stress and interacts functionally with multiple metals including Mn2+ and Cu2+, and not just Zn2+. Finally, we show direct inhibition of RNA polymerase II (Pol II) transcription initiation by thiolutin in vitro in support of classical studies that thiolutin can directly inhibit transcription in vitro. Inhibition requires both Mn2+ and appropriate reduction of thiolutin as excess DTT abrogates its effects. Pause prone, defective elongation can be observed in vitro if inhibition is bypassed. Thiolutin effects on Pol II occupancy in vivo are widespread but major effects are consistent with prior observations for Tor pathway inhibition and stress induction, suggesting that thiolutin use in vivo should be restricted to studies on its modes of action and not as an experimental tool.

Utilizing functional genomics in Saccharomyces cerevisiae to characterize food preservative compounds: A pilot study

Chemical preservatives are ubiquitously used to suppress the growth of or kill microorganisms across numerous industries, including the food industry. Utilizing yeast functional genomic techniques, genes and their functions can be observed at a genomic scale to elucidate how environmental stressors (e.g., chemical preservatives) impact microbial survival. These types of chemical genomics approaches can reveal genetic mutations that result in preservative resistance or sensitivity, assist in identification of preservative mechanism of action, and can be used to compare different preservatives for rational design of preservative mixtures. In this proof-of-concept study, we performed deletion and high-copy genetic expression screens to identify mutants that confer drug resistance to sodium benzoate, potassium sorbate, rosemary extract, and Natamax. By observing overlapping mutant genes between genetic screens, we were able to identify functional overlap between chemical preservatives and begin to explain mechanisms of action for these compounds.

Randomly barcoded transposon mutant libraries for gut commensals I: Strategies for efficient library construction

Randomly barcoded transposon mutant libraries are powerful tools for studying gene function and organization, assessing gene essentiality and pathways, discovering potential therapeutic targets, and understanding the physiology of gut bacteria and their interactions with the host. However, construction of high-quality libraries with uniform representation can be challenging. In this review, we survey various strategies for barcoded library construction, including transposition systems, methods of transposon delivery, optimal library size, and transconjugant selection schemes. We discuss the advantages and limitations of each approach, as well as factors to consider when selecting a strategy. In addition, we highlight experimental and computational advances in arraying condensed libraries from mutant pools. We focus on examples of successful library construction in gut bacteria and their application to gene function studies and drug discovery. Given the need for understanding gene function and organization in gut bacteria, we provide a comprehensive guide for researchers to construct randomly barcoded transposon mutant libraries.

From beer to breadboards: yeast as a force for biological innovation

The history of yeast Saccharomyces cerevisiae, aka brewer's or baker's yeast, is intertwined with our own. Initially domesticated 8,000 years ago to provide sustenance to our ancestors, for the past 150 years, yeast has served as a model research subject and a platform for technology. In this review, we highlight many ways in which yeast has served to catalyze the fields of functional genomics, genome editing, gene-environment interaction investigation, proteomics, and bioinformatics-emphasizing how yeast has served as a catalyst for innovation. Several possible futures for this model organism in synthetic biology, drug personalization, and multi-omics research are also presented.

John N, Zaidi F, Doyle F, Fasullo M. High-throughput screening of the Saccharomyces cerevisiae genome for 2-amino-3-methylimidazo [4,5-f] quinoline resistance identifies colon cancer-associated genes

Heterocyclic aromatic amines (HAAs) are potent carcinogenic agents found in charred meats and cigarette smoke. However, few eukaryotic resistance genes have been identified. We used Saccharomyces cerevisiae (budding yeast) to identify genes that confer resistance to 2-amino-3-methylimidazo[4,5-f] quinoline (IQ). CYP1A2 and NAT2 activate IQ to become a mutagenic nitrenium compound. Deletion libraries expressing human CYP1A2 and NAT2 or no human genes were exposed to either 400 or 800 µM IQ for 5 or 10 generations. DNA barcodes were sequenced using the Illumina HiSeq 2500 platform and statistical significance was determined for exactly matched barcodes. We identified 424 ORFs, including 337 genes of known function, in duplicate screens of the "humanized" collection for IQ resistance; resistance was further validated for a select group of 51 genes by growth curves, competitive growth, or trypan blue assays. Screens of the library not expressing human genes identified 143 ORFs conferring resistance to IQ per se. Ribosomal protein and protein modification genes were identified as IQ resistance genes in both the original and "humanized" libraries, while nitrogen metabolism, DNA repair, and growth control genes were also prominent in the "humanized" library. Protein complexes identified included the casein kinase 2 (CK2) and histone chaperone (HIR) complex. Among DNA Repair and checkpoint genes, we identified those that function in postreplication repair (RAD18, UBC13, REV7), base excision repair (NTG1), and checkpoint signaling (CHK1, PSY2). These studies underscore the role of ribosomal protein genes in conferring IQ resistance, and illuminate DNA repair pathways for conferring resistance to activated IQ.

How colonization bottlenecks, tissue niches, and transmission strategies shape protozoan infections

Protozoan pathogens such as Plasmodium spp., Leishmania spp., Toxoplasma gondii, and Trypanosoma spp. are often associated with high-mortality, acute and chronic diseases of global health concern. For transmission and immune evasion, protozoans have evolved diverse strategies to interact with a range of host tissue environments. These interactions are linked to disease pathology, yet our understanding of the association between parasite colonization and host homeostatic disruption is limited. Recently developed techniques for cellular barcoding have the potential to uncover the biology regulating parasite transmission, dissemination, and the stability of infection. Understanding bottlenecks to infection and the in vivo tissue niches that facilitate chronic infection and spread has the potential to reveal new aspects of parasite biology.

Resolving Deleterious and Near-Neutral Effects Requires Different Pooled Fitness Assay Designs

Pooled sequencing-based fitness assays are a powerful and widely used approach to quantifying fitness of thousands of genetic variants in parallel. Despite the throughput of such assays, they are prone to biases in fitness estimates, and errors in measurements are typically larger for deleterious fitness effects, relative to neutral effects. In practice, designing pooled fitness assays involves tradeoffs between the number of timepoints, the sequencing depth, and other parameters to gain as much information as possible within a feasible experiment. Here, we combined simulations and reanalysis of an existing experimental dataset to explore how assay parameters impact measurements of near-neutral and deleterious fitness effects using a standard fitness estimator. We found that sequencing multiple timepoints at relatively modest depth improved estimates of near-neutral fitness effects, but systematically biased measurements of deleterious effects. We showed that a fixed total number of reads, deeper sequencing at fewer timepoints improved resolution of deleterious fitness effects. Our results highlight a tradeoff between measurement of deleterious and near-neutral effect sizes for a fixed amount of data and suggest that fitness assay design should be tuned for fitness effects that are relevant to the specific biological question.

Fit-Seq2.0: An Improved Software for High-Throughput Fitness Measurements Using Pooled Competition Assays

The fitness of a genotype is defined as its lifetime reproductive success, with fitness itself being a composite trait likely dependent on many underlying phenotypes. Measuring fitness is important for understanding how alteration of different cellular components affects a cell's ability to reproduce. Here, we describe an improved approach, implemented in Python, for estimating fitness in high throughput via pooled competition assays.

Best Practices in Designing, Sequencing, and Identifying Random DNA Barcodes

Random DNA barcodes are a versatile tool for tracking cell lineages, with applications ranging from development to cancer to evolution. Here, we review and critically evaluate barcode designs as well as methods of barcode sequencing and initial processing of barcode data. We first demonstrate how various barcode design decisions affect data quality and propose a new design that balances all considerations that we are currently aware of. We then discuss various options for the preparation of barcode sequencing libraries, including inline indices and Unique Molecular Identifiers (UMIs). Finally, we test the performance of several established and new bioinformatic pipelines for the extraction of barcodes from raw sequencing reads and for error correction. We find that both alignment and regular expression-based approaches work well for barcode extraction, and that error-correction pipelines designed specifically for barcode data are superior to generic ones. Overall, this review will help researchers to approach their barcoding experiments in a deliberate and systematic way.

Functional analysis of the Aspergillus fumigatus kinome reveals a DYRK kinase involved in septal plugging is a novel antifungal drug target

More than 10 million people suffer from lung diseases caused by the pathogenic fungus Aspergillus fumigatus. The azole class of antifungals represent first line therapeutics for most of these infections however resistance is rising. Identification of novel antifungal targets that, when inhibited, synergise with the azoles will aid the development of agents that can improve therapeutic outcomes and supress the emergence of resistance. As part of the A. fumigatus genome-wide knockout program (COFUN), we have completed the generation of a library that consists of 120 genetically barcoded null mutants in genes that encode the protein kinase cohort of A. fumigatus. We have employed a competitive fitness profiling approach (Bar-Seq), to identify targets which when deleted result in hypersensitivity to the azoles and fitness defects in a murine host. The most promising candidate from our screen is a previously uncharacterised DYRK kinase orthologous to Yak1 of Candida albicans, a TOR signalling pathway kinase involved in modulation of stress responsive transcriptional regulators. Here we show that the orthologue YakA has been repurposed in A. fumigatus to regulate blocking of the septal pore upon exposure to stress via phosphorylation of the Woronin body tethering protein Lah. Loss of YakA function reduces the ability of A. fumigatus to penetrate solid media and impacts growth in murine lung tissue. We also show that 1-ethoxycarbonyl-beta-carboline (1-ECBC), a compound previously shown to inhibit Yak1 in C. albicans prevents stress mediated septal spore blocking and synergises with the azoles to inhibit A. fumigatus growth.

Large-scale genetic characterization of the model sulfate-reducing bacterium, Desulfovibrio vulgaris Hildenborough

Sulfate-reducing bacteria (SRB) are obligate anaerobes that can couple their growth to the reduction of sulfate. Despite the importance of SRB to global nutrient cycles and their damage to the petroleum industry, our molecular understanding of their physiology remains limited. To systematically provide new insights into SRB biology, we generated a randomly barcoded transposon mutant library in the model SRB Desulfovibrio vulgaris Hildenborough (DvH) and used this genome-wide resource to assay the importance of its genes under a range of metabolic and stress conditions. In addition to defining the essential gene set of DvH, we identified a conditional phenotype for 1,137 non-essential genes. Through examination of these conditional phenotypes, we were able to make a number of novel insights into our molecular understanding of DvH, including how this bacterium synthesizes vitamins. For example, we identified DVU0867 as an atypical L-aspartate decarboxylase required for the synthesis of pantothenic acid, provided the first experimental evidence that biotin synthesis in DvH occurs via a specialized acyl carrier protein and without methyl esters, and demonstrated that the uncharacterized dehydrogenase DVU0826:DVU0827 is necessary for the synthesis of pyridoxal phosphate. In addition, we used the mutant fitness data to identify genes involved in the assimilation of diverse nitrogen sources and gained insights into the mechanism of inhibition of chlorate and molybdate. Our large-scale fitness dataset and RB-TnSeq mutant library are community-wide resources that can be used to generate further testable hypotheses into the gene functions of this environmentally and industrially important group of bacteria.

High-throughput functional profiling of the human fungal pathogen Candida albicans genome

Candida albicans is a major fungal pathogen of humans. Although its genome has been sequenced more than two decades ago, there are still over 4300 uncharacterized C. albicans genes. We previously generated an ORFeome as well as a collection of destination vectors to facilitate overexpression of C. albicans ORFs. Here, we report the construction of ∼2500 overexpression mutants and their evaluation by in vitro spotting on rich medium and in a liquid pool experiment in rich medium, allowing the identification of genes whose overexpression has a fitness cost. The candidates were further validated at the individual strain level. This new resource allows large-scale screens in different growth conditions to be performed routinely. Altogether, based on the concept of identifying functionally related genes by cluster analysis, the availability of this overexpression mutant collection will facilitate the characterization of gene functions in C. albicans.

A universal sequencing read interpreter

Massively parallel DNA sequencing has led to the rapid growth of highly multiplexed experiments in biology. These experiments produce unique sequencing results that require specific analysis pipelines to decode highly structured reads. However, no versatile framework that interprets sequencing reads to extract their encoded information for downstream biological analysis has been developed. Here, we report INTERSTELLAR (interpretation, scalable transformation, and emulation of large-scale sequencing reads) that decodes data values encoded in theoretically any type of sequencing read and translates them into sequencing reads of another structure of choice. We demonstrated that INTERSTELLAR successfully extracted information from a range of short- and long-read sequencing reads and translated those of single-cell (sc)RNA-seq, scATAC-seq, and spatial transcriptomics to be analyzed by different software tools that have been developed for conceptually the same types of experiments. INTERSTELLAR will greatly facilitate the development of sequencing-based experiments and sharing of data analysis pipelines.

A multiplexed, three-dimensional pooling and next-generation sequencing strategy for creating barcoded mutant arrays: construction of a Schizosaccharomyces pombe transposon insertion library

Arrayed libraries of defined mutants have been used to elucidate gene function in the post-genomic era. Yeast haploid gene deletion libraries have pioneered this effort, but are costly to construct, do not reveal phenotypes that may occur with partial gene function and lack essential genes required for growth. We therefore devised an efficient method to construct a library of barcoded insertion mutants with a wider range of phenotypes that can be generalized to other organisms or collections of DNA samples. We developed a novel but simple three-dimensional pooling and multiplexed sequencing approach that leveraged sequence information to reduce the number of required sequencing reactions by orders of magnitude, and were able to identify the barcode sequences and DNA insertion sites of 4391 Schizosaccharomyces pombe insertion mutations with only 40 sequencing preparations. The insertion mutations are in the genes and untranslated regions of nonessential, essential and noncoding RNA genes, and produced a wider range of phenotypes compared to the cognate deletion mutants, including novel phenotypes. This mutant library represents both a proof of principle for an efficient method to produce novel mutant libraries and a valuable resource for the S. pombe research community.

Cellular barcoding of protozoan pathogens reveals the within-host population dynamics of Toxoplasma gondii host colonization

Cellular barcoding techniques are powerful tools to understand microbial pathogenesis. However, barcoding strategies have not been broadly applied to protozoan parasites, which have unique genomic structures and virulence strategies compared with viral and bacterial pathogens. Here, we present a CRISPR-based method to barcode protozoa, which we successfully apply to Toxoplasma gondii and Trypanosoma brucei. Using libraries of barcoded T. gondii, we evaluate shifts in the population structure from acute to chronic infection of mice. Contrary to expectation, most barcodes were present in the brain one month post-intraperitoneal infection in both inbred CBA/J and outbred Swiss mice. Although parasite cyst number and barcode diversity declined over time, barcodes representing a minor fraction of the inoculum could become a dominant population in the brain by three months post-infection. These data establish a cellular barcoding approach for protozoa and evidence that the blood-brain barrier is not a major bottleneck to colonization by T. gondii.

Distinct structural groups of histone H3 and H4 residues have divergent effects on chronological lifespan in Saccharomyces cerevisiae

We have performed a comprehensive analysis of the involvement of histone H3 and H4 residues in the regulation of chronological lifespan in yeast and identify four structural groups in the nucleosome that influence lifespan. We also identify residues where substitution with an epigenetic mimic extends lifespan, providing evidence that a simple epigenetic switch, without possible additional background modifications, causes longevity. Residues where substitution result in the most pronounced lifespan extension are all on the exposed face of the nucleosome, with the exception of H3E50, which is present on the lateral surface, between two DNA gyres. Other residues that have a more modest effect on lifespan extension are concentrated at the extremities of the H3-H4 dimer, suggesting a role in stabilizing the dimer in its nucleosome frame. Residues that reduce lifespan are buried in the histone handshake motif, suggesting that these mutations destabilize the octamer structure. All residues exposed on the nucleosome disk face and that cause lifespan extension are known to interact with Sir3. We find that substitution of H4K16 and H4H18 cause Sir3 to redistribute from telomeres and silent mating loci to secondary positions, often enriched for Rap1, Abf1 or Reb1 binding sites, whereas H3E50 does not. The redistribution of Sir3 in the genome can be reproduced by an equilibrium model based on primary and secondary binding sites with different affinities for Sir3. The redistributed Sir3 cause transcriptional repression at most of the new loci, including of genes where null mutants were previously shown to extend chronological lifespan. The transcriptomic profiles of H4K16 and H4H18 mutant strains are very similar, and compatible with a DNA replication stress response. This is distinct from the transcriptomic profile of H3E50, which matches strong induction of oxidative phosphorylation. We propose that the different groups of residues are involved in binding to heterochromatin proteins, in destabilizing the association of the nucleosome DNA, disrupting binding of the H3-H4 dimer in the nucleosome, or disrupting the structural stability of the octamer, each category impacting on chronological lifespan by a different mechanism.

An Integrated Approach Reveals DNA Damage and Proteotoxic Stress as Main Effects of Proton Radiation in S. cerevisiae

Proton radiotherapy (PRT) has the potential to reduce the normal tissue toxicity associated with conventional photon-based radiotherapy (X-ray therapy, XRT) because the active dose can be more directly targeted to a tumor. Although this dosimetric advantage of PRT is well known, the molecular mechanisms affected by PRT remain largely elusive. Here, we combined the molecular toolbox of the eukaryotic model Saccharomyces cerevisiae with a systems biology approach to investigate the physiological effects of PRT compared to XRT. Our data show that the DNA damage response and protein stress response are the major molecular mechanisms activated after both PRT and XRT. However, RNA-Seq revealed that PRT treatment evoked a stronger activation of genes involved in the response to proteotoxic stress, highlighting the molecular differences between PRT and XRT. Moreover, inhibition of the proteasome resulted in decreased survival in combination with PRT compared to XRT, not only further confirming that protons induced a stronger proteotoxic stress response, but also hinting at the potential of using proteasome inhibitors in combination with proton radiotherapy in clinical settings.

Barcoded Bulk QTL mapping reveals highly polygenic and epistatic architecture of complex traits in yeast

Mapping the genetic basis of complex traits is critical to uncovering the biological mechanisms that underlie disease and other phenotypes. Genome-wide association studies (GWAS) in humans and quantitative trait locus (QTL) mapping in model organisms can now explain much of the observed heritability in many traits, allowing us to predict phenotype from genotype. However, constraints on power due to statistical confounders in large GWAS and smaller sample sizes in QTL studies still limit our ability to resolve numerous small-effect variants, map them to causal genes, identify pleiotropic effects across multiple traits, and infer non-additive interactions between loci (epistasis). Here, we introduce barcoded bulk quantitative trait locus (BB-QTL) mapping, which allows us to construct, genotype, and phenotype 100,000 offspring of a budding yeast cross, two orders of magnitude larger than the previous state of the art. We use this panel to map the genetic basis of eighteen complex traits, finding that the genetic architecture of these traits involves hundreds of small-effect loci densely spaced throughout the genome, many with widespread pleiotropic effects across multiple traits. Epistasis plays a central role, with thousands of interactions that provide insight into genetic networks. By dramatically increasing sample size, BB-QTL mapping demonstrates the potential of natural variants in high-powered QTL studies to reveal the highly polygenic, pleiotropic, and epistatic architecture of complex traits.

MegaGate: A toxin-less gateway molecular cloning tool

Gateway cloning employs the use of the ccdb toxin and has low colony numbers, making it difficult to apply at scale to clone libraries of cDNA vectors. In this protocol, we describe MegaGate, a toxin-less Gateway technology capable of robust cDNA library cloning that is efficient, cheap, and scalable. MegaGate eliminates the ccdb toxin used in Gateway recombinase cloning and instead utilizes meganuclease-mediated digestion to eliminate background vectors during cloning and is 99.8% efficient with high colony numbers. For complete details on the use and execution of this protocol, please refer to Kramme et al. (2021).

An integrated pipeline for mammalian genetic screening

With the recent advancements in genome editing, next-generation sequencing (NGS), and scalable cloning techniques, scientists can now conduct genetic screens at unprecedented levels of scale and precision. With such a multitude of technologies, there is a need for a simple yet comprehensive pipeline to enable systematic mammalian genetic screening. In this study, we develop unique algorithms for target identification and a toxin-less Gateway cloning tool, termed MegaGate, for library cloning which, when combined with existing genetic perturbation methods and NGS-coupled readouts, enable versatile engineering of relevant mammalian cell lines. Our integrated pipeline for sequencing-based target ascertainment and modular perturbation screening (STAMPScreen) can thus be utilized for a host of cell state engineering applications.

Investigation of RNA metabolism through large-scale genetic interaction profiling in yeast

Gene deletion and gene expression alteration can lead to growth defects that are amplified or reduced when a second mutation is present in the same cells. We performed 154 genetic interaction mapping (GIM) screens with query mutants related with RNA metabolism and estimated the growth rates of about 700 000 double mutant Saccharomyces cerevisiae strains. The tested targets included the gene deletion collection and 900 strains in which essential genes were affected by mRNA destabilization (DAmP). To analyze the results, we developed RECAP, a strategy that validates genetic interaction profiles by comparison with gene co-citation frequency, and identified links between 1471 genes and 117 biological processes. In addition to these large-scale results, we validated both enhancement and suppression of slow growth measured for specific RNA-related pathways. Thus, negative genetic interactions identified a role for the OCA inositol polyphosphate hydrolase complex in mRNA translation initiation. By analysis of suppressors, we found that Puf4, a Pumilio family RNA binding protein, inhibits ribosomal protein Rpl9 function, by acting on a conserved UGUAcauUA motif located downstream the stop codon of the RPL9B mRNA. Altogether, the results and their analysis should represent a useful resource for discovery of gene function in yeast.

A genome-wide portrait of pervasive drug contaminants

Using a validated yeast chemogenomic platform, we characterized the genome-wide effects of several pharmaceutical contaminants, including three N-nitrosamines (NDMA, NDEA and NMBA), two related compounds (DMF and 4NQO) and several of their metabolites. A collection of 4800 non-essential homozygous diploid yeast deletion strains were screened in parallel and the strain abundance was quantified by barcode sequencing. These data were used to rank deletion strains representing genes required for resistance to the compounds to delineate affected cellular pathways and to visualize the global cellular effects of these toxins in an easy-to-use searchable database. Our analysis of the N-nitrosamine screens uncovered genes (via their corresponding homozygous deletion mutants) involved in several evolutionarily conserved pathways, including: arginine biosynthesis, mitochondrial genome integrity, vacuolar protein sorting and DNA damage repair. To investigate why NDMA, NDEA and DMF caused fitness defects in strains lacking genes of the arginine pathway, we tested several N-nitrosamine metabolites (methylamine, ethylamine and formamide), and found they also affected arginine pathway mutants. Notably, each of these metabolites has the potential to produce ammonium ions during their biotransformation. We directly tested the role of ammonium ions in N-nitrosamine toxicity by treatment with ammonium sulfate and we found that ammonium sulfate also caused a growth defect in arginine pathway deletion strains. Formaldehyde, a metabolite produced from NDMA, methylamine and formamide, and which is known to cross-link free amines, perturbed deletion strains involved in chromatin remodeling and DNA repair pathways. Finally, co-administration of N-nitrosamines with ascorbic or ferulic acid did not relieve N-nitrosamine toxicity. In conclusion, we used parallel deletion mutant analysis to characterize the genes and pathways affected by exposure to N-nitrosamines and related compounds, and provide the data in an accessible, queryable database.

A genome-scale yeast library with inducible expression of individual genes

The ability to switch a gene from off to on and monitor dynamic changes provides a powerful approach for probing gene function and elucidating causal regulatory relationships. Here, we developed and characterized YETI (Yeast Estradiol strains with Titratable Induction), a collection in which > 5,600 yeast genes are engineered for transcriptional inducibility with single-gene precision at their native loci and without plasmids. Each strain contains SGA screening markers and a unique barcode, enabling high-throughput genetics. We characterized YETI using growth phenotyping and BAR-seq screens, and we used a YETI allele to identify the regulon of Rof1, showing that it acts to repress transcription. We observed that strains with inducible essential genes that have low native expression can often grow without inducer. Analysis of data from eukaryotic and prokaryotic systems shows that native expression is a variable that can bias promoter-perturbing screens, including CRISPRi. We engineered a second expression system, Z3 EB42, that gives lower expression than Z3 EV, a feature enabling conditional activation and repression of lowly expressed essential genes that grow without inducer in the YETI library.

Gene Loss Predictably Drives Evolutionary Adaptation

Loss of gene function is common throughout evolution, even though it often leads to reduced fitness. In this study, we systematically evaluated how an organism adapts after deleting genes that are important for growth under oxidative stress. By evolving, sequencing, and phenotyping over 200 yeast lineages, we found that gene loss can enhance an organism’s capacity to evolve and adapt. Although gene loss often led to an immediate decrease in fitness, many mutants rapidly acquired suppressor mutations that restored fitness. Depending on the strain’s genotype, some ultimately even attained higher fitness levels than similarly adapted wild-type cells. Further, cells with deletions in different modules of the genetic network followed distinct and predictable mutational trajectories. Finally, losing highly connected genes increased evolvability by facilitating the emergence of a more diverse array of phenotypes after adaptation. Together, our findings show that loss of specific parts of a genetic network can facilitate adaptation by opening alternative evolutionary paths.

Spurious regulatory connections dictate the expression-fitness landscape of translation factors

During steady-state cell growth, individual enzymatic fluxes can be directly inferred from growth rate by mass conservation, but the inverse problem remains unsolved. Perturbing the flux and expression of a single enzyme could have pleiotropic effects that may or may not dominate the impact on cell fitness. Here, we quantitatively dissect the molecular and global responses to varied expression of translation termination factors (peptide release factors, RFs) in the bacterium Bacillus subtilis. While endogenous RF expression maximizes proliferation, deviations in expression lead to unexpected distal regulatory responses that dictate fitness reduction. Molecularly, RF depletion causes expression imbalance at specific operons, which activates master regulators and detrimentally overrides the transcriptome. Through these spurious connections, RF abundances are thus entrenched by focal points within the regulatory network, in one case located at a single stop codon. Such regulatory entrenchment suggests that predictive bottom-up models of expression-fitness landscapes will require near-exhaustive characterization of parts.

Inactivating histone deacetylase HDA promotes longevity by mobilizing trehalose metabolism

Histone acetylations are important epigenetic markers for transcriptional activation in response to metabolic changes and various stresses. Using the high-throughput SEquencing-Based Yeast replicative Lifespan screen method and the yeast knockout collection, we demonstrate that the HDA complex, a class-II histone deacetylase (HDAC), regulates aging through its target of acetylated H3K18 at storage carbohydrate genes. We find that, in addition to longer lifespan, disruption of HDA results in resistance to DNA damage and osmotic stresses. We show that these effects are due to increased promoter H3K18 acetylation and transcriptional activation in the trehalose metabolic pathway in the absence of HDA. Furthermore, we determine that the longevity effect of HDA is independent of the Cyc8-Tup1 repressor complex known to interact with HDA and coordinate transcriptional repression. Silencing the HDA homologs in C. elegans and Drosophila increases their lifespan and delays aging-associated physical declines in adult flies. Hence, we demonstrate that this HDAC controls an evolutionarily conserved longevity pathway.

Identification of Genetic Modifiers of TDP-43: Inflammatory Activation of Astrocytes for Neuroinflammation

Transactive response DNA-binding protein 43 (TDP-43) is a ubiquitously expressed DNA/RNA-binding protein linked to amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 has been implicated in numerous aspects of the mRNA life cycle, as well as in cell toxicity and neuroinflammation. In this study, we used the toxicity of the TDP-43 expression in Saccharomyces cerevisiae as an assay to identify TDP-43 genetic interactions. Specifically, we transformed human TDP-43 cDNAs of wild-type or disease-associated mutants (M337V and Q331K) en masse into 4653 homozygous diploid yeast deletion mutants and then used next-generation sequencing readouts of growth to identify yeast toxicity modifiers. Genetic interaction analysis provided a global view of TDP-43 pathways, some of which are known to be involved in cellular metabolic processes. Selected putative loci with the potential of genetic interactions with TDP-43 were assessed for associations with neurotoxicity and inflammatory activation of astrocytes. The pharmacological inhibition of succinate dehydrogenase flavoprotein subunit A (SDHA) and voltage-dependent anion-selective channel 3 (VDAC3) suppressed TDP-43-induced expression of proinflammatory cytokines in astrocytes, indicating the critical roles played by SDHA and VDAC3 in TDP-43 pathways during inflammatory activation of astrocytes and neuroinflammation. Thus, the findings of our TDP-43 genetic interaction screen provide a global landscape of TDP-43 pathways and may help improve our understanding of the roles of glia and neuroinflammation in ALS and FTD pathogenesis.

A genome-wide copper-sensitized screen identifies novel regulators of mitochondrial cytochrome c oxidase activity

Copper is essential for the activity and stability of cytochrome c oxidase (CcO), the terminal enzyme of the mitochondrial respiratory chain. Loss-of-function mutations in genes required for copper transport to CcO result in fatal human disorders. Despite the fundamental importance of copper in mitochondrial and organismal physiology, systematic identification of genes that regulate mitochondrial copper homeostasis is lacking. To discover these genes, we performed a genome-wide screen using a library of DNA-barcoded yeast deletion mutants grown in copper-supplemented media. Our screen recovered a number of genes known to be involved in cellular copper homeostasis as well as genes previously not linked to mitochondrial copper biology. These newly identified genes include the subunits of the adaptor protein 3 complex (AP-3) and components of the cellular pH-sensing pathway Rim20 and Rim21, both of which are known to affect vacuolar function. We find that AP-3 and Rim mutants exhibit decreased vacuolar acidity, which in turn perturbs mitochondrial copper homeostasis and CcO function. CcO activity of these mutants could be rescued by either restoring vacuolar pH or supplementing growth media with additional copper. Consistent with these genetic data, pharmacological inhibition of the vacuolar proton pump leads to decreased mitochondrial copper content and a concomitant decrease in CcO abundance and activity. Taken together, our study uncovered novel genetic regulators of mitochondrial copper homeostasis and provided a mechanism by which vacuolar pH impacts mitochondrial respiration through copper homeostasis.

Omics-Driven Biotechnology for Industrial Applications

Biomanufacturing is a key component of biotechnology that uses biological systems to produce bioproducts of commercial relevance, which are of great interest to the energy, material, pharmaceutical, food, and agriculture industries. Biotechnology-based approaches, such as synthetic biology and metabolic engineering are heavily reliant on "omics" driven systems biology to characterize and understand metabolic networks. Knowledge gained from systems biology experiments aid the development of synthetic biology tools and the advancement of metabolic engineering studies toward establishing robust industrial biomanufacturing platforms. In this review, we discuss recent advances in "omics" technologies, compare the pros and cons of the different "omics" technologies, and discuss the necessary requirements for carrying out multi-omics experiments. We highlight the influence of "omics" technologies on the production of biofuels and bioproducts by metabolic engineering. Finally, we discuss the application of "omics" technologies to agricultural and food biotechnology, and review the impact of "omics" on current COVID-19 research.

BaSDAS: a web-based pooled CRISPR-Cas9 knockout screening data analysis system

We developed the BaSDAS (Barcode-Seq Data Analysis System), a GUI-based pooled knockout screening data analysis system, to facilitate the analysis of pooled knockout screen data easily and effectively by researchers with limited bioinformatics skills. The BaSDAS supports the analysis of various pooled screening libraries, including yeast, human, and mouse libraries, and provides many useful statistical and visualization functions with a user-friendly web interface for convenience. We expect that BaSDAS will be a useful tool for the analysis of genome-wide screening data and will support the development of novel drugs based on functional genomics information.

Multiplexed Competition in a Synthetic Squid Light Organ Microbiome Using Barcode-Tagged Gene Deletions

Beneficial microbes play essential roles in the health and development of their hosts. However, the complexity of animal microbiomes and general genetic intractability of their symbionts have made it difficult to study the coevolved mechanisms for establishing and maintaining specificity at the microbe-animal host interface.

Beneficial symbioses between microbes and their eukaryotic hosts are ubiquitous and have widespread impacts on host health and development. The binary symbiosis between the bioluminescent bacterium Vibrio fischeri and its squid host Euprymna scolopes serves as a model system to study molecular mechanisms at the microbe-animal interface. To identify colonization factors in this system, our lab previously conducted a global transposon insertion sequencing (INSeq) screen and identified over 300 putative novel squid colonization factors in V. fischeri. To pursue mechanistic studies on these candidate genes, we present an approach to quickly generate barcode-tagged gene deletions and perform high-throughput squid competition experiments with detection of the proportion of each strain in the mixture by barcode sequencing (BarSeq). Our deletion approach improves on previous techniques based on splicing by overlap extension PCR (SOE-PCR) and tfoX-based natural transformation by incorporating a randomized barcode that results in unique DNA sequences within each deletion scar. Amplicon sequencing of the pool of barcoded strains before and after colonization faithfully reports on known colonization factors and provides increased sensitivity over colony counting methods. BarSeq enables rapid and sensitive characterization of the molecular factors involved in establishing the Vibrio-squid symbiosis and provides a valuable tool to interrogate the molecular dialogue at microbe-animal host interfaces.

IMPORTANCE Beneficial microbes play essential roles in the health and development of their hosts. However, the complexity of animal microbiomes and general genetic intractability of their symbionts have made it difficult to study the coevolved mechanisms for establishing and maintaining specificity at the microbe-animal host interface. Model symbioses are therefore invaluable for studying the mechanisms of beneficial microbe-host interactions. Here, we present a combined barcode-tagged deletion and BarSeq approach to interrogate the molecular dialogue that ensures specific and reproducible colonization of the Hawaiian bobtail squid by Vibrio fischeri. The ability to precisely manipulate the bacterial genome, combined with multiplex colonization assays, will accelerate the use of this valuable model system for mechanistic studies of how environmental microbes—both beneficial and pathogenic—colonize specific animal hosts.

Interrogation of kinase genetic interactions provides a global view of PAK1-mediated signal transduction pathways

Kinases are critical components of intracellular signaling pathways and have been extensively investigated with regard to their roles in cancer. p21-activated kinase-1 (PAK1) is a serine/threonine kinase that has been previously implicated in numerous biological processes, such as cell migration, cell cycle progression, cell motility, invasion, and angiogenesis, in glioma and other cancers. However, the signaling network linked to PAK1 is not fully defined. We previously reported a large-scale yeast genetic interaction screen using toxicity as a readout to identify candidate PAK1 genetic interactions. En masse transformation of the PAK1 gene into 4,653 homozygous diploid Saccharomyces cerevisiae yeast deletion mutants identified ∼400 candidates that suppressed yeast toxicity. Here we selected 19 candidate PAK1 genetic interactions that had human orthologs and were expressed in glioma for further examination in mammalian cells, brain slice cultures, and orthotopic glioma models. RNAi and pharmacological inhibition of potential PAK1 interactors confirmed that DPP4, KIF11, mTOR, PKM2, SGPP1, TTK, and YWHAE regulate PAK1-induced cell migration and revealed the importance of genes related to the mitotic spindle, proteolysis, autophagy, and metabolism in PAK1-mediated glioma cell migration, drug resistance, and proliferation. AKT1 was further identified as a downstream mediator of the PAK1-TTK genetic interaction. Taken together, these data provide a global view of PAK1-mediated signal transduction pathways and point to potential new drug targets for glioma therapy.

PhenoMIP: High-Throughput Phenotyping of Diverse Caenorhabditis elegans Populations via Molecular Inversion Probes

Whether generated within a lab setting or isolated from the wild, variant alleles continue to be an important resource for decoding gene function in model organisms such as Caenorhabditis elegans. With advances in massively parallel sequencing, multiple whole-genome sequenced (WGS) strain collections are now available to the research community. The Million Mutation Project (MMP) for instance, analyzed 2007 N2-derived, mutagenized strains. Individually, each strain averages ∼400 single nucleotide variants amounting to ∼80 protein-coding variants. The effects of these variants, however, remain largely uncharacterized and querying the breadth of these strains for phenotypic changes requires a method amenable to rapid and sensitive high-throughput analysis. Here we present a pooled competitive fitness approach to quantitatively phenotype subpopulations of sequenced collections via molecular inversion probes (PhenoMIP). We phenotyped the relative fitness of 217 mutant strains on multiple food sources and classified these into five categories. We also demonstrate on a subset of these strains, that their fitness defects can be genetically mapped. Overall, our results suggest that approximately 80% of MMP mutant strains may have a decreased fitness relative to the lab reference, N2. The costs of generating this form of analysis through WGS methods would be prohibitive while PhenoMIP analysis in this manner is accomplished at less than one-tenth of projected WGS costs. We propose methods for applying PhenoMIP to a broad range of population selection experiments in a cost-efficient manner that would be useful to the community at large.

High-Throughput Functional Genomics for Energy Production

Functional genomics remains a foundational field for establishing genotype-phenotype relationships that enable strain engineering. High-throughput (HTP) methods accelerate the Design-Build-Test-Learn cycle that currently drives synthetic biology towards a forward engineering future. Trackable mutagenesis techniques including transposon insertion sequencing and CRISPR-Cas-mediated genome editing allow for rapid fitness profiling of a collection, or library, of mutants to discover beneficial mutations. Due to the relative speed of these experiments compared to adaptive evolution experiments, iterative rounds of mutagenesis can be implemented for next-generation metabolic engineering efforts to design complex production and tolerance phenotypes. Additionally, the expansion of these mutagenesis techniques to novel bacteria are opening up industrial microbes that show promise for establishing a bio-based economy.

Leishmania differentiation requires ubiquitin conjugation mediated by a UBC2-UEV1 E2 complex

Post-translational modifications such as ubiquitination are important for orchestrating the cellular transformations that occur as the Leishmania parasite differentiates between its main morphological forms, the promastigote and amastigote. 2 E1 ubiquitin-activating (E1), 13 E2 ubiquitin-conjugating (E2), 79 E3 ubiquitin ligase (E3) and 20 deubiquitinating cysteine peptidase (DUB) genes can be identified in the Leishmania mexicana genome but, currently, little is known about the role of E1, E2 and E3 enzymes in this parasite. Bar-seq analysis of 23 E1, E2 and HECT/RBR E3 null mutants generated in promastigotes using CRISPR-Cas9 revealed numerous loss-of-fitness phenotypes in promastigote to amastigote differentiation and mammalian infection. The E2s UBC1/CDC34, UBC2 and UEV1 and the HECT E3 ligase HECT2 are required for the successful transformation from promastigote to amastigote and UBA1b, UBC9, UBC14, HECT7 and HECT11 are required for normal proliferation during mouse infection. Of all ubiquitination enzyme null mutants examined in the screen, Δubc2 and Δuev1 exhibited the most extreme loss-of-fitness during differentiation. Null mutants could not be generated for the E1 UBA1a or the E2s UBC3, UBC7, UBC12 and UBC13, suggesting these genes are essential in promastigotes. X-ray crystal structure analysis of UBC2 and UEV1, orthologues of human UBE2N and UBE2V1/UBE2V2 respectively, reveal a heterodimer with a highly conserved structure and interface. Furthermore, recombinant L. mexicana UBA1a can load ubiquitin onto UBC2, allowing UBC2-UEV1 to form K63-linked di-ubiquitin chains in vitro. Notably, UBC2 can cooperate in vitro with human E3s RNF8 and BIRC2 to form non-K63-linked polyubiquitin chains, showing that UBC2 can facilitate ubiquitination independent of UEV1, but association of UBC2 with UEV1 inhibits this ability. Our study demonstrates the dual essentiality of UBC2 and UEV1 in the differentiation and intracellular survival of L. mexicana and shows that the interaction between these two proteins is crucial for regulation of their ubiquitination activity and function.

The post-translational modification of proteins is key for allowing Leishmania parasites to transition between the different life cycle stages that exist in its insect vector and mammalian host. In particular, components of the ubiquitin system are important for the transformation of Leishmania from its insect (promastigote) to mammalian (amastigote) stage and normal infection in mice. However, little is known about the role of the enzymes that generate ubiquitin modifications in Leishmania. Here we characterise 28 enzymes of the ubiquitination pathway and show that many are required for life cycle progression or mouse infection by this parasite. Two proteins, UBC2 and UEV1, were selected for further study based on their importance in the promastigote to amastigote transition. We demonstrate that UBC2 and UEV1 form a heterodimer capable of carrying out ubiquitination and that the structural basis for this activity is conserved between Leishmania, Saccharomyces cerevisiae and humans. We also show that the interaction of UBC2 with UEV1 alters the nature of the ubiquitination activity performed by UBC2. Overall, we demonstrate the important role that ubiquitination enzymes play in the life cycle and infection process of Leishmania and explore the biochemistry underlying UBC2 and UEV1 function.

High-throughput analysis of adaptation using barcoded strains of Saccharomyces cerevisiae

BACKGROUND

Experimental evolution of microbes can be used to empirically address a wide range of questions about evolution and is increasingly employed to study complex phenomena ranging from genetic evolution to evolutionary rescue. Regardless of experimental aims, fitness assays are a central component of this type of research, and low-throughput often limits the scope and complexity of experimental evolution studies. We created an experimental evolution system in Saccharomyces cerevisiae that utilizes genetic barcoding to overcome this challenge.

RESULTS

We first confirm that barcode insertions do not alter fitness and that barcode sequencing can be used to efficiently detect fitness differences via pooled competition-based fitness assays. Next, we examine the effects of ploidy, chemical stress, and population bottleneck size on the evolutionary dynamics and fitness gains (adaptation) in a total of 76 experimentally evolving, asexual populations by conducting 1,216 fitness assays and analyzing 532 longitudinal-evolutionary samples collected from the evolving populations. In our analysis of these data we describe the strengths of this experimental evolution system and explore sources of error in our measurements of fitness and evolutionary dynamics.

CONCLUSIONS

Our experimental treatments generated distinct fitness effects and evolutionary dynamics, respectively quantified via multiplexed fitness assays and barcode lineage tracking. These findings demonstrate the utility of this new resource for designing and improving high-throughput studies of experimental evolution. The approach described here provides a framework for future studies employing experimental designs that require high-throughput multiplexed fitness measurements.

Multi-Faceted Systems Biology Approaches Present a Cellular Landscape of Phenolic Compound Inhibition in Saccharomyces cerevisiae

Synthetic biology has played a major role in engineering microbial cell factories to convert plant biomass (lignocellulose) to fuels and bioproducts by fermentation. However, the final product yield is limited by inhibition of microbial growth and fermentation by toxic phenolic compounds generated during lignocellulosic pre-treatment and hydrolysis. Advances in the development of systems biology technologies (genomics, transcriptomics, proteomics, metabolomics) have rapidly resulted in large datasets which are necessary to obtain a holistic understanding of complex biological processes underlying phenolic compound toxicity. Here, we review and compare different systems biology tools that have been utilized to identify molecular mechanisms that modulate phenolic compound toxicity in Saccharomyces cerevisiae. By focusing on and comparing functional genomics and transcriptomics approaches we identify common mechanisms potentially underlying phenolic toxicity. Additionally, we discuss possible ways by which integration of data obtained across multiple unbiased approaches can result in new avenues to develop yeast strains with a significant improvement in tolerance to phenolic fermentation inhibitors.

High-throughput mapping of the phage resistance landscape in E. coli

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage-specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

Measurements drive progress in directed evolution for precise engineering of biological systems

Precise engineering of biological systems requires quantitative, high-throughput measurements, exemplified by progress in directed evolution. New approaches allow high-throughput measurements of phenotypes and their corresponding genotypes. When integrated into directed evolution, these quantitative approaches enable the precise engineering of biological function. At the same time, the increasingly routine availability of large, high-quality data sets supports the integration of machine learning with directed evolution. Together, these advances herald striking capabilities for engineering biology.

A Signature-Tagged Mutagenesis (STM)-based murine-infectivity assay for Cryptococcus neoformans

Signature-tagged mutagenesis (STM) is a high-throughput genetic technique that can be used to investigate the function of genes by constructing a large number of mutant strains with unique DNA identification tags, pooling them, and screening them for a particular phenotypic trait. STM was first designed for the identification of genes that contribute to the virulence or infectivity of a pathogen in its host. Recently, this method has also been applied for the identification of mutants with specific phenotypes, such as antifungal drug resistance and proliferation. In the present study, we describe an STM method for the identification of genes contributing to the infectivity of Cryptococcus neoformans using a mutant library, in which each strain was tagged with a unique DNA sequence.

Investigation of Bar-seq as a method to study population dynamics of Saccharomyces cerevisiae deletion library during bioreactor cultivation

Background

Despite the latest advancements in metabolic engineering for genome editing and characterization of host performance, the successful development of robust cell factories used for industrial bioprocesses and accurate prediction of the behavior of microbial systems, especially when shifting from laboratory-scale to industrial conditions, remains challenging. To increase the probability of success of a scale-up process, data obtained from thoroughly performed studies mirroring cellular responses to typical large-scale stimuli may be used to derive crucial information to better understand potential implications of large-scale cultivation on strain performance. This study assesses the feasibility to employ a barcoded yeast deletion library to assess genome-wide strain fitness across a simulated industrial fermentation regime and aims to understand the genetic basis of changes in strain physiology during industrial fermentation, and the corresponding roles these genes play in strain performance.

Results

We find that mutant population diversity is maintained through multiple seed trains, enabling large scale fermentation selective pressures to act upon the community. We identify specific deletion mutants that were enriched in all processes tested in this study, independent of the cultivation conditions, which include MCK1, RIM11, MRK1, and YGK3 that all encode homologues of mammalian glycogen synthase kinase 3 (GSK-3). Ecological analysis of beta diversity between all samples revealed significant population divergence over time and showed feed specific consequences of population structure. Further, we show that significant changes in the population diversity during fed-batch cultivations reflect the presence of significant stresses. Our observations indicate that, for this yeast deletion collection, the selection of the feeding scheme which affects the accumulation of the fermentative by-product ethanol impacts the diversity of the mutant pool to a higher degree as compared to the pH of the culture broth. The mutants that were lost during the time of most extreme population selection suggest that specific biological processes may be required to cope with these specific stresses.

Conclusions

Our results demonstrate the feasibility of Bar-seq to assess fermentation associated stresses in yeast populations under industrial conditions and to understand critical stages of a scale-up process where variability emerges, and selection pressure gets imposed. Overall our work highlights a promising avenue to identify genetic loci and biological stress responses required for fitness under industrial conditions.

Model-driven analysis of mutant fitness experiments improves genome-scale metabolic models of Zymomonas mobilis ZM4

Genome-scale metabolic models have been utilized extensively in the study and engineering of the organisms they describe. Here we present the analysis of a published dataset from pooled transposon mutant fitness experiments as an approach for improving the accuracy and gene-reaction associations of a metabolic model for Zymomonas mobilis ZM4, an industrially relevant ethanologenic organism with extremely high glycolytic flux and low biomass yield. Gene essentiality predictions made by the draft model were compared to data from individual pooled mutant experiments to identify areas of the model requiring deeper validation. Subsequent experiments showed that some of the discrepancies between the model and dataset were caused by polar effects, mis-mapped barcodes, or mutants carrying both wild-type and transposon disrupted gene copies-highlighting potential limitations inherent to data from individual mutants in these high-throughput datasets. Therefore, we analyzed correlations in fitness scores across all 492 experiments in the dataset in the context of functionally related metabolic reaction modules identified within the model via flux coupling analysis. These correlations were used to identify candidate genes for a reaction in histidine biosynthesis lacking an annotated gene and highlight metabolic modules with poorly correlated gene fitness scores. Additional genes for reactions involved in biotin, ubiquinone, and pyridoxine biosynthesis in Z. mobilis were identified and confirmed using mutant complementation experiments. These discovered genes, were incorporated into the final model, iZM4_478, which contains 747 metabolic and transport reactions (of which 612 have gene-protein-reaction associations), 478 genes, and 616 unique metabolites, making it one of the most complete models of Z. mobilis ZM4 to date. The methods of analysis that we applied here with the Z. mobilis transposon mutant dataset, could easily be utilized to improve future genome-scale metabolic reconstructions for organisms where these, or similar, high-throughput datasets are available.

Bar-seq strategies for the LeishGEdit toolbox

The number of fully sequenced genomes increases steadily but the function of many genes remains unstudied. To accelerate dissection of gene function in Leishmania spp. and other kinetoplastids we previously developed a streamlined pipeline for CRISPR-Cas9 gene editing, which we termed LeishGEdit. To facilitate high-throughput mutant screens we have adapted this pipeline by barcoding mutants with unique 17-nucleotide barcodes, allowing loss-of-function screens in mixed populations. Here we present primer design and analysis tools that facilitate these bar-seq strategies. We have developed a standalone easy-to-use pipeline to design CRISPR primers suitable for the LeishGEdit toolbox for any given genome and have generated a list of 14,995 barcodes. Barcodes and oligo sequences are now accessible through our website www.leishgedit.net allowing researchers to pursue bar-seq experiments in all currently available TriTrypDB genomes (release 41). This will streamline CRISPR bar-seq assays in kinetoplastids, enabling pooled mutant screens across the community.

A decade of advances in transposon-insertion sequencing

It has been 10 years since the introduction of modern transposon-insertion sequencing (TIS) methods, which combine genome-wide transposon mutagenesis with high-throughput sequencing to estimate the fitness contribution or essentiality of each genetic component in a bacterial genome. Four TIS variations were published in 2009: transposon sequencing (Tn-Seq), transposon-directed insertion site sequencing (TraDIS), insertion sequencing (INSeq) and high-throughput insertion tracking by deep sequencing (HITS). TIS has since become an important tool for molecular microbiologists, being one of the few genome-wide techniques that directly links phenotype to genotype and ultimately can assign gene function. In this Review, we discuss the recent applications of TIS to answer overarching biological questions. We explore emerging and multidisciplinary methods that build on TIS, with an eye towards future applications.

In this Review, several experts discuss progress in the decade since the development of transposon-based approaches for bacterial genetic screens. They describe how advances in both experimental technologies and analytical strategies are resulting in insights into diverse biological processes.

Essential roles for deubiquitination in Leishmania life cycle progression

The parasitic protozoan Leishmania requires proteasomal, autophagic and lysosomal proteolytic pathways to enact the extensive cellular remodelling that occurs during its life cycle. The proteasome is essential for parasite proliferation, yet little is known about the requirement for ubiquitination/deubiquitination processes in growth and differentiation. Activity-based protein profiling of L. mexicana C12, C19 and C65 deubiquitinating cysteine peptidases (DUBs) revealed DUB activity remains relatively constant during differentiation of procyclic promastigote to amastigote. However, when life cycle phenotyping (bar-seq) was performed on a pool including 15 barcoded DUB null mutants created in promastigotes using CRISPR-Cas9, significant loss of fitness was observed during differentiation and intracellular infection. DUBs 4, 7, and 13 are required for successful transformation from metacyclic promastigote to amastigote and DUBs 3, 5, 6, 8, 10, 11 and 14 are required for normal amastigote proliferation in mice. DUBs 1, 2, 12 and 16 are essential for promastigote viability and the essential role of DUB2 in establishing infection was demonstrated using DiCre inducible gene deletion in vitro and in vivo. DUB2 is found in the nucleus and interacts with nuclear proteins associated with transcription/chromatin dynamics, mRNA splicing and mRNA capping. DUB2 has broad linkage specificity, cleaving all the di-ubiquitin chains except for Lys27 and Met1. Our study demonstrates the crucial role that DUBs play in differentiation and intracellular survival of Leishmania and that amastigotes are exquisitely sensitive to disruption of ubiquitination homeostasis.

Yeast-Based Genetic Interaction Analysis of Human Kinome

Kinases are critical intracellular signaling proteins. To better understand kinase-mediated signal transduction, a large-scale human-yeast genetic interaction screen was performed. Among 597 human kinase genes tested, 28 displayed strong toxicity in yeast when overexpressed. En masse transformation of these toxic kinase genes into 4653 homozygous diploid yeast deletion mutants followed by barcode sequencing identified yeast toxicity modifiers and thus their human orthologs. Subsequent network analyses and functional grouping revealed that the 28 kinases and their 676 interaction partners (corresponding to a total of 969 genetic interactions) are enriched in cell death and survival (34%), small-molecule biochemistry (18%) and molecular transport (11%), among others. In the subnetwork analyses, a few kinases were commonly associated with glioma, cell migration and cell death/survival. Our analysis enabled the creation of a first draft of the kinase genetic interactome network and identified multiple drug targets for inflammatory diseases and cancer, in which deregulated kinase signaling plays a pathogenic role.

Genetic interaction profiles of regulatory kinases differ between environmental conditions and cellular states

Cell growth and quiescence in eukaryotic cells is controlled by an evolutionarily conserved network of signaling pathways. Signal transduction networks operate to modulate a wide range of cellular processes and physiological properties when cells exit proliferative growth and initiate a quiescent state. How signaling networks function to respond to diverse signals that result in cell cycle exit and establishment of a quiescent state is poorly understood. Here, we studied the function of signaling pathways in quiescent cells using global genetic interaction mapping in the model eukaryotic cell, Saccharomyces cerevisiae (budding yeast). We performed pooled analysis of genotypes using molecular barcode sequencing (Bar-seq) to test the role of ~4,000 gene deletion mutants and ~12,000 pairwise interactions between all non-essential genes and the protein kinase genes TOR1, RIM15, and PHO85 in three different nutrient-restricted conditions in both proliferative and quiescent cells. We detect up to 10-fold more genetic interactions in quiescent cells than proliferative cells. We find that both individual gene effects and genetic interaction profiles vary depending on the specific pro-quiescence signal. The master regulator of quiescence, RIM15, shows distinct genetic interaction profiles in response to different starvation signals. However, vacuole-related functions show consistent genetic interactions with RIM15 in response to different starvation signals, suggesting that RIM15 integrates diverse signals to maintain protein homeostasis in quiescent cells. Our study expands genome-wide genetic interaction profiling to additional conditions, and phenotypes, and highlights the conditional dependence of epistasis.

Competition experiments in a soil microcosm reveal the impact of genetic and biotic factors on natural yeast populations

The ability to measure microbial fitness directly in natural conditions and in interaction with other microbes is a challenge that needs to be overcome if we want to gain a better understanding of microbial fitness determinants in nature. Here we investigate the influence of the natural microbial community on the relative fitness of the North American populations SpB, SpC and SpC* of the wild yeast Saccharomyces paradoxus using DNA barcodes and a soil microcosm derived from soil associated with oak trees. We find that variation in fitness among these genetically distinct groups is influenced by the microbial community. Altering the microbial community load and diversity with an irradiation treatment significantly diminishes the magnitude of fitness differences among populations. Our findings suggest that microbial interactions could affect the evolution of yeast lineages in nature by modulating variation in fitness.