Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis

Synthetic lethality and the minimal genome size problem

Gene knockout studies suggest that ~300 genes in a bacterial genome and ~1,100 genes in a yeast genome cannot be deleted without loss of viability. These single-gene knockout experiments do not account for negative genetic interactions, when two or more genes can each be deleted without effect, but their joint deletion is lethal. Thus, large-scale single-gene deletion studies underestimate the size of a minimal gene set compatible with cell survival. In yeast Saccharomyces cerevisiae, the viability of all possible deletions of gene pairs (2-tuples), and of some deletions of gene triplets (3-tuples), has been experimentally tested. To estimate the size of a yeast minimal genome from that data, we first established that finding the size of a minimal gene set is equivalent to finding the minimum vertex cover in the lethality (hyper)graph, where the vertices are genes and (hyper)edges connect k-tuples of genes whose joint deletion is lethal. Using the Lovász-Johnson-Chvatal greedy approximation algorithm, we computed the minimum vertex cover of the synthetic-lethal 2-tuples graph to be 1,723 genes. We next simulated the genetic interactions in 3-tuples, extrapolating from the existing triplet sample, and again estimated minimum vertex covers. The size of a minimal gene set in yeast rapidly approaches the size of the entire genome even when considering only synthetic lethalities in k-tuples with small k. In contrast, several studies reported successful experimental reductions of yeast and bacterial genomes by simultaneous deletions of hundreds of genes, without eliciting synthetic lethality. We discuss possible reasons for this apparent contradiction.IMPORTANCEHow can we estimate the smallest number of genes sufficient for a unicellular organism to survive on a rich medium? One approach is to remove genes one at a time and count how many of such deletion strains are unable to grow. However, the single-gene knockout data are insufficient, because joint gene deletions may result in negative genetic interactions, also known as synthetic lethality. We used a technique from graph theory to estimate the size of minimal yeast genome from partial data on synthetic lethality. The number of potential synthetic lethal interactions grows very fast when multiple genes are deleted, revealing a paradoxical contrast with the experimental reductions of yeast genome by ~100 genes, and of bacterial genomes by several hundreds of genes.

macroH2A1 drives nucleosome dephasing and genome instability in histone humanized yeast

In addition to replicative histones, eukaryotic genomes encode a repertoire of non-replicative variant histones, providing additional layers of structural and epigenetic regulation. Here, we systematically replace individual replicative human histones with non-replicative human variant histones using a histone replacement system in yeast. We show that variants H2A.J, TsH2B, and H3.5 complement their respective replicative counterparts. However, macroH2A1 fails to complement, and its overexpression is toxic in yeast, negatively interacting with yeast's native histones and kinetochore genes. To isolate yeast with macroH2A1 chromatin, we uncouple the effects of its macro and histone fold domains, revealing that both domains suffice to override native nucleosome positioning. Furthermore, both uncoupled constructs of macroH2A1 exhibit lower nucleosome occupancy, decreased short-range chromatin interactions (<20 kb), disrupted centromeric clustering, and increased chromosome instability. Our observations demonstrate that lack of a canonical histone H2A dramatically alters chromatin organization in yeast, leading to genome instability and substantial fitness defects.

Histone deacetylase Hos2 regulates protein expression noise by potentially modulating the protein translation machinery

Non-genetic variations derived from expression noise at transcript or protein levels can result in cell-to-cell heterogeneity within an isogenic population. Although cells have developed strategies to reduce noise in some cellular functions, this heterogeneity can also facilitate varying levels of regulation and provide evolutionary benefits in specific environments. Despite several general characteristics of cellular noise having been revealed, the detailed molecular pathways underlying noise regulation remain elusive. Here, we established a dual-fluorescent reporter system in Saccharomyces cerevisiae and performed experimental evolution to search for mutations that increase expression noise. By analyzing evolved cells using bulk segregant analysis coupled with whole-genome sequencing, we identified the histone deacetylase Hos2 as a negative noise regulator. A hos2 mutant down-regulated multiple ribosomal protein genes and exhibited partially compromised protein translation, indicating that Hos2 may regulate protein expression noise by modulating the translation machinery. Treating cells with translation inhibitors or introducing mutations into several Hos2-regulated ribosomal protein genes-RPS9A, RPS28B and RPL42A-enhanced protein expression noise. Our study provides an effective strategy for identifying noise regulators and also sheds light on how cells regulate non-genetic variation through protein translation.

Evidence for the Role of the Mitochondrial ABC Transporter MDL1 in the Uptake of Clozapine and Related Molecules into the Yeast Saccharomyces cerevisiae

Clozapine is an antipsychotic drug whose accumulation in white cells can sometimes prove toxic; understanding the transporters and alleles responsible is thus highly desirable. We used a strategy in which a yeast (Saccharomyces cerevisiae) CRISPR-Cas9 knock-out library was exposed to cytotoxic concentrations of clozapine to determine those transporters whose absence made it more resistant; we also recognised the structural similarity of the fluorescent dye safranin O (also known as safranin T) to clozapine, allowing it to be used as a surrogate marker. Strains lacking the mitochondrial ABC transporter MDL1 (encoded by YLR188W) showed substantial resistance to clozapine. MDL1 overexpression also conferred extra sensitivity to clozapine and admitted a massive increase in the cellular and mitochondrial uptake of safranin O, as determined using flow cytometry and microscopically. Yeast lacking mitochondria showed no such unusual accumulation. Mitochondrial MDL1 is thus the main means of accumulation of clozapine in S. cerevisiae. The closest human homologue of S. cerevisiae MDL1 is ABCB10.

Natural variation in yeast reveals multiple paths for acquiring higher stress resistance

BACKGROUND

Organisms frequently experience environmental stresses that occur in predictable patterns and combinations. For wild Saccharomyces cerevisiae yeast growing in natural environments, cells may experience high osmotic stress when they first enter broken fruit, followed by high ethanol levels during fermentation, and then finally high levels of oxidative stress resulting from respiration of ethanol. Yeast have adapted to these patterns by evolving sophisticated "cross protection" mechanisms, where mild 'primary' doses of one stress can enhance tolerance to severe doses of a different 'secondary' stress. For example, in many yeast strains, mild osmotic or mild ethanol stresses cross protect against severe oxidative stress, which likely reflects an anticipatory response important for high fitness in nature.

RESULTS

During the course of genetic mapping studies aimed at understanding the mechanisms underlying natural variation in ethanol-induced cross protection against H2O2, we found that a key H2O2 scavenging enzyme, cytosolic catalase T (Ctt1p), was absolutely essential for cross protection in a wild oak strain. This suggested the absence of other compensatory mechanisms for acquiring H2O2 resistance in that strain background under those conditions. In this study, we found surprising heterogeneity across diverse yeast strains in whether CTT1 function was fully necessary for acquired H2O2 resistance. Some strains exhibited partial dispensability of CTT1 when ethanol and/or salt were used as mild stressors, suggesting that compensatory peroxidases may play a role in acquired stress resistance in certain genetic backgrounds. We leveraged global transcriptional responses to ethanol and salt stresses in strains with different levels of CTT1 dispensability, allowing us to identify possible regulators of these alternative peroxidases and acquired stress resistance in general.

CONCLUSIONS

Ultimately, this study highlights how superficially similar traits can have different underlying molecular foundations and provides a framework for understanding the diversity and regulation of stress defense mechanisms.

High-throughput classification of S. cerevisiae tetrads using deep learning

