Introduction

Brazil has been at the frontline of the bioethanol industry since the 1970s and the sugarcane-based production process employed in Brazil constitute one of the most cost-efficient ethanologenic processes (Walker 2010). Having produced about 27 billion liters of ethanol in 2015, Brazil is the second largest ethanol producer in the world, only exceeded by the USA (Renewable Fuels Association 2016). The production process employed in Brazil is unlike other ethanol processes and includes yeast recycling, acid wash of the yeast biomass, and non-aseptic conditions, which altogether result in a continuously changing and highly stressful environment (Della-Bianca et al. 2013). Despite this, the Brazilian bioethanol process shows a high efficiency with yields up to 93% of theoretic stoichiometric conversion (Amorim et al. 2011).

Foam formation is considered one of the main drawbacks of industrial fermentation processes, especially on large-scale processes, and Brazilian ethanol fermentations are no exception. Foam consists of gas bubble dispersions on liquid, solid/liquid, and solid systems (Eisner et al. 2007). A number of features greatly influence its occurrence in bioreactors, including medium composition, gas introduction/formation, and strain-specific characteristics of the production organism (Junker 2007). During ethanol fermentation by Saccharomyces cerevisiae strains, one of the main factors affecting foaming is the yeast itself. Foam stability is strongly influenced by the proteins adsorbing the gas/liquid interface (Prins and Van Riet 1987), where mannoproteins from the yeast cell wall have shown to increase the foam stability to a great extent (Blasco et al. 2011). In the Brazilian ethanol setup, large volumes of biomass are used (about 10% w v−1), in a cell recycling system, and screening strategies have focused in selecting yeast with non-foaming phenotypes (Basso et al. 2008). However, since the process is non-sterile, the appearance of wild yeast strains with foaming phenotypes occurs regularly as contaminants resulting in increased foam formation. Having a yeast strain with a foaming phenotype is highly deleterious for the process economics, since a larger volume of the vessel is taken up by the foam, and it increases overall production costs, due to consumption of antifoam agents (AFA). This is especially important for production of low-value products such as ethanol (Basso et al. 2008).

Different chemical antifoams have different impacts on microbial physiology and cell growth (Routledge and Bill 2012), and the choice of the chemical used should take this into consideration. Several studies have considered the impact of AFA on cell physiology mainly as a result of the physicochemical effects on medium, such as reduction of the liquid surface tension, and decreased oxygen transfer rate (Morao et al. 1999; Garcia-Ochoa and Gomez 2009; Routledge et al. 2011). Still, some studies demonstrate a more direct impact on microbial physiology by the addition of AFA. One widely used AFA based on polypropylene glycol was shown to have toxic effects on Bacillus thuringiensis during sewage sludge fermentation for biopesticide production, by affecting the transport of nutrients and oxygen through the cell walls (Vidyarthi et al. 2000). Another AFA based on silicone polymers was shown to stimulate glycerol production, at the expense of ethanol, in S. cerevisiae steady-state chemostat cultures, under conditions of low oxygen and excess glucose (Grosz and Stephanopoulos 1990). These studies demonstrate that selection of AFA for a particular process, should take into consideration the physiological impact it might have on the microbial cell culture in a context dependent manner, rather than only the cost of the chemical.

This study has investigated the impact of AFA on the physiology and transcriptome of the industrial ethanol S. cerevisiae strain CAT-1. We show that while commercially available AFA used in the yeast fermentation literature are innocuous at concentrations relevant for defoaming, industrial AFA used at Brazilian ethanol plants have severe negative impacts on the production performance of S. cerevisiae. Effects of AFA treatments were investigated in both defined and in molasses-based media, and transcriptome analysis suggested that alterations in lipid metabolism, as a consequence of AFA exposure, might account for a decreased stress tolerance in Brazilian ethanol fermentations.

Materials and methods

Strain and media

The industrial S. cerevisiae strain CAT-1 (CBMAI 0957, Coleção Brasileira de Microrganismos de Ambiente e Indústria) was used for all fermentation assays in this study. Single colonies were picked from plates and used for inoculation of precultures in yeast extract peptone dextrose (YPD) (containing 10 g l−1 yeast extract, 20 g l−1 peptone, and 20 g l−1 glucose) or yeast nitrogen base (YNB) (BD Difco™ cat. no. 291940) media and incubated at 32 °C and 200 rpm overnight. Sugarcane must was prepared by diluting sugarcane molasses in tap water to 20 °Brix, followed by centrifugation to remove solids. All media were autoclaved for sterilization. Water-soluble AFA were diluted in sterile demineralized water while water-insoluble AFA were diluted in 99% ethanol and added to the fermentation media. Ethanol was added to appropriate media to compensate for the ethanol contribution from ethanol-dissolved AFA. Industrial AFA samples (denoted as Ind_B and Ind_Z) are composed of two products, and they were used in all cultivations as a mixture of an antifoam product (AF_B and AF_Z, respectively) and a dispersant product (D_B and D_Z, respectively) at 1:1 proportion (Table 1).

