Physicochemical properties of ammonia and pH-matched solutions
Ammonia solutions were prepared from liquid ammonia to molar concentrations of 0.01, 0.025, 0.05, 0.1, 0.25, 0.5 and 1 M, equating to pH values between 8 and 11 (Fig. 1a). Positive control (PC) was unamended yeast media. The percentage abundance of NH3 to NH4+ at these concentrations and pH values is given in Fig. 1b. From concentrations of 0.05 M ammonia and higher, NH3 accounts for over 50% of the total NH3/NH4+ in solution. Solutions of increasing ammonia concentration showed decreased oxygen concentrations; however, this decrease was non-significant from pH-matched counterparts (Fig. 1c). Alterations to water activity with increasing ammonia and pH were also statistically non-significant (Fig. 1d).
Ammonia sets a distinct molarity threshold for viability and growth
Bacteria were grown in ammonia solutions utilising two system types: a ‘closed-air’ system – to determine growth limits in a perpetually ammonia exposed environments, and an ‘open-air’ system – to determine growth limits in an environment permitting ammonia dispersal. Fig. 2a shows the final colony forming units of H. meridiana after incubation in the closed-air system for 72 h. Colonies were evident following incubation in 0.01 M, 0.025 M and 0.05 M ammonia, but colony number decreased as NH3 increased relative to NH4+ (Fig. 1b). A lower number of viable colonies compared to positive control was observed in 0.01 M ammonia (p < 0.01). No significant difference in colony number was observed when H. meridiana was incubated in 0.025 M compared to positive control (p = 0.056) or 0.01 M ammonia (p = 0.323). This concentration of ammonia has a pH within the optimum range for H. meridiana (pH 9) (Fig. 1a). Lower viable cell numbers were observed when incubated in 0.05 M compared to positive control (p < 0.0001) and 0.025 M ammonia (p < 0.01). No viable colonies were present following incubation in 0.1 M ammonia or ammonia solutions of higher concentration (Fig. 2a).
Fig. 2b. shows the growth dynamics of H. meridiana Slthf1 over 48 h in increasing concentrations of ammonia in an open-air system. Lag phase duration, doubling time (Td) and final OD600 are presented in Fig. 2c-e, respectively. At concentrations exceeding 0.01 M ammonia, lag phase was greater with higher ammonia concentrations; incubation in 0.1 M and 0.25 M ammonia prolongs lag phase time by 6-fold and 25-fold compared to the positive control, respectively (Fig. 2c). The lag phase time in 0.1 M ammonia was higher than that in 0.05 M ammonia (p < 0.05). Likewise, the lag phase time in 0.25 M ammonia was higher than that in 0.1 M ammonia (p < 0.0001) (Fig. 2c). Td was also higher with increasing ammonia concentration, with higher Td compared to the positive control observed in 0.05 M (p < 0.05), 0.1 M (p ≤ 0.0001) and 0.25 M (p < 0.0001) solutions (Fig. 2d). The Td was not significantly altered between positive control and 0.01 M (p = 0.965) or 0.1 M and 0.25 M solutions (p = 0.9972). Final OD600 after 48 h growth showed no statistically significant difference from positive control in 0.01 M (p = 0.932), 0.025 M (p = 0.330), 0.1 M (p = 0.248) and 0.25 M solutions (p = 0.135), but higher OD600 values were observed in 0.05 M compared to the positive control (p < 0.01) (Fig. 2e). Cell density was lower in 0.5 and 1 M solutions compared to the positive control (0.5 M, p < 0.01; 1 M, p < 0.001) (Fig. 2e), where no distinguishable growth occurred within 48 hours (Fig. 2b).
The prolonged lag phase in cells inoculated into 0.1 M and 0.25 M ammonia indicate bactericidal reduction to the cell population. Alternatively, ammonia may act bacteriostatically, the effects of which cease once ammonia has dispersed. To assess this, temporal cell viability was compared with ammonia concentration over time in 0.1 M (Fig. 2f) and 0.25 M ammonia solutions (Fig. 2g). After 4 h exposure, viable cells in 0.1 M ammonia show a 1000-fold decrease from 0 h (t=28.28, df=4, p < 0.0001) (Fig. 2f). In 0.25 M ammonia, viable cells at 4 h are reduced 10-fold from 0 h (t=6.499, df=4, p < 0.01) (Fig. 2g). Bactericidal reduction to cell populations within 0.1 M ammonia solutions cease between 4 h and 8 h, where a small but non-significant increase in cell viability is observed t=2.421, df=4, p = 0.0727) (Fig. 2f). Likewise, bactericidal effects are absent between 4 h and 8 h in 0.25 M ammonia solutions where cell numbers stabilise (t=0.8934, df=4, p = 0.442) and increase from 8 h to 16 h (t=11.76, df=4, p < 0.001) (Fig. 2g). The increase in cell viability within 0.1 and 0.25 M ammonia solutions at 8 h and 16 h, respectively, aligns with diminished ammonia levels to ≤ 0.05 M at 4 h and 13 hr, respectively (Fig. 2f, 2g). Thus, lag phase extension in the open-air system reflects two events: (1) an immediate bactericidal effect and (2) growth initiation after ammonia levels drop to sub-bactericidal levels (≤ 0.05 M) in surviving populations.
Ammonia toxicity is independent from pH toxicity
Ammonia is a weak base that raises solution pH with increasing concentration. To delimitate ammonia toxicity from pH toxicity, growth experiments were repeated in NaOH solutions pH-matched to ammonia solutions. Fig. 3a presents the final CFU/mL of H. meridiana in solutions following 72 h closed-air system incubation. No significant differences were observed between cells grown in positive control and pH-matched solutions up to pH 10.78 (equivalent to 1 M ammonia). Fig. 3b presents the growth curve of H. meridiana grown under pH-matched NaOH solutions in an open-air system, with extrapolated parameters of lag phase (Fig. 3c), Td (Fig. 3d) and final OD600 (Fig. 3e) presented in log2 fold changes (log2FC) from growth in ammonia solutions of the same pH. Significantly lower lag phases were observed in cells grown in pH 8.96, 9.38, 9.73 and 10.18 compared to those grown at 0.025 M (t=4.399, df=3, p < 0.05), 0.05 M (t=11.26, df=3, p < 0.01), 0.1 M (t=7.357, df=3, p < 0.001) and 0.25 M ammonia (t=31.72, df=3, p < 0.0001) (Fig. 3c). The Td compared between cells grown in ammonia and pH-matched counterparts was non-significantly different except for those grown in pH 9.38 (t=3.542, df=3, p < 0.05), which showed a lower Td (Fig. 3d). Final OD600 at 48 h was lower in pH 8.96 (t=4.882, df=3, p < 0.05), pH 9.38 (t=6.810, df=3, p < 0.01), pH 9.73 (t=4.123, df=3, p < 0.05) and pH 10.18 (t=7.206, df=3, p < 0.01) NaOH-solutions compared to growth in ammonia counterparts, while higher OD600 was observed in pH 10.49 (t=16.09, df=3, p < 0.001) and pH 10.78 solutions (t=16.74, df=3, p < 0.001) compared to ammonia counterparts (Fig. 3e). These findings confirm that lag phase extension and absent growth in 0.5 and 1 M ammonia cannot be attributed to pH increases.
To broaden the pH comparison beyond NaOH, lag phase analysis was used to explore whether the prolonged lag phase in 0.25 M ammonia (Fig. 2c) could be replicated in KOH, Na2SiO3, K2CO3, and Na2CO3 at pH 10.18 (Fig. 3f). No differences in lag phase duration were observed between H. meridiana grown in Na2SiO3 (p = 0.779) or Na2CO3 (p = 0.996) compared to NaOH. Longer lag phases were seen in KOH (p < 0.05) and K2CO3 (p < 0.0001) compared to NaOH, possibly reflecting reduced adaptation to potassium in H. meridiana. All high-pH solutions showed significantly shorter lag phases than 0.25 M ammonia (p < 0.001), confirming the pH-independent toxicity of ammonia.
Ammonia toxicity exerts distinct changes to bacterial morphology
Specific bactericidal effects of ammonia on cells can be observed in Fig. 4. Under control growth conditions (unaltered yeast media) (Fig. 4a-b), H. meridiana exhibited ribosome-rich cytoplasm, visible nucleoids and clear division of outermembrane, periplasm and inner membrane. Cells exhibited irregular, undulating outermembrane, enlarged periplasmic spaces and PHA-like granules. Upon 2 h treatment of 1 M ammonia (Fig. 4c-d), cells showed intracellular aggregation and loss of ribosomes with some showing cytoplasmic loss and lysis. The periplasmic space volume decreased, and detachment of the inner membrane from the outermembrane was observed. Cells showed expansion of electrolucent cavities, possibly surrounding nucleoids. Condensed material within the electrolucent cavities exhibited splayed morphologies, possibly indicating disruption to DNA supercoiling. Cells exposed to a NaOH solutions pH-matched to 1 M ammonia (10.78 pH) also showed some cell lysis events (Fig. 4e-f). However, there were few morphological differences from cells grown in control media. Differences include uniform outermembrane, and fewer PHA-like granules.
Metabolic pathways altered in response to ammonia and high pH exposure
The survival of cells in up to 0.25 M ammonia in an open-air system suggests adaptations enabling tolerance to ammonia. Thus, untargeted metabolomics was performed on H. meridiana in unamended yeast media (hereafter denoted ‘control’) and yeast media with 0.25 M ammonia at pH 10.18 (hereafter denoted ‘0.25 M ammonia’) to determine variations in metabolism that may account for the differences observed. To delimitate high pH adaptations from ammonia adaptations, metabolomics was also performed following exposure to yeast media adjusted to pH 10.18 with NaOH (hereafter denoted ‘NaOH pH 10.18’). Growth kinetics and sampling points in these conditions are indicated in Fig. 5a. Principal component analysis (PCA) (Fig. 5b) separated 0.25 M ammonia and NaOH pH 10.18 conditions from the control, with overlap between 0.25 M ammonia and NaOH pH 10.18 conditions indicating metabolic similarity. Univariate volcano analysis identified 23 features significantly altered in 0.25 M ammonia/control (Fig. 5c), 10 features significantly altered in 0.25 M ammonia/NaOH pH 10.18 (Fig. 5d), and 28 features significantly altered in NaOH pH 10.18/control (Fig. 5e). The volcano analysis dataset for each comparison group is shown in Supplementary Table 2-4.
The most significantly altered metabolites between these conditions were identified by ANOVA (Supplementary Table 5), depicted as a heatmap (Fig. 5f). Overall, high similarity was found between 0.25 M ammonia and NaOH pH 10.18 exposed samples; both conditions generally exhibited higher levels of unsaturated phospholipid and lower levels of linoleic acid and derivatives. Intermediates that feed glycerophospholipid biosynthesis, CMP-sialic acid (F = 2021.9, p < 0.0001) and glycerol-3-phosphate (F = 2022.2, p < 0.0001), were more abundant compared to the control. Amino acid pathway intermediates indole-3-ethanol (F = 2327.8, p < 0.0001), N, N-dimethylglycine (F = 30.626, p < 0.001) and N-acetylglutamate (F = 25.653, p < 0.01) were also altered compared to control. There were reduced levels of N-acetylserotonin (F = 18876, p < 0.0001) and N-acetyl-L-aspartic acid (F = 24.091, p < 0.01) that could suggest reactions with acetyl donors were less favourable at high pH in both conditions.
Ammonia exposure elicits a unique metabolomic response
Despite general similarities between cells exposed to 0.25 M ammonia and NaOH pH 10.18, ANOVA indicated five metabolites significantly altered only in 0.25 M ammonia exposed H. meridiana. Box plots of these metabolites are presented in Fig. 6a-e. Samples exposed to 0.25 M ammonia showed significantly higher levels of unidentified metabolites, metabolite i at m/z=215.084 (Fig. 6a) and metabolite ii at m/z=157.0973 (Fig. 6b), compared to NaOH pH 10.18 exposed samples (metabolite i, p < 0.0001; metabolite ii, p < 0.0001) and control samples (metabolite i, F = 2112.3, p < 0.0001; metabolite ii, F = 4865.6, p < 0.0001). Library annotations matched metabolite i to nitrogen-containing heterocyclic compound atrazine (89.6% Q score) and metabolite ii to hydrocarbon 2,7-dimethylnaphthalene (96.5% Q score). Structures of atrazine and 2,7-dimethylnaphthalene are depicted in Fig. 6a and Fig. 6b, respectively. Multivariate analysis also revealed exposure to 0.25 M ammonia increased the relative abundance of pantothenate (F = 28.839, p < 0.001) (Fig. 6c), and amino acids D-allo-isoleucine (F = 28.686, p < 0.001) (Fig. 6d), and alanine (F = 25.391, p < 0.01) (Fig. 6e), compared to those in NaOH pH 10.18 and control conditions.