L.
McDonnell
,
D.
Yusufu
and
A.
Mills
*
School of Chemistry and Chemical Engineering, Queen's University Belfast, Stranmillis Road, Belfast, BT9 5AG, UK. E-mail: andrew.mills@qub.ac.uk
First published on 25th September 2024
The bacterial load (BL), or total viable count, of aerobes can be measured using micro-respirometry, %O2-μR, in which the consumption of dissolved O2 is monitored with respect to incubation time, t. In %O2-μR the ‘bioreactor’ often comprises a canonical plastic tube with a small %O2 sensor; it is simple, fast and accurate and used in automated, commercial instruments for measuring BL. Here we show that it is also possible to measure BL using a new form of micro-respirometry, %CO2-μR, in which the production of CO2 in the growth medium is monitored. In %CO2-μR, the ‘bioreactor’ is the same as that used in %O2-μR, but with a small 3D printed, colour-based %CO2 indicator set in its base and its apparent absorbance, A′, is measured at any t, as it is related to the %CO2 dissolved in the inoculated growth medium. Under aerobic conditions, different inoculations of the facultative anaerobe, E. coli, of different concentrations (101–108 colony forming units (CFU) per mL) are used to generate a series of A′ vs. t profiles, and a straight-line calibration curve. Statistical comparative analysis of the results generated in the above %CO2-μR study, to those generated for the same system but using a commercial %O2-μR system, is used to demonstrate method equivalence. A study of the same system, under anaerobic conditions, using %CO2-μR, shows that %CO2-μR is suitable for measuring the BL of anaerobes. The potential of %CO2-μR for measuring the bacterial load of CO2-generating aerobes and anaerobes is discussed briefly.
Given the importance of making BL measurements in research, the food industry, the environment, and healthcare, there is an increasing need for a rapid, inexpensive method that can be automated to allow the high throughput of samples. In the case of aerobe BL measurements, one very popular alternative to PCM is O2 micro-respirometry, which involves monitoring the consumption of O2 due to bacterial growth and respiration, and which herein is referred to as %O2-μR.1,10–12 Note, although the rate of O2 consumption, as measured using %O2-μR, is most often used as a measure of the rate of microbial metabolism, the rate of generation of a product of metabolism, such as CO2, can also be used for this purpose.
In %O2-μR, the decrease in the concentration of O2 (units: %O2, where 100% = a partial O2 pressure of 1 atm), due to the respiration of a dispersion of an aerobe in an initially air-saturated growth medium, is monitored with respect to incubation time, t, and a value for the threshold time, TT, taken to reach a set %O2 value (the threshold level) is measured. For each bacterial species under study, the system is first calibrated by measuring the value of TT with respect to initial bacterial load, for a series of different loads, and the results used to produce a calibration plot of log(CFU per mL) vs. TT, which is usually a straight line. This calibration line is used subsequently to calculate the log(CFU mL−1) value, i.e., BL, in any subsequent sample, from its measured value for TT.1,11,12
In %O2-μR, an O2-sensitive luminescent dye is usually employed to monitor the %O2 in the growth medium via its luminescence lifetime, τ, and the Stern–Volmer equation,
τo/τ = 1 + Ksv%O2 | (1) |
Clearly, %O2-μR cannot be used to evaluate anaerobe BL levels. However, given most aerobes, and many anaerobes, generate CO2 as they respire,13 it should be possible to measure the BL levels of both these types of microbial species using a micro-respirometry method (%CO2-μR) based on the increase in the dissolved CO2 concentration (units: %CO2) in the growth medium with respect to t. %CO2-μR requires a %CO2 indicator that works in the growth medium and, although many colour and luminescence-based %CO2 indicators exist, most are based on ink films that are not stable in a typical growth medium which contains salts, either due to dye leaching or ion-exchange, which usually leads to the permanent protonation of the dye.14 In addition, %CO2-μR would require the absorbance of the indicator to be measured using a UV/Vis spectrophotometer, which would contribute significantly to its cost.
The above two apparent barriers to the widespread use of %CO2-μR for measuring the BL of aerobes and anaerobes have been addressed separately in recent work carried out by this group. Thus, this group has reported a method of encapsulating CO2-sensitive pigment particles in an extrudable polymer film, like low density polyethylene, LDPE, which can then be 3D printed as a thin plastic film %CO2 indicator that can measure the %CO2 in fresh and seawater over long a long time periods (>3 months).15,16 The LDPE makes these %CO2 indicators very stable, even in very salty aqueous medium, since it acts as a gas-permeable, ion-impermeable, membrane that covers the pigment particles that are CO2-sensitive, in much the same way as the polymer membrane, often polypropylene, protects the integral pH electrode in the Severinghaus CO2 electrode.17 This group has also demonstrated that expensive UV/Vis spectrophotometry, for measuring the absorbance, A, of a %CO2 indicator, can be readily replaced by inexpensive digital photography, coupled with digital colour analysis, DCA, which provides a measure of the indicator's apparent absorbance, A′, which is related directly to A.14,18 In this paper a 3D printed %CO2 indicator, analysed using digital photography + DCA, is used to create, for the first time, a simple %CO2-μR method for measuring BL of both aerobes and anaerobes.
The experimental conditions used to produce the filament and 3D printed CO2-sensitive indicator were the same as those reported previously for a 3D printed O2-sensitive indicator.20 As shown in Fig. 1, the 3D-printed XB/LDPE (CO2-sensitive) indicator was blue in air, and yellow in the presence of 100% CO2.
In this work, all digital photography was carried out using a Canon EOS 700D digital camera (Canon Inc., London, UK) to record images of the 3D printed %CO2 indicator through the glass front window of the incubator, lit using a D65 daylight lamp (HiraLite 14 W full spectrum Daylight lamp, Amazon, UK). Before any photographic images were taken, the custom white balance feature was activated on the camera using a white card, so that the camera locked into the colour temperature of the D65 daylight lamp from the light reflected from the card.
When comparing the different abilities of %O2-μR for measuring BL (an established reference–technique like PCM) and that of %CO2-μR (a new method) an Altman–Bland (A–B) analysis was used.22 Thus, in this work, the difference, d, was plotted against the mean, log(CFU)Av, where,
d = log(CFU mL−1)%O2-μR − log(CFU mL−1)%CO2-μR | (2) |
log(CFU)Av = (log(CFU mL−1)%O2-μR + log(CFU mL−1)%CO2-μR)/2 | (3) |
All PCM and μR BL assays were carried out in triplicate and the average value taken, and in all measurements the CI was ≤±0.1 log(CFU mL−1).
R = (A′0 − A′)/(A′ − A′∞) = α%CO2 | (4) |
The sensitivity of the XB/LDPE indicator towards %CO2 was measured by placing the indicator in 9 mL of nutrient broth (NB, see S1 in the ESI†) at 30 °C, sparging the solution with different known %CO2/N2 blends, and photographing the indicator each time as illustrated in Fig. 2(a). Each image, associated with a known value of %CO2, was analysed using DCA to obtain the corresponding value of A′, which in turn was used to calculate a value for R, using eqn (4).
![]() | ||
Fig. 2 (a) Photographic images of the XB/LDPE indicator in NB at 30 °C, saturated with different known levels of %CO2; (b) plot of R vs. %CO2, where R was calculated using eqn (4) and the DCA calculated values of A′ derived from the photos in (a). |
The subsequent straight line plot R vs. %CO2, arising from this analysis of the photos in Fig. 2(a), is shown in Fig. 2(b), with m, i.e., α, = 0.73 ± 0.04%−1, which indicates a halfway colour transition value (from blue to yellow) (= 1/α) of 1.38 ± 0.08%CO2.
In other work, the 90% response (0 to 5%), and recovery (5 to 0%), times of the 3D-printed XB/LDPE indicator were calculated by measuring A′ vs. time, t, as the atmosphere surrounding the indicator was changed from air to 5%CO2 and back again in a continuous cycle. The results of this work are shown in S4, Fig. S2 in the ESI,† and values for the 90% response and recovery times in the variation in A′ of 8 and 14 min, respectively, were calculated. Reassuringly, as we shall see, in a typical μR experiment the period over which the dissolved %CO2 starts to change rapidly is ca. 60 min, and so considerably longer that the 90% response time of the %CO2 indicator.
An important aspect of this work is the method of production of the %CO2 indicator, i.e., 3D printing, as it is inexpensive, quick, easy to use, and scalable. In order to test how reproducible the 3D printing process was in the production of the XB/LDPE %CO2 indicators, twenty were printed and the values of A′0, A′∞, and α, exhibited by each indicator, were measured in NB at 30 °C.
The results of this work are illustrated in Fig. 3, and show that the measured values of A′0, A′∞, and α for the twenty 3D printed XB/LDPE %CO2 indicators are very consistent, with an average sensitivity (α) of 0.74 ± 0.02%−1, which equates to a % variation of only ca. 2.3%. This % variation in sensitivity value is only slightly higher than that that (1.7%) reported by others for an array of %O2 indicators produced using a commercial ink made for %O2-μR.23
Glucose (C6H12O6) + 6O2 → 6CO2 + 6H2O | (5) |
As noted earlier, the BL of an aerobe is routinely measured using %O2-μR, by monitoring the consumption of O2, but given respiration reaction (5) it follows that it should also be possible to measure the BL for a sample by monitoring the generation of CO2, i.e. using %CO2-μR, instead of %O2-μR. Thus, in a typical %CO2-μR experiment, 1 mL of 104 CFU mL−1 of the E. coli bacteria, made by diluting the stock dispersion of the E. coli of known concentration (see S2 in the ESI†), was used to inoculate 9 mL of NB in a 15 mL Falcon™ tube. The inoculated solution was inverted to mix the two solutions, and then incubated at 30 °C. Digital photography was then used to measure the variation in A′ of the %CO2 indicator, and, at the same time, the Pyroscience O2 meter was used to measure the variation in the lifetime, τ, of the Oculer %O2 indicator which was also present, with respect to incubation time, t. Fig. 4(a) shows a typical collection of photographic images, from which the values of A′ vs. t, illustrated in Fig. 4(b), were calculated using DCA. The latter values were then used to calculate the %CO2 (using eqn (4) and α = 0.74%−1) vs. t data points illustrated in Fig. 4(c). Similarly, the simultaneously measured values of τ vs. t in Fig. 4(b), combined with eqn (1) and Ksv = 0.09%O2−1, allowed the %O2vs. t data points in Fig. 4(c) to be calculated.
A notable feature of the results in Fig. 4(b) is the initial changes in the apparent absorbance, A′, and luminescence lifetime, τ, values in the first 1.5 h of the incubation period. This feature is due to the inoculated Falcon™ tube bioreactor, with the indicators set in its base, see Fig. S1 in ESI,† warming up when first placed inside the incubator. Another, and initially puzzling, feature of the %O2 and %CO2vs. t profiles illustrated in Fig. 4(b), is the apparent faster O2 consumption kinetics compared to that of CO2 generation, so that, by the time all 21% of the dissolved O2 has been consumed (ca. 5.8 h), only 0.2%CO2 has been generated. Interestingly, similar results have been reported by Arain et al. in a μR study in which the variation of %O2 and %CO2 (measured indirectly using a pH indicator) were measured with respect to t.24 It appears likely that this delay is due to the combination of two kinetic features namely, (i) the lag between the consumption of O2 compared with that CO2 production inside the cell, because mechanistically the former is a precursor to the latter via a complex series of reactions, and (ii) the faster kinetics of O2 diffusion into the cell compared to those for CO2 diffusion out of the cell, where it is detected.
The reproducibility of the %CO2-μR method was tested by repeating the same experiment as above ten times and using the results to generate the ten A′ vs. t plots shown in S5, Fig. S3 in the ESI.† These results were used to calculate an average value of TT = 7.43 ± 0.18 h, i.e., a % variance of ±2.4%, which indicates that the %CO2-μR methodology for one inoculum at least is very reproducible. To test the %CO2-μR method more fully, the same set up as above was used to generate a series of A′ vs. t profiles for different inoculum concentrations of E. coli under aerobic conditions, the results of which are illustrated in Fig. 5(a), with the original photographs, from which the values of A′ were calculated, shown in S6, Fig. S4 in the ESI.† These results were then used to construct the usual straight line calibration plot of log(CFU mL−1), vs. TT, which is illustrated in Fig. 6(b) and, as noted before for %O2-μR, can then be used to determine the BL load in any subsequent samples.
![]() | ||
Fig. 5 (a) A′ vs. incubation time, t, profiles, calculated using DCA from the photos in Fig. S4 in the ESI,† for a series of 1 mL, 10-fold dilutions of the E. coli stock culture, covering 108 to 101 CFU mL−1, in units of ten from left to right, used to inoculate 9 mL of NB under aerobic conditions. The red dashed line highlights the value of TT, for each inoculum, for when A′ = 0.7; (b) initial inoculum BL, log(CFU mL−1), vs. TT, derived from the data in (a). A line of best fit analysis of the data (dashed line) reveals m ± Δm = −0.68 ± 0.02 log(CFU mL−1) h−1 and c ± Δc = 9.4 ± 0.2 log(CFU mL−1). |
![]() | ||
Fig. 6 BL values of 14 different E. coli samples were measured using %CO2-μR and %O2-μR systems and the results plotted in (a) as log(CFU mL−1)%CO2-μRvs. log(CFU mL−1)%O2-μR and (b) as d vs. log(CFU)Av, as defined by eqn (2) and (3). ![]() ![]() |
The above results suggest that the XB/LDPE %CO2 indicator can be used in μR to measure BL under aerobic conditions and is a first reported example of %CO2-μR in which the %CO2 was measured directly using a %CO2 sensor.
The simplest and most common method employed to compare the results, generated by two methods that measure the same quantity, is to plot them against each other, since ideally, a straight line should be generated with the gradient, m = 1, the intercept, c = 0 and the square of the correlation coefficient, r2 = 1.22 A rough guide to the degree of the equivalence of the two methods is provided by how far the actual values of m, c and r, deviate from the ideal, namely, m = 1, c = 0 and r = 1, and whether these ideal values of m and c fall within the limits of agreement described by m ± Δm and c ± Δc, respectively.
The plot of log(CFU mL−1)%CO2-μRvs. log(CFU mL−1)%O2-μR, for the 14 samples is illustrated in Fig. 6(a) and reveals values for m, c and r, of 0.993 ± 0.06, −0.0013 ± 0.30 and 0.9941, respectively. These results suggest that %CO2-μR and %O2-μR, can be considered equivalent, since the ideal values of 1 and 0 lie with the limits of agreement (e.g., 0.993 + 0.06 to 0.993 − 0.06 for m) for m and c, respectively.
The above straight line plot analysis of the data is not sufficiently rigorous to be sure the two methods are equivalent and so it is usual to also construct an A–B plot of the data. The A–B plot of the data in Fig. 6(a), illustrated in Fig. 6(b), shows the data points are randomly distributed about (= 0.03) and within the limits of agreement (0.03 ± 0.46). A Shapiro–Wilk normality test on the differences showed the data were normally distributed (P = 0.548). In addition, in Fig. 6(b), it can be seen that the grey shaded area, defined by the confidence intervals of
, (0.03 ± 0.06), encapsulates the black broken line of equality (d = 0), which indicates that there is little, or no, significant bias between the two methods.25 Thus, the A–B plot illustrated in Fig. 6(b) indicates that the new %CO2-μR method is equivalent to the reference %O2-μR method for measuring BL.
The latter idea was tested using E. coli and the same growth NB as before plus sodium sulfite, Na2SO3, as the latter is very effective in removing dissolved O2 from aqueous solution.28 Thus, in a typical experiment, a 1 mL inoculum of 104 CFU mL−1E. coli was added to 9 mL of the NB, containing 50 mg of Na2SO3, in a Falcon™ tube with both a XB/LDPE %CO2 and Oculer %O2 indicator set in its base, see Fig. S1 in the ESI.† The inoculated Falcon™ tube was then quickly stoppered, incubated at 30 °C and the apparent absorbance, A′, of the %CO2 indicator, and τ of the %O2 indicator monitored with incubation time, t.
The results of this work in the form of plots of τ (for the %O2 indicator) and A′ (derived from photographs of the %CO2 indicator using DCA) vs. t are shown in Fig. 7 and show, as expected, that the concentration of dissolved O2, %O2, remained approximately 0 throughout the 22 h run (since τ ≈ τo), whereas the concentration of %CO2 increased (and so A′ decreased) with increasing t, due to the fermentation-driven growth of the E. coli. A comparison of the two A′ vs. t profiles due to the generation of CO2, illustrated in Fig. 4(b) and 7, under aerobic and anaerobic conditions, respectively, shows that the TT values, for an otherwise identical μR experiment, are 7.5 and 12 h, respectively. This feature is consistent with the well-known observation that, under the same conditions, E. coli's growth kinetics due to respiration in air are faster than that associated with fermentation in the absence of air.
The results illustrated in Fig. 7 suggest that it should be possible to use %CO2-μR to measure the BL of CO2-generating facultative and non-facultative anaerobes. To test this proposition, the same experimental conditions as used above were employed to generate a series of A′ vs. t profiles for different concentrations of E. coli under anaerobic conditions and the results are illustrated in Fig. 8(a). Note, the original photographs, from which the values of A′ were calculated, are given in S6, Fig. S5 in the ESI.†
![]() | ||
Fig. 8 (a) A′ vs. t profiles, calculated using DCA and the photos illustrated in Fig. S5 in the ESI,† for a series of 1 mL, 10-fold dilutions of the E. coli stock culture, from 108–101 CFU mL−1, in units of ten from left to right, used to inoculate 9 mL of NB under anaerobic conditions. The red dashed line identifies the TT value, for each inoculum, for when A′ = 0.7. (b) Log(CFU mL−1), vs. TT, based on the data in (a). A line of best fit analysis of the data (dashed line) reveals m ± Δm = −0.52 ± 0.01 log(CFU mL−1) h−1 and 10.4 ± 0.2 log(CFU mL−1), respectively. |
The above results, based on E. coli dispersions under anaerobic conditions, show that the XB/LDPE %CO2 indicator can be used in μR for assessing the BL of CO2-generating anaerobes as it is able to generate a good straight calibration line plot of log(CFU mL−1) vs. TT, as illustrated in Fig. 8(b). As with aerobic μR for measuring BL, such a calibration plot can then be used to determine a value for the log(CFU mL−1) in any subsequent sample containing the anaerobe under test, from its measured TT value recorded under the same anaerobic conditions.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4an01016g |
This journal is © The Royal Society of Chemistry 2024 |