RNA-seq analysis of Pichia anomala reveals important mechanisms required for survival at low pH
© Fletcher et al. 2015
Received: 12 June 2015
Accepted: 31 August 2015
Published: 16 September 2015
The product yield and titers of biological processes involving the conversion of biomass to desirable chemicals can be limited by environmental stresses encountered by the microbial hosts used for the bioconversion. One of these main stresses is growth inhibition due to exposure to low pH conditions. In order to circumvent this problem, understanding the biological mechanisms involved in acid stress response and tolerance is essential. Characterisation of wild yeasts that have a natural ability to resist such harsh conditions will pave the way to understand the biological basis underlying acid stress resistance. Pichia anomala possesses a unique ability to adapt to and tolerate a number of environmental stresses particularly low pH stress giving it the advantage to outcompete other microorganisms under such conditions. However, the genetic basis of this resistance has not been previously studied.
To this end, we isolated an acid resistant strain of P. anomala, performed a gross phenotypic characterisation at low pH and also performed a whole genome and total RNA sequencing. By integrating the RNA-seq data with the genome sequencing data, we found that several genes associated with different biological processes including proton efflux, the electron transfer chain and oxidative phosphorylation were highly expressed in P. anomala cells grown in low pH media. We therefore present data supporting the notion that a high expression of proton pumps in the plasma membrane coupled with an increase in mitochondrial ATP production enables the high level of acid stress tolerance of P. anomala.
Our findings provide insight into the molecular and genetic basis of low pH tolerance in P. anomala which was previously unknown. Ultimately, this is a step towards developing non-conventional yeasts such as P. anomala for the production of industrially relevant chemicals under low pH conditions.
KeywordsPichia anomala Low pH Stress Biochemicals RNA-seq Genome sequencing ATP production
The conversion of lignocellulosic feedstocks into next-generation fuels and other industrially relevant chemicals has taken central stage in recent years. This involves several upstream and downstream processes which can be summarized as pre-treatment, saccharification, fermentation and product recovery . Microorganisms are generally used as cell factories for the biological conversion of feedstocks to the final products. In order to reduce costs and improve the efficiency of this process, it is preferable to use micro-organisms that can survive the harsh environment created by the pre-treatment step and also be tolerant to the toxicity of the final product. When using lignocellulosic materials as feedstock, an acid pre-treatment step is often required to effectively breakdown the lignin fibers trapping the carbohydrates that need to be hydrolysed before fermentation . There is also the release of weak acids and other toxic compounds in the process resulting in severe consequences on the growth and performance of the microbial host. Furthermore, in the production of biofuels and organic acids from lignocellulose or other feedstocks, the principal stress encountered by the cell is low pH stress which results in growth inhibition and low productivity. Therefore, not only must the micro-organism used be able to survive the acid pre-treatment stage, but also the low pH conditions required for the production and product recovery.
Fungi have been used in the food industry for several centuries and species such as Saccharomyces cerevisiae and Aspergillus niger have now been established as the most widely used microbial hosts for the production of biochemicals such as biofuels, organic acids, among others from biomass [3, 4]. Despite the wide use of S. cerevisiae in industrial fermentations, its use is restricted by its limited resistance to the different stress conditions encountered under an industrial setting, particularly acid stress.
Currently, strains of S. cerevisiae are being engineered to overcome the stress constraints posed by low pH environments typical for an industrial setting . As an alternative approach, however, isolation and utilization of microbes that can survive under these harsh acidic conditions can prove useful. Furthermore, studying the biological mechanisms behind the tolerance of superior yeasts will be important in revealing novel stress resistance mechanisms. Therefore, these strains can be developed as novel cell factories since they already possess the tolerance phenotype. One example is Pichia anomala (also known as Wickerhamomyces anomalus or Hansenula anomala), a yeast that is widespread in nature and has been isolated from different environments including the skin of fruits and contaminated wine . Pichia anomala is a robust yeast and possesses innate stress adaptation mechanisms making it resistant to a number of environmental stresses and it has also been described as having ‘potential biotechnology applications’ such as in bioremediation, antimicrobial and bioethanol production [6, 7]. Some studies have mentioned P. anomala’s resistance to low pH [7, 8]. A number of acid stress response mechanisms have been identified in S. cerevisiae and other fungi which include maintenance of the cell wall structure, metal metabolism and proton efflux by the membrane ATPases [9, 10].
While the advent of new sequencing technologies has revolutionized our knowledge on model organisms at the genome scale, these technologies also facilitate the study of novel interesting species, which have potential medical or industrial application. The draft genome of P. anomala was recently published and annotated to identify genes that code for proteins with antimicrobial activities . However, to date, the genetic basis of the strain’s resistance to low pH remains elusive as there is no information available on the subject.
Here, we first studied the physiological properties of a strain of P. anomala, isolated from low pH cultures. Its identity was confirmed based on the DNA sequence of its Internal Transcribed Sequence (ITS) region. We then, used whole genome sequencing and RNA-seq analysis to identify biological mechanisms that can be linked to the survival of this yeast at low pH.
Identification of P. anomala
We set up a series of shake flask cultures at low pH. Initially, the cultures were inoculated with Saccharomyces cerevisiae CEN.PK113-7D, but during the course of an evolution experiment we identified rapid appearance of five yeast strains that we then isolated from the low pH (pH 3.0) cultures. We amplified by PCR the conserved Internal Transcribed Sequence (ITS) region of the yeasts and that of S. cerevisiae. All five isolated strains produced a band of size ~650 bp. These PCR band sizes were compared to that produced by S. cerevisiae which had a size of ~840 bp (Additional file 1: Figure S1). A BLASTN (http://www.blast.ncbi.nlm.nih.gov/Blast) search of the sequenced PCR products revealed that the isolated yeast strains were Wickerhamomyces anomalus (also known as P. anomala) with a 100 % identity and coverage. Thus, we realized that this yeast clearly had a strong competitive advantage at low pH over S. cerevisiae and we therefore undertook a detailed physiological characterisation of this yeast.
Genome sequencing and assembly
P. anomala genome statistics
Number of contigs >500 bp
Total assembly size
Genome fraction (%)
# predicted genes (unique)
# predicted genes (≥300 bp)
# predicted genes ≥1500 bp)
# predicted genes (≥3000 bp)
Defining the transcriptome of Pichia anomala at low pH
Genes involved in proton efflux (H+-ATPases)
In Fig. 3b, the genes that mapped to the hydrogen-exporting ATPase activity GO terms are shown and these consist of the plasma membrane P-type H+-ATPase (FPKM = 60163.70), the alpha subunit of the F1 complex (FPKM = 23823.70), subunit C of the F0/V0 complex (FPKM = 3863.77) and subunits C and D of the V1 complex of the vacuolar H+-ATPase with FPKM values of 1448.06 and 1416.50 respectively (Additional file 8: Table S2).
Oxidative phosphorylation and the electron transport chain
Among the highly expressed genes, the following GO terms were also over-represented—oxidoreductase activity (GO:0016491), ubiquinol cytochrome-c activity (GO:0008121), electron transfer activity (GO:0009055) and succinate dehydrogenase activity (GO:0000104) alongside the terms mapped to the H+-ATPase activity (Fig. 3b). The high expression of genes associated with these GO terms reflect an active production of ATP through mitochondrial respiration in the P. anomala strain. In addition, some genes were mapped to the protein serine/threonine kinase activity GO term (GO:0004674) that might play a role in phosphorylation of specific serine and threonine residues in the H+-ATPase or associated proteins [24, 25].
We also performed a homology search to describe the P. anomala electron transport (respiratory) chain based on the S. cerevisiae oxidative phosphorylation chain. Using the KEGG database , we gathered all the protein components of S. cerevisiae electron transport chain and used them as queries to perform a blastp in JGI genome browser for P.anomala (http://www.genome.jgi.doe.gov/pages/blast-query.jsf?db=Wican1). Indeed, all the homologs of the S. cerevisiae respiratory chain were also found in P. anomala with an average sequence identity of 62.38 % (Additional file 10: Table S4). Several of the homologs were highly expressed in the P. anomala strain we isolated (Additional file 8: Table S2). It is worth noting, however, that unlike most yeasts, S. cerevisiae has several peripheral membrane NADH dehydrogenases instead of Complex I in the respiratory chain . Therefore, in order to determine the presence of Complex I in P. anomala, we used the Complex I subunits from Pichia pastoris available in the KEGG database  for a blastp search against P. anomala. Interestingly, we found 17 P. pastoris Complex I subunits present in the P. anomala genome with an average identity of 61.81 % (Additional file 10: Table S4). These subunits were not found in the S. cerevisiae respiratory chain. To determine if these subunits were expressed, we checked their FPKM values in our RNA-seq data and found that 12 of the identified subunits were highly expressed in the P. anomala strain having FPKM values in the range of 34–180. The most highly expressed subunits were the acyl carrier protein (ACP) of Complex I, NADH:ubiquinone oxidoreductase, F subunit and NDUFA5/B13 subunit with an average FPKM of 160 which confirms the active presence of Complex I in P. anomala (Additional file 8: Table S2).
Here, we studied low pH stress responses in an acid resistant yeast, P. anomala and have demonstrated how whole genome sequencing and RNA-seq can be used to characterise low pH stress responses in a non-model yeast strain. Yeasts have become a more preferred microbial cell factory for different industrial applications. Currently, S. cerevisiae which has been very well-characterised is mostly used as the main workhorse for several biological processes. However, the industrial use of this strain is challenged by its limited ability to deal with different stress conditions imposed during industrial fermentations. One of these stresses is low pH or acid stress. Apart from improving the tolerance of S. cerevisiae to stress, an alternative approach will be to isolate and characterise novel yeast strains that have a natural resistance to stresses such as acid stress  in order to identify novel tolerance mechanisms in these strains.
One well-established mechanism by which cells deal with acid stress resulting from intracellular accumulation of protons is by pumping the excess protons out of the cytoplasm into the extracellular environment or into the vacuole . In yeasts and other micro-organisms, this process is mainly carried out by H+-ATPases (proton pumps). Several components of ATPases were highly expressed with the membrane H+-ATPase being the most highly expressed gene as shown in the transcriptomic data presented in this work. Plasma membrane H+-ATPases belong to a large family of molecular pumps involved in ion transport across the membrane  and have already been shown to be activated in yeast when the external environment is acidified . The H+-ATPases pump out excess protons which diffuse into the cytoplasm in an ATP dependent process thus maintaining the internal pH of the cell and has been shown to be activated during acid stress showing its relevance in pH control . Our data confirm previous work by Fredlund et al.  who first reported the unique low pH resistance phenotype of P. anomala. To the best of our knowledge, the P. anomala H+-ATPase has not been previously reported to play a role in acid stress tolerance since prior to this present study there has not been any genetic and transcriptomic study on acid stress responses displayed by P. anomala. Subunits of the vacuolar H+-ATPase (V-ATPase) were also highly expressed even though the level of expression as measured by the FPKM values was not as high as the plasma membrane H+-ATPase. The V-ATPase works in tandem with the plasma membrane H+-ATPase to de-acidify the cytosol by pumping excess protons from the cytosol into the vacuole . Therefore, the expression of P. anomala V-ATPase in this study shows its contribution to increasing low pH tolerance in the yeast. The high expression of H+-ATPases provides an obvious explanation to acid stress tolerance. However, a closer look at the structure of the P. anomala plasma membrane H+-ATPase shows a high level of homology to the S. cerevisiae plasma membrane H+-ATPase (Pma1p) even though P. anomala is more resistant to low pH than S. cerevisiae. Therefore, increased expression of the P. anomala H+-ATPases alone cannot explain the basis of its high resistance to low pH.
Expulsion of protons by H+-ATPases is an energy intensive process which drains the cell of most of its intracellular ATP . Interestingly, we also observed a high expression of succinate dehydrogenase and several oxidoreductases in our work, which was an indication of an active energy generation mechanism by the cell through the respiratory chain and may provide a further explanation to the low pH tolerance of P. anomala. Succinate dehydrogenase is a key component and a functional member of the electron transport chain  which is closely linked to the oxidative phosphorylation pathway where the oxidoreductases mediate the conversion of mitochondrial ADP to ATP . In effect, ATP production by the mitochondria will be increased by the electron transport chain and the ATP produced is supplied to the cell to meet the ATP requirements to exclude protons out of the cell by the plasma membrane H+-ATPase. There is usually a high energy demand by microbial cells under environmental stress  and it has also been mentioned by de Kok et al.  that the H+/ATP stoichiometry of the plasma membrane H+-ATPase determines the ATP requirement of the cell for maintenance and H+ ion homeostasis. Therefore, increasing ATP supply within the cell is beneficial in alleviating the stress and this has indeed been reported by Zhou et al.  where acid tolerance of Candida glabrata was improved by supplying the cell with excess amounts of ATP. The structure of the respiratory chain in P. anomala in comparison with that of S. cerevisiae further highlights the ability of P. anomala to efficiently generate energy from respiration. Although core proteins were conserved in the respiratory chain of both yeasts, S. cerevisiae lacks most subunits of Complex I as a result of an evolutionary adaptation to fermentative anaerobic growth . Conversely, the presence of Complex I in P. anomala will enable the yeast to produce substantial amounts of energy predominantly through the electron transport chain since this Complex I generates 40 % of the transmembrane electrochemical gradient needed to synthesise ATP from ADP in the mitochondria .
In the presence of high extracellular amounts of glucose, S. cerevisiae, ferments the substrate to yield ethanol. Even though ATP is produced in this pathway, the yield is not high enough to meet the energy demand of the cell under acid stress unlike in P. anomala, which respires actively to generate much more ATP as shown in our study. Similarly, Chen et al.  observed a down-regulation of genes involved in the TCA cycle and in respiratory metabolism during acid stress in S. cerevisiae even though pyruvate decarboxylase activity was increased signifying a direction of carbon flux toward ethanol production while respiration was turned off. Furthermore, other reports characterising the effect of low pH on S. cerevisiae at both genomic and transcriptomic levels do not show an activation of mitochondrial ATP production. Instead, other mechanisms including activation of the protein kinase pathway which plays a role in maintaining cell shape and integrity  were identified. Also, upregulation of metal metabolism has been shown to be important in acid stress resistance in S. cerevisiae . Under conditions of acid stress, priority is given to dealing with proton accumulation in the cell by extruding the excess protons into the extracellular environment through the plasma membrane ATPase. Data presented in this study show that even though both S. cerevisiae and P. anomala consume glucose at the same rate under low pH conditions, S. cerevisiae has a lower biomass yield on glucose but a higher ethanol and glycerol production rate compared to P. anomala. This implies that P. anomala, through respiration resulting from its Crabtree negative nature, is able to generate a larger pool of ATP to deal with the stress and for maintenance (evidenced from the high biomass yield) whereas S. cerevisiae uses its limited supply of ATP (from glycolysis and fermentation) to provide the chemical energy needed by the plasma membrane ATPase to transport the excess protons out of the cell at the expense of cell maintenance reflected in a low biomass yield.
In summary, although the plasma membrane H+-ATPase of both P. anomala and S. cerevisiae show a high level of homology and both strains consume glucose at the same rate during acid stress, P. anomala proved to be more resistant to low pH than S. cerevisiae. Our findings suggest that the mere presence of a proton pump does not necessarily lead to tolerance to acid stress but an improved supply of chemical energy (ATP, in this case) is needed to power the proton pump. Yet, other factors that can remarkably influence the activity of Pma1p include the lipid composition of the membrane  and the regulatory mechanism underlying Pma1p expression. Thus, in vivo validation of Pma1p activity is needed to back transcriptional data. While this study characterises tolerance of P. anomala to inorganic acid stress, it will still be interesting to investigate transcriptional responses of this yeast to stress posed by different organic acids since tolerance to low pH does not necessarily imply tolerance to organic acids at low pH particularly because apart from the stress stemming from the low pH environment, accumulation of the anion form of an organic acid in the cytoplasm will also be toxic to the cell.
Isolation and identification of the yeast
Five strains of yeasts were isolated from different low pH (pH 3.0) cultures and were identified as P. anomala by a PCR amplification of the Internal Transcribed Sequence (ITS) region of the yeast genome which includes the conserved sequences of 5.8S and 28S ribosomal DNA . The ITS 1 (TCCGTAGGTGAACCTGCGG) and ITS 4 (TCCTCCGCTTATTGATATGC) primers  were used for the PCR amplification. The yeast genomic DNA was extracted using the LiOAc-SDS method as described by Looke et al. . The PCR products were sequenced by Eurofins Genomics (Ebersberg, Germany) and the species were identified through a BLASTN search (http://www.blast.ncbi.nlm.nih.gov/Blast.cgi).
Growth rate determination
Pure cultures of the yeast strains were obtained by isolating individual colonies on solid YPD (10 g/L yeast extract, 20 g/L peptone, 22 g/L glucose) media plates. Fresh minimal medium [7.5 g/L(NH4)2SO4, 14.4 g/L KH2PO4, 0.5 g/L MgSO4·7H2O, 20 g/L glucose, 1 mL/L vitamin solution and 1 mL/L trace elements solution) with pH adjusted to 3.0 with 1 M HCl was inoculated with overnight cultures of the strains to an initial OD600 of 0.05. The cultures were incubated with shaking (200 rpm) at 30 °C and the cell density (OD600) of the cultures were measured every 2 h with a Thermo Scientific Genesys20 spectrophotometer (Waltham, USA). In a parallel experiment, the growth rate of S. cerevisiae CEN.PK 113-7D was obtained for comparison with the P. anomala strains under the same conditions. The growth rate was determined from the slope obtained from a semi-logarithmic plot of OD600 against time during the exponential phase of the cultures. The top-performing P. anomala strain out of the five isolated strains was used for the subsequent experiments and deposited in our strain collection as P. anomala EF01.
Quantification of glucose and fermentation metabolites
Aliquots of the supernatants of the S. cerevisiae cell cultures and that of the top-performing P. anomala were taken every 2 h, kept at 4 °C and were used to determine the concentration of substrates and other fermentation metabolites by HPLC analysis. Using an RI and UV detectors, an Aminex HPX-87H column (Biorad, Irvine, CA, USA) was used to separate and quantify glucose, ethanol and glycerol in the supernatants. The column was eluted with 5 mM H2SO4 eluent at a flow rate of 0.6 mL/min.
Sample and library preparation
Genomic DNA was extracted from top-performing P. anomala strain using the Qiagen Gentra Puregene kit (Hilden, Germany). To obtain the total RNA of the P. anomala strain, 10 mL cultures were grown for ~8 h in minimal medium (pH 3.0) until they reached the mid-log growth phase (OD = 1). The total RNA was then isolated from these cultures in three biological replicates using the Qiagen RNAeasy Mini Kit (Hilden, Germany).
Sequencing and data preprocessing
The Illumina sequencing platform (HiSeq 2000, 100-bp paired-end) at the Theragen Bio Institute, South Korea was used to sequence the genomic DNA libraries while the Multiplex Illumina HiSeq High Output mode (2 × 100 bp paired-end) was used to sequence the RNA libraries. The sequences were pre-processed by Trimmomatic  in order to remove adaptors, leading low quality or N bases (below quality 3), quality trimming (4-base wide sliding window and quality cut off of 15), and removing reads below 36 bases long.
De novo assembly and gene catalog assessment
For the de novo assembly we first ran Velvet  using VelvetOptimizer with the following parameters: –exp_cov auto (automatic calculation of expected coverage), –scaffolding (scaffolding of contigs with paired-end reads) and –min_- contig_lgth 200 (mimimun contig length = 200). The optimal k-mer length was determined by adjusting the k-mer length from 31 to 75 bp in 4-bp increments and using the k-mer for which the N50 and the maximum contig length reached the highest value (Additional file 3: Table S1). The resulting contigs were then re-assembled with CAP3 v10/15/07 using standard parameters. To compare, we also ran ABySS  and SPAdes . For the ABySS we used the optimized suggested k-mer size of 37 based on work done by Haridas et al.  while keeping the other parameters as default. For the SPAdes we used the k-size range between k21 and k65 with the increment size of 4. To check if the yeast was already sequenced we blasted the longest scaffolds against JGI fungi sequenced genome databases. All the scaffolds were matched with average identity score of 98 % to the Wickerhamomyces anomalus NRRL Y-366-8, also known as P. anomala and Hansenula anomala in the JGI with the genome assembly size of 14.15 Mbp and 6423 gene models. For the rest of analysis we used W. anomalus NRRL Y-366-8 as our genome of reference. To compare the quality of our de novo assemblies we used the QUAST program  with W. anomalus as the reference genome and chose –gene-finder and –eukaryote parameters to predict the genes.
Among the assemblers that we used ABySS output had the highest genome coverage (>80 %) and the coverage of the Velvet-assembled genome was less than 50 %.
For the RNA-seq analysis we used the Tuxedo pipeline (Additional file 11: Figure S7) . However, we used Cuffnorm instead of Cuffdiff to count the reads (new extension of Cufflinks) as the former normalized the read counts to the sample librAdditionalry size.
Go enrichment analysis
As P. anomala is not supported by an annotation package for GO enrichment analysis of the each of the four gene categories we used the GOstats and Category R packages which are specifically designed for hypergeometric testing with unsupported model organisms . Therefore, we created a custom data structure of the GO ontologies and then each selected gene category was used as input to perform the hypergeometric test. The output of the GOstats was visualized using REVIGO .
Sequence and structural bioinformatics analysis
The orthologs of the P. anomala plasma membrane H+-ATPase were detected using the BLASTp program . Jalview was used for the pairwise alignment and visualization . The secondary structures of the S. cerevisiae and P. anomala H+-ATPases were predicted as described in [48, 49]. For the prediction and alignment of the structures we used the RaptorX server .
The experiments were designed by JN, VS and EF. EF carried out the experiments. SK extracted the genomic DNA, performed all the quality checks and had them sequenced. The bioinformatics analyses and interpretation of the data were carried out by AF. AF prepared and assembled the figures and tables. The manuscript was written by AF and EF. All authors read and approved the final manuscript.
The authors would like to acknowledge Samsung and the Novo Nordisk Foundation for funding this project. We would also like to thank the Science for Life Laboratory (SciLifeLabs), Sweden and Theragen Bio Institute, South Korea for providing assistance with the RNA-seq and the whole genome sequencing respectively.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- van Zyl WH, Lynd LR, den Haan R, McBride JE. Consolidated bioprocessing for bioethanol production using Saccharomyces cerevisiae. Adv Biochem Eng Biotechnol. 2007;108:205–35.Google Scholar
- Jung YH, Kim IJ, Kim HK, Kim KH. Dilute acid pretreatment of lignocellulose for whole slurry ethanol fermentation. Bioresour Technol. 2013;132:109–14.View ArticleGoogle Scholar
- Schrickx JM, Raedts MJ, Stouthamer H, van Verseveld HW. Organic acid production by Aspergillus niger in recycling culture analyzed by capillary electrophoresis. Anal Biochem. 1995;231:175–81.View ArticleGoogle Scholar
- Zhou YJ, Buijs N, Siewers V, Nielsen J. Fatty acid-derived biofuels and chemicals production in Saccharomyces cerevisiae. Front Bioeng Biotechnol. 2014;2(September):1–6.Google Scholar
- Tanaka K, Ishii Y, Ogawa J, Shima J. Enhancement of acetic acid tolerance in Saccharomyces cerevisiae by overexpression of the Haa1 gene, encoding a transcriptional activator. Appl Environ Microbiol. 2012;78:8161–3.View ArticleGoogle Scholar
- Walker GM. Pichia anomala: cell physiology and biotechnology relative to other yeasts. Antonie van Leeuwenhoek Int J Gen Mol Microbiol. 2011;99:25–34.View ArticleGoogle Scholar
- Kostas ET, White D, Du C, Cook DJ. Selection of yeast strains for bioethanol production from UK seaweeds. J Appl Phycol 2015.Google Scholar
- Fredlund E, Druvefors U, Boysen M, Lingsten K, Schnurer J. Physiological characteristics of the biocontrol yeast J121. FEMS Yeast Res. 2002;2:395–402.Google Scholar
- Brandao RL, Rosa JCC, Nicoli JR, Almeida MVS, do Carmo AP, Queiros HT, Castro IM. Investigating acid stress response in different Saccharomyces strains. J Mycol. 2014;2014:9.Google Scholar
- Kawahata M, Masaki K, Fujii T, Iefuji H. Yeast genes involved in response to lactic acid and acetic acid: acidic conditions caused by the organic acids in Saccharomyces cerevisiae cultures induce expression of intracellular metal metabolism genes regulated by Aft1p. FEMS Yeast Res. 2006;6:924–36.View ArticleGoogle Scholar
- Schneider J, Rupp O, Trost E, Jaenicke S, Passoth V, Goesmann A, Tauch A, Brinkrolf K. Genome sequence of Wickerhamomyces anomalus DSM 6766 reveals genetic basis of biotechnologically important antimicrobial activities. FEMS Yeast Res. 2012;12:382–6.View ArticleGoogle Scholar
- Zha Y, Hossain AH, Tobola F, Sedee N, Havekes M, Punt PJ. Pichia anomala 29X: a resistant strain for lignocellulosic biomass hydrolysate fermentation. FEMS Yeast Res. 2013;13:609–17.View ArticleGoogle Scholar
- Fredlund E, Blank LM, Schnürer J, Sauer U, Passoth V. Oxygen- and glucose-dependent regulation of central carbon metabolism in Pichia anomala. Appl Environ Microbiol. 2004;70:5905–11.View ArticleGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol I. ABySS : a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.View ArticleGoogle Scholar
- Bankevich A, Nurk S, Antipov D, Gurevich A, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin A, Sirotkin A, Vyahhi N, Tesler G, Alekseyev M, Pevzner P. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.View ArticleGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.View ArticleGoogle Scholar
- Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29:1072–5.View ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7:562–78.View ArticleGoogle Scholar
- Supek F, Bošnjak M, Škunca N, Šmuc T. Revigo summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.View ArticleGoogle Scholar
- Fernandes AR, Sá-Correia I. The activity of plasma membrane H+-ATPase is strongly stimulated during Saccharomyces cerevisiae adaptation to growth under high copper stress, accompanying intracellular acidification. Yeast. 2001;18:511–2.View ArticleGoogle Scholar
- Viegas C, Almeida PF, Cavaco M, Sá-Correia I. The H+-ATPase in the plasma membrane of Saccharomyces cerevisiae is activated during growth latency in octanoic acid-supplemented medium accompanying the decrease in intracellular pH and cell viability. Appl Environ Microbiol. 1998;64:779–83.Google Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.View ArticleGoogle Scholar
- Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2-A multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25:1189–91.View ArticleGoogle Scholar
- Portillo F, Eraso P, Serrano R. Analysis of the regulatory domain of yeast plasma membrane H+-ATPase by directed mutagenesis and intragenic suppression. FEBS Lett. 1991;287:71–4.View ArticleGoogle Scholar
- Goossens A, de La Fuente N, Forment J, Serrano R, Portillo F. Regulation of yeast H(+)-ATPase by protein kinases belonging to a family dedicated to activation of plasma membrane transporters. Mol Cell Biol. 2000;20:7654–61.View ArticleGoogle Scholar
- Okuda S, Yamada T, Hamajima M, Itoh M, Katayama T, Bork P, Goto S, Kanehisa M. KEGG Atlas mapping for global analysis of metabolic pathways. Nucleic Acids Res. 2008;36((Web Server issue)):423–6.View ArticleGoogle Scholar
- Yamashita T, Nakamaru-Ogiso E, Miyoshi H, Matsuno-Yagi A, Yagi T. Roles of bound quinone in the single subunit NADH-quinone oxidoreductase (Ndi1) from Saccharomyces cerevisiae. J Biol Chem. 2007;282:6012–20.View ArticleGoogle Scholar
- Ruyters S, Mukherjee V, Verstrepen KJ, Thevelein JM, Willems KA, Lievens B. Assessing the potential of wild yeasts for bioethanol production. J Ind Microbiol Biotechnol. 2014;42:39–48.View ArticleGoogle Scholar
- Piper P, Calderon CO, Hatzixanthis K, Mollapour M. Weak acid adaptation: the stress response that confers yeasts with resistance to organic acid food preservatives. Microbiology. 2001;147:2635–42.View ArticleGoogle Scholar
- Morth JP, Pedersen BP, Buch-Pedersen MJ, Andersen JP, Vilsen B, Palmgren MG, Nissen P. A structural overview of the plasma membrane Na+, K+-ATPase and H+-ATPase ion pumps. Nat Rev Mol Cell Biol. 2011;12:60–70.View ArticleGoogle Scholar
- Eraso P, Gancedo C. Activation of yeast plasma membrane ATPase by acid pH during growth. FEBS Lett. 1987;224:187–92.View ArticleGoogle Scholar
- Martínez-Muñoz GA, Kane P. Vacuolar and plasma membrane proton pumps collaborate to achieve cytosolic pH homeostasis in yeast. J Biol Chem. 2008;283:20309–19.View ArticleGoogle Scholar
- Ullah A, Chandrasekaran G, Brul S, Smits GJ. Yeast adaptation to weak acids prevents futile energy expenditure. Front Microbiol. 2013;4:142.View ArticleGoogle Scholar
- Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Electron-transport chains and their proton pumps. In: Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P, editors. Molecular biology of the cell, vol. 4. New York: Garland Science; 2002.Google Scholar
- Zhou J, Liu L, Chen J. Improved ATP supply enhances acid tolerance of Candida glabrata during pyruvic acid production. J Appl Microbiol. 2011;110:44–53.View ArticleGoogle Scholar
- de Kok S, Yilmaz D, Daran JM, Pronk JT, van Maris AJA. In vivo analysis of Saccharomyces cerevisiae plasma membrane ATPase Pma1p isoforms with increased in vitro H+/ATP stoichiometry. Int J Gen Mol Microbiol. 2012;102:401–6.Google Scholar
- Marcet-Houben M, Marceddu G, Gabaldón T. Phylogenomics of the oxidative phosphorylation in fungi reveals extensive gene duplication followed by functional divergence. BMC Evol Biol. 2009;9:295.View ArticleGoogle Scholar
- Papa S, Petruzzella V, Scacco S. Electron transport, structure, redox-coupled protonmotive activity and pathological disorders of respiratory chain complexes. In: Lajtha A, Dienel G, Gibson G, editors. Handbook of neurochemistry and molecular neurobiology. Berlin Heidelberg: Springer; 2007. p. 93–118.View ArticleGoogle Scholar
- Chen AKL, Gelling C, Rogers PL, Dawes IW, Rosche B. Response of Saccharomyces cerevisiae to stress-free acidification. J Microbiol. 2009;47:1–8.View ArticleGoogle Scholar
- Claret S, Gatti X, Thoraval D, Crouzet M. The Rgd1p Rho GTPase-activating protein and the Mid2p cell wall sensor are required at low pH for protein kinase C pathway activation and cell survival in Saccharomyces cerevisiae. Eukaryot Cell. 2005;4:1375–86.View ArticleGoogle Scholar
- Aguilera F, Peinado R, Millán C, Ortega JM, Mauricio JC. Relationship between ethanol tolerance, H+-ATPase activity and the lipid composition of the plasma membrane in different wine yeast strains. Int J Food Microbiol. 2006;110:34–42.View ArticleGoogle Scholar
- Fujita S, Senda Y, Nakaguchi S, Hashimoto T. Multiplex PCR using internal transcribed spacer 1 and 2 regions for rapid detection and identification of yeast strains. J Clin Microbiol. 2001;39:3617–22.View ArticleGoogle Scholar
- Löoke M, Kristjuahan K, Kristjuhan A. Extraction of genomic DNA from yeasts for PCR-based applications. Biotechniques. 2011;50:325–8.Google Scholar
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.View ArticleGoogle Scholar
- Haridas S, Breuill C, Bohlmann J, Hsiang T. A biologist’s guide to de novo genome assembly using next-generation sequence data: a test with fungal genomes. J Microbiol Methods. 2011;86:368–75.View ArticleGoogle Scholar
- Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–8.View ArticleGoogle Scholar
- Altschup SF, Science C, Pennsylvania T, University S, Park U. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.View ArticleGoogle Scholar
- Cole C, Barber JD, Barton GJ. The Jpred 3 secondary structure prediction server. Nucleic Acids Res. 2008;36:W197–201.View ArticleGoogle Scholar
- Cuff J, Barton GJ. Application of multiple sequence alignment profiles to improve protein secondary structure prediction. Proteins. 2000;40:502–11.View ArticleGoogle Scholar
- Peng J, Xu J. RaptorX: exploiting structure information for protein alignment by statistical inference. Proteins. 2011;79(Suppl 10):161–71.View ArticleGoogle Scholar