Table 1 Commercial and industrial antifoam agents (AFA) used in this study

Batch fermentation

All fermentations were inoculated from a preculture to an initial OD600 of 0.1 and incubated at 32 °C. Small-scale fermentations were conducted in 96-well Costar flat-bottom microtiter plates with a working volume of 100 μl YPD medium and linear shaking on “fast” using a Biotek Synergy Mx plate reader (Winooski, VT, USA). Growth was monitored by measuring OD600 every 15 min. Shake flask fermentations were conducted in 500-ml Erlenmeyer flasks in YNB medium, with 200 rpm orbital shaking and samples for biomass and metabolite quantification were collected every 1.5 h. All fermentations were conducted in biological triplicates.

Yeast propagation and fed-batch fermentation with cell recycling

To propagate sufficient yeast biomass for inoculation of molasses fermentations, a YPD preculture was incubated overnight at 32 °C without shaking. YPD cultures were subsequently fed in discrete portions with propagation medium (sugarcane must diluted to 10 °Brix and enriched with yeast extract (5 g l−1) and peptone (10 g l−1)) to a final volume of 2 l during the course of 3 days.

Fed-batch fermentations representing the typical industrial process were carried out in 50-ml falcon tubes. In the first cycle, cells from the propagation culture were added to each tube in an amount corresponding to 8% (w v−1) of the final volume. Cells were fed with 25 ml sugarcane must in three equal-sized portions with 1.5 h intervals. Cultures were incubated for 7 h at 32 °C without agitation and left at room temperature overnight. In the following day, cells were separated from the fermentation wine by centrifugation (3220 rcf, 5 min) and resuspended in wine (30% wet weight w w−1) to simulate the industrial centrifuge efficiency. Cells were further diluted in demineralized water (1:1) before addition of 1 M sulphuric acid to a final pH of 2.5. After incubation in acid at room temperature for 1 h, feeding of sugarcane must was initiated, restarting the process. Industrial AFA were administered as in the industry, with the dispersant added during the acid wash and antifoam after 1 h of fermentation.

Metabolite concentrations in the fermentation wine were determined with high-performance liquid chromatography (HPLC) and biomass as wet weight. Moreover, CO2 formation was followed during the fermentation by measuring the weight loss of the fermentation tubes and assuming all loss of weight could be attributed to CO2 production. Yields were calculated on a C-mole basis relative to total reducing sugars, i.e., sucrose, glucose, and fructose. The wine returned, and hence the addition of extracellular metabolites from the previous fermentation, was taken into account and corrected for in the yield calculations.

Metabolite quantification

Metabolite concentrations were determined by high-performance liquid chromatography (HLPC). Sugars (sucrose, glucose, and fructose) were separated on an Aminex HPX-87H ion-exclusion column (Hercules, CA, USA), isocratically eluted at 25 °C with 5 mM H2SO4 as the mobile phase at a flow rate of 0.6 ml min−1. Glycerol, ethanol, and organic acids were determined similarly, but with a higher flow rate (0.7 ml min−1) and temperature (65 °C). Detection was performed refractrometrically. For the sugarcane fermentations, ethanol was measured by distillation of 10 ml fermentation wine and subsequent density determination using an Anton Paar electronic density meter (Graz, Austria) (Basso et al. 2008).

Determination of cell viability

Samples from sugarcane fermentations were diluted 1000 times in demineralized water to a final concentration of approximately 500 cells ml−1. Cells were dyed with propidium iodide (PI) (215 nmol ml−1) and incubated for 5–10 min in the dark before applying the samlpes on a flow cytometer Accuri C6 (BD, Franklin Lakes, NJ, USA) according to the manufacturer’s recommendations. Threshold of FSC-H gate was set to 200,000 to avoid large particles. Quantification of viability was done by separating two populations generated from the histogram of the PI fluorophore detector FL3-H.

RNA extraction, sequencing, and quantification

