Yeast Propagation Control: Low Frequency Electrochemical Impedance Spectroscopy as an Alternative for Cell Counting Chambers in Brewery Applications

: Electrochemical impedance spectroscopy is a powerful tool in life science for cell and pathogen detection, as well as for cell counting. The measurement principles and techniques using impedance spectroscopy are highly diverse. Di ﬀ erences can be found in used frequency range ( β or α regime), analyzed quantities, like charge transfer resistance, dielectric permittivity of double layer capacitance and in o ﬀ - or online usage. In recent contributions, applications of low-frequency impedance spectroscopy in the α regime were tested for determination of cell counts and metabolic burden in Escherichia coli and Saccharomyces cerevisiae . The established easy to use methods showed reasonable potential in the lab scale, especially for S. cerevisiae . However, until now, measurements for cell counts in food science are generally based on Thoma cell counting chambers. These microscopic cell counting methods decelerate an easy and quick prediction of yeast viability, as they are labor intensive and result in a time delayed response signal. In this contribution we tested our developed method using low frequency impedance spectroscopy locally at an industrial brewery propagation site and compared results to classic cell counting procedures.


Introduction
Cultivation technology using microorganisms, like bacteria and yeasts, is not only applied in pharmaceutical technology, but also in branches like food and bulk chemical production and highly emerging in waste-to-value concepts [1]. In all those applications, process monitoring (process analytical technology-PAT) is of utmost importance to control the system and react timely upon metabolic changes. pH, dissolved oxygen (dO 2 ) and off-gas analytics are widely applied and can be regarded as state of the art in today's industrial cultivation process to aim for product quality and safety. The most important parameter in bioprocesses, the biomass, is generally analyzed using offline methods like cell chamber counting, cell dry weight or optical density (OD 600 ). As control and analytical tools aiming for biomass determination are often cost intensive, they are used for high value products like in the pharmaceutical biotech, by default. In contrast, classical bulk food products-such as yeast and beer-are low value products with costs of only some Euro/kg. These products are therefore produced in uncontrolled environments and often result in batch to batch variations and In this study, we applied the developed impedance measurement technique in the α-regime in a laboratory lab scale run using media and yeast from Stiegl brewery (Salzburg, Austria). In a second step we tested our system at the brewery on site. Samples were taken from the propagation reactor and put to a self-built test setup in order to maintain sterility of the propagation reactor. The samples were temperature-controlled and measured. Different state-of-the-art methods were applied for determination of the corresponding total biomass-dry cell weight (DCW), optical density (OD 610 ) and offline flow cytometry with propidium iodide for cell physiology evaluation. For the on-site measurement offline determination using cell counting chambers were performed. With this knowledge, we were able to correlate the total biomass to the extracted C DL and show the straight-forward application of the system for cell number determination even in industrial-scale applications.

Expression Host and Cultivation
Lab cultivations were performed using the S. cerevisiae strain, supplied by Brauerei Stiegl (Salzburg, Austria). Batch propagation was performed in a stainless-steel Sartorius Biostat Cplus bioreactor (Sartorius, Göttingen, Germany) with a 10 L working volume. We added 1 L of supplied yeast suspension (Stiegl) to an autoclaved reactor with 4 L of wort. The stirrer was operated at 1400 rpm stirrer speeds with an aeration of 2 vvm pressurized air. Temperature was set to 12 • C throughout the whole process.
Industrial propagation was performed at Brauerei Stiegl in Salzburg using a stainless-steel propagator with a volume of 1500 L. Process aeration was alternated with pumping of the broth. Temperature in the propagator was uncontrolled, Figure 1a. In order to maintain sterility intact, samples were taken through a sampling valve and put into a custom-made cooling reactor with included impedance measurement setup. The so taken samples were analyzed using the atline measurement setup given in detail in Reference [38].
Sensors 2020, 16, x FOR PEER REVIEW 3 of 12 In this study, we applied the developed impedance measurement technique in the α-regime in a laboratory lab scale run using media and yeast from Stiegl brewery (Salzburg, Austria). In a second step we tested our system at the brewery on site. Samples were taken from the propagation reactor and put to a self-built test setup in order to maintain sterility of the propagation reactor. The samples were temperature-controlled and measured. Different state-of-the-art methods were applied for determination of the corresponding total biomass-dry cell weight (DCW), optical density (OD610) and offline flow cytometry with propidium iodide for cell physiology evaluation. For the on-site measurement offline determination using cell counting chambers were performed. With this knowledge, we were able to correlate the total biomass to the extracted CDL and show the straightforward application of the system for cell number determination even in industrial-scale applications.