Meiotic crossovers play a vital role in proper chromosome segregation and evolution of most sexually reproducing organisms. Meiotic recombination can be visually observed in Saccharomyces cerevisiae tetrads using linked spore-autonomous fluorescent markers placed at defined intervals within the genome, which allows for analysis of meiotic segregation without the need for tetrad dissection. To automate the analysis, we developed a deep learning-based image recognition and classification pipeline for high-throughput tetrad detection and meiotic crossover classification. As a proof of concept, we analyzed a large image data set from wild-type and selected gene knock-out mutants to quantify crossover frequency, interference, chromosome missegregation, and gene conversion events. The deep learning-based method has the potential to accelerate the discovery of new genes involved in meiotic recombination in S. cerevisiae such as the underlying factors controlling crossover frequency and interference.

Ngk1 kinase-mediated N-acetylglucosamine metabolism promotes UDP-GlcNAc biosynthesis in Saccharomyces cerevisiae

N-acetylglucosamine (GlcNAc) is an important structural component of the cell wall chitin, N-glycans, glycolipids, and GPI-anchors in eukaryotes. GlcNAc kinase phosphorylates GlcNAc into GlcNAc-6-phosphate, a precursor of uridine diphosphate N-acetylglucosamine (UDP-GlcNAc) that serves as a substrate for glycan synthesis. Although GlcNAc kinase is found widely in organisms ranging from microorganisms to mammals, it has never been found in the model yeast Saccharomyces cerevisiae. Here, we demonstrate the presence of GlcNAc metabolism for UDP-GlcNAc biosynthesis in S. cerevisiae through Ngk1, a GlcNAc kinase we discovered previously. The overexpression or deletion of Ngk1 in the presence of GlcNAc affected the amount of both UDP-GlcNAc and chitin, suggesting that GlcNAc metabolism via Ngk1 promotes UDP-GlcNAc synthesis. Our data suggest that the Ngk1-mediated GlcNAc metabolism compensates for the hexosamine pathway, a known pathway for UDP-GlcNAc synthesis.

CRISPRi functional genomics in bacteria and its application to medical and industrial research

SUMMARYFunctional genomics is the use of systematic gene perturbation approaches to determine the contributions of genes under conditions of interest. Although functional genomic strategies have been used in bacteria for decades, recent studies have taken advantage of CRISPR (clustered regularly interspaced short palindromic repeats) technologies, such as CRISPRi (CRISPR interference), that are capable of precisely modulating expression of all genes in the genome. Here, we discuss and review the use of CRISPRi and related technologies for bacterial functional genomics. We discuss the strengths and weaknesses of CRISPRi as well as design considerations for CRISPRi genetic screens. We also review examples of how CRISPRi screens have defined relevant genetic targets for medical and industrial applications. Finally, we outline a few of the many possible directions that could be pursued using CRISPR-based functional genomics in bacteria. Our view is that the most exciting screens and discoveries are yet to come.

Search for protein kinase(s) related to cell growth or viability maintenance in the presence of ethanol in budding and fission yeasts

Alcohol fermentation comprises two phases: phase 1, alcohol fermentation occurs while yeast cells proliferate; phase 2, growth stops and alcohol fermentation continues. We categorized genes related to proliferation in low ethanol (phase 1) and viability in high ethanol (phase 2) as Alcohol Growth Ability (AGA) and Alcohol Viability (ALV), respectively. Although genes required for phase 1 are examined in budding yeast, those for phase 2 are unknown. We set conditions for ALV screening, searched for protein kinases (PKs) related to ALV in budding yeast, and expanded two screenings to fission yeast. Bub1 kinase was important for proliferation in low ethanol but not for viability in high ethanol, suggesting that the important PKs differ between the two phases. It was indeed the case. Further, 3 common PKs were identified as AGA in both yeasts, suggesting that the important cellular mechanism in phase 1 is conserved in both yeasts, at least partially.

A dynamic network model predicts the phenotypes of multicellular clusters from cellular properties

Cell division without cell separation produces multicellular clusters in budding yeast. Two fundamental characteristics of these clusters are their size (the number of cells per cluster) and cellular composition: the fractions of cells with different phenotypes. Using cells as nodes and links between mother and daughter cells as edges, we model cluster growth and breakage by varying three parameters: the cell division rate, the rate at which intercellular connections break, and the kissing number (the maximum number of connections to one cell). We find that the kissing number sets the maximum possible cluster size. Below this limit, the ratio of the cell division rate to the connection breaking rate determines the cluster size. If links have a constant probability of breaking per unit time, the probability that a link survives decreases exponentially with its age. Modeling this behavior recapitulates experimental data. We then use this framework to examine synthetic, differentiating clusters with two cell types, faster-growing germ cells and their somatic derivatives. The fraction of clusters that contain both cell types increases as either of two parameters increase: the kissing number and difference between the growth rate of germ and somatic cells. In a population of clusters, the variation in cellular composition is inversely correlated (r2 = 0.87) with the average fraction of somatic cells in clusters. Our results show how a small number of cellular features can control the phenotypes of multicellular clusters that were potentially the ancestors of more complex forms of multicellular development, organization, and reproduction.

Proline-directed yeast and human MAP kinases phosphorylate the Dot1p/DOT1L histone H3K79 methyltransferase

Disruptor of telomeric silencing 1 (Dot1p) is an exquisitely conserved histone methyltransferase and is the sole enzyme responsible for H3K79 methylation in the budding yeast Saccharomyces cerevisiae. It has been shown to be highly phosphorylated in vivo; however, the upstream kinases that act on Dot1p are almost entirely unknown in yeast and all other eukaryotes. Here, we used in vitro and in vivo kinase discovery approaches to show that mitogen-activated protein kinase HOG1 (Hog1p) is a bona fide kinase of the Dot1p methyltransferase. In vitro kinase assays showed that Hog1p phosphorylates Dot1p at multiple sites, including at several proline-adjacent sites that are consistent with known Hog1p substrate preferences. The activity of Hog1p was specifically enhanced at these proline-adjacent sites on Dot1p upon Hog1p activation by the osmostress-responsive MAP kinase kinase PBS2 (Pbs2p). Genomic deletion of HOG1 reduced phosphorylation at specific sites on Dot1p in vivo, providing further evidence for Hog1p kinase activity on Dot1p in budding yeast cells. Phenotypic analysis of knockout and phosphosite mutant yeast strains revealed the importance of Hog1p-catalysed phosphorylation of Dot1p for cellular responses to ultraviolet-induced DNA damage. In mammalian systems, this kinase-substrate relationship was found to be conserved: human DOT1L (the ortholog of yeast Dot1p) can be phosphorylated by the proline-directed kinase p38β (also known as MAPK11; the ortholog of yeast Hog1p) at multiple sites in vitro. Taken together, our findings establish Hog1p and p38β as newly identified upstream kinases of the Dot1p/DOT1L H3K79 methyltransferase enzymes in eukaryotes.

Elucidation and engineering mitochondrial respiratory-related genes for improving bioethanol production at high temperature in Saccharomyces cerevisiae