Samples for transcriptome analysis were collected from cycle 5 of the sugarcane fermentations, after 4 h of incubation, spun down, resuspended in RNAlater (Qiagen, Hilden, Germany, cat. no. 76106) and incubated for 1 h, at 4 °C. Subsequently, cells were pelleted and stored at −80 °C until further analysis. Biomass samples were lysed on a FastPrep instrument (Lysing Matrix Y from MP Biomedicals, Hampton, NH, USA, cat. no. 116960050) and total RNA extracted according to the manufacturer’s recommendations (TRIzol Plus RNA purification kit, Invitrogen, Carlsbad, CA, USA, cat. no. 12183-555). cDNA libraries were prepared using Illumina TruSeq kit (Illumina, San Diego, CA, USA, cat. no. RS-122-2001) with enzymatic rRNA depletion. Sequencing was done using Illumina HiSeq 2000 technology yielding an average of 26 M, 99-nt long paired-ends reads per sample, with an average insert distance of 300 nt’s. Reads in fastq format were quality assessed using FastQC (v0.10.1) (Andrews 2012) and trimmed according to varying GC-content at the 5′ and 3′ ends, and a quality cut-off (Q < 20). Trimmed sequence reads were mapped to the S. cerevisiae s288C reference genome (version EF4, Ensembl) (Engel et al. 2014) using Tophat (v2.0.9) (Kim et al. 2013), and transcripts were quantified on a count basis using HTSeq (v0.6.0) (Anders et al. 2015) with “intersection-strict” mode for overlapping reads. Data normalization and calculation of gene level statistics were conducted using the R package DESeq2 (Love et al. 2014), and significantly differentially expressed genes were identified as having Benjamini-Hochberg-corrected p values less than 0.01. Transcripts with overall counts less than 15 among all samples were removed for the further analysis. The transcriptome data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through the GEO Series accession number GSE103004.

Gene-set analysis

Gene sets with amiGO gene ontology (GO) terms were downloaded from Ensembl (Yates et al. 2016), and gene-metabolite associations for reporter metabolite analysis were extracted from the yeast 7 genome-scale metabolic model (Aung et al. 2013). Gene set analysis was performed with the normalized transcript counts using the R package Piano (Varemo et al. 2013), and gene sets with Benjamini-Hochberg corrected p-values less than 0.01 were considered significant.

Linear models

Linear models were applied to test whether mRNA levels were a linear function of cell viability. Modeling was performed using the R implementation of an ANOVA F test with p < 0.01 as criterium to reject the null-hypothesis, i.e., that gene expression and viability were uncorrelated. Enrichment of GO terms among genes with significant positive or negative correlations was tested with all expressed transcripts as background, using the hypergeometric test implemented in amiGO (Carbon et al. 2009).

Results

The differential effects of commercial and industrial antifoam agents on yeast physiology under defined laboratory conditions

In this study, a total of seven AFA were utilized and a description of each is provided in Table 1. Three commercial AFA were selected as commonly used in the fermentation literature for yeast physiological studies (Brochado et al. 2010; Routledge et al. 2011; Basso et al. 2011), while four industrial AFA were acquired as samples from two Brazilian ethanol plants and are currently being used in the Brazilian ethanol industry. As two AFA are used in conjunction to control foam in the industrial process (Basso et al. 2008; Della-Bianca et al. 2013), the four industrial AFA represent two combination treatments for defoaming (Ind_B and Ind_Z). These combination treatments consist of one dispersant (D_B or D_Z) which is added into the yeast slurry prior feeding, to prevent foam formation, and an antifoam (AF_B or AF_Z) which is added when foam formation occurs during fermentation (Supplementary Fig. S1).

A physiological characterization of the industrial yeast S. cerevisiae strain CAT-1 was performed initially in a microplate reader setup, in order to assess the effect of the AFA on growth rate at different concentrations. Secondly, a more detailed investigation was done using shake-flask batch cultivations in defined medium (YNB) to study the impact of industrial AFA on major physiological parameters, such as product yields.