Expression Host and Cultivation
Lab cultivations were performed using the S. cerevisiae strain, supplied by Brauerei Stiegl (Salzburg, Austria). Batch propagation was performed in a stainless-steel Sartorius Biostat Cplus bioreactor (Sartorius, Göttingen, Germany) with a 10 L working volume. We added 1 L of supplied yeast suspension (Stiegl) to an autoclaved reactor with 4 L of wort. The stirrer was operated at 1400 rpm stirrer speeds with an aeration of 2 vvm pressurized air. Temperature was set to 12°C throughout the whole process.
Industrial propagation was performed at Brauerei Stiegl in Salzburg using a stainless-steel propagator with a volume of 1,500 L. Process aeration was alternated with pumping of the broth. Temperature in the propagator was uncontrolled, Figure 1a. In order to maintain sterility intact, samples were taken through a sampling valve and put into a custom-made cooling reactor with included impedance measurement setup. The so taken samples were analyzed using the atline measurement setup given in detail in Reference [38].  [38].

Process Analytics
For the lab scale experiment, we determined biomass by cell dry weight (CDW) measurements and OD600. Following methodological procedure was performed: 2 mL of the cultivation broth were centrifuged at 4500 g, subsequently washed with 0.9% NaCl solution and centrifuged again. The cells were put into a furnace at 105 °C for 48 h and afterwards the dry pellet was evaluated gravimetrically

Process Analytics
For the lab scale experiment, we determined biomass by cell dry weight (CDW) measurements and OD 600 . Following methodological procedure was performed: 2 mL of the cultivation broth were centrifuged at 4500 g, subsequently washed with 0.9% NaCl solution and centrifuged again. The cells were put into a furnace at 105 • C for 48 h and afterwards the dry pellet was evaluated gravimetrically using an analytical scale. These measurements were performed in triplicates and the mean error for CDW was always below 3%. OD 600 was measured at a UV/VIS photometer Genisys 20 (Thermo Scientific, Waltham, MA, US). For determination of sugar in the fermentation broth we used a Supelco C-610H HPLC column (Supelco, Bellefonte, PA, USA) on an Ultimate 3000 HPLC system (Thermo Scientific, Waltham, MA, US) with 0.1% H 3 PO 4 as a running buffer at 0.5 mL/min. Ethanol in the cultivation broth was analyzed using an Aminex HPLC column (Biorad, Hercules, CA, USA) on an Agilent 1100 System (Agilent Systems, Santa Clara, CA, USA) with 40 mM H 2 SO 4 as a running buffer at 0.6 mL/min. The cultivation off-gas was analyzed by following gas sensors: IR for CO 2 and ZrO 2 -based for O 2 (Blue Sens Gas analytics, Herten, Germany). Details on the flow cytometry method can be found elsewhere [39]. A CytoSense flow cytometer (CytoBuoy, Woerden, Netherlands) with two forward scatter (FSC), one sideward scatter (SSC) and two fluorescence channels (green, red) was used for analysis of the yeast cells. The implemented laser had a wavelength of 488 nm. The configuration of the filter set was 515-562 ± 5 nm for the green fluorescence channel (FL-green, used for fluorescein diacetate) and 605-720 ± 5 nm for the red fluorescence channel (FL-red, used for propidium iodide) Industrial cell counting measurements were performed using Thoma cell measurements. The number of cells in the suspension can be determined by Equation (1) With 400 × 10 4 being the factor of Thoma chamber. Differentiation between viable and dead cells was not performed. A sketch of the used Thoma chamber is given in Supplementary Figure S1.