Industrial manufacturing of bioproducts, especially bioethanol, can benefit from high-temperature fermentation, which requires the use of thermotolerant yeast strains. Mitochondrial activity in yeast is closely related to its overall metabolism. However, the mitochondrial respiratory changes in response to adaptive thermotolerance are still poorly understood and have been rarely utilized for developing thermotolerant yeast cell factories. Here, adaptive evolution and transcriptional sequencing, as well as whole-genome-level gene knockout, were used to obtain a thermotolerant strain of Saccharomyces cerevisiae. Furthermore, thermotolerance and bioethanol production efficiency of the engineered strain were examined. Physiological evaluation showed the boosted fermentation capacity and suppressed mitochondrial respiratory activity in the thermotolerant strain. The improved fermentation produced an increased supply of adenosine triphosphate required for more active energy-consuming pathways. Transcriptome analysis revealed significant changes in the expression of the genes involved in the mitochondrial respiratory chain. Evaluation of mitochondria-associated gene knockout confirmed that ADK1, DOC1, or MET7 were the key factors for the adaptive evolution of thermotolerance in the engineered yeast strain. Intriguingly, overexpression of DOC1 with TEF1 promoter regulation led to a 10.1% increase in ethanol production at 42 °C. The relationships between thermotolerance, mitochondrial activity, and respiration were explored, and a thermotolerant yeast strain was developed by altering the expression of mitochondrial respiration-related genes. This study provides a better understanding on the physiological mechanism of adaptive evolution of thermotolerance in yeast.

N2-Acetylornithine deacetylase functions as a Cys-Gly dipeptidase in the cytosolic glutathione degradation pathway in Arabidopsis thaliana

Glutathione (GSH) is required for various physiological processes in plants, including redox regulation and detoxification of harmful compounds. GSH also functions as a repository for assimilated sulfur and is actively catabolized in plants. In Arabidopsis, GSH is mainly degraded initially by cytosolic enzymes, γ-glutamyl cyclotransferase, and γ-glutamyl peptidase, which release cysteinylglycine (Cys-Gly). However, the subsequent enzyme responsible for catabolizing this dipeptide has not been identified to date. In the present study, we identified At4g17830 as a Cys-Gly dipeptidase, namely cysteinylglycine peptidase 1 (CGP1). CGP1 complemented the phenotype of the yeast mutant that cannot degrade Cys-Gly. The Arabidopsis cgp1 mutant had lower Cys-Gly degradation activity than the wild type and showed perturbed concentrations of thiol compounds. Recombinant CGP1 showed reasonable Cys-Gly degradation activity in vitro. Metabolomic analysis revealed that cgp1 exhibited signs of severe sulfur deficiency, such as elevated accumulation of O-acetylserine (OAS) and the decrease in sulfur-containing metabolites. Morphological changes observed in cgp1, including longer primary roots of germinating seeds, were also likely associated with sulfur starvation. Notably, At4g17830 has previously been reported to encode an N2-acetylornithine deacetylase (NAOD) that functions in the ornithine biosynthesis. The cgp1 mutant did not show a decrease in ornithine content, whereas the analysis of CGP1 structure did not rule out the possibility that CGP1 has Cys-Gly dipeptidase and NAOD activities. Therefore, we propose that CGP1 is a Cys-Gly dipeptidase that functions in the cytosolic GSH degradation pathway and may play dual roles in GSH and ornithine metabolism.

Cytoduction preserves genetic diversity following plasmid transfer into pooled yeast libraries

Much of our understanding of functional genomics derives from insights gained from large strain libraries including the yeast deletion collection, the GFP and TAP-tagged libraries, QTL mapping populations, among others [1-5]. A limitation of these libraries is that it is not easy to introduce reporters or make genetic perturbations to all strains in these collections. Tools such as Synthetic Genetic Arrays allow for the genetic manipulation of these libraries but are labor intensive and require specialized equipment for high throughput pinning [6]. Manipulating a diverse library en mass without losing diversity remains challenging. Ultimately, this limitation stems from the inefficiency of transformation, which is the standard method for genetic manipulation in yeast. Here, we develop a method that uses cytoduction (mating without nuclear fusion) to transfer plasmids directionally from a "Donor" to a diverse pool of "Recipient" strains. Because cytoduction uses mating, it is a natural process and is orders-of-magnitude more efficient than transformation, enabling the introduction of plasmids into high-diversity libraries with minimal impact on the diversity of the population.

Screening and In Silico Analyses of the Yeast Saccharomyces cerevisiae Σ1278b Bank Mutants Using Citral as a Natural Antimicrobial

The antimicrobial function of citral, one of the main compounds of the essential oils (EO) of the Citrus genus, and widely used by the food industry toward spoilage yeast, was previously proven. In this study, the possible mode of action of citral against yeast cells was evaluated by using a global deletome approach. Firstly, the suitability of Saccharomyces cerevisiae Σ1278b to serve as model yeast was assessed by determining its sensitivity to citral (MIC = 0.5 μL/mL). Subsequently, the complete library of Σ1278b haploid mutants deleted in 4019 non-essential genes was screened to identify potential molecular targets of citral. Finally, the deleted genes in the 590 mutants showing increased citral resistance was analyzed with an in-silico approach (Gene Ontology). The significantly enriched GO Terms were "cytoplasm", "vacuole", and "mitochondrion" (cellular components); "catalytic activity" (molecular function); "pseudohyphal growth" (biological process). For molecular function, resistant mutants were grouped into thiosulfate sulfur transferase activity, transferase activity, and oxidoreductase activity; for cellular components, resistant mutants were grouped as: cytoplasm, intracellular organelle, membrane-bounded organelle, mitochondrion, organelle membrane, and vacuole; and finally, with regard to biological process, deleted genes were grouped as: pseudohyphal growth, mitochondrion organization, lipid metabolic process, DNA recombination and repair, and proteolysis. Interestingly, many identified genes were associated with the cellular response to oxidative stress and ROS scavenging. These findings have important implications for the development of citral-based antimicrobials and the elucidation of its mechanism of action.

A U2 snRNP-specific protein, U2A', is involved in stress response and drug resistance in Cryptococcus deneoformans

The spliceosome, a large complex containing five conserved small ribonucleoprotein particles (snRNPs) U1, U2, U4, U5 and U6, plays important roles in precursor messenger RNA splicing. However, the function and mechanism of the spliceosomal snRNPs have not been thoroughly studied in the pathogenic yeast Cryptococcus deneoformans. In this study, we identified a U2A' homologous protein as a component of the cryptococcal U2 snRNP, which was encoded by the LEA1 gene. Using the "suicide" CRISPR-Cas9 tool, we deleted the LEA1 gene in C. deneoformans JEC21 strain and obtained the disruption mutant lea1Δ. The mutant showed a hypersensitivity to 0.03 % sodium dodecyl sulfate, as well as disordered chitin distribution in cell wall observed with Calcofluor White staining, which collectively illustrated the function of U2A' in maintenance of cell wall integrity. Further examination showed that lea1Δ displayed a decreased tolerance to lower or elevated temperatures, osmotic pressure and oxidative stress. The lea1Δ still exhibited susceptibility to geneticin and 5-flucytosine, and increased resistance to ketoconazole. Even, the mutant had a reduced capsule, and the virulence of lea1Δ in the Galleria mellonella model was decreased. Our results indicate that the U2A'-mediated RNA-processing has a particular role in the processing of gene products involved in response to stresses and virulence.

Development of a telomere vector-based approach to overcome limitations caused by lethal phenotypes in the study of essential genes in Magnaporthe oryzae

Reverse genetic approaches are common tools in genomics for elucidating gene functions, involving techniques such as gene deletion followed by screening for aberrant phenotypes. If the generation of gene deletion mutants fails, the question arises whether the failure stems from technical issues or because the gene of interest (GOI) is essential, meaning that the deletion causes lethality. In this report, we introduce a novel method for assessing gene essentiality using the phytopathogenic ascomycete Magnaporthe oryzae. The method is based on the observation that telomere vectors are lost in transformants during cultivation without selection pressure. We tested the hypothesis that essential genes can be identified in deletion mutants co-transformed with a telomere vector. The M. oryzae gene MoPKC, described in literature as essential, was chosen as GOI. Using CRISPR/Cas9 technology transformants with deleted GOI were generated and backed up by a telomere vector carrying a copy of the GOI and conferring fenhexamid resistance. Transformants in which the GOI deletion in the genome was not successful lost the telomere vector on media without fenhexamid. In contrast, transformants with confirmed GOI deletion retained the telomere vector even in absence of fenhexamid selection. In the latter case, the maintenance of the telomere indicates that the GOI is essential for the surveillance of the fungi, as it would have been lost otherwise. The method presented here allows to test for essentiality of genes when no mutants can be obtained from gene deletion approaches, thereby expanding the toolbox for studying gene function in ascomycetes.