Microplate cultivations in YPD medium with 10 g l−1 glucose were used to estimate the effect of different concentrations of AFA on the maximum specific growth rate (μ max). In the case of commercial AFA, concentrations up to 100 times the manufacturer’s recommendation were tested. The presence of Antifoam C and P2000 decreased growth rates as compared to the control (no AFA added) in concentrations of 750 and 7500 mg l−1, respectively, while Antifoam 204 displayed no effect on growth rate up to 375 mg l−1 (Fig. 1a) (higher concentrations of Antifoam 204 did not emulsify and hence were not tested). Assessing the industrial AFA revealed that both treatments start to reduce growth rate at 60 mg l−1 (Fig. 1b), which corresponds to the average concentration used in the industrial setting (Pecege-Esalq 2012). To further investigate how metabolite production and growth characteristics were influenced by the two industrial AFA, shake-flask cultivations were carried out in YNB containing 10 g l−1 glucose. Under these conditions, while yield coefficients virtually remained unchanged by the presence of industrial AFA, specific metabolite production and consumption rates were negatively affected in view of the reduced growth rates (Fig. 1c).

Fig. 1
figure 1

The effect of commercial and industrial AFA on growth rate and product yields of S. cerevisiae CAT-1 cultures. a Maximum specific growth rate (h−1) in YPD medium in microplate reader cultures in the presence of different concentrations of three commercial AFA samples. b Maximum specific growth rate (h−1) in YPD medium in microplate reader cultures in the presence of different concentrations of two industrial AFA samples (products Ind_B and Ind_Z are composed of an antifoam and a dispersant, at equal amounts). c Product yields during growth of S. cerevisiae CAT-1 in shake-flask cultures on YNB with 10 g l−1 glucose as sole carbon source, in the presence of AFA Ind_B (30 mg l−1 AF_B and 30 mg l−1 D_B) and Ind_Z (30 mg l−1 AF_Z and 30 mg l−1 D_Z). Yield on biomass (Y sx), on ethanol (Y se), on glycerol (Y sg) and on acetate (Y sa). Values represent the average of triplicate experiments ± standard deviations

Industrial AFA impair cell viability and fermentation kinetics at molasses-based fermentation with cell recycle

After this initial evaluation, the influence of industrial AFA on the fermentation performance of S. cerevisiae was investigated at the physiological and at the transcriptional levels, in a setup that mimicked as far as possible the industrial ethanol process employed in Brazil. The two industrial defoaming treatments currently used in the Brazilian bioethanol industry (Ind_B and Ind_Z) were assessed and benchmarked against a control experiment without the addition of AFA (Table 1). The fermentations were conducted as laboratory-scale replications of the industrial production process, as described elsewhere (Basso et al. 2014; Raghavendran et al. 2017) (Supplementary Fig. S1). In brief, this included yeast cell recycle, acid wash, and fed-batch cultivations in sugarcane must for five consecutive fermentation cycles. For all fermentations, the loss of CO2 was monitored continuously and fermentation yields were determined. The two components of the industrial defoaming treatments were added separately as done in industry, with dispersants (D_B or D_Z) added during the acid wash, while antifoams (AF_B or AF_Z) were added after 1 h of fermentation, where foam issues are more prominent in industrial fermentations (H.V. Amorim, pers. comm.).

When simulating the industrial fermentation conditions, the overall fermentation yields on ethanol, CO2, glycerol, and acetate showed no significant differences among the treatments. Likewise, residual sugar levels were quite similar among conditions (Table 2). The biomass gain along the cultivation cycles was virtually absent (data not shown), which is in agreement with observations under industrial conditions, especially when a molasses-based must is used (Amorim et al. 2011). The carbon balance for all fermentations reached 87–90% of the total carbon input; however, the missing carbon can be attributed to ethanol evaporation and an underestimation of the CO2 formation, since CO2 loss was estimated by weighing the flasks, and could not be properly monitored during feeding and antifoam addition periods.

Table 2 Physiological parameters during yeast fed-batch fermentation of sugarcane must. Yield coefficients and residual sugars in the 5th fermentation cycle of a fed-batch culture of strain S. cerevisiae CAT-1 with no AFA added (control) and with the addition of two industrial AFA (Ind_B and Ind_Z). All yields were calculated based on total consumed sugars (sucrose, glucose, and fructose)

In the first cycles, both industrial AFA treatments (Ind_B and Ind_Z) showed cumulative CO2 production levels comparable to the control, although a reduction was observed in the later cycles (cycle 4 and cycle 5). Treatment Ind_Z imposed the greatest CO2 reduction, showing cumulative production levels corresponding to 52% of the control, while Ind_B treated cells were producing 77% of the CO2 in the control (Fig. 2a). The CO2 evolution over time showed that the CO2 reduction was mainly apparent in the beginning of the fermentations where an increased lag phase was observed for the AFA treatments (Fig. 2b). Ethanol was not quantified during fermentations (only at the end of each cycle), but it can be assumed that the CO2 formation has a 1:1 stoichiometry to ethanol production. Based on this assumption, ethanol productivity was reduced in the same proportion as the CO2 levels, and hence resulted in a reduced ethanol yield after 7 h of fermentation (in the industrial process, the fermentations are typically terminated after 7–9 h). The final CO2 emission (after ~ 24 h) reached the same level for all treatments emphasizing that the total yield is not affected.