Impedance Measurements
Analysis of viable cell concentrations is of utmost importance for industrial applications. The majority of these systems rely on β-dispersion (10 7 -10 4 Hz) and are strongly dependent on process parameters (e.g., stirring, temperature, pH, salt and substrate concentration, etc.) and the cultivation phase (exponential growth phase, starvation phase, etc.) [12,33]. For our impedance measurements we exploited a different physical phenomenon (α-dispersion), which inherits valuable information regarding the biomass concentration. The "α-dispersion effect", at frequencies below 1 to 10 kHz is most likely a result of deformation of ionic species around the cell membranes. The dielectric response was therefore proportional to the ionic charge gathered around the membrane of adsorbed cells on the stainless-steel electrode [20,21]. Impedance measurements were recorded in the range of 10 6 to 10 −1 Hz with amplitudes up to 100 mV using the N4L PSM 1735 frequency analyzer (Newton4th Ltd., Leicester, UK). No DC voltage was applied during the measurement. Because measurements in this frequency range are largely determined by the double-layer region between the electrode and the media, rather minor interferences with the process parameters (aeration and stirring) were to be expected. For measurement an online probe presented in a recent publication was used for the lab scale cultivation [38].

Results and Discussion
Yeast propagation is a crucial process step in industrial breweries. Differences in the cell numbers after propagation directly affect the quality of the produced beer and result in batch to batch variations and in the worst-case scenario to batch loss. A reliable online or even inline method for viable cell number and metabolic state determination would revolutionize this step. Therefore, we tested the proposed measurement technique based on α-dispersion [37] in real food technological environment in a small scale fermenter (10 L) in the lab, as well as in a large scale propagation tank with a volume of 1500 L on the production site.