Yeast poly(A)-binding protein (Pab1) controls translation initiation in vivo primarily by blocking mRNA decapping and decay

Poly(A)-binding protein (Pab1 in yeast) is involved in mRNA decay and translation initiation, but its molecular functions are incompletely understood. We found that auxin-induced degradation of Pab1 reduced bulk mRNA and polysome abundance in a manner suppressed by deleting the catalytic subunit of decapping enzyme (dcp2Δ), demonstrating that enhanced decapping/degradation is the major driver of reduced mRNA abundance and protein synthesis at limiting Pab1 levels. An increased median poly(A) tail length conferred by Pab1 depletion was also nullified by dcp2Δ, suggesting that mRNA isoforms with shorter tails are preferentially decapped/degraded at limiting Pab1. In contrast to findings on mammalian cells, the translational efficiencies (TEs) of many mRNAs were altered by Pab1 depletion; however, these changes were broadly diminished by dcp2∆, suggesting that reduced mRNA abundance is a major driver of translational reprogramming at limiting Pab1. Thus, assembly of the closed-loop mRNP via PABP-eIF4G interaction appears to be dispensable for normal translation of most yeast mRNAs in vivo. Interestingly, histone mRNAs and proteins are preferentially diminished on Pab1 depletion dependent on Dcp2, accompanied by activation of internal cryptic promoters in the manner expected for reduced nucleosome occupancies, revealing a new layer of post-transcriptional control of histone gene expression.

A protein network refinement method based on module discovery and biological information

BACKGROUND

The identification of essential proteins can help in understanding the minimum requirements for cell survival and development to discover drug targets and prevent disease. Nowadays, node ranking methods are a common way to identify essential proteins, but the poor data quality of the underlying PIN has somewhat hindered the identification accuracy of essential proteins for these methods in the PIN. Therefore, researchers constructed refinement networks by considering certain biological properties of interacting protein pairs to improve the performance of node ranking methods in the PIN. Studies show that proteins in a complex are more likely to be essential than proteins not present in the complex. However, the modularity is usually ignored for the refinement methods of the PINs.

METHODS

Based on this, we proposed a network refinement method based on module discovery and biological information. The idea is, first, to extract the maximal connected subgraph in the PIN, and to divide it into different modules by using Fast-unfolding algorithm; then, to detect critical modules according to the orthologous information, subcellular localization information and topology information within each module; finally, to construct a more refined network (CM-PIN) by using the identified critical modules.

RESULTS

To evaluate the effectiveness of the proposed method, we used 12 typical node ranking methods (LAC, DC, DMNC, NC, TP, LID, CC, BC, PR, LR, PeC, WDC) to compare the overall performance of the CM-PIN with those on the S-PIN, D-PIN and RD-PIN. The experimental results showed that the CM-PIN was optimal in terms of the identification number of essential proteins, precision-recall curve, Jackknifing method and other criteria, and can help to identify essential proteins more accurately.

Transcriptome-wide analysis of the function of Ded1 in translation preinitiation complex assembly in a reconstituted in vitro system

We have developed a deep sequencing-based approach, Rec-Seq, that allows simultaneous monitoring of ribosomal 48S preinitiation complex (PIC) formation on every mRNA in the translatome in an in vitro reconstituted system. Rec-Seq isolates key early steps in translation initiation in the absence of all other cellular components and processes. Using this approach, we show that the DEAD-box ATPase Ded1 promotes 48S PIC formation on the start codons of >1000 native mRNAs, most of which have long, structured 5'-untranslated regions (5'UTRs). Remarkably, initiation measured in Rec-Seq was enhanced by Ded1 for most mRNAs previously shown to be highly Ded1-dependent by ribosome profiling of ded1 mutants in vivo, demonstrating that the core translation functions of the factor are recapitulated in the purified system. Our data do not support a model in which Ded1acts by reducing initiation at alternative start codons in 5'UTRs and instead indicate it functions by directly promoting mRNA recruitment to the 43S PIC and scanning to locate the main start codon. We also provide evidence that eIF4A, another essential DEAD-box initiation factor, is required for efficient PIC assembly on almost all mRNAs, regardless of their structural complexity, in contrast to the preferential stimulation by Ded1 of initiation on mRNAs with long, structured 5'UTRs.

Restricted glycolysis is a primary cause of the reduced growth rate of zinc-deficient yeast cells

Zinc is required for many critical processes, including intermediary metabolism. In Saccharomyces cerevisiae, the Zap1 activator regulates the transcription of ∼80 genes in response to Zn supply. Some Zap1-regulated genes are Zn transporters that maintain Zn homeostasis, while others mediate adaptive responses that enhance fitness. One adaptive response gene encodes the 2-cysteine peroxiredoxin Tsa1, which is critical to Zn-deficient (ZnD) growth. Depending on its redox state, Tsa1 can function as a peroxidase, a protein chaperone, or a regulatory redox sensor. In a screen for possible Tsa1 regulatory targets, we identified a mutation (cdc19S492A) that partially suppressed the tsa1Δ growth defect. The cdc19S492A mutation reduced activity of its protein product, pyruvate kinase isozyme 1 (Pyk1), implicating Tsa1 in adapting glycolysis to ZnD conditions. Glycolysis requires activity of the Zn-dependent enzyme fructose-bisphosphate aldolase 1, which was substantially decreased in ZnD cells. We hypothesized that in ZnD tsa1Δ cells, the loss of a compensatory Tsa1 regulatory function causes depletion of glycolytic intermediates and restricts dependent amino acid synthesis pathways, and that the decreased activity of Pyk1S492A counteracted this depletion by slowing the irreversible conversion of phosphoenolpyruvate to pyruvate. In support of this model, supplementing ZnD tsa1Δ cells with aromatic amino acids improved their growth. Phosphoenolpyruvate supplementation, in contrast, had a much greater effect on growth rate of WT and tsa1Δ ZnD cells, indicating that inefficient glycolysis is a major factor limiting yeast growth. Surprisingly however, this restriction was not primarily due to low fructose-bisphosphate aldolase 1 activity, but instead occurs earlier in glycolysis.

Cytoplasmic redox imbalance in the thioredoxin system activates Hsf1 and results in hyperaccumulation of the sequestrase Hsp42 with misfolded proteins

Cells employ multiple systems to maintain homeostasis when experiencing environmental stress. For example, the folding of nascent polypeptides is exquisitely sensitive to proteotoxic stressors including heat, pH, and oxidative stress, and is safeguarded by a network of protein chaperones that concentrate potentially toxic misfolded proteins into transient assemblies to promote folding or degradation. The redox environment itself is buffered by both cytosolic and organellar thioredoxin and glutathione pathways. How these systems are linked is poorly understood. Here, we determine that specific disruption of the cytosolic thioredoxin system resulted in constitutive activation of the heat shock response in Saccharomyces cerevisiae and accumulation of the sequestrase Hsp42 into an exaggerated and persistent juxtanuclear quality control (JUNQ) compartment. Terminally misfolded proteins also accumulated in this compartment in thioredoxin reductase (TRR1)-deficient cells, despite apparently normal formation and dissolution of transient cytoplasmic quality control (CytoQ) bodies during heat shock. Notably, cells lacking TRR1 and HSP42 exhibited severe synthetic slow growth exacerbated by oxidative stress, signifying a critical role for Hsp42 under redox-challenged conditions. Finally, we demonstrated that Hsp42 localization patterns in trr1∆ cells mimic those observed in chronically aging and glucose-starved cells, linking nutrient depletion and redox imbalance with management of misfolded proteins via a process of long-term sequestration.

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.