Fig. 2
figure 2

The influence of industrial AFA on S. cerevisae CAT-1 fed-batch fermentation of sugarcane must. a Accumulated CO2 loss (in g l−1) over five fermentation cycles, b Accumulated CO2 loss (in g l−1) in the course of the 5th cycle, and c yeast viability at the end of each fermentation cycle, measured as the percentage of viable cells in a population of viable and non-viable cells. Figure legend: white fill (control, no AFA added), gray fill (Ind_B, composed of 150 mg l−1 AF_B and 150 mg l−1 D_B), black fill (Ind_Z, composed of 150 mg l−1 AF_Z and 150 mg l−1 D_Z). Values represent the average of triplicate experiments ± standard deviations

For all treatments, the viability declined for each fermentation cycle as expected (Basso et al. 2008), highlighting the stressful nature of the process (Fig. 2c). The viability profile showed a similar pattern as the CO2 emission, with the industrial AFA treatments reducing viability compared to the control. Viability results explain the underlying reason for the observed reduction in CO2 levels, as fewer living cells were present to produce CO2. This trend of reduced viability was consistent for all industrial AFA treatment experiments performed and moreover, proved to be concentration dependent (data not shown).

Transcriptome analysis reveals that AFA induce additional stress responses in yeast cells during industrial conditions

In order to better understand the physiological perturbations caused by the AFA treatments, samples were collected for transcriptome analysis 4 h after initiation of the last fermentation cycle (cycle 5), which was 1 h after the feeding had ended. RNA sequencing was performed using an Illumina HiSeq instrument with the paired-end method. Considering that the quality of the published genome assembly of the strain CAT-1 was not suitable as a reference for read mapping (Babrzadeh et al. 2012), reads were aligned to the genome of the laboratory strain S. cerevisiae S288C. For each sample, an average of 78.3% ± 2.1 of the sequencing reads were successfully mapped to the reference genome. The relatively low mapping frequency was anticipated due to the distant relation between CAT-1 and the reference strain S288C (Babrzadeh et al. 2012). Reads mapping to open reading frames in S288C were quantified, and a total of 6357 transcripts were identified as being expressed in at least one of the samples. Quality and reproducibility of the RNA-seq samples were assessed by a pairwise comparison of transcript levels for each replicate and showed a Pearson correlation coefficient of at least 0.94 between all samples (Supplementary Fig. S2).

To assess the global effects of the treatments on gene expression, a principal component analysis (PCA) was conducted based on all measured transcript levels (Fig. 3a and Supplementary Fig. S3). The PCA showed that the replicates of each treatment clustered together as expected, except for one of the control treatments which was deviating on the first PC, although was comparable on the second PC. All samples were included in the analysis since the results proved to be similar to excluding the deviating control sample. The relative effects of the treatments were in agreement with the observed physiological changes with the Ind_B treatment resulting in only minor changes, while the Ind_Z treated cells showed more drastic changes in expression levels, as seen by its relative distance to the control in the first two PC dimensions. A differential expression analysis further confirmed the relative pertubations caused by the AFA treatments as a total of 297 genes in Ind_B and 439 genes in Ind_Z were significantly differentially expressed relative to the control (adjusted p value < 0.01). Among these genes, a total of 69 were up-regulated and 42 were down-regulated in both species (Fig. 3).

Fig. 3
figure 3

Overview of transcriptional pertubations caused by the AFA treatments. a Principal component analysis (PCA) of normalized transcript levels in the three replicates of each experimental condition. Samples were collected 4 h after initiation of the last fermentation cycle. b Significantly differentially expressed up-regulated and c down-regulated genes caused by the two AFA treatments relative to the control. Figure legend: white fill (control, no AFA added), gray fill (Ind_B, composed of 150 mg l−1 AF_B and 150 mg l−1 D_B), black fill (Ind_Z, composed of 150 mg l−1 AF_Z and 150 mg l−1 D_Z)