Labscale Propagation
For the feasibility measurement a lab scale propagation was performed using raw materials and yeasts supplied by Stiegl brewery. We used constant temperature of 12 • C for microbial growth and sampled regularly in intervals of about 2 h. Dissolved oxygen was always higher than 70%, indicating complete aerobic growth. Table 1 shows results of process analytics. The cultivation reveals well-known behavior for S. cerevisiae indicating, beside aerobic growth, ethanol fermentation based on the limits in aerobic capacity of the cells, so called Crabtree effect [40]. So, highest amounts of 1.7% ethanol were produced within the lab scale cultivation, which are believed to show no inhibitory effects on the yeast growth [41]. Viability was also monitored via flow cytometry using propidium iodide as marker for cell death. Within the run almost no dead cells were found (0.1% in the last sample-Supplementary Figure S2). In Figure 2 the raw data for impedance measurement during the lab scale cultivation are presented, right after inoculation and at elevated batch times. a presents the Nyquist plot, while b shows the Bode plot with absolute value of impedance.
Sensors 2020, 16, x FOR PEER REVIEW 5 of 12 sampled regularly in intervals of about 2 h. Dissolved oxygen was always higher than 70%, indicating complete aerobic growth. Table 1 shows results of process analytics. .02 The cultivation reveals well-known behavior for S. cerevisiae indicating, beside aerobic growth, ethanol fermentation based on the limits in aerobic capacity of the cells, so called Crabtree effect [40]. So, highest amounts of 1.7% ethanol were produced within the lab scale cultivation, which are believed to show no inhibitory effects on the yeast growth [41]. Viability was also monitored via flow cytometry using propidium iodide as marker for cell death. Within the run almost no dead cells were found (0.1% in the last sample-Supplementary Figure S2). In Figure 2 the raw data for impedance measurement during the lab scale cultivation are presented, right after inoculation and at elevated batch times. a presents the Nyquist plot, while b shows the Bode plot with absolute value of impedance.  As clearly visible in the raw data, a small shoulder is visible for the media contribution at high frequencies, while at lower frequencies the contribution of the double layer is predominant. In general, such data are fitted using complex circuit models.
As already given in recent publications [37,38] the data can be fitted using constant phase elements (CPE) in combination with a resistance (R), so called R-CPE element, and resulted in an excellent fit of the raw data. Using CPE n and Q fitting parameters are used to calculate the sample capacitance (C) according to Equation (2). [42], with C being the sample capacitance and R the real part of the R-CPE element.
However, we experienced problems in correct fitting of the R element since the error was tremendous high and this makes it impossible to determine the accurate capacitance in this case via Equation (2). These problems were also experienced in recent studies using low cell densities [37,38]. Therefore, we fitted our data using in a straightforward way: an offset resistance for cable and contact resistance and an ideal capacitor equivalent circuit element for the double layer capacitance were applied, see Equation (3).
Sensors 2020, 16, x FOR PEER REVIEW 6 of 12 As clearly visible in the raw data, a small shoulder is visible for the media contribution at high frequencies, while at lower frequencies the contribution of the double layer is predominant. In general, such data are fitted using complex circuit models.
As already given in recent publications [37,38] the data can be fitted using constant phase elements (CPE) in combination with a resistance (R), so called R-CPE element, and resulted in an excellent fit of the raw data. Using CPE n and Q fitting parameters are used to calculate the sample capacitance (C) according to Equation (2). [42], with C being the sample capacitance and R the real part of the R-CPE element.
However, we experienced problems in correct fitting of the R element since the error was tremendous high and this makes it impossible to determine the accurate capacitance in this case via Equation (2). These problems were also experienced in recent studies using low cell densities [37,38]. Therefore, we fitted our data using in a straightforward way: an offset resistance for cable and contact resistance and an ideal capacitor equivalent circuit element for the double layer capacitance were applied, see Equation (3).
The equivalent circuit shown in Equation (3) was fitted by complex non-linear least square fitting (CNLS) using the software ZView (Scribner, Southern Pines, NC, USA) using a resistance for the cable/setup offset and an ideal capacitanceThe results for the fit are displayed in Figure 2c. The fit with the ideal capacitance is not ideal, but give C values with an error of about 3%. The capacitance values for the entire cultivation are given in Figure 3a. As the raw data showed certain fluctuations in the signal, we smoothed the signal using a polynomial smoothing procedure using the OriginPro 2016G (Northampton, MA, USA). This results in the impedance signal given in Figure 3b. Certain plateaus are visible in the impedance data. Such effects were already observed in recent publications and indicate a change in metabolism of the cells affecting ionic cloud and membrane composition. First plateau at about 25 h resembles the depletion of the monosaccharide glucose, which is the favored carbon source for most heterotroph microorganisms. A switch to disaccharides in metabolism is clearly indicated. The main disaccharide in wort is maltose, but also different polysaccharides are present, which cannot be separated by our respective HPLC method and not metabolized. Upon carbon source depletion at about 45 h also the impedance signal stagnates, indicating the direct response to cell number in the cultivation broth, even when complex media, like wort, is used. This is also indicated by reduction of ethanol in the broth, since S. cerevisiae starts feeding on alcohol as last available carbon source. Figure 4a shows the linear fit of the impedance signal extracted using the equivalent circuit in Equation (3). The fit quality is high except for deviations in the beginning, when only a low number of cells are present in the broth and upon metabolism change as discussed before. We predicted the CDW based on the fit in Figure 4a.

R offset C DL
The equivalent circuit shown in Equation (3) was fitted by complex non-linear least square fitting (CNLS) using the software ZView (Scribner, Southern Pines, NC, USA) using a resistance for the cable/setup offset and an ideal capacitanceThe results for the fit are displayed in Figure 2c. The fit with the ideal capacitance is not ideal, but give C values with an error of about 3%. The capacitance values for the entire cultivation are given in Figure 3a. As the raw data showed certain fluctuations in the signal, we smoothed the signal using a polynomial smoothing procedure using the OriginPro 2016G (Northampton, MA, USA). This results in the impedance signal given in Figure 3b. Certain plateaus are visible in the impedance data. Such effects were already observed in recent publications and indicate a change in metabolism of the cells affecting ionic cloud and membrane composition. First plateau at about 25 h resembles the depletion of the monosaccharide glucose, which is the favored carbon source for most heterotroph microorganisms. A switch to disaccharides in metabolism is clearly indicated. The main disaccharide in wort is maltose, but also different polysaccharides are present, which cannot be separated by our respective HPLC method and not metabolized. Upon carbon source depletion at about 45 h also the impedance signal stagnates, indicating the direct response to cell number in the cultivation broth, even when complex media, like wort, is used. This is also indicated by reduction of ethanol in the broth, since S. cerevisiae starts feeding on alcohol as last available carbon source. Figure 4a shows the linear fit of the impedance signal extracted using the equivalent circuit in Equation (3). The fit quality is high except for deviations in the beginning, when only a low number of cells are present in the broth and upon metabolism change as discussed before. We predicted the CDW based on the fit in Figure 4a.
As process analytics in Stiegl brewery rely on cell counting methods, we also evaluated our cell numbers in the broth using flow cytometry as explained in the material part. We did a correlation of the CDW to the cell numbers determined by flow cytometry. Linear correlation is given in Figure 4b and used for calculation of the cells/mL in Figure 4c. Maximum cell numbers are about 100 million cells/mL during the lab scale cultivation. The results are given in Figure 4c. Except in the beginning the fit quality is high enough and gives a reasonable approximation of the CDW in the broth.  As process analytics in Stiegl brewery rely on cell counting methods, we also evaluated our cell numbers in the broth using flow cytometry as explained in the material part. We did a correlation of the CDW to the cell numbers determined by flow cytometry. Linear correlation is given in Figure 4b and used for calculation of the cells/mL in Figure 4c. Maximum cell numbers are about 100 million cells/mL during the lab scale cultivation. The results are given in Figure 4c. Except in the beginning the fit quality is high enough and gives a reasonable approximation of the CDW in the broth.

