Genome-wide identification of resistance genes and transcriptome regulation in yeast to accommodate ammonium toxicity
BMC Genomics volume 23, Article number: 514 (2022)
Ammonium is an important raw material for biomolecules and life activities, and the toxicity of ammonium is also an important ecological and agricultural issue. Ammonium toxicity in yeast has only recently been discovered, and information on its mechanism is limited. In recent years, environmental pollution caused by nitrogen-containing wastewater has been increasing. In addition, the use of yeast in bioreactors to produce nitrogen-containing compounds has been developed. Therefore, research on resistance mechanisms that allow yeast to grow under conditions of high concentrations of ammonium has become more and more important.
To further understand the resistance mechanism of yeast to grow under high concentration of ammonium, we used NH4Cl to screen a yeast non-essential gene-deletion library. We identified 61 NH4Cl-sensitive deletion mutants from approximately 4200 mutants in the library, then 34 of them were confirmed by drop test analysis. Enrichment analysis of these 34 genes showed that biosynthesis metabolism, mitophagy, MAPK signaling, and other pathways may play important roles in NH4Cl resistance. Transcriptome analysis under NH4Cl stress revealed 451 significantly upregulated genes and 835 significantly downregulated genes. The genes are mainly enriched in: nitrogen compound metabolic process, cell wall, MAPK signaling pathway, mitophagy, and glycine, serine and threonine metabolism.
Our results present a broad view of biological pathways involved in the response to NH4Cl stress, and thereby advance our understanding of the resistance genes and cellular transcriptional regulation under high concentration of ammonium.
Ammonium is a building block for amino acids, nucleic acids, polysaccharides, and other important cellular structural components in all organisms, and at the same time is also an important raw material for the synthesis of secondary metabolites such as alkaloids and polyamines . Ammonium is a paradoxical nutrient ion. First, ammonium is the key nitrogen source for all life forms, including bacteria, fungi, protists, plants, and animals. Whereas the ammonium used by plants mainly comes from natural nitrogen fixation, in heterotrophic cells, it is mainly derived from cellular metabolites [2, 3]. Second, high concentrations of ammonium are toxic to organisms [4,5,6]. The downstream molecular events of NH4+ as a nutrient component have been extensively studied; it is mainly related to biosynthesis, carbon metabolism, energy metabolism, primary nitrogen metabolism, plant hormones, cell wall stability, and signaling pathways . However, the mechanism underlying ammonium toxicity, and how to inhibit it, are areas of intense current research [8,9,10].
In plants, excessive NH4+ can inhibit photosynthesis, block the growth of plant roots and leaves , and promote the accumulation of reactive oxygen species . When ammonium is the primary nitrogen source, it usually significantly inhibits plant growth, manifesting as a phenotype of severely inhibited root growth and chlorosis, which is a well-known manifestation of ammonium salt toxicity . For vertebrates, high concentrations of ammonium may displace K+ and depolarize neurons, leading to the activation of NMDA-type glutamate receptor, which leads to excessive Ca2+ influx and subsequent cell death in the central nervous system, ultimately leading to convulsions, coma, and even death . In addition, ammonium salts are also important environmental pollutants in the ecosystem. Excessive use of ammonium fertilizer will not only cause soil acidification, resulting in the loss of soil nutrients and the aggravation of the harm of heavy metal pollution, but also lead to secondary salinization of the soil, resulting in large-scale reductions in crop yield and deterioration of quality .
The harm of ammonium to animals and plants has long been known, and extensively studied; however, the discovery of ammonium toxicity in Saccharomyces cerevisiae is relatively recent . The reason for this is that yeast cells have extremely high tolerance to ammonium under normal culture conditions . In Saccharomyces cerevisiae, ammonium uptake is facilitated by ammonium permeases (Mep1, Mep2, and Mep3), which transport NH4+ ions and conduct ATP-dependent export of protons to maintain intracellular charge balance and pH. This process can put a great burden on the energy metabolism of cells, especially when there is an excess of ammonium salts outside the cells . With growing interest in the biosynthesis of nitrogen-containing compounds using yeast bioreactors, the mechanisms of ammonia resistance in yeast cells have also received increased attention . Therefore, it is particularly important to fully understand the mechanism of action of ammonium metabolism in yeast cells and to find and discover resistance genes and pathways. This could aid in the development of ways to increase cellular ammonia resistance while reducing intracellular energy expenditure.
To gain a comprehensive and in-depth understanding of the causes and mechanisms underlying ammonium toxicity in yeast cells, we performed experiments using two genome-wide strategies. First, we used the yeast gene deletion library (SGA collection) to screen NH4Cl resistance-related mutants. Next, we analyzed transcriptome changes in yeast cells at concentrations of ammonium that significantly inhibited growth. Finally, the present study reveals the key genes and cellular pathways underlying resistance to NH4Cl in the yeast genome, which provides an experimental basis for further cultivating high ammonium resistant strains and designing schemes to reduce ammonium toxicity.
Genome-wide screen for NH4Cl-sensitive mutants
To understand the specific mechanism of yeast under NH4Cl stress, we used the SGA-V2 library to screen for mutants that are sensitive to NH4Cl. First, we diluted the control strain his∆ in plates with a gradient of concentrations of NH4Cl. The growth of the strain began to be inhibited under a concentration of 500 mM NH4Cl (Fig. S1). This is consistent with previous reports that yeast have high NH4Cl tolerance . After determining the concentration range of NH4Cl inhibition on yeast, one plate (SGA-V2–3) from the library was randomly selected for a preliminary experiment. The NH4Cl concentration was chosen with reference to our previous criteria ; that is, the colony size of the control strain was reduced by nearly half, and some mutants showed significant growth inhibition. This criterion was met when the NH4Cl concentration reached 800 mM (Fig. 1), so we chose 800 mM NH4Cl as the treatment concentration for the genome-wide screen. About 4200 genes in the yeast gene deletion collection were screened. Based on the criteria described in the method, a total of 61 mutants were considered sensitive (Table S1).
Drop test verifies high-throughput screening results
To verify the reliability of the high-throughput screening results, we performed drop test experiments on each of the 61 deletion strains. Among these deletion strains, 34 remained sensitive to NH4Cl in the drop test experiments (Fig. 2). Growth of these 34 deletion strains was inhibited, compared with control strains, in the presence of 800 mM NH4Cl. Information on these 34 mutants, together with their quantitative fitness scores in genome-wide screening, are summarized in Table 1. Based on the results of the drop test, we classified the degree of NH4Cl inhibition of the mutants. Colony growth only in the first column was defined as “serious” inhibition, colony growth only in the first two columns was defined as “obvious” inhibition, “moderate” referred to colony growth in the first three columns, and “slight” referred to the first four columns with colony growth. Seven of the strains showed a relatively serious inhibitory effect, including sla1Δ, swf1Δ, yme1Δ, slt2Δ, hog1Δ, vps51Δ, and ure2Δ. In the other 27 mutants, 10 strains were obviously inhibited (dos2Δ, rvs167Δ, vrp1Δ, gtr1Δ, yor008c-aΔ, ser1Δ, ser2Δ, aim29Δ, gyp1Δ, and aim44Δ); 11 strains were moderately inhibited (stb5Δ, tps1Δ, gtr2Δ, meh1Δ, uth1Δ, pat1Δ, erg6Δ, erg2Δ, ubp15Δ, pfa4Δ, and mnn10Δ); and 6 strains were slightly inhibited (inp53Δ, sgf29Δ, lcb4Δ, rtt103Δ, ecm8Δ, and snf1Δ) (Fig. 2). The deleted genes in the various mutants were distributed in a variety of different cellular pathways, indicating that cellular mechanisms involved in NH4Cl cellular resistance involve not single, but multiple, pathways.
Since low concentrations of NH4Cl have no obvious inhibitory effect on yeast growth, we used a high concentration of 800 mM NH4Cl for high-throughput screening. This concentration resulted in an obvious inhibitory effect on wild-type yeast. To verify the sensitivity of the mutants we screened to NH4Cl, we performed gradient NH4Cl experiments on the 34 mutants. A concentration gradient of 0, 100, 300, and 500 mM NH4Cl was set up for the drop test, in accordance with the previous serial dilution method. At a concentration of 100 mM, the growth of yme1Δ and ure2Δ was severely inhibited, indicating that these two genes are essential for cellular resistance to NH4Cl toxicity (Fig. 3). In addition, rvs167Δ, slt2Δ, vps51Δ, gtr1Δ, ser1Δ, aim29Δ, erg2Δ, erg6Δ, gyp1Δ, and aim44Δ were significantly inhibited at a concentration of 300 mM NH4Cl. At a concentration of 500 mM NH4Cl, most of the 34 mutants showed growth inhibition, while the growth of the control strain was still relatively normal, indicating that the 34 genes obtained by high-throughput screening played different roles in the resistance of cells to NH4Cl. Deletion of these genes makes cells sensitive to different degrees of NH4Cl stress. Previous studies have revealed a set of multidrug resistant genes whose deletions are associated with sensitivity to multiple compounds with diverse modes of action that are utilized by cells in a wide range of stress responses [20,21,22]. Among the 34 genes corresponding to sensitive deletion mutants, 13 were previously thought to contribute to general multidrug resistance, including ERG2, ERG6, GTR1, PAT1, RVS167, SLA1, SLT2, SNF1, URE2, VPS51, VRP1, YKR007W, and YME1.
Classification of genes related to NH4Cl resistance
To further understand the function of the screened genes, we conducted a bioinformatics analysis of these genes. First, to understand the localization of these genes, we analyzed and mapped all of the genes for the cellular locations of their encoded proteins (Fig. 4A). Most of the gene products were located in the cytoplasm, nucleus, and endoplasmic reticulum (ER). In addition, MNN10 and VPS51 were located in the Golgi apparatus; GTR1, GTR2, and MEH1 were located in the vacuole; and UTH1 and YME1 were located in the mitochondria. Further, several of the identified genes encoded products localized to more than one cellular location: SGF29, SER2, and SLT2 were distributed in both the cytoplasm and the nucleus; YOR008C-A was located in both the ER and vacuole; LCB4 was located in both the ER and Golgi apparatus; and ERG6 was distributed in both the ER and mitochondria. GO and KEGG analysis of the 34 genes found that they were mainly enriched in categories related to Gtr1-Gtr2 GTPase complex, cytoskeleton, MAPK activity, TOR signaling, mitophagy, and metabolic pathways (Fig. 4B and C, Table S2). In addition, functional classification of the 34 resistance genes was performed according to the functional description of the Saccharomyces genome database (https://www.yeastgenome.org/) (Fig. 4D). The genes for five of the deletion mutants were unidentified in terms of function. Of the remaining 29, the genes could be divided into seven groups: lipid metabolism- and ammonia metabolism-related genes, post-translational modification, mitophagy, Gtr1-Gtr2 GTPase complex, DNA and RNA synthesis-related genes, and cell resistance- and endocytosis-related genes. Overall, this analysis revealed the complexity of the machinery involved in resistance to NH4Cl toxicity. NH4Cl may influence cells by affecting both growth and metabolism, and the integrity of these pathways plays an important role in cells’ resistance to ammonium toxicity.
Transcriptome analysis in response to NH4Cl stress
In order to further understand the mechanism of NH4Cl toxicity on cells, we cultured wild-type strain BY4741 with different concentrations of NH4Cl to find suitable processing conditions for transcriptome sequencing (Fig. S2). Under the condition of 800 mM NH4Cl treatment, the cells grew for 8 hours and entered the logarithmic growth phase. At this time, the effect of NH4Cl on the cells was already reflected at the transcription level. Therefore, transcriptome sequencing analysis was performed using this treatment condition, and untreated samples grown at the same time were used as a control.
We analyzed the expression levels of 5966 genes in total, of which there were 451 significantly upregulated genes and 835 significantly downregulated genes. Then, we performed GO and KEGG analysis on the 451 significantly upregulated genes. The main enrichment pathways were cellular nitrogen compound metabolic process, cell wall, MAPK signaling pathway, oxidative phosphorylation, and TCA cycle (Fig. 5A and B). GO and KEGG analysis of 835 significantly downregulated genes found that these genes were mainly enriched in ubiquitin-mediated proteolysis, protein phosphorylation, lipid metabolic process, mitophagy, and glycine, serine, and threonine metabolism (Fig. 5C and D). It has been demonstrated that under stress, cells tend to optimize cellular resources for stress adaptation, leading to the massive expression of genes involved in stress adaptation accompanied by inhibition of expression of genes involved in proliferation and cell cycle progression .
To better compare our screening and transcriptome sequencing analysis results, we performed transcriptomic analysis of the expression of the 34 deleted genes. The significantly upregulated genes were TPS1, UTH1, and GYP1. The significantly downregulated genes were SER1, SER2, ERG6, AIM44, and YOR008C-A. The genes for which expression levels were not significantly downregulated were ECM8, SLA1, SGF29, PAT1, MNN10, RTT103, RVS167, DOS2, SWF1, GTR2, SLT2, STB5, MEH1, VPS51, AIM29, HOG1, VRP1, GTR1, ERG2, UBP15, URE2, PFA4, INP53, LCB4, YME1, and SNF1 (Fig. 6A). Among the 34 genes screened, most showed no obvious expression changes in the transcriptome, a lack of such correlation was also previously reported .
Next, a correlation estimate was performed for all genes identified in both the mutant screening and transcriptomic analysis. As shown (Fig. 6B), the genes with significantly upregulated expression and score ≤ − 0.2 were TPS1 and TPS2, and the genes with significantly downregulated expression and score ≤ − 0.2 were BNA1, ICS2, AIM44, TOM70, ERG6, BUG1, PUN1, SMI1, YOR008C-A, SER1, and SER2. The genes with significantly downregulated expression and score ≥ 0.2 were GFD2, COQ10, MVP1, RKM3, and PRR1.
Compared with other organisms, the tolerance of yeast to ammonium is unimaginable, and the discovery of ammonium toxicity in yeast occurred significantly later than that of other stresses. It was not until the discovery of the relationship between NH4+ and K+ that the major mystery of how yeast could accomplish this was solved . However, more mechanisms, especially in high concentrations of ammonium under normal K+ culture conditions, have not been reported. We aimed to learn more about the mechanisms, through (1) screening of yeast deletion mutant libraries, to identify resistance genes, and (2) transcriptome sequence analysis, to identify transcriptional changes, under high-concentration ammonium treatment. In our results, we found many shared genes and pathways that were involved in resistance to multiple stresses, and also found many genes and pathways that were specific to ammonium resistance.
Carbon and nitrogen metabolism maintains cellular ammonium resistance
One of the ways in which cells respond to excess intracellular ammonium is to synthesize ammonium into amino acids and exclude the excess to the extracellular space [15, 16]. This process consumes a large amount of carbon skeleton, which is needed for processes such as glycolysis and the TCA cycle. This can be confirmed in the enrichment results of upregulated genes in our transcriptome analysis. Studies by others have come to similar conclusions. For example, In a genome-wide high-throughput screening study of multiple N-nitrosamine compounds using a library of yeast deletion mutants, Joseph Uche Ogbede et al. found that ammonium sulfate caused growth defects in mutants of the arginine biosynthesis pathway . Xueping Tian et al. discovered overall upregulation of starch synthesis and degradation pathways, as well as glycolysis and the TCA cycle, which could provide abundant carbon skeleton material for excess NH4+ assimilation and avoid carbon deficiency, a mechanism that might also be an important feature of the duckweed response to NH4+ toxicity . The TCA cycle serves as a mitochondria-based hub for the final steps in carbon skeleton oxidative catabolism for carbohydrates, amino acids, and fatty acids . NH4+ can be effectively detoxified into amino acids to maintain the availability of the carbon skeleton by utilizing the TCA cycle . Chandran et al. discovered that upregulated genes are mainly enriched in metabolic processes for diverse amino acids as well as nitrogen compounds, by exploring the genomic responses in rice roots dealing with 0.5 mM (NH4)2SO4 . Rui Wang et al. found that a complex physiological and genetic regulatory network of processes including nitrogen metabolism, carbon metabolism, abiotic stress response, and secondary metabolism at the root and leaf levels was involved in NH4+ resistance in Myriophyllum aquaticum . The same response has been observed in mudskippers under NH4+ stress. Xinxin You et al. discovered that reducing the catabolism of protein and amino acids could be an effective way to slow down internal ammonia accumulation . It is reasonable that genes in various other cellular metabolic processes involving nitrogen could be downregulated. Similar to their results, our upregulated genes were also enriched in these pathways, which indicates that it is possible to improve tolerance to NH4+ toxicity in yeast by increasing the expression of the TCA cycle, glycolysis, starch and sucrose metabolism (Fig. 5B).
In the GO enrichment analysis of the 34 genes corresponding to sensitive deletion mutants, URE2, GTR1, and GTR2 were enriched in several pathways at the highest enrichment level. According to our experimental results, ure2Δ is one of two mutants showing significant growth inhibition at 100 mM NH4Cl, and URE2 is not a multidrug resistant gene. The growth of both gtr1Δ and gtr2Δ was reduced by more than 40% compared with the wild type in the screening results. Ure2 is involved in the inhibition of nitrogen catabolism . When an optimal nitrogen source is available, Ure2 acts as a transcriptional corepressor and downregulates the expression of many genes involved in nitrogen utilization. The Gtr1-Gtr2 GTPase complex is composed of Gtr1 and Gtr2, and is important for sensing the presence of amino acids in the medium by activating TORC1 . Gtr1 is a subunit of the EGO complex, which is responsible for activating TORC1 in response to the utilization of amino acids. The combination of EGO complex and TOR can positively regulate microautophagy . In contrast, gtr1Δ reduces the activity of TORC1, which increases the expression level of nitrogen transporters and ammonium consumption but reduces amino acid consumption . GTR2 encodes a Ras-like small GTPase that plays a role in regulating nutrition-responsive TORC1 kinase signal transduction, exocytosis sorting of endosomes, and epigenetic control of gene expression . Deletion of GTR1 or GTR2 may render cells incapable of activating TORC1, resulting in cell toxicity due to the inability to degrade ammonia in vivo under high ammonia stress.
Osmotic stress is one of the causes of ammonium toxicity
TPS1 is one of the few genes in our screen for which deletion affects cell growth, but which is significantly upregulated in the transcriptome. TPS1 encodes trehalose 6-phosphate synthase and catalyzes the first step in trehalose biosynthesis. In Saccharomyces cerevisiae, trehalose is a major reserve carbohydrate involved in responses to thermal, osmotic, oxidative, and ethanol stresses . Genes in the trehalose metabolic pathway, including TPS1, TPS2, and NTH1, were all upregulated. TPS2 is responsible for reaction catalysis of step 2 in trehalose synthesis. The NTH1 gene product might contribute to trehalose mobilization in S. cerevisiae under saline stress conditions . NTH1, as well as trehalose biosynthesis genes, was upregulated under saline stress conditions . Cells may be able to reduce the toxicity caused by NH4Cl through saline stress by increasing the expression of these genes (TPS1, TPS2, and NTH1).
In addition, Hog1 is an important regulator of transcription in conditions of osmotic stress in yeast . Fang Li et al. found that Hog1 and Slt2 were downregulated under osmotic and cell wall stresses, respectively . Hog1-mediated transcriptional control of ERG genes accounts for most of the downregulation of sterol levels upon osmotic stress . In our transcriptome results, the expression of HOG1, SLT2, ERG2, and ERG6 were all downregulated, and erg2Δ and erg6Δ were also more sensitive than other mutants, showing severe growth inhibition at 300 mM NH4Cl. These results suggest that NH4Cl may also be toxic to cells through osmotic stress.
Mitophagy plays an important role in ammonium resistance
Among the 34 sensitive mutants, KEGG and functional classification were all enriched in the mitophagy pathway, including yme1Δ, slt2Δ, and hog1Δ, of which yme1Δ showed severe growth inhibition at 100 mM NH4Cl, and slt2Δ was significantly inhibited at 300 mM NH4Cl. The mitochondrial protease Yme1 plays an important role in the ability of cells lacking tafazzin function to maintain mitochondrial structural integrity, mitochondrial quality control, and mitochondrial autophagy . Yme1 is an ATP-dependent protease located on the inner mitochondrial membrane that is required for the growth of yeast lacking a complete mitochondrial genome . The processing of Atg32 by Yme1 is an important regulatory mechanism of mitophagy . Slt2 and Hog1 are both required for mitophagy. Slt2 is a MAPK of the cell wall integrity pathway and is necessary for the degradation of mitochondria. Hog1, another member of the MAPK family, is a kinase that regulates and is regulated by Sch9p and is independent of the PKA and TOR pathways in response to stress . Hog1 plays a role in Atg32 phosphorylation, and this is necessary for mitophagy [48, 49].
Mitochondria are important organelles that provide cellular energy and the carbon skeleton building blocks for the synthesis of macromolecular substances. Therefore, the removal of damaged mitochondria through mitophagy is essential for maintaining proper cell function . When cells are exposed to ammonium, higher demands are placed on mitochondrial function. This is because, first, every time NH4+ enters the cell, H+ is transported outside the cell, resulting in an increase in intracellular pH . Ammonium uptake requires the consumption of large amounts of ATP for maintaining intracellular pH by using the plasma membrane bound H+-ATPase .. In addition, studies have found that ineffective transmembrane ammonium cycling in some species can also lead to energy consumption, because a large fraction of intracellular ammonium leaks out of the cell through the membrane. This ineffective transport causes great energy loss, resulting in adverse symptoms [52, 53]. Second, the assimilation of ammonium into amino acids requires a large amount of carbon skeleton, which leads to harmful effects due to insufficient carbon building blocks in the cell . The high concentration of ammonium increases the synthesis rate of glutamic acid, alanine, and glycine, and accelerates the consumption of glucose and glutamine. The increase in the concentration of ammonium salt not only increases the energy requirements of the cells, but also affects the TCA cycle. These all place high demands on mitochondrial function, and the impairment of mitophagy leads to the sensitivity of cells to ammonium toxicity.
In this study, we analyzed the resistance mechanism that allow yeast to grow under high concentration of ammonium using yeast genome-wide screening. Our results showed that, out of nearly 4200 mutants, 34 mutants were identified and confirmed by drop test as being vulnerable to NH4Cl. Furthermore, functional enrichment analysis indicated that these 34 genes were mainly involved in lipid metabolism and ammonia metabolism, post-translational modification, mitophagy, Gtr1-Gtr2 GTPase complex, DNA and RNA synthesis-related genes, cell resistance, and endocytosis-related genes. Transcriptome analysis further supported the accuracy of our screening results and demonstrated that cells significantly upregulated carbon and nitrogen metabolism, TCA cycle and other stress adaptation pathways, and downregulated cell growth-related pathways under NH4Cl stress. These results can provide us with a clearer understanding of the resistance mechanism that allow yeast to grow under NH4Cl.
Materials and methods
Genome-wide screen to identify gene deletion mutants sensitive to NH4Cl
The gene deletion library (SGA-V2) was kindly provided by Prof. Charlie Boone, University of Toronto, Canada [54, 55]. The library of non-essential haploid deletion strains containing about 4200 mutants was started from S. cerevisiae strain BY4741 (MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0). Each mutant was constructed by replacing the corresponding ORF with a KanMX cassette, and the strains in the library were arranged in 384 format. Mutant his3Δ::KanR in this library was designated as the control strain and added as a border around four edges of each plate . The pinning steps for library handling were performed using a SINGER ROTOR HDA Robot (Singer Instruments, UK). The deletion library was cloned on YPD agar plates (with G418 added) and grown at 30 °C for 2 days. The 384 strains in each plate were then transferred to a new agar plate, and each colony was repeated four times to finally form an array of 1536 colonies. The entire library in 1536 array was then cloned onto a solid plate with or without 800 mM NH4Cl and grown at 30 °C for 2 days. The growth status of the colonies was photographed by PhenoBooth (Singer Instruments, UK). Images were analyzed using SGAtools (http://sgatools.ccbr.utoronto.ca/) to evaluate the growth of the colonies , and compare the growth of each mutant with or without NH4Cl. First, images of plates with colonies were processed to give quantified colony sizes for the screen. Next, the colony sizes were normalized and filtered within plates, taking into account position effects and other confounding factors. Ratios of normalized colony sizes from NH4Cl-treated and untreated mutants were used as a measure of sensitivity. A score > 0 represented positive interaction, i.e., increased colony size, while a score < 0 represented negative interaction, i.e., decreased colony size. According to previous studies , a score of < − 0.2 generally indicates a strong effect, so we used scores < − 0.2 and P < 0.05 as thresholds to identify significantly affected mutants. Each experiment was repeated three times.
Yeast strains were grown in YPD + G418 medium which contained 1% yeast extract, 2% peptone, 2% glucose, and 200 mg/mL G418. The 800 mM NH4Cl medium consisted of 1% yeast extract, 2% peptone, 2% glucose, 200 mg/mL G418, and 800 mM NH4Cl. All the strains which were mentioned in this article were cultivated in an incubator at 30 °C.
For the drop test, we picked up the selected mutant from the SGA-V2 library and cultured it overnight in 3 mL YPD + G418 medium at 30 °C. Then, the strains were diluted to OD600 = 0.1, and cultured at 30 °C until the mid-logarithmic stage. The culture was then continuously diluted tenfold in sterile water in 96-well plates and the dilutions were spotted on plates with or without NH4Cl. Images were taken after culturing at 30 °C for 2 days, and the growth of yeast colonies was observed. Shown are representative drop tests from three independent replicate assays.
Bioinformatics enrichment analysis and functional annotation
The GO Term Finder in the Saccharomyces Genome Database (https://www.yeastgenome.org/) was used to analyze the enrichment of GO terms in 34 genes . Then, we used the gene-list analysis tool Metascape (http://metascape.org) with the background of genes corresponding to SGA-V2 library, to conduct KEGG pathway analysis, with a chosen P-value of less than 0.05 . The functional annotation results were mapped and clustered using Cytoscape (version 3.8.0) .
Transcriptional RNA sequence analysis
Yeast cells treated with 0 mM and 800 mM NH4Cl for 8 hours were used for transcriptome sequencing, and four sample replicates were set for each treatment. Total RNA was extracted using Trizol reagent (Thermofisher, 15,596,018), using the Bioanalyzer 2100 and RNA 6000 Nano LabChip Kits (Agilent, CA, USA, 5067–1511) to analyze the total RNA quantity and purity, the sequencing cDNA library was constructed using high-quality RNA with an RNA integrity number (RIN) > 7.0. The library was sequenced using the Illumina Novaseq™ 6000 sequence platform. Cutadapt (https://cutadapt.readthedocs.io/en/stable/, version: cutadapt-1.9) was used to remove the Illumina adapter contamination and for trimming the reads and clipping the low-quality bases to get high-quality clean reads. Then, sequence quality was verified using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/, version: 0.11.9). After that, approximately 6 Gb of clean reads were produced . Differentially expressed genes (DEGs) analysis was performed by DESeq2 software, between two different groups. The genes with the parameter of false discovery rate (FDR) below 0.05 and absolute fold change ≥1.5 were considered differentially expressed genes . Then ClueGO in Cytoscape was used to perform the enrichment analysis of GO function and KEGG pathway for the significantly upregulated genes and the significantly downregulated genes.
Availability of data and materials
All the sequencing data generated in this study have been deposited in the Sequence Read Archive (SRA) database under accession number PRJNA795291 (https://dataview.ncbi.nlm.nih.gov/object/PRJNA795291?reviewer=3v8dgbe27ar7v1vbh4uscmhmvv).
Bagh K, Hiraoki T, Thorpe TA, Vogel HJ. Nitrogen-15 NMR studies of nitrogen metabolism in Picea glauca buds. Plant Physiol Biochem. 2004;42(10):803–9.
Liu Y, von Wirén N. Ammonium as a signal for physiological and morphological responses in plants. J Exp Bot. 2017;68(10):2581–92.
Dimski DS. Ammonia metabolism and the urea cycle: function and clinical implications. J Vet Intern Med. 1994;8(2):73–8.
Bittsánszky A, Pilinszky K, Gyulai G, Komives T. Overcoming ammonium toxicity. Plant Sci. 2015;231:184–90.
Müller T, Walter B, Wirtz A, Burkovski A. Ammonium toxicity in bacteria. Curr Microbiol. 2006;52(5):400–6.
Randall DJ, Tsui TKN. Ammonia toxicity in fish. Mar Pollut Bull. 2002;45(1–12):17–23.
Ariz I, Asensio AC, Zamarreño AM, García-Mina JM, Aparicio-Tejo PM, Moran JF. Changes in the C/N balance caused by increasing external ammonium concentrations are driven by carbon and energy availabilities during ammonium nutrition in pea plants: the key roles of asparagine synthetase and anaplerotic enzymes. Physiol Plant. 2013;148(4):522–37.
Hachiya T, Inaba J, Wakazaki M, Sato M, Toyooka K, Miyagi A, et al. Excessive ammonium assimilation by plastidic glutamine synthetase causes ammonium toxicity in Arabidopsis thaliana. Nat Commun. 2021;12(1):4944.
Wang P. Friend or foe: how nitrate antagonizes ammonium toxicity. Plant Physiol. 2021;186(1):210–1.
Lu Q, Han P, Chen F, Liu T, Li J, Leng L, et al. A novel approach of using zeolite for ammonium toxicity mitigation and value-added Spirulina cultivation in wastewater. Bioresour Technol. 2019;280:127–35.
Coskun D, Britto DT, Li M, Becker A, Kronzucker HJ. Rapid ammonia gas transport accounts for futile transmembrane cycling under NH3/NH4+ toxicity in plant roots. Plant Physiol. 2013;163(4):1859–67.
Patterson K, Cakmak T, Cooper A, Lager I, Rasmusson AG, Escobar MA. Distinct signalling pathways and transcriptome response signatures differentiate ammonium- and nitrate-supplied plants. Plant Cell Environ. 2010;33(9):1486–501.
Britto DT, Kronzucker HJ. NH4+ toxicity in higher plants: a critical review. J Plant Physiol. 2002;159(6):567–84.
Ahmed M, Rauf M, Mukhtar Z, Saeed NA. Excessive use of nitrogenous fertilizers: an unawareness causing serious threats to environment and human health. Environ Sci Pollut Res Int. 2017;24(35):26983–7.
Hess DC, Lu W, Rabinowitz JD, Botstein D. Ammonium toxicity and potassium limitation in yeast. PLoS Biol. 2006;4(11):e351.
Reisser C, Dick C, Kruglyak L, Botstein D, Schacherer J, Hess DC. Genetic basis of ammonium toxicity resistance in a sake strain of yeast: a mendelian case. G3. 2013;3(4):733–40.
Cueto-Rojas HF, Milne N, van Helmond W, Pieterse MM, van Maris AJA, Daran JM, et al. Membrane potential independent transport of NH3 in the absence of ammonium permeases in Saccharomyces cerevisiae. BMC Syst Biol. 2017;11(1):49.
Turner TL, Kim H, Kong II, Liu JJ, Zhang GC, Jin YS. Engineering and evolution of Saccharomyces cerevisiae to produce biofuels and chemicals. Adv Biochem Eng Biotechnol. 2018;162:175–215.
Jin X, Zhang J, An T, Zhao H, Fu W, Li D, et al. A genome-wide screen in Saccharomyces cerevisiae reveals a critical role for oxidative phosphorylation in cellular tolerance to lithium hexafluorophosphate. Cells. 2021;10(4):888.
Parsons AB, Brost RL, Ding H, Li Z, Zhang C, Sheikh B, et al. Integration of chemical-genetic and genetic interaction data links bioactive compounds to cellular target pathways. Nat Biotechnol. 2004;22(1):62–9.
Parsons AB, Lopez A, Givoni IE, Williams DE, Gray CA, Porter J, et al. Exploring the mode-of-action of bioactive compounds by chemical-genetic profiling in yeast. Cell. 2006;126(3):611–25.
Hillenmeyer ME, Fung E, Wildenhain J, Pierce SE, Hoon S, Lee W, et al. The chemical genomic portrait of yeast: uncovering a phenotype for all genes. Science. 2008;320(5874):362–5.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Giaever G, Chu AM, Ni L, Connelly C, Riles L, Véronneau S, et al. Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002;418(6896):387–91.
Birrell GW, Brown JA, Wu HI, Giaever G, Chu AM, Davis RW, et al. Transcriptional response of Saccharomyces cerevisiae to DNA-damaging agents does not identify the genes that protect against these agents. Proc Natl Acad Sci U S A. 2002;99(13):8778–83.
Ogbede JU, Giaever G, Nislow C. A genome-wide portrait of pervasive drug contaminants. Sci Rep. 2021;11(1):12487.
Tian X, Fang Y, Jin Y, Yi Z, Li J, Du A, et al. Ammonium detoxification mechanism of ammonium-tolerant duckweed (Landoltia punctata) revealed by carbon and nitrogen metabolism under ammonium stress. Environ Pollut. 2021;277:116834.
Akram M. Citric acid cycle and role of its intermediates in metabolism. Cell Biochem Biophys. 2014;68(3):475–8.
Vega-Mas I, Cukier C, Coleto I, González-Murua C, Limami AM, González-Moro MB, et al. Isotopic labelling reveals the efficient adaptation of wheat root TCA cycle flux modes to match carbon demand under ammonium nutrition. Sci Rep. 2019;9(1):8925.
Chandran AKN, Priatama RA, Kumar V, Xuan Y, Je BI, Kim CM, et al. Genome-wide transcriptome analysis of expression in rice seedling roots in response to supplemental nitrogen. J Plant Physiol. 2016;200:62–75.
Wang R, Xu S, Sun H, Feng S, Jiang C, Zhou S, et al. Complex regulatory network allows Myriophyllum aquaticum to thrive under high-concentration ammonia toxicity. Sci Rep. 2019;9(1):4801.
You X, Chen J, Bian C, Yi Y, Ruan Z, Li J, et al. Transcriptomic evidence of adaptive tolerance to high environmental ammonia in mudskippers. Genomics. 2018;110(6):404–13.
Coschigano PW, Magasanik B. The URE2 gene product of Saccharomyces cerevisiae plays an important role in the cellular response to the nitrogen source and has homology to glutathione s-transferases. Mol Cell Biol. 1991;11(2):822–32.
Valbuena N, Guan K-L, Moreno S. The Vam6 and Gtr1-Gtr2 pathway activates TORC1 in response to amino acids in fission yeast. J Cell Sci. 2012;125(Pt 8):1920–8.
Dubouloz F, Deloche O, Wanke V, Cameroni E, De Virgilio C. The TOR and EGO protein complexes orchestrate microautophagy in yeast. Mol Cell. 2005;19(1):15–26.
Molinet J, Salinas F, Guillamon JM, Martinez C. GTR1 affects nitrogen consumption and TORC1 activity in Saccharomyces cerevisiae under fermentation conditions. Front Genet. 2020;11:519.
Eskandari R, Ratnayake L, Lakin-Thomas PL. Shared components of the FRQ-less oscillator and TOR pathway maintain rhythmicity in Neurospora. J Biol Rhythm. 2021;36(4):329–45.
Pereira MD, Eleutherio ECA, Panek AD. Acquisition of tolerance against oxidative damage in Saccharomyces cerevisiae. BMC Microbiol. 2001;1(1):11.
Novo MT, Beltran G, Rozès N, Guillamón JM, Mas A. Effect of nitrogen limitation and surplus upon trehalose metabolism in wine yeast. Appl Microbiol Biotechnol. 2005;66(5):560–6.
Hirasawa T, Nakakura Y, Yoshikawa K, Ashitani K, Nagahisa K, Furusawa C, et al. Comparative analysis of transcriptional responses to saline stress in the laboratory and brewing strains of Saccharomyces cerevisiae with DNA microarray. Appl Microbiol Biotechnol. 2006;70(3):346–57.
Nadal-Ribelles M, Conde N, Flores O, González-Vallinas J, Eyras E, Orozco M, et al. Hog1 bypasses stress-mediated down-regulation of transcription by RNA polymerase II redistribution and chromatin remodeling. Genome Biol. 2012;13(11):R106.
Li F, Wang Z-L, Zhang L-B, Ying S-H, Feng M-G. The role of three calcineurin subunits and a related transcription factor (Crz1) in conidiation, multistress tolerance and virulence in Beauveria bassiana. Appl Microbiol Biotechnol. 2015;99(2):827–40.
Montañés FM, Pascual-Ahuir A, Proft M. Repression of ergosterol biosynthesis is essential for stress resistance and is mediated by the Hog1 MAP kinase and the Mot3 and Rox1 transcription factors. Mol Microbiol. 2011;79(4):1008–23.
Gaspard GJ, McMaster CR. The mitochondrial quality control protein Yme1 is necessary to prevent defective mitophagy in a yeast model of Barth syndrome. J Biol Chem. 2015;290(14):9284–98.
Kominsky DJ, Brownson MP, Updike DL, Thorsness PE. Genetic and biochemical basis for viability of yeast lacking mitochondrial genomes. Genetics. 2002;162(4):1595–604.
Wang K, Jin M, Liu X, Klionsky DJ. Proteolytic processing of Atg32 by the mitochondrial i-AAA protease Yme1 regulates mitophagy. Autophagy. 2013;9(11):1828–36.
Pascual-Ahuir A, Proft M. The Sch9 kinase is a chromatin-associated transcriptional activator of osmostress-responsive genes. EMBO J. 2007;26(13):3098–108.
Leskoske KL, Roelants FM, Emmerstorfer-Augustin A, Augustin CM, Si EP, Hill JM, et al. Phosphorylation by the stress-activated MAPK Slt2 down-regulates the yeast TOR complex 2. Genes Dev. 2018;32(23–24):1576–90.
Mao K, Wang K, Zhao M, Xu T, Klionsky DJ. Two MAPK-signaling pathways are required for mitophagy in Saccharomyces cerevisiae. J Cell Biol. 2011;193(4):755–67.
Ding W-X, Yin X-M. Mitophagy: mechanisms, pathophysiological roles, and analysis. Biol Chem. 2012;393(7):547–64.
Krupa SV. Effects of atmospheric ammonia (NH3) on terrestrial vegetation: a review. Environ Pollut. 2003;124(2):179–221.
Britto DT, Siddiqi MY, Glass AD, Kronzucker HJ. Futile transmembrane NH4+ cycling: a cellular hypothesis to explain ammonium toxicity in plants. Proc Natl Acad Sci U S A. 2001;98(7):4255–8.
Kronzucker HJ, Britto DT, Davenport RJ, Tester M. Ammonium toxicity and the real cost of transport. Trends Plant Sci. 2001;6(8):335–7.
Tong AH, Evangelista M, Parsons AB, Xu H, Bader GD, Pagé N, et al. Systematic genetic analysis with ordered arrays of yeast deletion mutants. Science. 2001;294(5550):2364–8.
van Leeuwen J, Pons C, Mellor JC, Yamaguchi TN, Friesen H, Koschwanez J, et al. Exploring genetic suppression interactions on a global scale. Science. 2016;354(6312):aag0839.
Tong AHY, Boone C. 16 high-throughput strain construction and systematic synthetic lethal screening in Saccharomyces cerevisiae. Methods Microbiol. 2007;36:369–707.
Wagih O, Usaj M, Baryshnikova A, VanderSluis B, Kuzmin E, Costanzo M, et al. SGAtools: one-stop analysis and visualization of array-based genetic interaction screens. Nucleic Acids Res. 2013;41:W591–6.
Chen X, Ji B, Hao X, Li X, Eisele F, Nyström T, et al. FMN reduces amyloid-β toxicity in yeast by regulating redox status and cellular metabolism. Nat Commun. 2020;11(1):867.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10–2.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
We thank Charles Boone (Toronto University, Toronto, ON, Canada) for providing us with the SGA-V2 collection.
Open access funding provided by University of Gothenburg. This research was funded by grants from the National Natural Science Foundation of China (32000387) to X-LC, Scientific Research Foundation of Zhejiang A&F University (2021LFR053) to X-JJ, as well as the Swedish Cancer Fund (Cancerfonden) [19 0069] and the Swedish Research Council (Vetenskapsrådet) [VR 2019–03604] to B-DL.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Drop test of his∆ in plates with a gradient of concentrations of NH4C1.
Data from three independent genome-wide screening experiments.
GO and KEGG-enriched gene list for the 34 NH4Cl-sensitive mutants.
Growth curve of BY4741 at different concentrations of NH4C1 treatment.
Transcriptome sequencing data under 800 mM vs. 0 mM NH4Cl treatment.
About this article
Cite this article
Fu, W., Cao, X., An, T. et al. Genome-wide identification of resistance genes and transcriptome regulation in yeast to accommodate ammonium toxicity. BMC Genomics 23, 514 (2022). https://doi.org/10.1186/s12864-022-08742-y
- Saccharomyces cerevisiae
- Genome-wide screen