Yeast zinc cluster transcription factors involved in the switch from fermentation to respiration show interdependency for DNA binding revealing a novel type of DNA recognition

In budding yeast, fermentation is the most important pathway for energy production. Under low-glucose conditions, ethanol is used for synthesis of this sugar requiring a shift to respiration. This process is controlled by the transcriptional regulators Cat8, Sip4, Rds2 and Ert1. We characterized Gsm1 (glucose starvation modulator 1), a paralog of Rds2 and Ert1. Genome-wide analysis showed that Gsm1 has a DNA binding profile highly similar to Rds2. Binding of Gsm1 and Rds2 is interdependent at the gluconeogenic gene FBP1. However, Rds2 is required for Gsm1 to bind at other promoters but not the reverse. Gsm1 and Rds2 also bind to DNA independently of each other. Western blot analysis revealed that Rds2 controls expression of Gsm1. In addition, we showed that the DNA binding domains of Gsm1 and Rds2 bind cooperatively in vitro to the FBP1 promoter. In contrast, at the HAP4 gene, Ert1 cooperates with Rds2 for DNA binding. Mutational analysis suggests that Gsm1/Rds2 and Ert1/Rds2 bind to short common DNA stretches, revealing a novel mode of binding for this class of factors. Two-point mutations in a HAP4 site convert it to a Gsm1 binding site. Thus, Rds2 controls binding of Gsm1 at many promoters by two different mechanisms: regulation of Gsm1 levels and increased DNA binding by formation of heterodimers.

In vivo CRISPR screening directly targeting testicular cells

CRISPR-Cas9 short guide RNA (sgRNA) library screening is a powerful approach to understand the molecular mechanisms of biological phenomena. However, its in vivo application is currently limited. Here, we developed our previously established in vitro revival screening method into an in vivo one to identify factors involved in spermatogenesis integrity by utilizing sperm capacitation as an indicator. By introducing an sgRNA library into testicular cells, we successfully pinpointed the retinal degeneration 3 (Rd3) gene as a significant factor in spermatogenesis. Single-cell RNA sequencing (scRNA-seq) analysis highlighted the high expression of Rd3 in round spermatids, and proteomics analysis indicated that Rd3 interacts with mitochondria. To search for cell-type-specific signaling pathways based on scRNA-seq and proteomics analyses, we developed a computational tool, Hub-Explorer. Through this, we discovered that Rd3 modulates oxidative stress by regulating mitochondrial distribution upon ciliogenesis induction. Collectively, our screening system provides a valuable in vivo approach to decipher molecular mechanisms in biological processes.

Plasma membrane damage limits replicative lifespan in yeast and induces premature senescence in human fibroblasts

Plasma membrane damage (PMD) occurs in all cell types due to environmental perturbation and cell-autonomous activities. However, cellular outcomes of PMD remain largely unknown except for recovery or death. In this study, using budding yeast and normal human fibroblasts, we found that cellular senescence-stable cell cycle arrest contributing to organismal aging-is the long-term outcome of PMD. Our genetic screening using budding yeast unexpectedly identified a close genetic association between PMD response and replicative lifespan regulations. Furthermore, PMD limits replicative lifespan in budding yeast; upregulation of membrane repair factors ESCRT-III (SNF7) and AAA-ATPase (VPS4) extends it. In normal human fibroblasts, PMD induces premature senescence via the Ca2+-p53 axis but not the major senescence pathway, DNA damage response pathway. Transient upregulation of ESCRT-III (CHMP4B) suppressed PMD-dependent senescence. Together with mRNA sequencing results, our study highlights an underappreciated but ubiquitous senescent cell subtype: PMD-dependent senescent cells.

BnaA4.BOR2 contributes the tolerance of rapeseed to boron deficiency by improving the transport of boron from root to shoot

Boron (B) is essential for plant growth. However, the molecular mechanism of B transport in rapeseed (Brassica napus L.) is unknown well. Here, we report that B transporter BnaA4.BOR2 is involved in the transport of B from root to shoot and its distribution in shoot cell wall and flower in rapeseed. The results of GUS staining and in-situ PCR analysis showed that BnaA4.BOR2 is mainly expressed in cortex and endodermis of root tip meristem zone and endodermis of mature zone. BnaA4.BOR2 was mainly localized in plasma membrane and showed B transport activity in yeast. Overexpression of Bna4.BOR2 could rescue the phenotype of Arabidopsis mutant bor2-2 under low-B condition. Furthermore, knockout of BnaA4.BOR2 could significantly enhance the sensitivity of rapeseed mutants to B deficiency, including inhibition of root elongation and biomass decrease of roots and shoots. The B concentration in xylem sap of BnaA4.BOR2 mutants was significantly decreased under B deficiency, which resulted in significantly lower B concentrations in shoot cell wall at seedling stage and flower organ at reproductive stage compared to that of wild-type QY10. The growth of BnaA4.BOR2 mutants were severely inhibited, exhibiting a typical B-deficient phenotype of "flowering without seed setting", leading to a sharp decrease in seed yield in B deficient soil. Taken together, these results indicate that BnaA4.BOR2 is critical for rapeseed growth and seed yield production under low B level, which is mainly expressed in cortex and endodermis, and contributed to the transport of B from roots to shoots and its distribution in shoot.

Amino acid metabolism and MAP kinase signaling pathway play opposite roles in the regulation of ethanol production during fermentation of sugarcane molasses in budding yeast

Sugarcane molasses is one of the main raw materials for bioethanol production, and Saccharomyces cerevisiae is the major biofuel-producing organism. In this study, a batch fermentation model has been used to examine ethanol titers of deletion mutants for all yeast nonessential genes in this yeast genome. A total of 42 genes are identified to be involved in ethanol production during fermentation of sugarcane molasses. Deletion mutants of seventeen genes show increased ethanol titers, while deletion mutants for twenty-five genes exhibit reduced ethanol titers. Two MAP kinases Hog1 and Kss1 controlling the high osmolarity and glycerol (HOG) signaling and the filamentous growth, respectively, are negatively involved in the regulation of ethanol production. In addition, twelve genes involved in amino acid metabolism are crucial for ethanol production during fermentation. Our findings provide novel targets and strategies for genetically engineering industrial yeast strains to improve ethanol titer during fermentation of sugarcane molasses.

A CUG codon-adapted anchor-away toolkit for functional analysis of genes in Candida albicans