To identify the most apparent transcriptional responses to the AFA treatments and to group the affected genes based on function, gene set analysis was conducted with the normalized transcript counts for each treatment in comparison to the untreated cells (Supplementary Data 1). Gene sets included gene ontology (GO) terms and metabolites connected to genes through the corresponding enzymatic reaction and identified using the reporter feature algorithm (Patil and Nielsen 2005). In addition to these analyses, linear models were applied between the gene expression of each gene in a sample relative to the cell viability observed in the sample, in order to identify gene expression patterns that could be explained by the physiological response. The main transcriptional responses identified are summarized in Fig. 4.

Fig. 4
figure 4

Schematic representation of the affected parts of yeast metabolism upon exposure to industrial AFA. Color key indicate the transcript fold changes in the AFA-treated cells relative to the control and arrows can represent multiple enzymatic reactions. Samples were collected 4 h after initiation of the last fermentation cycle

The most significantly enriched GO term among the up-regulated genes in Ind_B was transmembrane transport (GO:0055085), and some of these transporters were related to sulfur metabolism, i.e., S-methylmethionine transport (GO:0015806) and cysteine transport (GO:0042883). In addition, other sulfur assimilation processes like sulfonate dioxygenase activity (GO:0000907) and sulfur compound catabolic process (GO:0044273) were enriched. For Ind_Z other GO terms related to sulfur metabolism was enriched among up-regulated genes, including cysteine-and methionine biosynthetic process (GO:0019344/0009086) and sulfate assimilation (GO:0000103). The main up-regulated genes within the mentioned GO terms were JLP1 in the Ind_B treatment and MET3 in Ind_Z (Fig. 4). JLP1 encodes an Fe(II)dependent sulfonate/alpha-ketoglutarate dioxygenase which has been shown to be induced under sulfur starvation (Zhang et al. 2001) and is believed to make use of aliphatic sulfonates such as taurine, cysteate, and isethionate as alternative sulfur sources (Hogan et al. 1999). MET3 catalyzes the primary step in the intracellular sulfate activation, where sulfate is reduced to sulfide which is involved in cysteine and methionine metabolism (Gierest et al. 1985). Interestingly, both up-regulated genes propose that sulfur starvation is taking place in both treatments, but also that the cells attempt to cope with the stress using different strategies: Ind_B-treated cells, by recycling intracellular sulfur and Ind_Z-treated cells by increasing the assimilation of extracellular sulfate. Investigating the significantly overrepresented reporter metabolites within these treatments revealed sulfur-containing compounds as well. For Ind_B, sulfate, S-methylmethionine, taurine, and its degradation product aminoacetaldehyde were seen, in agreement with the JLP1 up-regulation. For Ind_Z, only sulfate was significantly overrepresented. Specifically investigating the metabolism of cysteine, methionine and the methionine salvage pathway did not indicate any marked changes in the gene expression controlling these processes. In addition, it was noted that the high-affinity cysteine transporter YCT1 was strongly up-regulated in Ind_B-treated cells (Kaur and Bachhawat 2007). All in all, both industrial AFA treatments induced sulfur starvation responses, although the effects seemed to be more pronounced in Ind_B.

Beside an up-regulation of sulfate assimilation, the Ind_Z treatment generally contained a larger number of significantly enriched GO terms among up-regulated genes compared to Ind_B, indicating a more diverse and global response caused by the treatment. These GO terms included plasma membrane organization (GO:0007009), cellular amino acid biosynthetic process (GO:0008652), and response to stress (GO:0006950). One of these stress response genes was the third most up-regulated gene in Ind_Z, HSP12, which encodes a membrane-bound protein involved in maintaining membrane organization (Welker et al. 2010). Other stress-responsive genes were among the most up-regulated genes in Ind_Z (> 1.5 log2 fold changes) and included HSP30, BTN2, DDR2, OPI10, and CWP2.

Among the down-regulated genes in both Ind_B and Ind_Z, GO terms were strongly enriched for the biosynthesis of lipids and sterols, in particular ergosterol (GO:0006696). Reporter metabolites revealed an ergosterol intermediate, ergosta-5,7,24(28)-trienol, to be underrepresented in both treatments. Further, the fatty acid synthesis pathway was down-regulated including ACC1, FAS1, FAS2, and OLE1, while the elongases ELO1, ELO2, and ELO3 were largely unchanged (Fig. 4). The activation of fatty acids for beta-oxidation by FAA4 was also down-regulated in both treatments. Reporter metabolite analysis showed that the lipids, palmitoyl-CoA in Ind_B, and glycerol-3-phosphocholine in Ind_Z were underrepresented. Palmitoyl-CoA is a thioester, and hence, its underrepresentation is in agreement with the indications that the cells were experiencing sulfur starvation and that fatty acid synthesis and beta-oxidation were down-regulated.

