- 1Chemistry and Biochemistry Department, California State Polytechnic University, Pomona, CA, United States
- 2Blue Marble Institute of Science, Seattle, WA, United States
Introduction: Stringent cleaning procedures during spacecraft assembly are critical to maintaining the integrity of life-detection missions. To ensure cleanliness, NASA spacecraft are assembled in cleanroom facilities, where floors are routinely cleansed with Kleenol 30 (K30), an alkaline detergent.
Methods: Through metabolomic and cultivation approaches, we show that cultures of spacecraft-associated Acinetobacter tolerate up to 1% v/v K30 and are fully inhibited at ≥2%; in comparison, NASA cleanrooms are cleansed with ~0.8-1.6% K30.
Results: For A. johnsonii 2P08AA (isolated from a cleanroom floor), cultivations with 0.1% v/v K30 yield (1) no changes in cell density at late-log phase, (2) modest decreases in growth rate (~17%), (3) negligible lag phase times, (4) limited changes in the intracellular metabolome, and (5) increases in extracellular sugar acids, monosaccharides, organic acids, and fatty acids. For A. radioresistens 50v1 (isolated from a spacecraft surface), cultivations yield (1) ~50% survivals, (2) no changes in growth rate, (3) ~70% decreases in the lag phase time, (4) differential changes in intracellular amino acids, compatible solutes, nucleotide-related metabolites, dicarboxylic acids, and saturated fatty acids, and (5) substantial yet differential impacts to extracellular sugar acids, monosaccharides, and organic acids.
Discussion: These combined results suggest that (1) K30 manifests strain-dependent impacts on the intracellular metabolomes, cultivation kinetics, and survivals, (2) K30 influences extracellular trace element acquisition in both strains, and (3) K30 is better tolerated by the floor-associated strain. Hence, this work lends support towards the hypothesis that repeated cleansing during spacecraft assembly serve as selective pressures that promote tolerances towards the cleaning conditions.
1. Introduction
Maintaining low biological contamination in cleanroom-type facilities are critical components to spacecraft assembly (Rummel, 1992; NASA, 2011; Frick et al., 2014), crewed spacecraft exploration (Spry et al., 2020), delivery of healthcare facilities (Shams et al., 2016; Rutala and Weber, 2019), and the manufacturing of pharmaceuticals (Sandle, 2015). Common to these confined environments and endeavors are the routine cleansing of non-critical hard surfaces and support equipment with ≥70% v/v ethyl alcohol, 70% v/v isopropyl alcohol, 3% v/v hydrogen peroxide, and/or disposable disinfectant wipes containing quaternary ammonium compounds (Agalloco and Carleton, 2007; Rutala et al., 2008; Wong et al., 2013; Frick et al., 2014; Checinska Sielaff et al., 2019).
For pharmaceutical facilities, and to lesser extent healthcare facilities, the floors are additionally cleansed with disinfectants and/or detergents such as benzalkonium chloride (and other quaternary ammonium compounds), hydrogen peroxide, sodium hypochlorite, and other alternatives (Murtough et al., 2001; Rutala et al., 2008; Sandle, 2012; Eissa et al., 2014; Han et al., 2021; Tembo et al., 2022). While debated in efficacy (Daschner et al., 1982; Meyer and Cookson, 2010; Suleyman et al., 2018), the surface cleaning agents are often rotated in these facilities to minimize potential microbial resistance and/or selection toward the cleaning conditions (Murtough et al., 2001; Sandle, 2012).
For NASA cleanrooms (ISO class 8) where Mars and Europa spacecraft are assembled, the floors are routinely cleansed with ~0.8-1.6% v/v Kleenol 30 (Mogul et al., 2018; Danko et al., 2021), which is an alkaline detergent formulation. Due to the proprietary nature of Kleenol 30, scant information is available regarding the precise chemical composition. Available safety data sheets indicate a composition containing sodium dodecyl benzene sulfonate (1% w/w), polyethylene glycol mono-nonylphenyl ether (1–5% w/w), sodium metasilicate (1–5% w/w), ethylenediaminetetraacetic acid (EDTA; 2% w/w), sodium metasilicate (1–5% w/w), potassium hydroxide (KOH; ~2% v/v; presumably prepared from concentrated KOH), and 2-butoxyethanol (10–15% w/w). When considering the mechanism of cleansing for Kleenol 30, the combined chemical ingredients imply adjustments toward alkalinity (KOH), sequestering of transition and alkaline earth metals through chelation by the metal-binding reagents (EDTA and sodium metasilicate), and sequestering, emulsification, and removal of organics (non-polar, polar, and charged organics) by detergent action (polyethylene glycol mono-nonylphenyl ether, dodecyl benzene sulfonate, and sodium metasilicate).
Yet, despite these robust cleaning practices, spacecraft assembly facilities harbor a persistent yet low abundance microbial bioburden (Danko et al., 2021; Hendrickson et al., 2021). Further, molecular measures reveal the presence of a core microbiome across the spacecraft assembly facilities, International Space Station (Checinska et al., 2015; Checinska Sielaff et al., 2019), and clinical and operating facilities (Shams et al., 2016; Ellingson et al., 2020; Perry-Dow et al., 2022). In spacecraft assembly facilities, the microbiomes are reasonably diverse though generally low in total counts with measures of ~103–104 16S rRNA m−2 (intact cells), ~104–105 16S rRNA m−2 (total cells), 1–40 OTU m−2, ~102–104 colony forming units m−2, and ~ 101 spores m−2 (La Duc et al., 2012; Mahnert et al., 2015; Moissl-Eichinger et al., 2015).
Among the dominant members of the core microbiome across spacecraft and healthcare facilities are the Acinetobacter (La Duc et al., 2012; Mora et al., 2016; Hendrickson et al., 2021). The Acinetobacter are a Gram-negative bacterial genus associated with desiccation tolerance (McCoy et al., 2012; Farrow et al., 2018), radiation tolerance (La Duc et al., 2003; McCoy et al., 2012), bioemulsification (Pirog et al., 2021), biofilm formation (Yeom et al., 2013; Pompilio et al., 2021), and multi-drug resistance (Peleg et al., 2008).
In Mogul et al. (2018), we showed that Acinetobacter radioresistens 50v1, which was isolated from the surface of the pre-flight Mars Odyssey orbiter, could be cultivated in the presence of 1% v/v Kleenol 30, a floor detergent, while utilizing ethanol, a surface cleaner, as a sole carbon source. Additionally, cultivations of A. radioresistens 50v1 with Kleenol 30 were shown to yield tri, penta, and octaethylene glycols, which is suggestive of partial biodegradation of polyethylene glycol mono-nonylphenyl ether (a component of the Kleenol 30 formulation) via scission of the ether linkage (White, 1993) – to yield the mixed polyethylene glycols (e.g., tri, penta, and octaethylene glycol).
In addition, studies show that clinical strains of Acinetobacter baumanii tolerate benzalkonium chloride, which is a biocide and surface cleaner, through biofilm formation (Rajamohan et al., 2009) and increased expression of efflux pump genes (Fernández-Cuenca et al., 2015; Srinivasan et al., 2015). Further, survivability studies on cleanroom-associated Acinetobacter reveal extreme tolerances toward aqueous hydrogen peroxide, a surface and floor disinfectant; where the tolerances are perhaps the highest among Gram-negative bacteria (McCoy et al., 2012; Derecho et al., 2014; Mogul et al., 2018). These combined survival characteristics are likely key contributors to the Acinetobacter being among the dominant members of the spacecraft and healthcare microbiomes.
Therefore, to obtain molecular and quantitative insights into the tolerances of the cleanroom-associated Acinetobacter, we measured the impacts of Kleenol 30 on the survivals, cultivation kinetics, and intracellular and extracellular metabolomes. By design, we profiled Acinetobacter strains isolated from differing sub-locations within the NASA cleanrooms for spacecraft assembly (e.g., floor and spacecraft surface) which were subjected to differing cleaning regimes (e.g., Kleenol 30 for the floors vs. isopropyl alcohol and/or ethyl alcohol for spacecraft surfaces). Hence, comparisons across these sub-environments lend support toward the hypothesis that repeated cleansing under the conditions of spacecraft assembly serve as selective pressures that favor or promote biochemical tolerances toward the local cleaning conditions.
2. Materials and methods
2.1. Materials and conditions
Bacterial strains of A. johnsonii 2P08AA and A. radioresistens 50v1 were obtained through the Planetary Protection Culture Collection at the Jet Propulsion Laboratory. Stocks of ethanol (ENG Scientific) and Kleenol 30 (Mission Laboratories, Los Angeles, CA; Clovis Janitorial) were sterile filtered and stored as 200 μL aliquots at 4°C. Stock solutions of 25 mM Fe2+ were prepared using ferrous ammonium sulfate (Fe(NH4)2(SO4)2·6H2O; EM Science) in ultrapure water (18 MΩ−1), followed by sterile filtration (0.2 μm syringe filter, VWR), and storage as 200 μL aliquots at 4°C (with no visible precipitation over long term storage).
Concentrated minimal media (5x MM) were prepared using 30.0 g Na2HPO4·7H2O (Sigma), 15.0 g KH2PO4 (EM Science), 2.50 g NaCl (Fisher Scientific), 5.00 g NH4Cl (EM Science) per liter of ultrapure water. Low-osmolarity media (0.2x M9) were prepared by adding 0.4929 g MgSO4·7H2O (EM Science) and 0.0147 g CaCl2·2H2O (EM Science) to 200.0 mL 5x MM and dilution to 1.00 L using ultrapure water (to yield 1x M9), followed by an additional 5-fold dilution (200 mL 1x M9 in 1.00 L) to yield 0.2x M9. Final concentrations for ions in 0.2x M9 were 20 μM Ca2+, 400 μM Mg2+, 4.4 mM K+, 10.7 mM Na+, and 3.7 mM NH4+, along with 5.4 mM Cl−, 4.5 mM H2PO4−, 4.4 mM HPO42−, and 0.4 mM SO42−.
Lysogeny broth media (LB) were prepared using 10.0 g tryptone (VWR), 5.0 g yeast extract (Amresco), 5.0 g NaCl (Fisher Scientific) and 1 mL of 1 M NaOH (Sigma-Aldrich) per liter of ultrapure water. Agar plates were prepared using 1.0 L LB and 15 g agar (Amresco).
Cultivations were performed in 15 mL screw cap test tubes (Pyrex, borosilicate glass; 13 × 100 mm), which were washed with tap water, rinsed with distilled water, and autoclaved between experiments. All screw cap test tubes were capped tightly and wrapped with parafilm during cultivation to prevent loss of ethanol, which served as a sole supplied carbon source. Cultivations were performed with mild agitation (200 rpm) using a New Brunswick Scientific Innova 4,200 incubator. Cell densities were monitored during cultivation by following agar plate counts and changes in optical density (OD) at 600 nm (Spectronic 20 Genesys). All media were autoclaved at 121°C and 15 psi for 45 min.
2.2. Cultivations with Kleenol 30
Glycerol stocks of A. radioresistens 50v1 and A. johnsonii 2P08AA were separately streaked onto LB agar plates and incubated at 32°C for ~24 h. Isolated colonies were inoculated into 2.0 mL 0.2x M9 containing 25 μM Fe2+ and 0–2.0% v/v Kleenol 30 (or 0–40.0 μL). Pre-cultures were initiated by addition of 20.0 μL 90% v/v ethanol to yield final concentrations of 150 mM or 1.0% v/v ethanol. Pre-cultures were grown to late log phase (OD ~0.6 at ~12 h). Target cultures were prepared using fresh media (2.00 mL) which were inoculated with 20.0 μL (1:100 dilution) of the respective pre-cultures and initiated by addition of 1.0% v/v ethanol.
Temporal changes in OD were followed for 0–15 h. Growth kinetics were characterized by regression analysis (Microsoft Excel) using a modified version of the Gompertz equation (Begot et al., 1996), which describes a non-linear bacterial growth model, and yields the parameters of growth rate (k), lag time (L), and an estimate of the maximum change in relative biomass ( ). All regressions were minimized by least squares analysis. Control cultures containing no inoculate showed no growth (OD ≤ 0.002), as did cultures containing inoculate but no ethanol (OD ≤ 0.002), which indicated negligible biological contamination and accumulation of abiotic particles during cultivation.
Survivals of late-log phase cultures (OD ~0.4) were assessed by plating onto LB agar plates. In control experiments, plate counts using 0.2x M9 agar plates supplemented with ethanol (just prior to use) yielded irreproducible results, when compared to LB agar plate, likely due to variances in adsorption of ethanol (under our conditions). For the plate count assays, therefore, aliquots (100 μL) of the cultures (in 0.2x M9/Fe) were transferred to 2.5 mL microcentrifuge tubes, decimally diluted by 106-fold using 0.2x M9, and spread (20 μL) onto LB agar plates using sterile plastic cell spreaders. All plates were sealed with parafilm and incubated at 32°C for 24 h. Plates with ≤300 colonies were enumerated and expressed as colony forming using per mL of the parent culture (cfu mL−1).
2.3. Metabolomics of Acinetobacter cultivations in Kleenol 30
Cultures of A. radioresistens 50v1 and A. johnsonii 2P08AA were prepared as described and harvested during early stationary phase (OD ~0.6–0.7, 7–8 h) by centrifugation (3,500 rpm) at 4°C for 15 min (Beckman Coulter Allegra 21R). After centrifugation, the cell pellets and supernatant fractions were separated and, respectively, treated. The supernatants were saved as 500 μL aliquots, dried using a DNA 110 Savant DNA SpeedVac, and stored at −80°C. The pelleted cells were washed twice with 0.2x M9 and stored at −80°C.
Samples were characterized by the West Coast Metabolomics Center using gas chromatography and time-of-flight mass spectrometry (GC-TOF/MS). In brief, dried cells and culture broth were separately resuspended in 2 mL of pre-chilled (−20°C) and degassed extraction solvent (acetonitrile:isopropanol:water, 3:3:2), vortexed for 30 s, shaken for 5 min at 4°C, clarified by centrifugation (~12,000 g), and the resulting supernatant evaporated to dryness. Samples were derivatized by resuspension in 10 μL of 40 mg/mL methoxyamine hydrochloride in pyridine (30°C, 1.5 h) followed by addition of 41 μL N-methyl-N-(trimethylsilyl) trifluoroacetamide (80°C, 30 min); fatty acid methyl esters (e.g., C8-C30) were additionally added to serve as retention index markers (Barupal et al., 2019).
Samples were transferred to crimp top vials and separated on an Agilent 6,890 Gas Chromatograph equipped with a Gerstel automatic liner exchange system (ALEX), multipurpose sample (MPS2) dual rail, Gerstel CIS cold injection system (Gerstel, Muehlheim, Germany), and built-in gas purifier (Airgas, Radnor PA). Chromatographic separation was afforded using a 30 m (0.25 mm i.d.) Rtx-5Sil MS column (0.25 μm 95% dimethyl 5% diphenyl polysiloxane film) with an additional 10 m integrated guard column (Restek, Bellefonte PA). Carrier gas was 99.9999% pure Helium with a constant flow of 1 mL/min.
Sample volumes of 0.5 μL were injected with a 10 μL s−1 injection speed on a spitless injector with a purge time of 25 s. Temperature profile included 1 min at 50°C for the oven temperature, an increase of 20°C min−1 to 330°C across 14 min, with a final hold for 5 min at 330°C. A temperature of 280°C was used for the transfer line between the gas chromatograph and Leco Pegasus IV time of flight mass spectrometer (single mass analyzer). The measured mass range was 85–500 Da, scan rate was at 17 spectra s−1, electron energy was 70 eV, ion source temperature was 250°C, and detector voltage was 1850 V. Spectra were acquired using the Leco ChromaTOF software vs. 2.32 (St. Joseph, MI) and processed using the BinBase database system (Fiehn et al., 2005, 2010; Trigg et al., 2019), which quantifies (signal-to-noise ratio of 5:1) and matches mass peaks to the Fiehn mass spectral library using retention index (RI) and reference mass spectral information (which are internally compiled by the West Coast Metabolomic Center). Metabolomic data (Study ID ST002380, DatatrackID: 3552) were stored at the NIH Common Fund’s National Metabolomics Data Repository (Sud et al., 2016).
2.4. Statistical analyses
Comparisons of the cultivation kinetic parameters were conducted using unpaired Student’s t-tests and one- and two-way ANOVA analyses (Microsoft Excel), where statistical relevance was accepted at p < 0.05, normal distribution of the data was supported by Shapiro–Wilk tests (Past 4.10; Hammer et al., 2001), and variances were assumed as equal. Metabolomic data were compared using univariate, multivariate, and visual approaches.
Discrete changes in the metabolomes (e.g., trends for a single metabolite) were identified using unpaired Student’s t-tests (Microsoft Excel) with corrections for multiple testing using a false discovery rate (FDR) of FDR ≤ 0.20 (Benjamini and Hochberg, 1995). For the 700 metabolites in the total study (whole cell and extracellular metabolomes from 2 bacterial strains treated with and without Kleenol 30), normal distributions (Shapiro–Wilk Tests) were indicated for ~92% (646) of the metabolites. To account for potential underestimations of normality given the sample size (n = 3), metabolites exhibiting thresholds for normality of p > 0.03 (Shapiro–Wilk Tests) and significance of p < 0.05 (Student’s t-tests) were carried forward in the univariate and visual assessments of the data.
Broad structural and metabolic trends were obtained by visualizing the changes in the metabolomes (p < 0.05) using MetaMapp1 and Cytoscape 3.9.12 (Barupal et al., 2012), which constitute a statistical organizational approach to yield visual maps of metabolites arranged by known structural patterns and metabolic pathways. Confirmation of broad changes in the metabolomes were obtained using ChemRich3 (Barupal and Fiehn, 2017), a statistical enrichment tool that compares groups of metabolites based on structural and biochemical classes. For this study, the standard and user-defined classifications in ChemRich included amino acids, monosaccharides, sugar acids, organic acids, fatty acids, lipid-related metabolites, nucleotide-related metabolites, and compatible solutes. Multivariate tests were conducted using canonical correspondence analyses (Past 4.10; Hammer et al., 2001) to correlate the metabolomic trends to differing cultivation growth parameters and conditions.
3. Results
3.1. Cultivation and survival in the presence of Kleenol 30
In Figure 1A we show that the cleanroom-associated Acinetobacter tested in this study, A. johnsonii 2P08AA and A. radioresistens 50v1, readily grow in the presence of ≤1% v/v Kleenol 30 (K30) under low-osmolarity aqueous conditions with ethanol serving as the sole supplied carbon source. Cultivations on 2% v/v K30 did not yield measurable cell densities.
Figure 1. Impacts of 0–1.0% v/v Kleenol 30 (K30) on (A) the survival and (B–D) cultivation kinetics of Acinetobacter johnsonii 2P08AA (red) and Acinetobacter radioresistens 50v1 (blue). Survival (N/N⊝) is expressed as the ratio of plate counts (cfu mL−1) from cultures grown in the absence of K30 (N⊝) and presence of K30 (N). Fitted regressions of growth curves yielded changes in (B) growth rates (h−1), (C) lag time (h), and (D) maximum relative biomass (log (N/N0)), which is assumed to represent a ratio of the total biomass (N) at stationary phase and the biomass at the start of the culture (N0). Error bars represent the standard errors (n = 4–5, growth curves; n = 3–5, plate counts). Univariate tests are represented as asterisks (*) for t-tests with p < 0.05, hashtags (#) for t-tests with p ≥ 0.05, and asterisks (*) with an underlying line for one or two-way ANOVA with p < 0.05.
In 0.1% v/v K30, cultures of the floor-associated A. johnsonii 2P08AA show no apparent loss in survival, as plate counts (~late log phase) effectively show no difference (p > 0.05) in the presence (2.1 ± 0.5 × 107 cfu mL−1) and absence (2.6 ± 0.1 × 107 cfu mL−1) of the detergent. These trends indicate quantitative survival of the 2P08AA strain in the presence of 0.1% v/v K30.
For the spacecraft surface-associated A. radioresistens 50v1, in contrast, survivability in 0.1% v/v K30 readily decreases (p = 0.006) to 46 ± 6% (propagated error) – as calculated by comparison of cultures cultivated in the presence (4.0 ± 0.4 × 107 cfu mL−1) and absence (8.8 ± 0.7 × 107 cfu mL−1) of 0.1% v/v K30. These trends are indicative of ~50% survival in the presence of 0.1% v/v K30.
In 1.0% v/v K30, survivabilities for the 2P08AA strain (26 ± 12%) and 50v1 strain (22 ± 8%) decrease to similar values – as indicated by comparison of cultures cultivated in the presence (2P08AA, 0.7 ± 0.3 × 107 cfu mL−1; 50v1, 1.9 ± 1.1 × 107 cfu mL−1) and absence (2P08AA, 2.6 ± 0.1 × 107 cfu mL−1; 50v1, 8.8 ± 0.7 × 107 cfu mL−1) of 1.0% v/v K30. These trends are indicative of ~70–80% inhibition in the presence of 1.0% v/v K30 across both strains.
3.2. Impacts of Kleenol 30 on the cultivation kinetics
Displayed in Figures 1B–D and Figure 2 are the cultivation kinetic parameters and associated growth curves for A. johnsonii 2P08AA and A. radioresistens 50v1. Optical measurements (Figure 2, empty blue circles) for cultures were converted to the ratiometric changes in relative biomass ( ) by assuming a direct relationship between optical transmittance and cell density for the initial culture media ( ) and at the time of measurement ( ), as outlined in Begot et al. (1996). Changes in relative biomass (Figure 2, filled red circles, fitted line) were expressed as a log function ( ), plotted over time, and fit to a modified version of the Gompertz equation (Eq. 1), as detailed in Begot et al. (1996) and Mogul et al. (2018). Minimized regressions yielded the parameters of growth rate (k; cell divisions h−1), estimated time in lag phase (L; h), and log of the maximum increase in relative biomass ( ; unitless dimension).
Figure 2. Measured and fitted growth curves for A. johnsonii 2P08AA (left panel) and A. radioresistens 50v1 (right panel) cultured in the presence of 0–1.0% v/v Kleenol 30 at 28°C in 0.2x M9 media containing 25 μM Fe2+ and 150 mM ethanol (the sole supplied carbon source). Optical density measurements (blue empty circles; right-hand y-axis) were converted to log(N/N0) (red filled circles; left-hand y-axis) as described and fit by non-linear regression to yield the parameters of time in lag phase, growth rates, and maximum ratiometric and logarithmic change in biomass (log (N/N0)) at stationary phase.
To account for optical scattering, the maximum increases in relative biomass ( ) were broadly interpreted as maximum increases in intact cells, cellular aggregation, and/or intra- and extracellular polymeric substances at stationary phase. Reported in this study are the averaged values from the minimized regressions along with standard errors (n = 3–5 biological replicates).
Comparisons of the growth rates (k) show ~18% decreases for cultivations of the 2P08AA strain in 0.1% v/v K30 (p = 0.029) – as calculated by comparison of the difference in growth rates (Figure 1B) between the absence (0.50 ± 0.01 h−1) and presence (0.41 ± 0.02 h−1) of 0.1% v/v K30. The 50v1 strain, in contrast, shows no change in growth rate (p = 0.05) in 0.1% v/v K30 – as calculated by comparison of rates (Figure 1B) in the absence (0.39 ± 0.01 h−1) and presence (0.40 ± 0.02 h−1) of 0.1% v/v K30. When cultivated in 1.0% v/v K30, growth rates for the 2P08AA strain (0.22 ± 0.03 h−1) and 50v1 strain (0.35 ± 0.01 h−1) decrease by ~56% (p = 0.001) and ~ 13% (p = 0.043), respectively, as calculated by comparison rates (Figure 1B) in the absence (2P08AA, 0.50 ± 0.01 h−1; and 50v1, 0.39 ± 0.01 h−1) of the detergent.
Comparisons by one factor ANOVA analyses confirm that cultivations in the presence of 0.1–1.0% v/v K30 significantly impact the respective growth rates for the 2P08AA (p = 0.038, f = 4.256) and 50v1 (p = 0.0002, f = 4.459) strains. Additionally, two factor ANOVA analyses confirm that cultivations in the presence of K30 differentially impact the growth rates of the 2P08AA and 50v1 strains (p = 0.015, f = 5.318). These combined trends for exponential phase behavior are indicative of higher tolerances toward K30 by the 50v1 strain at 0.1 and 1.0% K30. These trends are opposite to those from plate counts at late-log phase, which indicate a higher tolerance for the 2P08AA strain.
Comparisons of the lag times (L) in 0.1 and 1.0% v/v K30 show decreases for both strains during cultivation (Figure 1C), which is unexpected since the addition of detergents was presumed to inhibit growth and increase the time in lag phase. For the 2P08AA strain, the native lag time of 0.77 ± 0.23 h−1 (absence of K30) reduces to negligible values (~0 h) in 0.1 and 1.0% v/v K30. For the 50v1 strain, the native lag time of 2.61 ± 0.24 h−1 is ~3-fold longer than 2P08AA strain and decreases by ~70% in 0.1% v/v K30 (0.76 ± 0.15 h−1) and to negligible values in 1.0% v/v K30.
Comparisons by one factor ANOVA analyses confirm that cultivations in the presence of 0.1–1.0% v/v K30 significantly impact the respective lag times for the 2P08AA (p = 0.00008, f = 4.256) and 50v1 (p = 0.015, f = 4.459) strains. Additionally, two factor ANOVA analyses confirm that cultivations in the presence of K30 differentially impact the lag times for the 2P08AA and 50v1 strains (p = 0.006, f = 5.318). These combined trends indicate that K30 induces accelerated entry into the exponential phase in a concentration dependent manner, where the lag phase for the 2P08AA strain is effectively eliminated at ~10-fold lower K30 abundances when compared to the 50v1 strain – under these conditions.
Comparisons of the maximum change in relative biomass ( ) for cultivations in 0.1% v/v K30 show unexpected increases for both strains (Figure 1D). For the 2P08AA strain, cultivations in 0.1% v/v K30 yield ~40% increases (p = 0.026) in maximum apparent biomass when comparing values obtained in absence (2.30 ± 0.02) and presence of (2.44 ± 0.04) of 0.1% v/v K30 and accounting for the log transformation. Similarly, for the 50v1 strain, cultivations in 0.1% v/v K30 yield ~55% increases (p = 0.048) in maximum apparent biomass when comparing values obtained in absence (2.28 ± 0.01) and presence of (2.47 ± 0.04) of 0.1% v/v K30 and accounting for the log transformation. When cultivated in 1.0% v/v K30, the 2P08AA strain (1.53 ± 0.18) and 50v1 strain (2.07 ± 0.03) exhibit decreases of ~80% and ~ 40% in maximum apparent biomass, respectively, which is consistent with the decreased growth rates under these conditions (~56 and ~ 13%, respectively).
Comparisons by one factor ANOVA analyses confirm that cultivations in the presence of 0.1–1.0% v/v K30 significantly impact the respective maximum changes in relative biomass for the 2P08AA (p = 0.000003, f = 4.256) and 50v1 strains (p = 0.000004, f = 4.459). Additionally, two factor ANOVA analyses confirm that cultivations in the presence of K30 differentially impact the maximum changes in relative biomass for the 2P08AA and 50v1 strains (p = 0.007, f = 5.318). Combined, these trends suggest that 0.1% v/v K30 induces changes in apparent biomass in stationary phase.
3.3. Impacts of Kleenol 30 on the intracellular metabolomes
3.3.1. Metabolomic considerations
Metabolomes from whole cells of A. johnsonii 2P08AA and A. radioresistens 50v1 were compared after cultivation in the absence and presence of 0.1% v/v Kleenol 30 (early stationary phase). Per strain, whole cell extracts yielded a total of 119 intracellular metabolites and free metabolites from the membrane (e.g., free fatty acids and monoacylglycerols). Metabolites from both strains showing important changes (p < 0.05) due to cultivation with K30 are listed in Table 1, as are changes (fold-changes) in the respective abundances; changes retaining significance (p < 0.05, FDR ≤ 0.20) after correction for multiple testing are underlined.
Table 1. Numbered list of metabolites from whole cell extracts of Acinetobacter johnsonii 2P08AA and Acinetobacter radioresistens 50v1, where the magnitude of decreases (↓) or increases (↑) in the relative abundances after cultivation in 0.1% v/v Kleenol 30 (p < 0.05) are listed in differing columns and presented in heat-map type coloring to visually represent the magnitude of decreasing (↓: light to dark blue) or increasing (↑: light green to dark red) abundances; the column label of ID# represents the metabolites numbered in Figure 3, CID is the associated PubChem Compound ID, RI is the associated retention index, m/z is the mass to charge ratio for the associated derivatized metabolite, fold-changes represent the ratio of abundances obtained in presence and absence of the K30, asterisks (*) indicate relevance after correction for multiple testing (FDR ≤ 0.20), dashes indicate that the metabolite was not detected, and gray boxes indicate that the option is not applicable.
Displayed in Figure 3 are metabolomic maps for the 2P08AA and 50v1 strains, which readily show the changes resulting from cultivation in K30. The metabolomic maps (MetaMapp graphs) in Figure 3 are organized via networks of KEGG reactant pairs (black edges or arrows) and Tanimoto chemical similarity scores (gray edges or arrows). Detected metabolites are presented as nodes (circles). Important changes are highlighted as red nodes (increases in relative abundance), blue nodes (decreases in relative abundance), or yellow nodes (no change). Node sizes signify the relative degree of change (or fold-change).
Figure 3. Maps showing the impacts of cultivation in 0.1% v/v K30 on the intracellular and extracellular metabolomes of A. johnsonii 2P08AA (A–C) and A. radioresistens 50v1 (B,D). Metabolites are represented as nodes and organized by structural and metabolic connections (MetaMapp and Cytoscape 3.9.1). Important changes in abundance (p < 0.05) are represented as red (increases), blue (decreases), and yellow (no change) nodes, where node sizes represent the degree of change. Relevant metabolites that show important (black text; p < 0.05) and significant changes (red text; p < 0.05; FDR ≤ 0.20), as listed in Tables 1, 2, are displayed using the associated identity numbers (ID#). Shaded circles represent biochemical classes (red text) that show significant changes (ChemRich); biochemical classes showing no change are listed for comparison purposes (black text).
For the 2P08AA strain (Figure 3A), cultivations in 0.1% v/v K30 yield no statistically discernable impacts to the whole cell metabolome. No changes are observed after accounting for multiple testing (FDR ≤ 0.20, n = 119). No changes are observed when parsing the data using ChemRich. While potential increases (p < 0.05) are observed in gluconate, analyses in ChemRich provide no support for broad changes in sugar acid content. Despite the potential (p < 0.05) decreases in arabinose and increases in sorbitol, analyses in ChemRich show no support for broad changes in monosaccharide composition. Similarly, ChemRich analyses provide no support for broad changes in lipid-related metabolites despite the potential (p < 0.05) decreases in phosphoethanolamine and increases in ethanolamine and behenic acid (22:0; long chain saturated fatty acid).
In contrast, the 50v1 strain displays multiple changes in the intracellular metabolome due to cultivation in 0.1% v/v K30 (Table 1). After correction for multiple testing (n = 119), significant decreases are observed for lysine, homoserine, ornithine, mannitol, adenosine-5′-monophosphate (5’-AMP), 5′-methylthioadenosine (5′-MTA), and 2,5-dihydroxypyrazine. The whole cell metabolomic maps for the 50v1 strain in Figure 3B were parsed using MetaMapp and ChemRich and expand on the univariate analyses to reveal differential changes in amino acids, compatible solutes, nucleotide-related metabolites, dicarboxylic acids, sugar acids, and saturated fatty acids.
Impacts to amino acids (p = 9.1×10−6; Q = 6.7 × 10−5; ChemRich) are supported by decreases in the relative abundances of standard amino acids (glycine, lysine, proline, serine, and valine), non-standard amino acids (β-alanine, homoserine, ornithine, and citrulline), and a peptide (alanine–alanine). Dimerization of glycine under the analytical conditions is suggested by detection of 2,5-dihydroxypyrazine (Haffenden and Yaylayan, 2005) with the observed decreases in 2,5-dihydroxypyrazine in cultures with K30 being generally consistent with the decreases in glycine.
Impacts to compatible solutes (p = 9.6×10−6; Q = 6.7 × 10−5; ChemRich) are supported by decreases in mannitol (Empadinhas and da Costa, 2008; Tam et al., 2022), as well as glycerate and 3-phosphoglycerate, which are involved in the biosynthesis of compatible solutes (Franceus et al., 2017). Impacts to nucleotide-related metabolites (p = 5.2×10−4; Q = 2.4 × 10−3; ChemRich) are supported by decreases in 5’-AMP, MTA, pyrophosphate, thymidine, and uracil, along with increases in cytidine-5′-monophosphate (5’-CMP). Impacts to dicarboxylic acids (p = 4.8×10−3; Q = 1.7 × 10−2; ChemRich) are supported by increases in fumarate, malate, and tartrate. No impacts to saturated fatty acids are observed (p = 0.11; Q = 0.31; ChemRich) despite the decreases in isopentadecanoic acid (branched saturated fatty acid; 13-methyl-14:0, or 15:0 iso) and increases in stearic acid (saturated fatty; 18:0).
3.3.2. Metabolic considerations
Cross referencing of the metabolomic changes to the KEGG database yields probable impacts to several intracellular pathways for the 50v1 strain. Provided in parentheses are the associated KEGG pathway or reaction identifiers.
Inhibition or reduction in Glycine, Serine, and Threonine Metabolism (map00260) is supported by decreases in intracellular glycerate, 3-phosphoglycerate, serine, glycine, and homoserine. Inhibition or reduction in Lysine Biosynthesis (map00300) and Valine, Leucine, and Isoleucine Biosynthesis (map00290) is supported by decreases in lysine and valine; however, no other changes are noted in the respective pathways.
Inhibition or reduction in proline and arginine biosynthesis through Arginine and Proline Metabolism (map00330) and Arginine Biosynthesis (map00220) are supported by decreases in ornithine, citrulline, and proline. The lack of changes in abundances for urea indicate no measurable impact on the flux of the urea cycle – despite the decreases in ornithine and citrulline. The lack of observed arginine is suggestive of conversion of arginine to ornithine under the analytical conditions (Halket et al., 2005) or the limited presence of free arginine in the cell.
Inhibition or reduction in Peptidoglycan Biosynthesis (map00550) and/or peptide (protein) synthesis is suggested by decreases in alanine–alanine. Inhibition or reduction in β-alanine metabolism (map00410) is supported by decreases in β-alanine, which in turn is suggestive of reductions in arginine and uracil degradation. Differential changes in fatty acid metabolism (map01212) are suggested by the decreases in free isopentadecanoic acid (Reaction R02663) and increases in free stearic acid.
Inhibition or altered flux through glycolysis (map00010) is suggested through decreases in 3-phosphoglycerate. Activation or altered flux through components of the TCA cycle are suggested by increases in fumarate and malate. Inhibition or reduced flux of selected monosaccharides through Fructose and Mannose Metabolism (map00051) and Pentose and Glucuronate Interconversions (map00040) are, respectively, suggested by decreases in mannitol and arabinose.
Differential changes in Nucleotide Metabolism (map01232) are suggested by increases in 5’-CMP and 5′-MTA, along with decreases in 5’-AMP, pyrophosphate, thymidine, and uracil. Additionally, inhibition or reduction in Pyrimidine Metabolism (map00240) is suggested decreases by uracil and thymidine.
3.4. Impacts of Kleenol 30 on the extracellular metabolomes
Displayed in Figure 3 are biochemical maps that compare the extracellular metabolomes from A. johnsonii 2P08AA and A. radioresistens 50v1 (early stationary phase). Per strain, dried cultivation broth (cell free) yielded a total of 56 metabolites, which were processed and visualized as described above. Listed in Table 2 are the extracellular metabolomic changes in the 2P08AA and 50v1 strains after cultivation in the absence and presence 0.1% v/v K30. Pathway analyses were not conducted for extracellular metabolites.
Table 2. Numbered list of extracellular metabolites from cultures of A. johnsonii 2P08AA and A. radioresistens 50v1, where the magnitude of decreases (↓) or increases (↑) in the relative abundances after cultivation in 0.1% v/v Kleenol 30 (p < 0.05) are listed in differing columns and presented in heat-map type coloring to visually represent the magnitude of decreases (↓: light to dark blue) or increases (↑: light green to dark red); the ID# represents the metabolite number from Figure 3 and continues with the numbering order from Table 1, CID is the PubChem Compound ID, RI is the associated retention index, m/z is the mass to charge ratio for the associated derivatized metabolite, fold-changes represent the ratio of abundances obtained in presence and absence of the K30, asterisks (*) indicate relevance after correction for multiple testing (FDR ≤ 0.20), dashes indicate that the metabolite was not detected, gray boxes indicate that the option is not applicable.
When grouped by biochemical classes, cultivations of the 2P08AA strain in K30 yield differential impacts to extracellular sugar acids, monosaccharides, and organic acids (di- and monocarboxylic acids). Yet, univariate tests show no changes after correction for multiple testing. Impacts to sugar acids (p = 1.4×10−5; Q = 6.9×10−5; ChemRich) are supported by substantial increases in extracellular gluconolactone (22-fold) and glycerate (34-fold). Impacts to di- and monocarboxylic acids, which were grouped as organic acids (p = 3.0×10−4; Q = 7.6×10−4; ChemRich), are supported by increases in oxalate and decreases in adipic acid (hexanedioic acid) and 3,4-dihydroxybenzoate (protocatechuate). Impacts to extracellular monosaccharides (p = 7.6×10−3; Q = 1.3×10−2; ChemRich) are supported by substantial increases in extracellular fructose (4.5-fold), galactose (22-fold), and mannose (57-fold) and decrease (2.5-fold) in extracellular ribose (pentose).
In contrast, cultivations of the 50v1 strain in K30 – after correction for multiple testing (n = 56) – yield significant changes in the extracellular galactonate, gluconate, gluconolactone, glycerate, mannose, ribose, tartrate, 4-hydroxybutyrate, and α-ketoglutarate. Likewise, grouping by biochemical classes are reveals significant changes in the extracellular metabolome.
Impacts to sugar acids (p = 1.9×10−8; Q = 9.6×10−8; ChemRich) are supported by substantial increases in extracellular galactonate (36-fold), gluconate (46-fold), gluconolactone (8.8-fold), glycerate (9.3-fold), and 2-methylglycerate (1.9-fold). Impacts to extracellular organic acids (p = 7.8×10−2; Q = 0.16; ChemRich) are suggested by the significant increases in tartrate and significant decreases in α-ketoglutatate, along with the potential decreases in citramalate and 4-hydroxybutyrate. Impacts to extracellular monosaccharides (p = 9.6×10−2; Q = 0.16; ChemRich) are suggested by significant increases in extracellular mannose (20-fold) and significant decreases in ribose (2.5-fold).
3.5. Impacts of Kleenol 30 on cultivation and the intracellular metabolome
3.5.1. Canonical correlation analyses
In Figure 4, we use canonical correlation analyses (CCA) to, respectively, compare metabolomes from A. johnsonii 2P08AA (n = 3) and A. radioresistens 50v1 (n = 3) to differing quantitative descriptors characterizing changes in the cultivation media and growth kinetics in the absence and presence of 0.1% v/v K30. Correlations for the intracellular (or whole cell) metabolomes from late log-phase cultures were assessed against growth rates at exponential phase, plate counts at late-log phase, maximum relative biomass at stationary phase, and total detergent concentrations in the media, which were assumed to be 0.025% w/w in cultures containing 0.1% v/v K30.
Figure 4. Canonical correspondence analyses (CCA) that compare the impacts of cultivation in 0.1% v/v Kleenol 30 (K30) on (A) the intracellular metabolome of A. johnsonii 2P08AA, (B) the extracellular metabolome of A. johnsonii 2P08AA, (C) the intracellular metabolome of A. radioresistens 50v1, and (D) the extracellular metabolome of A. radioresistens 50v1 against the quantitative descriptors of growth rates at exponential phase (blue callout boxes, green vectors), plate counts at late-log phase (purple callout boxes, green vectors), maximum relative biomass at stationary phase (brown callout boxes, green vectors), and either the total concentrations of detergents (0.025% w/w; assumed) or chelators (0.02% w/w; assumed) in the cultivation media (black callout boxes, green vectors); dimension reduction of the triplicate metabolomic measures associated with 0% v/v K30 (red squares, pink circles, and red and pink diamonds) and 0.1% v/v K30 (black triangles) are provided, reduced terms are highlighted by the green and red circles and callout boxes, and metabolites are arrayed as text.
The CCA plots include the total array of detected metabolites (text), listed quantitative descriptors (vector lines), and transformed data after dimension reduction (glyphs). For the 2P08AA and 50v1 strains, stronger correlations along CCA Axis 1 (61.70, 69.46%) are observed when compared to CCA Axis 2 (24.72, 23.70%), respectively. Dimension reduction for the metabolomes from the 2P08AA and 50v1 strains (Figures 4A,B) respectively show clear separations between samples prepared in the absence and presence of 0.1% v/v K30. These trends suggest that the intracellular metabolomic compositions of both strains are impacted by K30 during cultivation. Described in the following sub-sections are comparisons to the quantitative descriptors. For selected metabolites, additional clarifications and/or fatty acid abbreviations are provided in parentheses.
3.5.2. Comparisons to survival
When considering changes in survival (Figure 4A), the vector describing plate counts for the 2P08AA strain (late log phase) overlaps with the reduced terms associated with 0% v/v K30 (filled red squares). This correlation suggests that survival trends for the 2P08AA strain (across 0 to 0.1% v/v K30) are associated with limited changes in the intracellular metabolome. This assessment is consistent with the associated univariate tests which indicate no changes in plate counts (p > 0.05) or metabolite abundances for the 2P08AA strain in 0.1% v/v K30 (FDR > 0.20).
From the CCA plot for the 2P08AA strain (Figure 4A), metabolites trending alongside the plate counts include trehalose, trehalose-6-phosphate, fructose-6-phosphate, and potentially glucose-6-phosphate. Univariate tests infer no changes in the abundances for these metabolites (p > 0.05 or FDR > 0.20). These trends suggest that the quantitative survival for the 2P08AA strain (in part) relates to lack of changes in central metabolism and compatible solute formation.
For the 50v1 strain (Figure 4B), the vector describing plate counts overlaps with the reduced terms associated with 0% v/v K30 (pink circles). This correlation suggests that survival trends for the 50v1 strain (across 0 to 0.1% v/v K30) are associated with limited changes in the intracellular metabolome. This assessment is contrary to the associated univariate tests that support ~50% reductions in survival (p < 0.05) and substantial changes in the metabolome (p < 0.05, FDR ≤ 0.20).
From the CCA plot for the 50v1 strain (Figure 4B), metabolites trending alongside the plate counts include 4-hydroxybenzoate, several amino acids (similar to those in Table 1), ribose-5-phosphate, N-acetylglucosamine, nicotinic acid, phenylacetic acid, bisphosphoglyercol, and potentially 5’-AMP and 5′-MTA. These trends suggest that survival of the 50v1 strain is related to changes in central metabolism, amino acid metabolism, nucleotide metabolism, and cell wall metabolism (e.g., metabolism of glycosoaminoglycans and peptidoglycans, which are based on N-acetylglucosamine).
Trends from the CCA plot for the 50v1 strain (Figure 4B) are additionally suggestive of changes in native benzene metabolism (e.g., phenylacetic acid and 4-hydroxybenzoate). However, univariate analyses yield no support (p > 0.05) for the apparent decreases in phenylacetic acid (~0.1-fold decreases) and 4-hydroxybenzoate (~0.9-fold decreases). Alternatively, described further in Section 3.5.3, phenylacetic acid and 4-hydroxybenzoate are potential biodegraded products of K30. Hence, trends for the 50v1 strain hint at a potential for biodegradation of K30 through benzene metabolism.
3.5.3. Comparisons to maximum relative biomass
When considering biomass, the vector describing the maximum relative biomass at stationary phase for the 2P08AA strain (Figure 4A) overlaps with the reduced terms associated with 0.1% v/v K30 (black triangles). This correlation suggests that trends in the maximum relative biomass for the 2P08AA strain (across 0 to 0.1% v/v K30) are associated with changes in the intracellular metabolome. This assessment is consistent with the associated univariate tests for the cultivation data that support increases in the maximum relative biomass (p < 0.05), yet counter to the associated univariate tests for the metabolomic data that indicate no change in the relative metabolite abundances after correction for multiple testing (FDR > 0.20).
From the CCA plot for the 2P08AA strain (Figure 4A), the metabolites trending alongside the maximum relative biomass include methionine, lignoceric acid (fatty acid; 23:0), oxoproline, and glucose-1-phosphate. Univariate tests indicate no change in the relative abundances for these metabolites (FDR > 0.20). These combined trends suggest that increases in maximum relative biomass relate to minimal or no changes in glucose metabolism, amino acid metabolism, and unsaturated fatty acid abundances.
For the 50v1 strain (Figure 4B), the vector describing the maximum relative biomass overlaps with the reduced terms associated with 0.1% v/v K30 (black triangles). This correlation suggests that trends in the maximum relative biomass for the 50v1 strain (across 0 to 0.1% v/v K30) are associated with changes in the intracellular metabolome. This assessment is consistent with the associated univariate tests on the cultivation data (p < 0.05) and metabolomic data (p < 0.05, FDR ≤ 0.20), which yield support for increases in the maximum relative biomass and differential changes in metabolite abundances. In contrast to the 2P08AA strain, these trends for 50v1 strain suggest that increases in the maximum relative biomass are associated with several intracellular metabolomic changes.
From the CCA plot for the 50v1 strain (Figure 4B), metabolites trending alongside the increases in maximum relative biomass include lignoceric acid (fatty acid; 23:0), tryptophan, xanthine, stearic acid (fatty acid; 18:0), α-ketoglutarate, phosphate, 4-hydroxybutanoate (4-hydroxybutyric acid), 4-aminobutanoate (4-aminobutyric acid), O-phosphoserine, 1-monopalmitin (1-palmitoylglycerol; palmitoyl = 16:0), behenic acid (fatty acid; 22:0), and 2-monoolein (2-oleoylglycerol; oleoyl = 18:1Δ9). These trends suggest that increases in the maximum relative biomass for the 50v1 strain are related (at the minimum) to changes in butanoate metabolism (e.g., 4-hydroxybutanoate and 4-aminobutanoate), amino acid metabolism (e.g., tryptophan, α-ketoglutarate, and O-phosphoserine) and phosphate metabolism (phosphate and O-phosphoserine).
The trends for the 50v1 strain (Figure 4B) also suggest a relation to changes in free fatty acids and monoacylglycerols; where the associated metabolites from the CCA plot (and respective univariate analyses in parentheses) include lignoceric acid (~1.2-fold increases; p > 0.05), stearic acid (~1.5-fold increases; p > 0.05), 1-monopalmitin (~2.2-fold increases; p > 0.05), behenic acid (~10-fold increases; p > 0.05), and 2-monoolein (~32-fold increases; p > 0.05). However, in comparison, ChemRich analyses yield no support for changes in the intracellular fatty acids.
Further, trends from the CCA plot for the 50v1 strain (Figure 4B) support a correlation between the maximum relative biomass and metabolism of polyhydroxyalkanoates (PHAs), which are intracellular polymers (e.g., poly-4-hydroxybutanoate and/or other co-polymers) used for carbon and energy storage (Khanna and Srivastava, 2005). Among the metabolites trending alongside the maximum relative biomass are α-ketoglutarate, glutamate, 4-aminobutanoate, and 4-hydroxybutanoate. These combined metabolites represent a stepwise pathway toward the synthesis of PHAs (KEGG map 00250 and map00650), where oxidized carbons are acquired through the TCA cycle (α-ketoglutarate) and sequentially shuttled through amino acid metabolism (glutamate) and butanoate metabolism (4-aminobutanoate) to yield a monomeric unit (4-hydroxybutanoate) commonly found in PHAs.
Together, these observations are relevant since the accumulation of PHAs can result in increases in relative optical densities (Slaninova et al., 2018), which is consistent with our growth curves and ensuing regression analyses that yield increases in the maximum relative biomass in the presence of 0.1% v/v K30. Hence, for the 50v1 strain, the increases in maximum relative biomass may relate to the increased synthesis of PHAs under the stresses imposed by K30.
3.5.4. Comparisons to growth rates
When considering growth rates (exponential phase), the vector for the 2P08AA strain (Figure 4A) overlaps with the reduced terms associated with 0% v/v K30 (red squares). This correlation suggests that growth rates trends for the 2P08AA strain (across 0 to 0.1% v/v K30) are associated with limited changes in the intracellular metabolome. This assessment is consistent with the associated univariate tests for growth rates that support moderate decreases of ~18% (p < 0.05) and the associated univariate tests for the metabolomes that imply no change in relative abundances (FDR > 0.20).
From the CCA plot of the 2P08AA strain, metabolites trending alongside the growth rates include palmitoleic acid (fatty acid; 16:1Δ9), 5′-MTA, ribose, several amino acids, oleic acid, N-acetylglucosamine, and 5’-AMP. These trends suggest that the decreases in growth rates are related to changes in amino acid metabolism, cell wall metabolism, and nucleotide metabolism. The trends also suggest relations to changes in free unsaturated fatty acid content through associations with palmitoleic acid (~0.7-fold decrease; p > 0.05) and oleic acid (~0.3-fold decrease; p > 0.05). However, ChemRich analyses yield no support for changes in unsaturated fatty acid metabolism.
For the 50v1 strain (Figure 4B), in contrast, the vector describing growth rates overlaps more closely with the reduced terms associated with 0.1% v/v K30 (black triangles). This correlation suggests that growth rate trends for the 50v1 strain (across 0 to 0.1% v/v K30) are associated with changes to the intracellular metabolome. In comparison, univariate tests on the cultivation data indicate no change in the growth rates (p > 0.05), while univariate tests on the metabolomic data support changes in the abundances for multiple metabolites (p < 0.05, FDR ≤ 0.20). These combined trends are suggestive of substantial post-exponential phase changes in the metabolomes, as growth rates for the 50v1 strain during exponential phase are unaltered by – and potentially not associated with – the multiple metabolomic changes observed during late-log phase.
3.5.5. Comparisons to detergent concentrations
When considering the K30 formulation (Figures 4A,B), the vectors describing detergent concentrations (~0–0.025% w/w) overlap with the reduced terms associated with 0.1% v/v K30 for both the 2P08AA and 50v1 strains (black triangles). These correlations suggest the metabolomes from both strains adjust in response to detergents in 0.1% v/v K30. We note the trends for detergent concentrations likely overlap with the maximum relative biomass term, which plots across the same respective quadrant.
For the 2P08AA strain, metabolites trending alongside the detergent concentrations include myristic acid (fatty acid; 14:0), pelargonic acid (fatty acid; 9:0), 4-hydroxybenzoate, phenylacetic acid, and behenic acid (fatty acid; 22:0). The trends for 4-hydroxybenzoate and phenylacetic acid lend support toward the potential for biodegradation of dodecyl benzene sulfonate and polyethylene glycol mono-nonylphenyl ether, the detergents from the K30 formulation. As described in Hashim et al. (1992), phenylacetic acid is a potential product of microbial degradation of dodecyl benzene sulfonate through oxidation (or metabolism) of the dodecyl group (to yield 4-sulfonylphenylacetic acid), followed by desulphonization (to yield phenylacetic acid). Similarly, we posit that biodegradation of polyethylene glycol mono-nonylphenyl ether could yield 4-hydroxybenzoate through scission of the polyethylene glycol units (to yield 4-nonylphenol), as described in Section 1 followed by oxidation or metabolism of the nonyl side group (to yield 4-hydroxybenzoate). In comparison, however, univariate tests yield no support for changes in the abundances for phenylacetic acid (~2.4-fold increases; p = 0.39) and 4-hydroxybenzoate (~1.5-fold increases; p = 0.11).
The trends for the 2P08AA strain are also suggestive of correlations between detergent concentrations and changes in free fatty acid content. Relevant metabolites from the CCA plot (and the associated univariate analyses in parentheses) include behenic acid (~15-fold increases; p < 0.05), myristic acid (~1.2-fold increases; p > 0.05), and pelargonic acid (~1.4-fold increases; p > 0.05). We note that behenic acid levels in fluvial biofilms increase after exposure to desiccation (Serra-Compte et al., 2018). We also speculate the increases in long chain saturated fatty acid content may promote aggregation with the detergents in the intracellular and/or membrane space, which could potentially serve as a survival mechanism and/or acquisition strategy during biodegradation. Analyses by ChemRich, however, yield no support for changes in fatty acid or lipid-related metabolites.
For the 50v1 strain (Figure 4B), metabolites that trend alongside detergent concentrations include tartrate, 1-monoolein (1-oleoylglycerol; oleoyl = 18:1Δ9), and octadecanol (fatty alcohol). These trends are suggestive of increases in an intracellular dicarboxylic acid (tartrate) and free lipid-related molecules (1-monoolein and octadecanol), which may, respectively, relate to metal retention in the cell and the aggregation of detergents in the intracellular or membrane space.
3.6. Impacts of Kleenol 30 on the cultivations and extracellular metabolomes
3.6.1. Canonical correlation analyses
In Figure 4, we use canonical correlation analyses (CCA) to, respectively, compare the extracellular metabolomes from A. johnsonii 2P08AA (n = 3) and A. radioresistens 50v1 (n = 3) from late log-phase cultures to the quantitative descriptors of growth rates at exponential phase, plate counts at late-log phase, maximum relative biomass at stationary phase, and total concentration of EDTA and sodium metasilicate in the media, which was assumed to be 0.02% w/w in cultures containing 0.1% v/v K30.
The CCA plots in Figures 4C,D, which compare results from 0 and 0.1% v/v K30, include the total array of detected metabolites (text), listed quantitative descriptors (lines), and transformed data after dimension reduction (glyphs). For the 2P08AA and 50v1 strains, stronger correlations along CCA Axis 1 (67.18, 82.53%) are observed when compared to CCA Axis 2 (32.81, 11.10%), respectively. Dimension reduction for the metabolomes from the 2P08AA and 50v1 strains (Figures 4C,D) respectively show clear separations between samples prepared in the absence and presence of 0.1% v/v K30. These trends suggest that the extracellular metabolomic compositions of both strains are impacted by K30 during cultivation. Described in the following sub-sections are comparisons to the quantitative descriptors. For selected metabolites, additional explanations and/or fatty acid abbreviations are provided in parentheses.
3.6.2. Comparisons to survival and growth rates
When considering survival and cultivations, the vectors for plate counts and growth rates for the 2P08AA strain (Figure 4C) overlap with the reduced terms associated with 0% v/v K30 (red diamonds). These correlations suggest that survival and growth rate trends for the 2P08AA strain (across 0 to 0.1% v/v K30) are associated with minimal changes in the extracellular metabolome. These assessments are roughly consistent with the associated univariate analyses that indicate quantitative survivals, ~18% decreases in growth rate, and a lack of changes in the extracellular metabolome (after correction for multiple testing).
From the CCA plot, metabolites from the 2P08AA strain that trend alongside the plate counts and growth rates include 2-isopropylmalic acid, pyrophosphate, phosphate, arachidic acid (fatty acid; 20:0), lauric acid (fatty acid; 12:0), capric acid (fatty acid; 10:0), 4-hydroxyphenyllactic acid, sorbitol, and isopentadecanoic acid (13-methyl-14:0). These trends for the 2P08AA strain suggest that survival (which was not impacted) and growth rates (which minimally decreased) relate to a lack of change in the abundances of extracellular saturated fatty acids, phosphate, pyrophosphate, and α-hydroxy acids (2-isopropylmalic acid and 4-hydroxyphenyllactic acid). However, the potential increases in extracellular sorbitol (~5-fold increase, p < 0.05; Table 1) are suggestive of efflux of a compatible solute.
For the 50v1 strain (Figure 4D), in contrast, the vector describing plate counts overlaps with the reduced terms associated with 0% v/v K30 (pink diamonds), while the vector describing growth rates moderately overlaps with the reduced terms associated with 0.1% v/v K30 (black triangles). These opposing correlations suggest that survival trends (across 0 to 0.1% v/v K30) are associated with minimal changes in the extracellular metabolome, while growth rate trends (across 0 to 0.1% v/v K30) are associated with more notable changes in the extracellular metabolome.
Metabolites trending alongside the plate counts for the 50v1 strain include 2-isopropylmalic acid, 2-hydroxyglutarate, succinic acid, and α-ketoglutarate. These results suggest that survival for the 50v1 strain is related to changes in extracellular α-hydroxy acids (2-isopropylmalic acid, 2-hydroxyglutarate, and α-ketoglutarate) and a dicarboxylic acid (succinate); while univariate tests support ~9-fold decreases in extracellular α-ketoglutarate. Metabolites trending alongside the growth rates for the 50v1 strain include tagatose and N-acetylglucosamine, which suggests that changes in extracellular monosaccharides relate to maintaining stable growth rates.
3.6.3. Comparisons to maximum relative biomass and chelator concentrations
When considering changes in the maximum relative biomass (stationary phase) and chelator concentrations (EDTA and sodium metasilicate from the K30 formulation), similar trends are observed for the 2P08AA and 50v1 strains (Figures 4C,D). The vectors describing chelator concentrations for both strains overlap with the reduced terms associated with 0.1% v/v K30 (black triangles). These correlations suggest the extracellular metabolomes from both strains adjust in response to the change in chelator concentration across 0 to 0.1% v/v K30 (e.g., the metal sequestering agents of EDTA and sodium metasilicate).
From the CCA plot for the 2P08AA strain, metabolites trending alongside changes in the maximum relative biomass and/or chelator concentration include threonate, fructose, oxalate, galactonate, gluconate, gluconolactone, mannose, galactose, and glycerate. These trends suggest that chelator concentrations influence the abundances of extracellular monosaccharides and sugar acids. While these trends are consistent with univariate tests (p < 0.05), no support for changes in the abundances of the associated metabolites (Table 2) are obtained after correction for multiple testing (FDR ≤ 0.20).
For the 50v1 strain, metabolites trending alongside changes in the maximum relative biomass and/or chelator concentrations include N-acetylglucosamine, tartrate, 2-methylglycerate, gluconolactone, glycerate, mannose galactonate, and gluconate. These results suggest that the chelator concentrations influence the abundances of extracellular monosaccharides, modified monosaccharides, and sugar acids. These assessments generally are consistent with univariate analyses.
4. Discussion
4.1. Spacecraft-associated Acinetobacter
Molecular insights into the tolerance of spacecraft-associated Acinetobacter toward Kleenol 30, a NASA cleanroom floor detergent, were obtained through cultivation, kinetic, and metabolomic approaches. Comparisons were conducted using A. johnsonii 2P08AA and A. radioresistens 50v1, which were isolated at, respectively, unique locations within the spacecraft assembly facilities, which were subjected to differing cleaning regimes.
The strain A. johnsonii 2P08AA was isolated from the floor of the Mars Phoenix lander assembly facility, which was routinely cleansed with Kleenol 30. The strain A. radioresistens 50v1 was isolated from the surface of the pre-flight Mars Odyssey orbiter, which was likely cleaned with isopropyl alcohol (2-propanol, isopropanol, IPA) and never exposed to Kleenol 30.
Hence, comparisons across these differing sub-environments provide initial insights into the impacts of repeated cleansing, under the conditions of spacecraft assembly, toward biochemical and cultivation selection.
4.2. Survival of the spacecraft-associated Acinetobacter against Kleenol 30
Cultivations containing 0–2.0% v/v Kleenol 30 (K30) were conducted in a low-osmolarity minimal media (0.2x M9, 25 μM Fe) containing 150 mM ethanol, which served as the sole supplied carbon source. Cultivation media containing 2.0% v/v K30 were bactericidal for the tested Acinetobacter, which suggests a biocide-type of activity for K30.
In contrast, in 1.0% v/v K30, appreciable cultivation tolerances were observed for the 2P08AA and 50v1 strains. Cultures containing 1.0% v/v K30 exhibited (1) ~20–30% survivals, (2) ~66 and 12% decreases in growth rates, and (3) eliminated lag phase times for the 2P08AA and 50v1 strains, respectively. These cultivation tolerances of <2% v/v K30 are significant since current NASA planetary protection practices utilize ~0.8-1.6% v/v K30 for floor cleansing. Our results, therefore, suggest that the spacecraft-associated Acinetobacter harbor a tolerance toward the current floor cleansing practices for NASA cleanrooms.
In 0.1% v/v K30, differential impacts were observed across the cultivation parameters. For the 2P08AA strain, cultures exhibited (1) no measurable changes in cell density at late-log phase (plate counts), (2) modest decreases in the growth rate (~18%), and (3) negligible times in the lag phase. These results indicate that K30 impacts the 2P08AA strain by accelerating growth out of the lag phase, slightly inhibiting growth during the exponential phase, and by yielding no changes cell density after cultivation to late-log phase.
For the 50v1 strain, cultures exhibited (1) ~50% survivals, (2) no changes in growth rate, and (3) ~70% decreases in the lag phase time. This is indicative of K30 impacting the 50v1 strain by mildly accelerating growth out of the lag phase, yielding no change in the growth rates during the exponential phase, and inhibiting cell density by ~50% when cultivated to late-log phase.
The differential cultivation trends observed in 0.1% v/v K30 are indicative of strain and growth-phase dependent responses to K30. In turn, these trends suggest that K30 reduces nutrient bioavailability in culture (e.g., alkaline earth and transition metals) and/or elicits specific biochemical responses (e.g., metabolism of K30 and/or enzymatic inhibition by K30) – which together manifest as acceleration out of the lag phase for both strains, slight inhibition in rates of cell division during the exponential phase for the 2P08AA strain, and moderate reductions in total cell density at late-log phase for the 50v1 strain.
4.3. Metabolomic profiling of the cultivation tolerances against Kleenol 30
To obtain metabolic insights into the tolerances toward Kleenol 30 (K30), we used cultures grown in 0.1% v/v K30 to obtain reliable biomasses at late-log phase for both strains (see comparison of 0.1 and 1.0% w/w K30 in Figures 1D). For A. johnsonii 2P08AA, cultivations with 0.1% v/v K30 exhibited (1) no significant and discreet changes in the abundances of intracellular metabolites, (2) several significant extracellular changes when considering broader classifications such as sugar acids, monosaccharides, organic acids (di- and monocarboxylic acids), and fatty acids, and (3) potential support for biodegradation and metabolism of both detergents in the K30 formulation.
The combined results for the 2P08AA strain suggest that the native and intracellular metabolic state of the 2P08AA strain readily accommodates cultivation in the presence of 0.1% v/v K30. The trends for the 2P08AA strain also suggest that survival is associated with adjustments in the extracellular metabolome and a potential for biodegradation of the Kleenol 30 detergents.
For A. radioresistens 50v1, cultivations in 0.1% v/v K30 exhibited (1) differential changes in the abundances of intracellular amino acids, compatible solutes, nucleotide-related metabolites, dicarboxylic acids, and saturated fatty acids, (2) substantial and differential changes in the abundances of extracellular sugar acids, monosaccharides, and organic acids, and (3) no clear statistical support for biodegradation and metabolism of the detergents from the K30 formulation (when excluding the potential for degradation via benzene metabolism).
These results for the 50v1 strain indicate that the intracellular and extracellular metabolomes of the 50v1 strain significantly adjust to accommodate cultivation in the presence of 0.1% v/v K30. These trends also suggest that survival is associated with intracellular changes to amino acid and peptide metabolism, nucleotide metabolism, central metabolism, fatty acid metabolism, and pyrimidine metabolism.
For the 50v1 strain, the changes in amino acid metabolism include decreases in lysine, proline, glycine, ornithine, homoserine, valine and others listed in Table 1. This is relevant since decreases in lysine and valine are observed during biofilm formation in A. baumanni 1,656–2 (Yeom et al., 2013), while decreases in amino acid metabolism are observed after treatments of A. baumanni AB5075 with mixed antibiotics (polymyxin B and rifampicin) (Zhao et al., 2021). In E.coli UTI189, decreases in intracellular lysine and valine are associated with biofilm formation (Lu et al., 2019). In Bifodobacterium bifudum, decreases in intracellular lysine, arginine, and proline are associated with biofilm formation (Sadiq et al., 2020). These observations suggest that the Acinetobacter undergo stress-like responses toward K30 that inhibit total cell mass (e.g., survival), which could include changes in early biofilm development.
Trends across nucleotide abundances and the cultivation terms are similarly suggestive of early biofilm-type development (Sarkar and Chakraborty, 2008; Castillo-Juarez et al., 2017) and/or quorum sensing related changes such as intracellular PHA synthesis (Kessler and Witholt, 2001; Xu et al., 2011). These trends include increases in the maximum relative biomass terms for the 2P08AA and 50v1 strains, decreases in 5′-MTA and 5’-AMP for the 50v1 strain, and substantial increases in extracellular monosaccharides (e.g., mannose, galactose, and fructose for the 2P08AA strain).
In comparison, biofilms of E. coli UTI189 show decreases in 5′-MTA and 5’-CMP, as well increases in metal-binding siderophores, when Fe3+ concentrations in culture decrease from 10 to 1 μM (Guo et al., 2021). Lowered 5′-MTA abundances are associated with biofilms of Pseudomonas aeruginosa PAO1 (Tielen et al., 2013). In Pseudomonas fluorescens PF08, exogenously added 5’-AMP inhibits biofilm formation (Wang et al., 2021). In Streptococcus pyogenes, the import of mannose induces lysozyme resistance via a quorum sensing pathway (Chang et al., 2015); while exogenously added galactose differentially inhibits and activates biofilm development in differing species of Streptococcus (Ryu et al., 2016).
4.4. Biochemical strategies against Kleenol 30
To obtain biochemical insights into the mechanisms of tolerance, the metabolomic trends were interpreted in context of the cultivation conditions. In this study, cultivations were conducted in low-osmolarity media (0.2x M9) supplemented with 25 μM Fe2+ (Fe(NH4)2(SO4)2), which slowly oxidized to Fe3+ under the aerobic cultivations in the capped glass tubes (per absorbance spectroscopy). However, in control experiments, cultures supplemented with [Fe(EDTA)]1− grew slower, when compared to Fe(NH4)2(SO4)2, while cultivations in acid-washed glassware yielded no growth. Therefore, under the cultivation conditions, trace elements and species such as Zn2+, Co2+, Cu2+, Mn2+, BO33− (borate; a B source), and MoO42− (molybdate; a Mo source) were likely obtained from trace contaminants in the media and/or the glassware.
These details are relevant since K30 contains EDTA and sodium metasilicate, which are chelators that sequester transition and alkaline earth metals. Hence, in culture, chelation by EDTA would likely restrict the bioavailability of transition metals and inhibit growth – as observed in control studies. Limited bioavailabilities for Ca2+ and Mg2+ are also expected since sodium metasilicate (Na2SiO3) is a water softening agent that functions through ion exchange via chelation of Ca2+ or Mg2+ (and release of Na+).
Under these considerations, therefore, our review of the extracellular metabolomes reveals changes in several biochemicals possessing metal-binding properties. Appraisal of the literature shows that sugar acids such as galactonate, gluconate, and gluconolactone form appreciably stable complexes with transition metals such as Fe3+, Zn2+, Co2+, Cu2+, and Mn2+ (Sawyer, 1964; Carper et al., 1989; Frutos et al., 1998; Gyurcsik and Nagy, 2000). For the alkaline earth metals of Ca2+ and Mg2+, stable complexes are formed with sugar acids such as gluconolactone and glycerate (Taga et al., 1978; Meehan et al., 1979; Tajmir-Riahi, 1990), monosaccharides such as mannose (Angyal, 1973), and organic acids such as oxalate (Chutipongtanate et al., 2013; Izatulina et al., 2018) and tartrate (Hawthorne et al., 1982). In addition, borate and Ca2+ form stable complexes with gluconate, glycerate, and tartrate (van Duin et al., 1987); while molybdate yields stable complexes with mannose and fructose (Spence and Kiang, 1963; Sauvage et al., 1996).
Biochemically, therefore, these observations suggest that the 2P08AA strain targets trace Ca2+ and Mg2+ acquisition through potential increases in extracellular glycerate, gluconolactone, and oxalate, but not via gluconate and galactonate, which bind transition and alkaline earth metals but were not detected. In contrast, the 50v1 strain potentially targets both transition and alkaline earth metals given the increases (p < 0.05, FDR ≤ 0.20) in extracellular gluconate, galactonate, glycerate, gluconolactone, and tartrate. In fact, increases in glycerate export for the 50v1 strain are supported by the increases in extracellular glycerate and decreases in intracellular glycerate. Control of intracellular Ca2+ retention is also suggested by increases in extra- and intracellular tartrate. Multivariate tests support these total trends and suggest that chelator concentrations (from the K30 formulation) influence the abundances of extracellular metal-binding biochemicals.
Thus, our total results suggest that up-regulations (per se) in extracellular trace element acquisition strategies are common survival features to the spacecraft-associated Acinetobacter. Future studies using concentrated and clarified cultivation media may yield spectral support for the trace metal chelation. Additional common survival features may include regulation of the bulk osmolarity, as suggested by decreases in intracellular compatible solutes for the 50v1 strain (e.g., glycerate, 3-phosphoglycerate, and mannitol), statistical associations (CCA plot) between survival and changes in intracellular compatible solutes (e.g., trehalose) for the 2P08AA strain, and increase in an extracellular compatible solute for the 2P08AA strain (e.g., sorbitol).
When considering intracellular changes, differing responses to 0.1% v/v K30 are observed by the 2P08AA and 50v1 strains. While limited intracellular responses are observed for the 2P08AA strain, substantial intracellular changes are observed for the 50v1 strain. Multivariate tests for the 50v1 strain support a potential for the synthesis of light scattering PHAs during stationary phase in response to 0.1% v/v K30. Relatedly, univariate tests suggest decreases in the export of α-ketoglutarate and 4-hydroxybutanoate (e.g., decreases in extracellular abundances), while multivariate tests show correlations between intracellular PHA synthesis (e.g., α-ketoglutarate, 4-aminobutanoate, and 4-hydroxybutanoate) and increases in the apparent maximum relative biomass.
For the 2P08AA strain, multivariate analyses highlight a potential for partial metabolism of the detergents (from the Kleenol 30 formulation) – which would serve as an additional biochemical strategy toward survival. While these exact trends were not observed for the 50v1 strain, multivariate tests do support correlations between survival and changes in the degradation of aromatic compounds (benzene metabolism). In support, GC–MS studies in Mogul et al. (2018) provided evidence for scission of polyethylene glycol mono-nonylphenyl ether by the 50v1 strain. Future targeted cultivation studies are required to confirm the extent of biodegradation of the detergent by these strains.
When considering impacts of K30 on the lag phase, we observe a differential concentration dependence, where the time in lag phase is eliminated in the presence of 0.1 v/v K30 for the 2P08AA strain, and in 1.0% v/v K30 for the 50v1 strain. We venture that metabolism of the alkyl sides chains on the detergents (via ω-oxidation of the terminal carbon and subsequent cycles of β-oxidation) may serve as activators for acceleration out of the lag phase, and that degradation of the aromatic detergent cores (via benzene metabolism) may sustain the 50v1 strain during exponential phase, where growth rates show no inhibition in the presence of 0.1% v/v K30. If so, carbon and energy acquisition from the detergents could serve as an added biochemical strategy toward survival. Alternatively, the K30 detergents could potentially disperse any low abundance cellular aggregates during the lag phase, thereby yielding higher OD values and the apparent increases in biomass. Under this assumption, the observed concentration dependence across the strains would be suggestive of the cellular aggregates from both strains having differing stabilities.
5. Conclusion
In conclusion, our combined cultivation and metabolomic results lend supports toward the hypothesis that repetitive cleansing with Kleenol 30 presents as selective pressures toward members of the spacecraft microbiome – where the phenotypic outcomes likely include cultivation tolerances and measurable biochemical responses. Fittingly, our described results indicate that Kleenol 30 is more readily tolerated by the floor-associated A. johnsonii 2P08AA, when compared to the spacecraft surface-associated A. radioresistens 50v1. However, both strains display reasonable tolerances toward 1% v/v Kleenol 30, which is higher than the current formulation (~0.8-1.6% v/v) used in NASA planetary protection practices. These results indicate that the spacecraft-associated Acinetobacter may harbor tolerances toward the floor cleaning conditions.
Such tolerances are likely important considerations for lower bioburden planetary missions, such as orbiter missions to Europa (Planetary Protection Category III), life detection missions to Mars (Planetary Protection Category IVb) or Europa (Planetary Protection Category IV), and investigations of Mars Special Regions (Planetary Protection Category IVc). In this context, our results suggest that higher floor cleanliness levels may be achieved with higher concentrations of Kleenol 30 (e.g., 2–3% v/v) or rotations with differing concentrations of Kleenol 30 (e.g., 1 and 3% v/v). Alternatively, rotation of differing floor cleansers could assist in achieving lower floor bioburdens (e.g., Kleenol 30 and a quaternary ammonium compound), with the caveat that the detergents and/or biocides retain compatibility with the spacecraft assembly procedures and spacecraft materials.
Data availability statement
Metabolomic data (Study ID ST002380, DatatrackID: 3552) are available through the Metabolomics Workbench (https://www.metabolomicsworkbench.org/) at the NIH Common Fund’s National Metabolomics Data Repository (NMDR) and can be accessed directly via the project DOI: http://dx.doi.org/10.21228/M8XH73.
Author contributions
RM, DM, SL, and BR contributed to analysis of the data, assisted in drafting of the report, approved of submission of the manuscript. RM analyzed the metabolomic data, performed the statistical analyses, prepared the manuscript, and is the corresponding author. DM and BR conducted the survival assays and cultivation kinetics. SL prepared cultures for metabolomics analysis. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by a grant (NNH18ZDA001N) from the Planetary Protection Research component of the National Aeronautics and Space Administration (NASA) Research Opportunities in Space and Earth Sciences program.
Acknowledgments
We acknowledge logistical support from the College of Science, Chemistry and Biochemistry Department, and Department of Biological Sciences at Cal Poly Pomona.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Agalloco, J. P., and Carleton, F. J. (2007). Validation of Pharmaceutical Processes. Boca Raton, Florida: CRC Press.
Angyal, S. (1973). “Complex formation between sugars and metal ions” in Carbohydrate Chemistry–VI. ed. W. M. Doane, Northern Regional Research Laboratory, Peoria, Illinois (Amsterdam, Netherlands: Elsevier), 131–146.
Barupal, D. K., and Fiehn, O. (2017). Chemical similarity enrichment analysis (ChemRICH) as alternative to biochemical pathway mapping for metabolomic datasets. Sci. Rep. 7:14567. doi: 10.1038/s41598-017-15231-w
Barupal, D. K., Haldiya, P. K., Wohlgemuth, G., Kind, T., Kothari, S. L., Pinkerton, K. E., et al. (2012). MetaMapp: mapping and visualizing metabolomic data by integrating information from biochemical pathways and chemical and mass spectral similarity. BMC Bioinform. 13, 1–15. doi: 10.1186/1471-2105-13-99
Barupal, D. K., Zhang, Y., Shen, T., Fan, S., Roberts, B. S., Fitzgerald, P., et al. (2019). A comprehensive plasma metabolomics dataset for a cohort of mouse knockouts within the international mouse phenotyping consortium. Meta 9:101. doi: 10.3390/metabo9050101
Begot, C., Desnier, I., Daudin, J. D., Labadie, J. C., and Lebert, A. (1996). Recommendations for calculating growth parameters by optical density measurements. J. Microbiol. Methods 25, 225–232. doi: 10.1016/0167-7012(95)00090-9
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Carper, W., Coffin, D., and Addis, J. (1989). Studies of magnesium and manganese interactions with gluconate and 1, 5-gluconolactone. Spectrochim. Acta Part A 45, 391–392. doi: 10.1016/0584-8539(89)80106-0
Castillo-Juarez, I., Lopez-Jacome, L. E., Soberón-Chávez, G., Tomás, M., Lee, J., Castañeda-Tamez, P., et al. (2017). Exploiting quorum sensing inhibition for the control of Pseudomonas aeruginosa and Acinetobacter baumannii biofilms. Curr. Top. Med. Chem. 17, 1915–1927. doi: 10.2174/1568026617666170105144104
Chang, J. C., Jimenez, J. C., and Federle, M. J. (2015). Induction of a quorum sensing pathway by environmental signals enhances group a streptococcal resistance to lysozyme. Mol. Microbiol. 97, 1097–1113. doi: 10.1111/mmi.13088
Checinska, A., Probst, A. J., Vaishampayan, P., White, J. R., Kumar, D., Stepanov, V. G., et al. (2015). Microbiomes of the dust particles collected from the international Space Station and spacecraft assembly facilities. Microbiome 3:50. doi: 10.1186/s40168-015-0116-3
Checinska Sielaff, A., Urbaniak, C., Mohan, G. B. M., Stepanov, V. G., Tran, Q., Wood, J. M., et al. (2019). Characterization of the total and viable bacterial and fungal communities associated with the international Space Station surfaces. Microbiome 7:50. doi: 10.1186/s40168-019-0666-x
Chutipongtanate, S., Sutthimethakorn, S., Chiangjong, W., and Thongboonkerd, V. (2013). Bacteria can promote calcium oxalate crystal growth and aggregation. J. Biol. Inorg. Chem. 18, 299–308. doi: 10.1007/s00775-012-0974-0
Danko, D. C., Sierra, M. A., Benardini, J. N., Guan, L., Wood, J. M., Singh, N., et al. (2021). A comprehensive metagenomics framework to characterize organisms relevant for planetary protection. Microbiome 9, 1–15. doi: 10.1186/s40168-021-01020-1
Daschner, F., Fahlberg, W. J., and Mallison, G. F. (1982). Floor disinfection in the United States. Infect. Control Hosp. Epidemiol. 3, 281–284. doi: 10.1017/S0195941700056289
Derecho, I., Mccoy, K. B., Vaishampayan, P., Venkateswaran, K., and Mogul, R. (2014). Characterization of hydrogen peroxide–resistant Acinetobacter species isolated during the Mars Phoenix spacecraft assembly. Astrobiology 14, 837–847. doi: 10.1089/ast.2014.1193
Eissa, M. E., Kamel, M., and Beshir, M. (2014). Study of antimicrobial power of amphoteric disinfectants of Tego series used in pharmaceutical industry. Egypt. J. Biol. Sci. G 6, 17–27. doi: 10.21608/eajbsg.2014.16626
Ellingson, K. D., Pogreba-Brown, K., Gerba, C. P., and Elliott, S. P. (2020). Impact of a novel antimicrobial surface coating on health care–associated infections and environmental bioburden at 2 urban hospitals. Clin. Infect. Dis. 71, 1807–1813. doi: 10.1093/cid/ciz1077
Empadinhas, N., and Da Costa, M. S. (2008). Osmoadaptation mechanisms in prokaryotes: distribution of compatible solutes. Int. Microbiol. 11, 151–161.
Farrow, J. M. 3rd, Wells, G., and Pesci, E. C. (2018). Desiccation tolerance in Acinetobacter baumannii is mediated by the two-component response regulator BfmR. PLoS One 13:e0205638. doi: 10.1371/journal.pone.0205638
Fernández-Cuenca, F., Tomás, M., Caballero-Moyano, F.-J., Bou, G., Martínez-Martínez, L., Vila, J., et al. (2015). Reduced susceptibility to biocides in Acinetobacter baumannii: association with resistance to antimicrobials, epidemiological behaviour, biological cost and effect on the expression of genes encoding porins and efflux pumps. J. Antimicrob. Chemother. 70, 3222–3229. doi: 10.1093/jac/dkv262
Fiehn, O., Garvey, W. T., Newman, J. W., Lok, K. H., Hoppel, C. L., and Adams, S. H. (2010). Plasma metabolomic profiles reflective of glucose homeostasis in non-diabetic and type 2 diabetic obese African-American women. PLoS One 5:e15234. doi: 10.1371/journal.pone.0015234
Fiehn, O., Wohlgemuth, G., and Scholz, M. (2005). “Setup and annotation of metabolomic experiments by integrating biological and mass spectrometric metadata” in International Workshop on Data Integration in the Life Sciences. ed. W. M. Doane, Northern Regional Research Laboratory, Peoria, Illinois (Berlin: Springer), 224–239.
Franceus, J., Pinel, D., and Desmet, T. (2017). Glucosylglycerate phosphorylase, an enzyme with novel specificity involved in compatible solute metabolism. Appl. Environ. Microbiol. 83:e01434. doi: 10.1128/AEM.01434-17
Frick, A., Mogul, R., Stabekis, P., Conley, C. A., and Ehrenfreund, P. (2014). Overview of current capabilities and research and technology developments for planetary protection. Adv. Space Res. 54, 221–240. doi: 10.1016/j.asr.2014.02.016
Frutos, A. A., Signorella, S., Sala, L. F., Escandar, G. M., Peregrin, J. M. S., and Moreno, V. (1998). Interaction of divalent metal ions with D-aldonic acids in the solid phase-II. Structural information fn1. Polyhedron 17, 3371–3378. doi: 10.1016/S0277-5387(98)00119-3
Guo, R., Luo, X., Liu, J., and Lu, H. (2021). Mass spectrometry based targeted metabolomics precisely characterized new functional metabolites that regulate biofilm formation in Escherichia coli. Anal. Chim. Acta 1145, 26–36. doi: 10.1016/j.aca.2020.12.021
Gyurcsik, B., and Nagy, L. (2000). Carbohydrates as ligands: coordination equilibria and structure of the metal complexes. Coord. Chem. Rev. 203, 81–149. doi: 10.1016/S0010-8545(99)00183-6
Haffenden, L. J., and Yaylayan, V. A. (2005). Mechanism of formation of redox-active hydroxylated benzenes and pyrazine in 13C-labeled glycine/D-glucose model systems. J. Agric. Food Chem 53, 9742–9746. doi: 10.1021/jf051610r
Halket, J. M., Waterman, D., Przyborowska, A. M., Patel, R. K., Fraser, P. D., and Bramley, P. M. (2005). Chemical derivatization and mass spectral libraries in metabolic profiling by GC/MS and LC/MS/MS. J. Exp. Biol. 56, 219–243. doi: 10.1093/jxb/eri069
Hammer, Ø., Harper, D. A., and Ryan, P. D. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4:9.
Han, Z., Pappas, E., Simmons, A., Fox, J., Donskey, C. J., and Deshpande, A. (2021). Environmental cleaning and disinfection of hospital rooms: a nationwide survey. Am. J. Infect. Cont. 49, 34–39. doi: 10.1016/j.ajic.2020.08.008
Hashim, M., Kulandai, J., and Hassan, R. (1992). Biodegradability of branched alkylbenzene sulphonates. J. Chem. Technol. Biotechnol. 54, 207–214. doi: 10.1002/jctb.280540302
Hawthorne, F., Borys, I., and Ferguson, R. (1982). Structure of calcium tartrate tetrahydrate. Acta Crystallogr. B 38, 2461–2463. doi: 10.1107/S0567740882009042
Hendrickson, R., Urbaniak, C., Minich, J. J., Aronson, H. S., Martino, C., Stepanauskas, R., et al. (2021). Clean room microbiome complexity impacts planetary protection bioburden. Microbiome 9:238. doi: 10.1186/s40168-021-01159-x
Izatulina, A. R., Gurzhiy, V. V., Krzhizhanovskaya, M. G., Kuz’mina, M. A., Leoni, M., and Frank-Kamenetskaya, O. V. (2018). Hydrated calcium oxalates: crystal structures, thermal stability, and phase evolution. Cryst. Growth Des. 18, 5465–5478. doi: 10.1021/acs.cgd.8b00826
Kessler, B., and Witholt, B. (2001). Factors involved in the regulatory network of polyhydroxyalkanoate metabolism. J. Biotechnol. 86, 97–104. doi: 10.1016/S0168-1656(00)00404-1
Khanna, S., and Srivastava, A. K. (2005). Recent advances in microbial polyhydroxyalkanoates. Process Biochem. 40, 607–619. doi: 10.1016/j.procbio.2004.01.053
La Duc, M. T., Nicholson, W., Kern, R., and Venkateswaran, K. (2003). Microbial characterization of the Mars odyssey spacecraft and its encapsulation facility. Environ. Microbiol. 5, 977–985. doi: 10.1046/j.1462-2920.2003.00496.x
La Duc, M. T., Vaishampayan, P., Nilsson, H. R., Torok, T., and Venkateswaran, K. (2012). Pyrosequencing-derived bacterial, archaeal, and fungal diversity of spacecraft hardware destined for Mars. Appl. Environ. Microbiol. 78, 5912–5922. doi: 10.1128/AEM.01435-12
Lu, H., Que, Y., Wu, X., Guan, T., and Guo, H. (2019). Metabolomics deciphered metabolic reprogramming required for biofilm formation. Sci. Rep. 9:13160. doi: 10.1038/s41598-019-49603-1
Mahnert, A., Vaishampayan, P., Probst, A. J., Auerbach, A., Moissl-Eichinger, C., Venkateswaran, K., et al. (2015). Cleanroom maintenance significantly reduces abundance but not diversity of indoor microbiomes. PLoS One 10:e0134848. doi: 10.1371/journal.pone.0134848
McCoy, K., Derecho, I., Wong, T., Tran, H., Huynh, T., La Duc, M., et al. (2012). Insights into the extremotolerance of Acinetobacter radioresistens 50v1, a gram-negative bacterium isolated from the Mars odyssey spacecraft. Astrobiology 12, 854–862. doi: 10.1089/ast.2012.0835
Meehan, E., Einspahr, H., and Bugg, C. E. (1979). Calcium binding to hydroxycarboxylic acids: crystal structure of calcium di-dl-glycerate dihydrate. Acta Crystallogr. B 35, 828–832. doi: 10.1107/S0567740879004994
Meyer, B., and Cookson, B. (2010). Does microbial resistance or adaptation to biocides create a hazard in infection prevention and control? J. Hosp. Infect. 76, 200–205. doi: 10.1016/j.jhin.2010.05.020
Mogul, R., Barding, G. A. Jr., Lalla, S., Lee, S., Madrid, S., Baki, R., et al. (2018). Metabolism and biodegradation of spacecraft cleaning reagents by strains of spacecraft-associated Acinetobacter. Astrobiology 18, 1517–1527. doi: 10.1089/ast.2017.1814
Moissl-Eichinger, C., Auerbach, A. K., Probst, A. J., Mahnert, A., Tom, L., Piceno, Y., et al. (2015). Quo vadis? Microbial profiling revealed strong effects of cleanroom maintenance and routes of contamination in indoor environments. Sci. Rep. 5:9156. doi: 10.1038/srep09156
Mora, M., Mahnert, A., Koskinen, K., Pausan, M. R., Oberauner-Wappis, L., Krause, R., et al. (2016). Microorganisms in confined habitats: microbial monitoring and control of intensive care units, operating rooms, cleanrooms and the international Space Station. Front. Microbiol. 7:1573. doi: 10.3389/fmicb.2016.01573
Murtough, S., Hiom, S., Palmer, M., and Russell, A. (2001). Biocide rotation in the healthcare setting: is there a case for policy implementation? J. Hosp. Infect. 48, 1–6. doi: 10.1053/jhin.2001.0950
NASA. (2011). Planetary Protection Provisions for Robotic Extraterrestrial Missions. Washington, DC: National Aeronautics and Space Administration.
Peleg, A. Y., Seifert, H., and Paterson, D. L. (2008). Acinetobacter baumannii: emergence of a successful pathogen. Clin. Microbiol. Rev. 21, 538–582. doi: 10.1128/CMR.00058-07
Perry-Dow, K. A., De Man, T. J., Halpin, A. L., Shams, A. M., Rose, L. J., and Noble-Wang, J. A. (2022). The effect of disinfectants on the microbial community on environmental healthcare surfaces using next generation sequencing. Am. J. Infect. Cont. 50, 54–60. doi: 10.1016/j.ajic.2021.08.027
Pirog, T., Lutsai, D., and Muchnyk, F. (2021). Biotechnological potential of the Acinetobacter genus bacteria. Mikrobiol. Zh. 83, 92–109. doi: 10.15407/microbiolj83.03.092
Pompilio, A., Scribano, D., Sarshar, M., Di Bonaventura, G., Palamara, A. T., and Ambrosi, C. (2021). Gram-negative bacteria holding together in a biofilm: the Acinetobacter baumannii way. Microorganisms 9:1353. doi: 10.3390/microorganisms9071353
Rajamohan, G., Srinivasan, V., and Gebreyes, W. A. (2009). Biocide-tolerant multidrug-resistant Acinetobacter baumannii clinical strains are associated with higher biofilm formation. J. Hosp. Infect. 73, 287–289. doi: 10.1016/j.jhin.2009.07.015
Rummel, J. D. (1992). Planetary protection policy (U.S.A.). Adv. Space Res. 12, 129–131. doi: 10.1016/0273-1177(92)90166-U
Rutala, W. A., and Weber, D. J. (2019). Best practices for disinfection of noncritical environmental surfaces and equipment in health care facilities: a bundle approach. Am. J. Infect. Cont. 47, A96–A105. doi: 10.1016/j.ajic.2019.01.014
Rutala, W. A., and Weber, D. J., Healthcare Infection Control Practices Advisory Committee. (2008). Guideline for Disinfection and Sterilization in Healthcare Facilities. Atlanta, Georgia: Centers for Disease Control and Prevention.
Ryu, E.-J., Sim, J., Sim, J., Lee, J., and Choi, B.-K. (2016). D-galactose as an autoinducer 2 inhibitor to control the biofilm formation of periodontopathogens. J. Microbiol. 54, 632–637. doi: 10.1007/s12275-016-6345-8
Sadiq, F. A., Yan, B., Zhao, J., Zhang, H., and Chen, W. (2020). Untargeted metabolomics reveals metabolic state of Bifidobacterium bifidum in the biofilm and planktonic states. LWT 118:108772. doi: 10.1016/j.lwt.2019.108772
Sandle, T. (2012). The CDC Handbook-a Guide to Cleaning and Disinfecting Clean Rooms. England: Grosvenor House Publishing.
Sandle, T. (2015). Pharmaceutical Microbiology: Essentials for Quality Assurance and Quality Control. Sawston, United Kingdom: Woodhead Publishing.
Sarkar, S., and Chakraborty, R. (2008). Quorum sensing in metal tolerance of Acinetobacter junii BB1A is associated with biofilm production. FEMS Microbiol. Lett. 282, 160–165. doi: 10.1111/j.1574-6968.2008.01080.x
Sauvage, J.-P., Verchère, J.-F., and Chapelle, S. (1996). A multinuclear NMR spectroscopy study of the tungstate and molybdate complexes of d-fructose and l-sorbose. Carbohydr. Res. 286, 67–76. doi: 10.1016/0008-6215(96)00048-1
Serra-Compte, A., Corcoll, N., Huerta, B., Rodríguez-Mozaz, S., Sabater, S., Barceló, D., et al. (2018). Fluvial biofilms exposed to desiccation and pharmaceutical pollution: new insights using metabolomics. Sci. Total Environ. 618, 1382–1388. doi: 10.1016/j.scitotenv.2017.09.258
Shams, A. M., Rose, L. J., Edwards, J. R., Cali, S., Harris, A. D., Jacob, J. T., et al. (2016). Assessment of the overall and multidrug-resistant organism bioburden on environmental surfaces in healthcare facilities. Infect. Control Hosp. Epidemiol. 37, 1426–1432. doi: 10.1017/ice.2016.198
Slaninova, E., Sedlacek, P., Mravec, F., Mullerova, L., Samek, O., Koller, M., et al. (2018). Light scattering on PHA granules protects bacterial cells against the harmful effects of UV radiation. App. Microbiol. Biotechnol. 102, 1923–1931. doi: 10.1007/s00253-018-8760-8
Spence, J., and Kiang, S.-C. (1963). Complexes of sugars with molybdate. J. Org. Chem. 28, 244–245. doi: 10.1021/jo01036a522
Spry, J. A., Siegel, B., Race, M., Kminek, G., Pratt, L., and Coustenis, A. (2020). Planetary Protection Requirements for Crewed Missions-an Update to the COSPAR Workshop Series: 2020. International Conference on Environmental Systems.
Srinivasan, V. B., Venkataramaiah, M., Mondal, A., and Rajamohan, G. (2015). Functional characterization of AbeD, an RND-type membrane transporter in antimicrobial resistance in Acinetobacter baumannii. PLoS One 10:e0141314. doi: 10.1371/journal.pone.0141314
Sud, M., Fahy, E., Cotter, D., Azam, K., Vadivelu, I., Burant, C., et al. (2016). Metabolomics workbench: an international repository for metabolomics data and metadata, metabolite standards, protocols, tutorials and training, and analysis tools. Nucleic Acids Res. 44, D463–D470. doi: 10.1093/nar/gkv1042
Suleyman, G., Alangaden, G., and Bardossy, A. C. (2018). The role of environmental contamination in the transmission of nosocomial pathogens and healthcare-associated infections. Cur. Infect. Dis. Rep. 20, 1–11. doi: 10.1007/s11908-018-0620-2
Taga, T., Ohashi, M., and Osaki, K. (1978). Interactions of calcium ions with carbohydrates: crystal structure of calcium DL-Glycerate Dihydrate. Bull. Chem. Soc. Japan 51, 1697–1700.
Tajmir-Riahi, H. (1990). Carbohydrate complexes with alkaline earth metal ions. Interaction of D-glucono-1,5-lactone with the mg(II), ca(II), Sr(II), and Ba(II) cations in the crystalline solid and aqueous solution. J. Inorg. Biochem. 39, 33–41. doi: 10.1016/0162-0134(90)80013-N
Tam, H.-K., König, P., Himpich, S., Ngu, N. D., Abele, R., Müller, V., et al. (2022). Unidirectional mannitol synthesis of Acinetobacter baumannii MtlD is facilitated by the helix–loop–helix-mediated dimer formation. Proc. Natl. Acad. Sci. U. S. A. 119:e2107994119. doi: 10.1073/pnas.2107994119
Tembo, G. M., Chaggar, G. K., Li, X., Teska, P. J., and Oliver, H. F. (2022). Evaluation of automated floor cleaning, disinfection, and application methods against Staphylococcus aureus. Am. J. Infect. Cont. 19:S0196-6553(22)00535-1. doi: 10.1016/j.ajic.2022.07.011
Tielen, P., Rosin, N., Meyer, A.-K., Dohnt, K., Haddad, I., Jänsch, L., et al. (2013). Regulatory and metabolic networks for the adaptation of Pseudomonas aeruginosa biofilms to urinary tract-like conditions. PLoS One 8:e71845. doi: 10.1371/journal.pone.0071845
Trigg, S. A., Mcelhany, P., Maher, M., Perez, D., Busch, D. S., and Nichols, K. M. (2019). Uncovering mechanisms of global ocean change effects on the Dungeness crab (Cancer magister) through metabolomics analysis. Sci. Rep. 9, 1–12. doi: 10.1038/s41598-019-46947-6
Van Duin, M., Peters, J. A., Kieboom, A. P., and Van Bekkum, H. (1987). Synergic coordination of calcium in borate-polyhydroxycarboxylate systems. Carbohydr. Res. 162, 65–78. doi: 10.1016/0008-6215(87)80201-X
Wang, Y., Wang, Y., Chen, J., Koseki, S., Yang, Q., Yu, H., et al. (2021). Screening and preservation application of quorum sensing inhibitors of Pseudomonas fluorescens and Shewanella baltica in seafood products. LWT 149:111749. doi: 10.1016/j.lwt.2021.111749
White, G. F. (1993). Bacterial biodegradation of ethoxylated surfactants. Pest. Sci. 37, 159–166. doi: 10.1002/ps.2780370209
Wong, W. C., Zalesak, S., Yoets, A., Capriotti, J., Smith, M. J., and Castro, V. A. (2013). Effectiveness of HEPA filter vacuum in removing transient microbial contaminants on cargo bags destined for the international Space Station. J. Occup. Environ. Hyg. 10, D71–D75. doi: 10.1080/15459624.2013.784179
Xu, C., Li, M., Huang, Y., Zhang, Z., Bian, Z., and Song, S. (2011). Involvement of quorum-sensing in biosynthesis of polyhydroxyalkanoates in Pseudomonas aeruginosa. Acta Microbiolog. Sin. 51, 769–775.
Yeom, J., Shin, J.-H., Yang, J.-Y., Kim, J., and Hwang, G.-S. (2013). 1H NMR-based metabolite profiling of planktonic and biofilm cells in Acinetobacter baumannii 1656-2. PLoS One 8:e57730. doi: 10.1371/journal.pone.0057730
Zhao, J., Han, M.-L., Zhu, Y., Lin, Y.-W., Wang, Y.-W., Lu, J., et al. (2021). Comparative metabolomics reveals key pathways associated with the synergistic activity of polymyxin B and rifampicin combination against multidrug-resistant Acinetobacter baumannii. Biochem. Pharmacol. 184:114400. doi: 10.1016/j.bcp.2020.114400
Keywords: Acinetobacter, detergent, survival, metabolomics, planetary protection, spacecraft, cleanrooms
Citation: Mogul R, Miller DR, Ramos B and Lalla SJ (2023) Metabolomic and cultivation insights into the tolerance of the spacecraft-associated Acinetobacter toward Kleenol 30, a cleanroom floor detergent. Front. Microbiol. 14:1090740. doi: 10.3389/fmicb.2023.1090740
Edited by:
Camilla Urbaniak, NASA Jet Propulsion Laboratory, United StatesReviewed by:
Markus Kalkum, Beckman Research Institute, United StatesJason M. Wood, University of Illinois at Chicago, United States
Copyright © 2023 Mogul, Miller, Ramos and Lalla. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rakesh Mogul, rmogul@cpp.edu