Promoter shutoff of essential genes in the diploid Candida albicans has often been insufficient to create tight, conditional null alleles due to leaky expression and has been a stumbling block in pathogenesis research. Moreover, homozygous deletion of non-essential genes has often been problematic due to the frequent aneuploidy in the mutant strains. Rapid, conditional depletion of essential genes by the anchor-away strategy has been successfully employed in Saccharomyces cerevisiae and other model organisms. Here, rapamycin mediates the dimerization of human FK506-binding protein (FKBP12) and FKBP12-rapamycin-binding (FRB) domain-containing target protein, resulting in relocalization to altered sub-cellular locations. In this work, we used the ribosomal protein Rpl13 as the anchor and took two nuclear proteins as targets to construct a set of mutants in a proof-of-principle approach. We first constructed a rapamycin-resistant C. albicans strain by introducing a dominant mutation in the CaTOR1 gene and a homozygous deletion of RBP1, the ortholog of FKBP12, a primary target of rapamycin. The FKBP12 and the FRB coding sequences were then CUG codon-adapted for C. albicans by site-directed mutagenesis. Anchor-away strains expressing the essential TBP1 gene or the non-essential SPT8 gene as FRB fusions were constructed. We found that rapamycin caused rapid cessation of growth of the TBP-AA strain within 15 minutes and the SPT8-AA strain phenocopied the constitutive filamentous phenotype of the spt8Δ/spt8Δ mutant. Thus, the anchor-away toolbox for C. albicans developed here can be employed for genome-wide analysis to identify gene function in a rapid and reliable manner, further accelerating anti-fungal drug development in C. albicans.

IMPORTANCE

Molecular genetic studies thus far have identified ~27% open-reading frames as being essential for the vegetative growth of Candida albicans in rich medium out of a total 6,198 haploid set of open reading frames. However, a major limitation has been to construct rapid conditional alleles of essential C. albicans genes with near quantitative depletion of encoded proteins. Here, we have developed a toolbox for rapid and conditional depletion of genes that would aid studies of gene function of both essential and non-essential genes.

The Vps13-like protein BLTP2 is pro-survival and regulates phosphatidylethanolamine levels in the plasma membrane to maintain its fluidity and function

Lipid transport proteins (LTPs) facilitate nonvesicular lipid exchange between cellular compartments and have critical roles in lipid homeostasis1. A new family of bridge-like LTPs (BLTPs) is thought to form lipid-transporting conduits between organelles2. One, BLTP2, is conserved across species but its function is not known. Here, we show that BLTP2 and its homolog directly regulate plasma membrane (PM) fluidity by increasing the phosphatidylethanolamine (PE) level in the PM. BLTP2 localizes to endoplasmic reticulum (ER)-PM contact sites34, 5, suggesting it transports PE from the ER to the PM. We find BLTP2 works in parallel with another pathway that regulates intracellular PE distribution and PM fluidity6, 7. BLTP2 expression correlates with breast cancer aggressiveness8-10. We found BLTP2 facilitates growth of a human cancer cell line and sustains its aggressiveness in an in vivo model of metastasis, suggesting maintenance of PM fluidity by BLTP2 may be critical for tumorigenesis in humans.

Transcriptome-wide analysis of the function of Ded1 in translation preinitiation complex assembly in a reconstituted in vitro system

We have developed a deep sequencing-based approach, Rec-Seq, that allows simultaneous monitoring of ribosomal 48S pre-initiation complex (PIC) formation on every mRNA in the translatome in an in vitro reconstituted system. Rec-Seq isolates key early steps in translation initiation in the absence of all other cellular components and processes. Using this approach we show that the DEAD-box ATPase Ded1 promotes 48S PIC formation on the start codons of >1000 native mRNAs, most of which have long, structured 5'-untranslated regions (5'UTRs). Remarkably, initiation measured in Rec-Seq was enhanced by Ded1 for most mRNAs previously shown to be highly Ded1-dependent by ribosome profiling of ded1 mutants in vivo, demonstrating that the core translation functions of the factor are recapitulated in the purified system. Our data do not support a model in which Ded1acts by reducing initiation at alternative start codons in 5'UTRs and instead indicate it functions by directly promoting mRNA recruitment to the 43S PIC and scanning to locate the main start codon. We also provide evidence that eIF4A, another essential DEAD-box initiation factor, is required for efficient PIC assembly on almost all mRNAs, regardless of their structural complexity, in contrast to the preferential stimulation by Ded1 of initiation on mRNAs with long, structured 5'UTRs.

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.

ECDEP: identifying essential proteins based on evolutionary community discovery and subcellular localization

BACKGROUND

In cellular activities, essential proteins play a vital role and are instrumental in comprehending fundamental biological necessities and identifying pathogenic genes. Current deep learning approaches for predicting essential proteins underutilize the potential of gene expression data and are inadequate for the exploration of dynamic networks with limited evaluation across diverse species.

RESULTS

We introduce ECDEP, an essential protein identification model based on evolutionary community discovery. ECDEP integrates temporal gene expression data with a protein-protein interaction (PPI) network and employs the 3-Sigma rule to eliminate outliers at each time point, constructing a dynamic network. Next, we utilize edge birth and death information to establish an interaction streaming source to feed into the evolutionary community discovery algorithm and then identify overlapping communities during the evolution of the dynamic network. SVM recursive feature elimination (RFE) is applied to extract the most informative communities, which are combined with subcellular localization data for classification predictions. We assess the performance of ECDEP by comparing it against ten centrality methods, four shallow machine learning methods with RFE, and two deep learning methods that incorporate multiple biological data sources on Saccharomyces. Cerevisiae (S. cerevisiae), Homo sapiens (H. sapiens), Mus musculus, and Caenorhabditis elegans. ECDEP achieves an AP value of 0.86 on the H. sapiens dataset and the contribution ratio of community features in classification reaches 0.54 on the S. cerevisiae (Krogan) dataset.

CONCLUSIONS

Our proposed method adeptly integrates network dynamics and yields outstanding results across various datasets. Furthermore, the incorporation of evolutionary community discovery algorithms amplifies the capacity of gene expression data in classification.

Large-scale genomic rearrangements boost SCRaMbLE in Saccharomyces cerevisiae

Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is a promising tool to study genomic rearrangements. However, the potential of SCRaMbLE to study genomic rearrangements is currently hindered, because a strain containing all 16 synthetic chromosomes is not yet available. Here, we construct SparLox83R, a yeast strain containing 83 loxPsym sites distributed across all 16 chromosomes. SCRaMbLE of SparLox83R produces versatile genome-wide genomic rearrangements, including inter-chromosomal events. Moreover, when combined with synthetic chromosomes, SCRaMbLE of hetero-diploids with SparLox83R leads to increased diversity of genomic rearrangements and relatively faster evolution of traits compared to hetero-diploids only with wild-type chromosomes. Analysis of the SCRaMbLEd strain with increased tolerance to nocodazole demonstrates that genomic rearrangements can perturb the transcriptome and 3D genome structure and consequently impact phenotypes. In summary, a genome with sparsely distributed loxPsym sites can serve as a powerful tool for studying the consequence of genomic rearrangements and accelerating strain engineering in Saccharomyces cerevisiae.

Phosphate starvation signaling increases mitochondrial membrane potential through respiration-independent mechanisms

Mitochondrial membrane potential directly powers many critical functions of mitochondria, including ATP production, mitochondrial protein import, and metabolite transport. Its loss is a cardinal feature of aging and mitochondrial diseases, and cells closely monitor membrane potential as an indicator of mitochondrial health. Given its central importance, it is logical that cells would modulate mitochondrial membrane potential in response to demand and environmental cues, but there has been little exploration of this question. We report that loss of the Sit4 protein phosphatase in yeast increases mitochondrial membrane potential, both by inducing the electron transport chain and the phosphate starvation response. Indeed, a similarly elevated mitochondrial membrane potential is also elicited simply by phosphate starvation or by abrogation of the Pho85-dependent phosphate sensing pathway. This enhanced membrane potential is primarily driven by an unexpected activity of the ADP/ATP carrier. We also demonstrate that this connection between phosphate limitation and enhancement of mitochondrial membrane potential is observed in primary and immortalized mammalian cells as well as in Drosophila. These data suggest that mitochondrial membrane potential is subject to environmental stimuli and intracellular signaling regulation and raise the possibility for therapeutic enhancement of mitochondrial function even in defective mitochondria.