Further, both treatments showed enrichment of the mitochondrion term (GO:0005739) among down-regulated genes as well as several other mitochondrial processes including mitochondrial inner membrane (GO:0005743), respiratory chain complex IV (GO:0005751), and for Ind_Z, also mitochondrial translation (GO:0032543) and mitochondrial ribosomes (GO:0005762/GO:0005763). We also observed that ferrocytochrome c was underrepresented in both treatments.

Although the main transcriptional perturbations were similar in the two industrial treatments, one opposite regulation was seen within cellular amino acid biosynthetic process (GO:0008652). In Ind_B, in particular, branched chain amino acids and lysine biosynthetic genes were down-regulated, while these pathways were unaffected or slightly up-regulated in the Ind_Z-treated cells (Fig. 4).

To determine the transcriptional responses to the physiological phenotypes, linear models were constructed to correlate viability and the expression levels of each gene in all replicates of each treatment. We found a total of 245 genes with significant positive correlation and 293 genes with significant inverse correlation to viability (p < 0.01), and these gene sets were investigated with an enrichment analysis using GO terms. The most enriched GO terms were seen among the genes with negative correlation to viability, meaning that they were higher expressed in the cultivation with low viability, and included the biosynthetic processes of hydrogen sulfide, cysteine, and methionine, in agreement with the observations in the gene set analysis above. Few GO terms were enriched among the genes with positive correlation to viability and they were related to mitochondrial organization and mitochondrial RNA metabolic processes. This observation agreed with the gene set analysis and confirmed that mitochondrial processes were repressed in the low viability/increased stress AFA treatments.

Discussion

During microbial fermentations, AFA are commonly added to prevent the detrimental effects of foam formation, such as loss of fermentor space, leading to increased production costs. AFA have previously been shown to alter growth characteristics of microorganisms, but only a limited number of studies have systematically investigated this. Current knowledge on the topic suggests that the biological effects of AFA exposure is highly dependent on AFA composition, organism of choice, and culture conditions (Routledge et al. 2014). Further, it has been shown that different AFA can lead to both increased and reduced growth rates of yeasts (Holmes et al. 2006) and that the yield of recombinant protein in Pichia pastoris and Eschericha coli can be increased upon AFA exposure (Koch et al. 1995; Routledge and Bill 2012). Despite it being documented that AFA can give rise to negative fermentation performance, no studies have investigated the impact of this under industrial relevant conditions or by using genome-wide transcriptional profiling.

In the present work, we have shown that the growth rate of industrial S. cerevisiae strain CAT-1 is reduced upon exposure to AFA that are used in the Brazilian ethanol industry. In microplate fermentations, at 60 mg l−1, corresponding to the average AFA concentration used in industrial settings (Pecege-Esalq 2012), we observed a 4–6% reduction in growth rates upon exposure to currently used industrial defoaming treatments (Ind_B and Ind_Z), and the response intensified with increasing concentrations. Microplate fermentations with addition of commercially available AFA, frequently used in the fermentation literature (Antifoam C, Sigma 204, and P2000), resulted in no change in growth rates at concentrations up to 375 mg l−1, approximately ten times the concentration of the industrial AFA tested here. Although the amount of AFA utilized in scientific publications is rarely stated, this correspond to concentration higher than what is commonly used for defoaming purposes in publications where this is explicit (Brochado et al. 2010) as well as the manufacturer’s recommendations. These findings suggest that industrially used AFA contain components more harmful to yeasts than what is used in commercially available AFA. One explanation for this might be that industrially used AFA are of varying quality due to efforts to keep production costs down.