Impedance Measurement Atline at Stiegl Brewery
After successful tests in the small scale we tested the impedance measurement technique on site at Stiegl brewery in a 1,500 L propagation reactor. To keep the sterile environment of the propagation system intact, broth was sampled and circulated inside the measurement apparatus built by Gruber GmbH (Salzburg, Austria). During propagation of yeast in the reactor temperature rose from about 8°C to 12°C in the end of the cultivation. The measurement system was cooled to 7.6 to 8°C for measurement. For inline measurement drifts in temperature during the run would have to be considered as the impedance signal is generally strongly dependent on the measurement temperature. Furthermore, aeration in the propagation system is very different to the lab scale setup.  As process analytics in Stiegl brewery rely on cell counting methods, we also evaluated our cell numbers in the broth using flow cytometry as explained in the material part. We did a correlation of the CDW to the cell numbers determined by flow cytometry. Linear correlation is given in Figure 4b and used for calculation of the cells/mL in Figure 4c. Maximum cell numbers are about 100 million cells/mL during the lab scale cultivation. The results are given in Figure 4c. Except in the beginning the fit quality is high enough and gives a reasonable approximation of the CDW in the broth.

Impedance Measurement Atline at Stiegl Brewery
After successful tests in the small scale we tested the impedance measurement technique on site at Stiegl brewery in a 1,500 L propagation reactor. To keep the sterile environment of the propagation system intact, broth was sampled and circulated inside the measurement apparatus built by Gruber GmbH (Salzburg, Austria). During propagation of yeast in the reactor temperature rose from about 8°C to 12°C in the end of the cultivation. The measurement system was cooled to 7.6 to 8°C for measurement. For inline measurement drifts in temperature during the run would have to be considered as the impedance signal is generally strongly dependent on the measurement temperature. Furthermore, aeration in the propagation system is very different to the lab scale setup.