Riboproteome remodeling during quiescence exit in Saccharomyces cerevisiae

The quiescent state is the prevalent mode of cellular life in most cells. Saccharomyces cerevisiae is a useful model for studying the molecular basis of the cell cycle, quiescence, and aging. Previous studies indicate that heterogeneous ribosomes show a specialized translation function to adjust the cellular proteome upon a specific stimulus. Using nano LC-MS/MS, we identified 69 of the 79 ribosomal proteins (RPs) that constitute the eukaryotic 80S ribosome during quiescence. Our study shows that the riboproteome is composed of 444 accessory proteins comprising cellular functions such as translation, protein folding, amino acid and glucose metabolism, cellular responses to oxidative stress, and protein degradation. Furthermore, the stoichiometry of both RPs and accessory proteins on ribosome particles is different depending on growth conditions and among monosome and polysome fractions. Deficiency of different RPs resulted in defects of translational capacity, suggesting that ribosome composition can result in changes in translational activity during quiescence.

Cross-feeding promotes heterogeneity within yeast cell populations

Cellular heterogeneity in cell populations of isogenic origin is driven by intrinsic factors such as stochastic gene expression, as well as external factors like nutrient availability and interactions with neighbouring cells. Heterogeneity promotes population fitness and thus has important implications in antimicrobial and anticancer treatments, where stress tolerance plays a significant role. Here, we study plasmid retention dynamics within a population of plasmid-complemented ura3∆0 yeast cells, and show that the exchange of complementary metabolites between plasmid-carrying prototrophs and plasmid-free auxotrophs allows the latter to survive and proliferate in selective environments. This process also affects plasmid copy number in plasmid-carrying prototrophs, further promoting cellular functional heterogeneity. Finally, we show that targeted genetic engineering can be used to suppress cross-feeding and reduce the frequency of plasmid-free auxotrophs, or to exploit it for intentional population diversification and division of labour in co-culture systems.

Hyperpolarisation of Mitochondrial Membranes Is a Critical Component of the Antifungal Mechanism of the Plant Defensin, Ppdef1

Plant defensins are a large family of small cationic proteins with diverse functions and mechanisms of action, most of which assert antifungal activity against a broad spectrum of fungi. The partial mechanism of action has been resolved for a small number of members of plant defensins, and studies have revealed that many act by more than one mechanism. The plant defensin Ppdef1 has a unique sequence and long loop 5 with fungicidal activity against a range of human fungal pathogens, but little is known about its mechanism of action. We screened the S. cerevisiae non-essential gene deletion library and identified the involvement of the mitochondria in the mechanism of action of Ppdef1. Further analysis revealed that the hyperpolarisation of the mitochondrial membrane potential (MMP) activates ROS production, vacuolar fusion and cell death and is an important step in the mechanism of action of Ppdef1, and it is likely that a similar mechanism acts in Trichophyton rubrum.

A peroxidase-derived ligand that induces Fusarium graminearum Ste2 receptor-dependent chemotropism

Introduction

The fungal G protein-coupled receptors Ste2 and Ste3 are vital in mediating directional hyphal growth of the agricultural pathogen Fusarium graminearum towards wheat plants. This chemotropism is induced by a catalytic product of peroxidases secreted by the wheat. Currently, the identity of this product, and the substrate it is generated from, are not known.

Methods and results

We provide evidence that a peroxidase substrate is derived from F. graminearum conidia and report a simple method to extract and purify the FgSte2-activating ligand for analyses by mass spectrometry. The mass spectra arising from t he ligand extract are characteristic of a 400 Da carbohydrate moiety. Consistent with this type of molecule, glycosidase treatment of F. graminearum conidia prior to peroxidase treatment significantly reduced the amount of ligand extracted. Interestingly, availability of the peroxidase substrate appears to depend on the presence of both FgSte2 and FgSte3, as knockout of one or the other reduces the chemotropism-inducing effect of the extracts.

Conclusions

While further characterization is necessary, identification of the F. graminearum-derived peroxidase substrate and the FgSte2-activating ligand will unearth deeper insights into the intricate mechanisms that underlie fungal pathogenesis in cereal crops, unveiling novel avenues for inhibitory interventions.

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.

Quantitative analysis of tRNA abundance and modifications by nanopore RNA sequencing

Transfer RNAs (tRNAs) play a central role in protein translation. Studying them has been difficult in part because a simple method to simultaneously quantify their abundance and chemical modifications is lacking. Here we introduce Nano-tRNAseq, a nanopore-based approach to sequence native tRNA populations that provides quantitative estimates of both tRNA abundances and modification dynamics in a single experiment. We show that default nanopore sequencing settings discard the vast majority of tRNA reads, leading to poor sequencing yields and biased representations of tRNA abundances based on their transcript length. Re-processing of raw nanopore current intensity signals leads to a 12-fold increase in the number of recovered tRNA reads and enables recapitulation of accurate tRNA abundances. We then apply Nano-tRNAseq to Saccharomyces cerevisiae tRNA populations, revealing crosstalks and interdependencies between different tRNA modification types within the same molecule and changes in tRNA populations in response to oxidative stress.

Integrating non-mammalian model organisms in the diagnosis of rare genetic diseases in humans

Next-generation sequencing technology has rapidly accelerated the discovery of genetic variants of interest in individuals with rare diseases. However, showing that these variants are causative of the disease in question is complex and may require functional studies. Use of non-mammalian model organisms - mainly fruitflies (Drosophila melanogaster), nematode worms (Caenorhabditis elegans) and zebrafish (Danio rerio) - enables the rapid and cost-effective assessment of the effects of gene variants, which can then be validated in mammalian model organisms such as mice and in human cells. By probing mechanisms of gene action and identifying interacting genes and proteins in vivo, recent studies in these non-mammalian model organisms have facilitated the diagnosis of numerous genetic diseases and have enabled the screening and identification of therapeutic options for patients. Studies in non-mammalian model organisms have also shown that the biological processes underlying rare diseases can provide insight into more common mechanisms of disease and the biological functions of genes. Here, we discuss the opportunities afforded by non-mammalian model organisms, focusing on flies, worms and fish, and provide examples of their use in the diagnosis of rare genetic diseases.

MEM-FET: Essential protein prediction using membership feature and machine learning approach

Proteins are played key roles in different functionalities in our daily life. All functional roles of a protein are a bit enhanced in interaction compared to individuals. Identification of essential proteins of an organism is a time consume and costly task during observation in the wet lab. The results of observation in wet lab always ensure high reliability and accuracy in the biological ground. Essential protein prediction using computational approaches is an alternative choice in research. It proves its significance rapidly in day-to-day life as well as reduces the experimental cost of wet lab effectively. Existing computational methods were implemented using Protein interaction networks (PPIN), Sequence, Gene Expression Dataset (GED), Gene Ontology (GO), Orthologous groups, and Subcellular localized datasets. Machine learning has diverse categories of features that enable to model and predict essential macromolecules of understudied organisms. A novel methodology MEM-FET (membership feature) is predicted based on features, that is, edge clustering coefficient, Average clustering coefficient, subcellular localization, and Gene Ontology within a compartment of common neighbors. The accuracy (ACC) values of the predicted true positive (TP) essential proteins are 0.79, 0.74, 0.78, and 0.71 for YHQ, YMIPS, YDIP, and YMBD datasets. An enriched set of essential proteins are also predicted using the MEM-FET algorithm. Ensemble ML also validated the proposed model with an accuracy of 60%. It has been predicted that MEM-FET algorithms outperform other existing algorithms with an ACC value of 80% for the yeast dataset.