To further investigate how the physiological perturbations observed in laboratory fermentations would translate into the industrial ethanol production process, we carried out laboratory simulations of the Brazilian sugarcane-based ethanol production process with recycling of the yeast biomass for 5 cycles (usually done for 200–230 days in industry, with 2 cycles per day). For these fermentations, we used AFA concentrations of five times the average concentration used in industrial settings, in order to make the effects more pronounced. Since the day to day variation in AFA usage in the industry can change up to tenfold, this can be considered as a plausible exposure level under industrial conditions. The physiological effects of the treatments showed a reduction in viability exerted by both treatments after five fermentation cycles. This led to a decrease in the CO2 evolution which can be used as a proxy for the ethanol production and hence suggests a decrease in ethanol productivity. Since the Brazilian ethanol production process is already extremely efficient with yields as high as 93% of the theoretical maximum (Amorim et al. 2011) and that more than half of the ethanol final costs are due to the costs of raw materials, any increase in efficiency would result in major economic gains (Della-Bianca et al. 2013). Taken together, the data suggests that industrial AFA is a potential stress factor to yeast cells in sugarcane fermentation, in addition to the commonly reported ones, such as thermal, osmotic and ethanolic stresses (Basso et al. 2008; Amorim et al. 2011; Della-Bianca et al. 2013).

At the transcriptional level, several of the affected cellular processes were seen in both treatments, despite Ind_Z giving rise to the most pronounced negative physiological phenotype. Interestingly, both treatments showed indications of sulfur starvation, which is surprising considering the sulfur availability from the sulfuric acid added for lowering the pH during recycling. Additionally, a down-regulation of mitochondrial processes was observed, which suggests that oxygen availability is reduced, as previously shown to be caused by AFA addition to culture medium (Koch et al. 1995). The Ind_Z treatment resulted in about half the viability of the control and showed a much more diverse and global transcriptional responses, including a stronger down-regulation of mitochondrial processes and an up-regulation of several stress-response genes which were not seen in Ind_B.

Lipid metabolism proved to be strongly affected in both treatments, in particular sterol and fatty acid synthesis, which was down-regulated with respect to the control. Previous work has shown that sterol abundance increased in Kluyveromyces bulgaricus and Saccharomyces uvarum under exposure to polyoxyalkylene glycol-oleic acid condensates in aerobic cultivations. The relative composition of sterols was unaffected in S. uvarum, while in K. bulgaricus, there was an increase in ergosterol ratio (Pawiroharsono et al. 1987). The authors further stated that the permeability to certain sugars relevant for molasses fermentation, such as glucose and sucrose, was affected by the presence of the AFA investigated. In the present study, under microaerobic conditions, we saw that the ergosterol biosynthetic pathway was down-regulated while only transcript levels of ERG24 and ERG4 were unchanged (Fig. 4). More recently, it was reported that S. cerevisiae grown under ethanol stress had elevated sterol levels while ERG genes were down-regulated (Lahtvee et al. 2016), and emphasizes that the overall sterol content cannot be inferred based on transcript levels alone. Another study showed that oleic acid gets incorporated into phospholipids and glycolipids, and that AFA increased the overall lipid content in the filamentous fungus Aspergillus niger (Nemec and Jernejc 2002). Thus, alterations in lipid biosynthesis may alter membrane fluidity (Lahtvee et al. 2016) and therefore increase sensitivity of yeast cells toward stress factors normally found in molasses-based fermentations (Della-Bianca et al. 2013).

Taken together, the observed reduction in growth rate and cell viability, as well as increased stress responses, could be explained by several factors. The repression of sterol metabolism suggests a change in sterol composition, likely as a response to an increased membrane permeability as previously reported for yeast upon AFA exposure (Pawiroharsono et al. 1987). In inhibitor-rich medium such as molasses (Della-Bianca et al. 2013), this increased permeability could be contributing to the reduced growth rate since the cells attempt to cope with inhibitors by inducing stress responses, e.g., to maintain membrane organization as seen in the up-regulation of HSP12 in this study. This requires an increased ATP demand which is then diverted away from growth processes. In turn, failing to efficiently remove inhibitors might contribute to the observed decrease in viability as well. Another factor is the oxygen transfer rate in the medium, which can be significantly increased or reduced upon AFA addition (Morao et al. 1999; Routledge et al. 2011). The observed down-regulation of mitochondrial processes suggests a reduction in oxygen availability and this could contribute to a reduced ATP generation via oxidative phosphorylation and hence hamper growth. In particular, since the scale-down molasses fermentations conducted in this study took place under microaerobic conditions, changes in dissolved oxygen concentrations due to the presence of AFA might have drastic effects on cell physiology.

Overall, this study has documented distinct negative effects associated with the use of industrial AFA currently used in the Brazilian ethanol fermentation industry. The perturbations were documented at physiological and transcriptional levels. Our results suggest that the decreased viability and induction of stress responses is strongly related to changes in lipid metabolism. These findings highlight the importance in developing less harmful defoaming strategies in industrial fermentations, in order to prevent compromising cellular performance of the production organism.