Impedance Measurement Atline at Stiegl Brewery
After successful tests in the small scale we tested the impedance measurement technique on site at Stiegl brewery in a 1500 L propagation reactor. To keep the sterile environment of the propagation system intact, broth was sampled and circulated inside the measurement apparatus built by Gruber GmbH (Salzburg, Austria). During propagation of yeast in the reactor temperature rose from about 8 • C to 12 • C in the end of the cultivation. The measurement system was cooled to 7.6 to 8 • C for measurement. For inline measurement drifts in temperature during the run would have to be considered as the impedance signal is generally strongly dependent on the measurement temperature. Furthermore, aeration in the propagation system is very different to the lab scale setup. The lab scale is an ideal stirred tank reactor with aeration through sparger and mechanic stirring. The industrial propagation relies on alternating pumping of the cells and aeration of the system. Therefore, differences in the metabolic state of the samples is to be expected during this process design as cells always react upon differences in the dissolved oxygen level and may switch between aerobic and anaerobic conditions. After sampling, cells were circulated in the built setup, compare to Figure 1b. In Figure 5a Chemosensors 2020, 8, 27 8 of 11 the impedance signal measured in the measurement setup is given. Shallow fluctuations are visible at each sample. The data points were averaged and compared to the measured cell concentrations.
The lab scale is an ideal stirred tank reactor with aeration through sparger and mechanic stirring. The industrial propagation relies on alternating pumping of the cells and aeration of the system. Therefore, differences in the metabolic state of the samples is to be expected during this process design as cells always react upon differences in the dissolved oxygen level and may switch between aerobic and anaerobic conditions. After sampling, cells were circulated in the built setup, compare to Figure 1b. In Figure 5a the impedance signal measured in the measurement setup is given. Shallow fluctuations are visible at each sample. The data points were averaged and compared to the measured cell concentrations. Cell numbers were analyzed by Stiegl in-house analytical methods using cell counting chambers. These values were compared to averaged capacitance value. Linear correlations are given in Figure 5b. A high linearity between the received cell numbers and the impedance signal is given. However, we see slight deviations to the lab scale system in the linear regression. While capacitance values extracted from the signal are very similar between the two measurement runs, absolute cell numbers differed. Differences in cell counting between cell chamber counting methods and flow cytometric analysis are also reported in literature for different applications [43,44]. Flow cytometry analyzes a high number of cells within a very short time and can distinguish between single, budding and dead cells. With cell chamber counting only a very limited number of cells can be analyzed at a time. Including staining would make it possible to distinguish further between living and dead cells. Summing up, cell counting chambers seems to underestimate the cell numbers in the cultivation broth as dedicated by the measured impedance signal. Therefore, calibration with the method of choice is crucial for exact determination of cell numbers used for further transfer of a certain cell number to the subsequent process step. Furthermore, differences in metabolic state between the labscale and the large-scale system could account for the signal differences. It is very likely that cells in the 1,500 L reactor are very homogenously active in metabolism and result in a very different signal, also based on the differences in aeration and mixing. Table 2 presents the calculated cell numbers based on the calibration for yeast propagation at the Stiegl brewery on site.  Cell numbers were analyzed by Stiegl in-house analytical methods using cell counting chambers. These values were compared to averaged capacitance value. Linear correlations are given in Figure 5b. A high linearity between the received cell numbers and the impedance signal is given. However, we see slight deviations to the lab scale system in the linear regression. While capacitance values extracted from the signal are very similar between the two measurement runs, absolute cell numbers differed. Differences in cell counting between cell chamber counting methods and flow cytometric analysis are also reported in literature for different applications [43,44]. Flow cytometry analyzes a high number of cells within a very short time and can distinguish between single, budding and dead cells. With cell chamber counting only a very limited number of cells can be analyzed at a time. Including staining would make it possible to distinguish further between living and dead cells. Summing up, cell counting chambers seems to underestimate the cell numbers in the cultivation broth as dedicated by the measured impedance signal. Therefore, calibration with the method of choice is crucial for exact determination of cell numbers used for further transfer of a certain cell number to the subsequent process step. Furthermore, differences in metabolic state between the lab-scale and the large-scale system could account for the signal differences. It is very likely that cells in the 1500 L reactor are very homogenously active in metabolism and result in a very different signal, also based on the differences in aeration and mixing. Table 2 presents the calculated cell numbers based on the calibration for yeast propagation at the Stiegl brewery on site. We want to highlight that we can do accurate measurement of a very shallow range of cell numbers in an industrial propagation for beer production in a large-scale system. Based on these results, this straightforward measurement principle might find further applications in different food industrial branches, using yeasts and different microorganism.

Conclusions and Outlook
Within this contribution, we transferred the established impedance measurement method in the α-dispersion regime to industrial large-scale application in the brewing industry. We tested propagation in a lab scale system using wort and yeasts from the industrial partner. With the at-line measurement system we were able to measure impedance of the propagation broth on site without risking sterility of the propagation system. The received correlation between cell counting (done at Stiegl) and the impedance signal can now be used for direct measurement of cell numbers during propagation. The next step would be the manufacturing of an inline measurement system in the propagation tank including temperature compensation of the system and direct visualization of cell numbers.