GetPrimers: A generalized PCR-based genetic targeting primer designer enabling easy and standardized targeted gene modification across multiple systems

Genetic targeting (e.g., gene knockout and tagging) based on polymerase chain reaction (PCR) is a simple yet powerful approach for studying gene functions. Although originally developed in classic budding and fission yeast models, the same principle applies to other eukaryotic systems with efficient homologous recombination. One-step PCR-based genetic targeting is conventionally used but the sizes of the homologous arms that it generates for recombination-mediated genetic targeting are usually limited. Alternatively, gene targeting can also be performed via fusion PCR, which can create homologous arms that are orders of magnitude larger, therefore substantially increasing the efficiency of recombination-mediated genetic targeting. Here, we present GetPrimers (https://www.evomicslab.org/app/getprimers/), a generalized computational framework and web tool to assist automatic targeting and verification primer design for both one-step PCR-based and fusion PCR-based genetic targeting experiments. Moreover, GetPrimers by design runs for any given genetic background of any species with full genome scalability. Therefore, GetPrimers is capable of empowering high-throughput functional genomic assays at multipopulation and multispecies levels. Comprehensive experimental validations have been performed for targeting and verification primers designed by GetPrimers across multiple organism systems and experimental setups. We anticipate GetPrimers to become a highly useful and popular tool to facilitate easy and standardized gene modification across multiple systems.

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

Humanization reveals pervasive incompatibility of yeast and human kinetochore components

Kinetochores assemble on centromeres to drive chromosome segregation in eukaryotic cells. Humans and budding yeast share most of the structural subunits of the kinetochore, whereas protein sequences have diverged considerably. The conserved centromeric histone H3 variant, CenH3 (CENP-A in humans and Cse4 in budding yeast), marks the site for kinetochore assembly in most species. A previous effort to complement Cse4 in yeast with human CENP-A was unsuccessful; however, co-complementation with the human core nucleosome was not attempted. Previously, our lab successfully humanized the core nucleosome in yeast; however, this severely affected cellular growth. We hypothesized that yeast Cse4 is incompatible with humanized nucleosomes and that the kinetochore represented a limiting factor for efficient histone humanization. Thus, we argued that including the human CENP-A or a Cse4-CENP-A chimera might improve histone humanization and facilitate kinetochore function in humanized yeast. The opposite was true: CENP-A expression reduced histone humanization efficiency, was toxic to yeast, and disrupted cell cycle progression and kinetochore function in wild-type (WT) cells. Suppressors of CENP-A toxicity included gene deletions of subunits of 3 conserved chromatin remodeling complexes, highlighting their role in CenH3 chromatin positioning. Finally, we attempted to complement the subunits of the NDC80 kinetochore complex, individually and in combination, without success, in contrast to a previous study indicating complementation by the human NDC80/HEC1 gene. Our results suggest that limited protein sequence similarity between yeast and human components in this very complex structure leads to failure of complementation.

Putative purine nucleoside interacting residues in the malaria parasite purine uptake transporter PfENT1 are critical for transporter function

Malaria remains a major public health threat for billions of people worldwide. Infection with obligate intracellular, unicellular parasites from the genus Plasmodium causes malaria. Plasmodium falciparum causes the deadliest form of human malaria. Plasmodium parasites are purine auxotrophic. They rely on purine import from the host red blood cell cytoplasm via equilibrative nucleoside transporters to supply substrates to the purine salvage pathway. We previously developed a high throughput screening assay to identify inhibitors of the P. falciparum Equilibrative Nucleoside Transporter Type 1 (PfENT1). Screening a small molecule library identified PfENT1 inhibitors that blocked proliferation of P. falciparum parasites in in vitro culture. The goal of the current work was to validate a high-resolution model of PfENT1 predicted by the AlphaFold protein structure prediction program. We superimposed the predicted PfENT1 structure on the human homologue structure, hENT1, and developed a structure-based sequence alignment. We mutated the residues in PfENT1 aligned with and flanking the residues in hENT1 that interact with the purine analog, nitrobenzylthioinosine (NBMPR). Mutation of the PfENT1 residues Q135, D287, and R291 that are predicted to form hydrogen bonds to purine nucleosides eliminated purine and pyrimidine transport function in various yeast-based growth and radiolabeled substrate uptake assays. Mutation of two flanking residues, W53 and S290, also resulted in inactive protein. Mutation of L50 that forms hydrophobic interactions with the purine nucleobase reduced transport function. Based on our results the AlphaFold predicted structure for PfENT1 may be useful in guiding medicinal chemistry efforts to improve the potency of our PfENT1 inhibitors.

Characterization of a novel RNAi yeast insecticide that silences mosquito 5-HT1 receptor genes

G protein-coupled receptors (GPCRs), which regulate numerous intracellular signaling cascades that mediate many essential physiological processes, are attractive yet underexploited insecticide targets. RNA interference (RNAi) technology could facilitate the custom design of environmentally safe pesticides that target GPCRs in select target pests yet are not toxic to non-target species. This study investigates the hypothesis that an RNAi yeast insecticide designed to silence mosquito serotonin receptor 1 (5-HTR1) genes can kill mosquitoes without harming non-target arthropods. 5-HTR.426, a Saccharomyces cerevisiae strain that expresses an shRNA targeting a site specifically conserved in mosquito 5-HTR1 genes, was generated. The yeast can be heat-inactivated and delivered to mosquito larvae as ready-to-use tablets or to adult mosquitoes using attractive targeted sugar baits (ATSBs). The results of laboratory and outdoor semi-field trials demonstrated that consumption of 5-HTR.426 yeast results in highly significant mortality rates in Aedes, Anopheles, and Culex mosquito larvae and adults. Yeast consumption resulted in significant 5-HTR1 silencing and severe neural defects in the mosquito brain but was not found to be toxic to non-target arthropods. These results indicate that RNAi insecticide technology can facilitate selective targeting of GPCRs in intended pests without impacting GPCR activity in non-targeted organisms. In future studies, scaled production of yeast expressing the 5-HTR.426 RNAi insecticide could facilitate field trials to further evaluate this promising new mosquito control intervention.

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

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

Genome-wide identification of resistance genes and cellular analysis of key gene knockout strain under 5-hydroxymethylfurfural stress in Saccharomyces cerevisiae

In bioethanol production, the main by-product, 5-hydroxymethylfurfural (HMF), significantly hinders microbial fermentation. Therefore, it is crucial to explore genes related to HMF tolerance in Saccharomyces cerevisiae for enhancing the tolerance of ethanol fermentation strains. A comprehensive analysis was conducted using genome-wide deletion library scanning and SGAtools, resulting in the identification of 294 genes associated with HMF tolerance in S. cerevisiae. Further KEGG and GO enrichment analysis revealed the involvement of genes OCA1 and SIW14 in the protein phosphorylation pathway, underscoring their role in HMF tolerance. Spot test validation and subcellular structure observation demonstrated that, following a 3-h treatment with 60 mM HMF, the SIW14 gene knockout strain exhibited a 12.68% increase in cells with abnormal endoplasmic reticulum (ER) and a 22.41% increase in the accumulation of reactive oxygen species compared to the BY4741 strain. These findings indicate that the SIW14 gene contributes to the protection of the ER structure within the cell and facilitates the clearance of reactive oxygen species, thereby confirming its significance as a key gene for HMF tolerance in S. cerevisiae.