Impact of Artificial Sputum Medium Formulation on Pseudomonas aeruginosa Secondary Metabolite Production.

In vitro culture media are being developed to understand how host site-specific nutrient profiles influence microbial pathogenicity and ecology. To mimic the cystic fibrosis (CF) lung environment, a variety of artificial sputum media (ASM) have been created. However, the composition of these ASM vary in the concentration of key nutrients, including amino acids, lipids, DNA, and mucin. In this work, we used feature-based molecular networking (FBMN) to perform comparative metabolomics of Pseudomonas aeruginosa, the predominant opportunistic pathogen infecting the lungs of people with CF, cultured in nine different ASM. We found that the concentration of aromatic amino acids and iron from mucin added to the media contributes to differences in the production of P. aeruginosa virulence-associated secondary metabolites. IMPORTANCE Different media formulations aiming to replicate in vivo infection environments contain different nutrients, which affects interpretation of experimental results. Inclusion of undefined components, such as commercial porcine gastric mucin (PGM), in an otherwise chemically defined medium can alter the nutrient content of the medium in unexpected ways and influence experimental outcomes.

INTRODUCTION

Cystic fibrosis (CF) is a genetic disease characterized by the accumulation of thick and viscous bronchial mucus, which provides a nutrient-rich environment for the growth of opportunistic pathogens (1). Pseudomonas aeruginosa, the primary bacterial opportunistic pathogen infecting the lungs of people with CF, contributes to lung function decline by growing at high density within biofilms that are recalcitrant to the innate immune response and antibiotic therapy (2–4). While the nutritional environments of most sites of infection are poorly defined, the composition of CF pulmonary mucus is known to be comprised of amino acids, mucin, DNA, lipids, and micronutrients, including metals (5–9). To study the microbial physiology of P. aeruginosa and other pathogens in vitro, a variety of artificial sputum media (ASM) have been created that aim to replicate the nutrient availability of the CF lung (9–20). The various compositions of ASM were assembled from direct measurement of nutrients from clinical samples, personal observations, literature precedent, and economic cost of materials.

Nine ASM are regularly used to investigate microbial physiology, biofilm development, antibiotic susceptibility, and interspecies interactions in the context of CF (see Table S1 in the supplemental material) (9–17). These ASM can be separated into two lineages based upon the foundational formulations (Fig. 1; Table S2), Soothill-derived ASM and the synthetic CF medium (SCFM) series. Soothill ASM was established in 1997 from literature concentrations for key nutrients identified in CF sputum, including mucin, DNA, sodium, potassium, chloride, and lipids (10). Egg yolk emulsion was used as a source of lecithin, and the chelator diethylenetriamine pentaacetate (DTPA) was used to bind free iron in the medium, as most iron is bound to host proteins in vivo. Soothill ASM was modified by the addition of 250 mg/liter of each of the 20 canonical amino acids to generate Romling ASM, as amino acids are present in CF sputum (12). The ASM published by Kirchner et al. in 2012 (Winstanley) is compositionally identical to Romling except for its buffering capacity (15). To reflect recent measurements of nutrients from CF sputum, Romling ASM was altered to create ASMDM ASM by adding bovine serum albumin (BSA), increasing mucin concentration, and reducing DNA levels (13). Derived from ASMDM ASM, Cordwell ASM has reduced buffering capacity, Casamino Acids (Cas AA) in place of individual amino acids, and ferritin instead of DTPA to replicate CF sputum iron availability (14). Cordwell ASM was modified via substitution of minimal essential medium (MEM) amino acids (MEM AA) for Cas AA, removal of BSA, and an increase in mucin concentration to make San Diego State University (SDSU) ASM (16).

FIG 1

Publication timeline of ASM formulations. Each formulation is plotted at the year it was first published. Open circles indicate media derived from the Soothill ASM recipe. Closed circles represent the SCFM ASM lineage. Nutrient alteration between formulations are indicated by arrows. +, added; −, removed; ↑, increased; ↓, decreased; →, different source compared to the preceding formulations. AA, amino acids; BSA, bovine serum albumin; Cas AA, Casamino Acids; MEM AA, MEM amino acids; DOPC, 1,2-dioleoyl-sn-glycero-3-phosphocholine; GlcNAc, N-acetylglucosamine, PABA, 4-aminobenzoic acid; NAD+, nicotinamide adenine dinucleotide; ADE, adenine; GUA, guanine; XAN, xanthine; HX, hypoxanthine. Recipes for ASM formulations are in detailed in Data Set S1.

Chemically defined SCFM1 was developed in 2007 from the average concentrations of ions, free amino acids, glucose, and lactate measured from CF sputum samples (9). SCFM1 was supplemented with DNA, mucin, N-acetylglucosamine (GlcNAc), and dioleoylphosphatidylcholine (DOPC) at concentrations similar to those measured from CF sputum to create SCFM2 (17). Interrogation of P. aeruginosa strain-specific genetic requirements for in vitro growth using sputum as a medium identified a set of genes located in the anabolic pathways for the biosynthesis of thiamine, nicotinamide adenine dinucleotide (NAD), purines, folate, branched-chain amino acids, and tryptophan. These results suggested that these metabolites are more available in CF sputum than SCFM2. Therefore, SCFM2 was further modified by the addition of p-aminobenzoic acid, NAD+, adenine, guanine, xanthine, and hypoxanthine to generate SCFM3 (17).

Although each of the nine ASM formulations is distinct in composition, experimental conclusions from different ASM claim to provide depictions of pathogen growth and physiology that are reflective of observations made from patient samples (12–14, 21–30). However, disparate nutrient concentration and bioavailability between ASM formulations likely affect microbial physiology and production of virulence factors. Therefore, we hypothesized that the compositional variations between ASM would be sufficient to promote differential production of P. aeruginosa secondary metabolites. These virulence-associated metabolites play important roles in mediating the interactions of P. aeruginosa with host cells and other microbes as well as acquisition of nutrients from the environment (31).

Herein, we describe the application of comparative metabolomics to P. aeruginosa PAO1 stationary biofilm cultures grown in nine ASM formulations. (9–17). PAO1 is a commonly used laboratory strain whose secondary metabolite production has been well documented (32, 33). We show that differential production of P. aeruginosa secondary metabolites between ASM formulations is driven in part due to variations in aromatic amino acid concentration as well as high concentrations of iron added to the media from the inclusion of commercial porcine gastric mucin (PGM). Importantly, due to their roles in virulence-associated phenotypes, altered secondary metabolite production by P. aeruginosa due to compositional differences likely influences the outcome of microbe-microbe interactions, virulence assays, and comparative assessment of experiments performed in different ASM. Additionally, the influence of crude complex nutrients, like mucin, on experimental outcomes may not be due to the nutrient itself but, rather, other constituents of the mixture.

RESULTS

ASM formulation affects P. aeruginosa PAO1 growth.

To assess the effect of ASM composition on growth, strain PAO1 was cultured in ASM (Fig. 2). Samples were grown and evaluated in the following three batches, with one plate per medium: batch 1, which includes LB, Soothill, SDSU, and Cordwell ASM; batch 2, which includes LB, Romling, Winstanley, ASMDM, and SCFM1 ASM; and batch 3, which includes LB, SCFM2, and SCFM3 ASM. Gross morphology of the biofilms was visualized using a stereo microscope. As expected, the cultures varied in structure and color depending upon the growth medium (Fig. 3A; see Fig. S1 and Table S3 in the supplemental material). Among ASM formulations, PAO1 formed apparent biofilm structures in all but Soothill ASM. The translucent but turbid cultures in Soothill ASM indicate that the lack of exogenous amino acids affects both secondary metabolite production and biofilm formation. Due to their nearly identical compositions, Romling and Winstanley ASM cultures grew as brown macrostructures suspended in solution. The nutrient alterations from Romling to ASMDM ASM strongly influenced the P. aeruginosa phenotype, resulting in a tan-green macrostructure at the air-liquid interface, with large open pores within the structure giving way to turbid growth toward the bottom of the well. Surprisingly, the composition of Cordwell ASM induced a phenotype similar to those observed in Romling and Winstanley ASM; PAO1 grew as a dense brown biofilm without clear structural features. Biofilms in SDSU ASM grew as a light-green structure with tendrils near the center of the well that reached the surface of the medium. For SCFM-series ASM, growth in SCFM1 was turbid and bright green, which was noticeably distinct from biofilms in SCFM2 and SCMF3, which were comparable in structure with film-like growth at the air-liquid interface, aggregates on the bottom of the well, and a marked blue color.

FIG 2

Experimental workflow. (1) P. aeruginosa PAO1 growth and secondary metabolite production in artificial sputum media (ASM) was evaluated in 24-well plates (12 control wells and 12 growth wells per medium). The media were inoculated in the following three batches: batch 1, LB, Soothill, SDSU, and Cordwell ASM; batch 2, LB, Romling, Winstanley, ASMDM, and SCFM1 ASM; and batch 3, LB, SCFM2, and SCFM3 ASM. All plates were covered and incubated statically at 37°C for 72 h. (2) Growth and control wells were photographed. (3) Samples from each batch were mechanically disrupted and aliquoted. (4) A subset of aliquots (n = 3) was used to measure growth (CFU per milliliter) and culture pH. (5) Aliquots of control and growth samples (n = 12) were chemically extracted for metabolomics analysis and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS). (6) Metabolomics data from growth samples were processed for feature-based molecular networking (FBMN) using MZmine 2.0. This data processing yielded two complementary summary files, a feature quantification table and MS/MS spectral summary. The feature quantification table provided molecular ion abundance across all samples. The MS/MS spectral summary provided a representative MS/MS spectrum for each molecular ion in the feature quantification table and was used to determine compound identity using molecular networking. (7) MetaboAnalyst v5.0 was used to analyze the feature quantification table using statistical methods. (8) The MS/MS spectral summary was analyzed via the Global Natural Products Social Molecular Networking (GNPS) online data analysis platform using the FBMN workflow to organize structurally related features into molecular families using MS/MS spectral similarity as a proxy for structural similarity. Spectral similarity was also used to annotate features as metabolites based on matches to the GNPS spectral libraries. Abundance data from the feature quantification table was overlaid onto the molecular network as pie charts to visually compare the relative amount of each feature present in samples from each ASM formulation.

FIG 3

P. aeruginosa PAO1 growth in ASM. (A) Representative phenotypes of PAO1 grown statically in 2 ml of each ASM formulation (n = 12) in 24-well plates for 72 h at 37°C. Photographs of replicate control and growth wells are in Fig. S1 in the supplemental material. (B) pH measurement of representative replicates (n = 3) after mechanical disruption of growth in ASM. (C) Representative replicates of PAO1 growth (n = 2 or 3) in ASM.

Importantly, gradients of pH have been shown to dictate microbial niche partitioning in in vitro studies (34). The phenotypic differences between ASM cultures were not driven by variations in pH. As these cultures were incubated statically, the lack of flow of nutrients in our system may have led to a buildup of metabolites that increase the pH (Fig. 3B). The alkalinity of samples after PAO1 growth and sample disruption suggests the nutritional profile of ASM or the culture conditions may not fully replicate the in vivo environment. Notably, only Soothill ASM cultures became acidic, which suggests that metabolites produced from amino acids may contribute to changes in pH.

PAO1 growth in ASM was enumerated as CFU per milliliter (Fig. 3C). While growth was robust, quantification was impacted by batch effects. In batch 1 (LB, Soothill, SDSU, and Cordwell ASM), P. aeruginosa grew to ∼1014 CFU/ml. However, growth in batch 2 (LB, Romling, Winstanley, and ASMDM ASM) was measured at ∼1010 CFU/ml. As biofilm disruption was affected by culture viscosity and aggregation, it is impossible to distinguish whether the difference in growth between batches is due to ASM composition or processing of samples.

Identification of secondary metabolites produced in ASM.

Although several P. aeruginosa secondary metabolites have authentic commercial standards, most must be identified using other methods, such as spectral matching to metabolite databases. Therefore, to quantify and identify the secondary metabolites produced by PAO1 in ASM, we performed feature-based molecular networking (FBMN) on untargeted liquid chromatography-tandem mass spectrometry (LC-MS/MS) metabolomics data (Fig. 2; Table S4). (35, 36). FBMN combines the quantitative assessment of feature finding from traditional metabolomics analysis with MS/MS molecular networking, which organizes metabolites into structurally related molecular families. Feature finding is a data reduction step in LC-MS analysis which creates a quantitative representation (feature table) of all molecular ions, or features, within a metabolomics data set. Each feature is defined by its precursor mass, retention time, and corresponding abundance. In data processing for FBMN, an additional step is applied to the LC-MS/MS data to provide an MS/MS spectral summary in which each feature is associated with an MS/MS spectrum. The MS/MS spectra in the spectral summary are then networked using the Global Natural Products Social Molecular Networking (GNPS) web-based platform (36). Based upon the assumption that molecules of similar structure will have similar MS/MS fragmentation patterns, molecular networking (MN) organizes the MS/MS data into molecular families based upon their spectral similarity. Within the molecular network, nodes (or clusters) represent a feature, and the thickness of the connecting lines (or edges) indicates the level of MS/MS spectral similarity between the features. The MS/MS spectra in the data set are simultaneously scored for their relatedness to the GNPS spectral libraries to provide initial identification of metabolites. The abundance of each feature in each sample group (e.g., individual ASM formulations) is overlaid onto the molecular network as a pie chart, wherein the different colors represent the different samples groups, and the size of a wedge indicates the relative abundance of that feature.

In total, 492 features were quantified and visualized in the molecular network (Fig. S2). Of those features, 80 (16%) had spectral matches to the GNPS libraries with 30 (6%) to phenazine, quinolone, rhamnolipid, and pyochelin reference spectra. Of the five phenazines produced by PAO1, 1-hydroxyphenazine (1-HP), pyocyanin (PYO), phenazine-1-carboxamide (PCN), and phenazine-1-carboxylic acid (PCA) were quantified (Fig. 4A) (37). From the pie chart representation of relative abundance, PYO levels reflect the blue color of SCFM2 and SCFM3, while 1-HP levels were highest in Romling, Winstanley, and ASMDM ASM.

FIG 4

Molecular annotation and relative production of PAO1 secondary metabolites using FBMN. Molecular families of P. aeruginosa secondary metabolites (as determined by MS/MS) identified in the feature-based molecular network are shown. Data includes extracts from PAO1 cultures in all ASM formulations (n = 12 biological replicates per formulation). Nodes represent distinct features (precursor mass and LC retention time). The width of the edges represent the similarity of the MS/MS fragmentation of the connected nodes, while the node color represents the relative feature abundance between ASM formulations. Black circles indicate a match between the data and the GNPS spectral libraries. Nodes annotated as compounds are labeled with an abbreviated compound name. Nodes that could not be assigned structures are labeled with the precursor mass of the feature. The adduct ions detected are indicated for all nodes. (A) Phenazine molecular family. 1-HP, 1-hydroxyphenazine; PYO, pyocyanin; PCN, phenazine-1-carboxamide; PCA, phenazine-1-carboxylic acid. (B) Quinolone molecular family. HHQ, 2-heptyl-4-quinolone; NHQ, 2-nonyl-4-quinolone; UHQ, 2-undecyl-4-quinolone; PQS, 2-heptyl-3-hydroxy-4-quinolone; HQNO, 2-heptyl-4-hydroxyquinoline N-oxide; NQNO, 2-nonyl-4-hydroxyquinoline N-oxide; UQNO, 2-undecyl-4-hydroxyquinoline N-oxide. All other nodes in the molecular family are annotated by alkyl chain length (e.g., C5, alkyl chain containing 5 carbons) and unsaturation (e.g., C9:1, alkyl chain containing 9 carbons and 1 double bond), followed by structural subclass assignment (-HQ, -PQS, and -QNO) to denote structural relatedness to HHQ, PQS, and HQNO, respectively. The structural subclass of nodes labeled with QNO/PQS could not be distinguished based upon MS/MS fragmentation pattern. Annotations with -OH designation indicate hydroxylation at unknown position within the molecule. (C) Rhamnolipid molecular family. Nodes corresponding to rhamnolipids are annotated by rhamnose units (Rha) followed by fatty acid chain lengths (e.g., Rha-C10-C10 is a rhamnolipid containing one rhamnose unit and two fatty acids containing 10 carbons). (D) Pyoverdine (PVD; left) and pyochelin (PCH; right) molecular families.

Of over 50 quinolones reported in literature, 43 were quantified from ASM cultures, including functionally characterized 2-heptyl-4-quinolone (HHQ), 2-heptyl-3-hydroxy-4-quinolone (Pseudomonas quinolone signal; PQS), and 2-heptyl-4-hydroxyquinoline N-oxide (HQNO) (Fig. 4B) (38, 39). As most of the structural variation of the quinolones was due to differences in alkyl chain length, the quinolones quantified were categorized into structural subclasses based upon their MS/MS spectral similarity to HHQ, PQS, and HQNO (38, 40). Twenty-three quinolones were categorized as HHQ-type, six were grouped as HQNO-type, and four were assigned as PQS-type quinolones. Twelve nodes could not be assigned to a specific structural subclass from the experimental data. The relative quantification displayed on the nodes demonstrated that PAO1 quinolone production was heterogeneous among ASM and did not align with substructure categories.

Structurally, P. aeruginosa rhamnolipids consist of 1 to 2 rhamnose units bound to 3‐(3‐hydroxyalkanoyloxy) alkanoic acids (HAAs), yielding mono- and di-rhamnolipids, respectively (41). Six rhamnolipids were identified, including three mono-rhamnolipids and three di-rhamnolipids with fatty acid HAAs of C10-C10, C10-C12:1, and C10-C12 (Fig. 4C). Unlike phenazine and quinolone levels, all rhamnolipids showed similar production patterns between ASM cultures.

Nodes in the molecular network were identified as pyochelin by a spectral match to the GNPS libraries (Fig. 4D). As illustrated in the pie chart visualization, pyochelin was primarily produced by PAO1 in SCFM1. Other nodes with similar coloration within the molecular network were prioritized for identification, as the production pattern suggested that these nodes corresponded to features related to pyochelin or under similar regulation. Manual assessment of the precursor masses, MS/MS spectra, and retention times led to the identification of pyoverdine E (succinamide isoform) and Fe3+-bound pyoverdine E (40). The identification of m/z 661.3262 could not be determined from the experimental data. Further analysis of the FBMN also identified ferribactin (succinamide isoform), the nonribosomal peptide cytoplasmic precursor to pyoverdines (Table S4) (42). Although pyoverdine biosynthesis was reported to be upregulated in P. aeruginosa SCFM1 cultures, the production of the pyoverdine molecular family by PAO1 is below the limit of detection when cultured in most media for the LC-MS/MS method used in this study (9).

PAO1 metabolome changes in response to ASM formulation.

To determine which molecular features differentiated the metabolome of PAO1 in different ASM, the feature table was analyzed using principal-component analysis. From the initial results (Fig. S3), the clustering of samples in the scores plot was due to batch effects, primarily driven by the abundance of benzethonium in batch 1 samples. Benzethonium is an ingredient of bactericidal surface disinfectants. As this compound was introduced during sample processing, it was removed from the data set, and the data were reanalyzed.

A principal-component model was created with five principal components. Extracts of replicate cultures grown in the same ASM clustered together within the two-dimensional space created by principal component 1 (PC1) and principal component 2 (PC2) of the scores plot, which explained 48.6% and 21.5% of the variation between samples, respectively (Fig. 5A). Four sample groups emerged from the principal-component analysis. Soothill ASM samples clustered within their own group. Samples from compositionally related media clustered together, with Romling, Winstanley, and ASMDM ASM samples in one group and samples from SCFM2 and SCFM3 in another. The fourth group consisted of samples from Cordwell, SDSU, and SCFM1 ASM. However, this grouping was influenced by the two-dimensional projection of the data as these media separate in principal component 3 (PC3). This result indicated that the metabolite profiles of replicates of PAO1 cultures in individual ASM are highly similar, and, unsurprisingly, nutrient composition of the ASM is likely the primary driver of differences in PAO1 secondary metabolite production.

FIG 5

Principal-component analysis and its loadings plots of untargeted metabolomics data of PAO1 grown in ASM formulations. (A) Principal-component analysis scores plot of PAO1 samples, colored by medium and drawn with 95% confidence ellipses. The three processing batches are indicated by shapes. Data points representing the 12 replicate samples per medium were closely clustered and distinct clusters were observed for metabolite profiles from cultures in (i) Romling, Winstanley, and ASMDM ASM; (ii) SCFM2 and SCFM3; (iii) Cordwell, SDSU, and SCFM1 ASM; and (iv) Soothill ASM. (B) Loadings plot from untargeted LC-MS/MS-based principal-component analysis of PAO1 in ASM formulations. Features are labeled by their identification number from the feature quantification table and were identified using the FBMN results (Fig. 4). Features of unknown structure are labeled with their precursor mass. Features along the x axis (loadings 1) were responsible for the separation observed along the horizontal axis (PC1) of panel A. Features along the y axis (loadings 2) were responsible for the separation observed along the vertical axis (PC2) of panel A.

Principal-component loadings estimate the extent to which each metabolite contributes to each of the principal components, with magnitude in the two-dimensional space signifying the contribution of a specific metabolite to sample clustering in the scores plot (Fig. 5B). The separation of sample groups in the scores plot is due to several highly abundant metabolites, including PYO, PCA, HHQ, HQNO, PQS, and a set of unknown metabolites, including feature 13 (m/z 336.1201; C18H19NO4 [M+Na]). PYO contributed positively to the separation of both PC1 and PC2, suggesting that its abundance influenced the clustering of SCFM2 and SCFM3 samples from the other ASM (Fig. 6A). PCA also contributed positively to PC2 but to a lesser extent than PYO. Lower production of PCA by PAO1 in Soothill and SCFM1 promotes separation of these samples along PC2 (Fig. 6B). Feature 13 and other related unknown metabolites were detected in smaller amounts in Romling, Winstanley, and ASMDM samples (Fig. 6C) contributing to separation along PC1.

FIG 6

Differential levels of features in ASM formulations identified in the principal-component analysis. Total normalized area under the curve (AUC) for metabolites identified from the principal-component analysis in samples from PAO1 grown in each ASM. Box plots represent the25 to 75th percentiles, with a line at the median. Error bars indicate the minimum to maximum. Individual sample values shown (n = 12 biological replicates per ASM formulation). (A) Pyocyanin (PYO). (B) Phenazine-1-carboxylic acid (PCA). (C) 2-heptyl-4-quinolone (HHQ). (D) 2-heptyl-4-hydroxyquinoline N-oxide (HQNO). (E) Pseudomonas quinolone signal (PQS). (F) A feature of unknown structure with a precursor mass of 336.1201.

While multiple quinolones promoted grouping of ASM samples in the scores plot, HHQ, HQNO, and PQS abundance had the largest effect. HHQ abundance was highest in Soothill, Romling, Winstanley, and ASMDM samples, which contributed to the separation of these ASM along both PC1 and PC2 (Fig. 6D). The amplitude of HQNO in the loading plot suggested that it contributed positively to sample separation along PC2, with a smaller effect along PC1. HQNO abundance was higher in Winstanley, Cordwell, SCFM2, and SCFM3 ASM cultures than in Soothill, Romling, ASMDM, and SCFM1 ASM cultures (Fig. 6E). The lack of tight clustering of SDSU samples in the scores plot is partially explained by the variability of HQNO measured. PQS contributed negatively to PC1 and positively to PC2, aligning with the position of Romling, Winstanley, and ASMDM samples in the left upper quadrant of the scores plot. Confirming the impact of PQS on the sample clustering within scores plot, PQS abundance was higher in Romling, Winstanley, and ASMDM cultures (Fig. 6F).

ASM formulation influences PAO1 secondary metabolite family production.

Principal-component analysis identified a small number of highly abundant phenazines and quionolones as differentially produced by PAO1 in different ASM. However, rhamnolipids and the siderophores pyochelin and pyoverdine were also identified in the molecular network. Therefore, to associate these metabolites with nutrient differences between ASM, we performed an analysis of variance (ANOVA) on the feature table. However, the number of features that had statistically different levels between ASM cultures were too numerous to investigate further (Data Set S2). Alternatively, we quantified the total levels of the phenazines, quinolones, rhamnolipids, and siderophores (Fig. 7).

FIG 7

Differential levels of P. aeruginosa secondary metabolite molecular families in ASM. Total normalized area under the curve (AUC) for each secondary metabolite molecular family in each ASM. Box plots represent the 25 to 75th percentiles, with a line at the median. Error bars indicate the minimum to maximum. Individual samples values shown (n = 12 biological replicates per ASM formulation). (A) Phenazines. (B) Quinolones. (C) Rhamnolipids. (D) Siderophores, including pyoverdine and pyochelin.

Total phenazine levels were lowest in Soothill ASM samples (Fig. 7A). Proportionally, PCA accounted for ∼90% of the total phenazines produced by PAO1 in Soothill ASM (Fig. S4). All other Soothill-derived ASM had increased levels of phenazines compared to Soothill ASM cultures, with higher proportions of 1-HP and PYO. SCFM1 cultures also had relatively low production of phenazines compared to the other ASM. However, SCFM2 and SCFM3 cultures had higher phenazine production than SCFM1 while maintaining the same percentage of each phenazine.

In all ASM, the quinolones produced by PAO1 were predominantly HHQ type (Fig. S5A), with HHQ, NHQ (C9-HQ), and db:NHQ (C9:1-HQ) as the most abundant forms (Fig. S5B). Soothill and the SCFM-series ASM cultures had similarly low levels of total quinolones (Fig. 7B). Romling, Winstanley, and ASMDM ASM samples had the highest levels of PQS-type and total quinolones (Fig. S5A). SDSU ASM samples had the lowest average abundance of quinolones, and the variation in total quinolone levels between samples reflects the variability in HQNO production (Fig. 6E). Proportionally, HHQ-type quinolones accounted for 49% (Cordwell ASM) to 76% (Soothill ASM) of total quinolones, followed by HQNO-type (15% [ASMDM ASM] to 41% [SCFM2]), and PQS-type (3% [SCFM1] to 17% [Romling ASM]).

Total rhamnolipid levels were similar in all ASM cultures, except for SCFM2 and SCFM3 (Fig. 7C). PAO1 in SCFM2 and SCFM3 produced approximately 2-fold less rhamnolipids than the SCFM1 cultures. The proportion of mono- to di-rhamnolipids produced by PAO1 was generally maintained across all cultures, ranging from 30% (SDSU ASM) to 44% (Romling ASM) mono-rhamnolipids (Fig. S6A). Similarly, the fatty acids incorporated into rhamnolipids were rather homogeneous across all cultures, with 62% (SDSU ASM) to 72% (Soothill ASM) incorporating C10-C10 HHAs (Fig. S6B). The other HHAs incorporated include C10-C12:1 (16% [Soothill ASM] to 22% [Romling ASM]) and C10-C12 (10% [ASMDM ASM] to 20% [SDSU ASM]).

In the FBMN, the pie chart representation of relative metabolite abundance indicated that pyochelin and pyoverdine levels were highest in SCFM1 cultures. Quantification of total siderophore levels across all ASM samples confirmed this observation (Fig. 7D; Fig. S7). As siderophore biosynthesis is primarily regulated by the bioavailability of environmental iron, this result suggested that SCFM1 was the only iron-limited ASM formulation (43).

Commercial PGM contains high levels of iron.

We theorized that PGM, the most complex additive of SCFM2 not present in SCFM1, contained a small molecule that was responsible for the differences in PAO1 secondary metabolite production between the cultures. To confirm the effect of PGM on secondary metabolite production, PAO1 was grown in SCFM1 and SCFM1 complemented with 0.5% PGM (PGM-SCFM1). Gross morphology of the biofilms was visualized using a stereomicroscope, and growth was enumerated as CFU per milliliter. Addition of PGM increased the opacity of the cultures and led to a log fold increase of PAO1 growth in PGM-SCFM1 compared to SCFM1 (Fig. S8). P. aeruginosa secondary metabolites were identified and quantified via FBMN (Table S6). The addition of PGM to SCFM1 was sufficient to alter the total levels of the phenazine, rhamnolipid, and siderophore secondary metabolite families (Fig. 8A). PAO1 produced higher levels of total phenazines in PGM-SCFM1 compared to SCFM1, driven by increase in abundance of 1-HP, PYO, and PCA (Fig. 8B). While the total quinolone levels were the same in both PGM-SCFM1 and SCFM1 cultures, the composition of the quinolone molecular family was altered, with an increase in HHQ and decrease in the C9 HQNO-like quinolone 2-nonyl-4-hydroxyquinoline N-oxide in PGM-SCFM1 samples. Both rhamnolipid and siderophore production by PAO1 was decreased in PGM-SCFM1 compared to SCFM1 cultures.

FIG 8

Effect of 0.5% (wt/vol) porcine gastric mucin (PGM) on PAO1 secondary metabolite production. (A) Fold-change difference in production of phenazine, quinolone, rhamnolipid, and siderophore secondary metabolite molecular families by PAO1 between SCFM1 and SCFM1 supplemented with 0.5% PGM (PGM-SCFM1; n = 4 biological replicates per medium). (B) Fold-change difference in production of phenazine and select quinolone metabolites by PAO1 between SCFM1 and PGM-SCFM1 (n = 4 biological replicates per medium). (C) Concentration of iron in 0.5% PGM and vehicle control (MOPS) (n = 6 technical replicates). Dashed line indicates the concentration of iron in SCFM1. 1-HP, 1-hydroxyphenazine; PYO, pyocyanin; PCN, phenazine-1-carboxamide; PCA, phenazine-1-carboxylic acid; HHQ, 2-heptyl-4-quinolone; HQNO, 2-heptyl-4-hydroxyquinoline N-oxide; NHQ, 2-nonyl-4-quinolone; NQNO, 2-nonyl-4-hydroxyquinoline N-oxide.

A mucin suspension was subjected to untargeted LC-MS/MS metabolomics analysis and inductively coupled plasma MS (ICP-MS) to quantify metals. Molecular networking analysis of PGM lipophilic extracts revealed the presence of amino acids, small peptides, and lipids (Data Set S5). ICP-MS analysis revealed that, in addition to other metals (Table S6), a 0.5% (wt/vol) solution of commercial PGM contains approximately 30 μM iron (Fig. 8C). Addition of PGM to SCFM1 increases the total added iron concentration approximately 10-fold.

DISCUSSION

To study microbial physiology in a CF-like nutrient environment in vitro, many ASM have been created with various levels of important nutrients, such as amino acids, mucin, DNA, and iron. (9–20). However, the effect of the different composition of ASM on microbial physiology is poorly understood. Therefore, we applied FBMN to perform comparative analysis of untargeted metabolomics data collected from PAO1 cultured in nine ASM formulations to determine how variations in nutrient composition affect PAO1 secondary metabolite production. Our data indicate a clear relationship between nutrient composition of ASM and virulence-associated metabolite levels produced by PAO1, driven by variation in the concentration of aromatic amino acids and high concentration of iron from the commercial PGM. Importantly, the differences in total and individual secondary metabolite levels produced by P. aeruginosa between ASM influences the ability to compare experimental conclusions, as these metabolites are associated with quorum sensing and virulence.

Phenazines are redox-active metabolites with important roles in electron cycling, oxidative stress, and iron acquisition (44). While 1-HP, PYO, and PCA have been detected in CF sputum samples, their detection has been spurious (29, 45, 46). The low levels of total phenazines produced by PAO1 in Soothill ASM may reflect a shift of chorismate utilization away from phenazine biosynthesis toward the production of aromatic amino acids, as these amino acids are not present in the medium (37). The relatively low level of phenazines produced by PAO1 in Soothill ASM may reflect a reduced virulence phenotype, as highly virulent P. aeruginosa strains have upregulated transcription of genes related to chorismate biosynthesis (47). Phenazine production by PAO1 was highest in SCFM2 and SCFM3. The blue phenotype of these cultures reflected high production of PYO. GlcNAc, a component of SCFM2 and SCFM3, has been shown to increase PYO production by P. aeruginosa, likely through a mechanism to sense and respond to the presence of peptidoglycan (48). Due to the relatively high levels of PYO produced, interactions of P. aeruginosa with lung epithelial cells or other microbes using SCFM2 or SCFM2 as a culture medium may result in increased killing due to the generation of superoxide radicals and increased biofilm formation (49, 50). 1-HP levels were highest in Romling, Winstanley, and ASMDM ASM samples. These ASM contain the highest concentration of total amino acids, which may serve as a preferential carbon source for P. aeruginosa strain PAO1 (51). Catabolism of amino acids to intermediates of the tricarboxylic acid (TCA) cycle may repress transcription of phzM, which encodes the methyltransferase required for PYO biosynthesis, and shift phenazine production toward accumulation of 1-HP (37, 52). 1-HP is an integral component of the interkingdom interaction between P. aeruginosa and the fungal CF pathogen Aspergillus fumigatus (40, 53). 1-HP generates reactive oxygen and nitrogen species and contributes to A. fumigatus iron starvation by chelating iron. As a result, A. fumigatus upregulates production of its siderophores to compete for environmental iron. Depending upon the ASM formulation used as culture medium, the 1-HP mediated interaction between P. aeruginosa and A. fumigatus may not be captured. In Soothill ASM, PCA accounted for ∼90% of the phenazines produced, and PCA levels were highest in ASM that contained both amino acids and mucins. The levels of PCA in ASM samples likely reflect an increase in iron uptake via the Feo system, as the primary function of PCA is the reduction of Fe3+ to Fe2+ (54, 55).

Quinolones produced by P. aeruginosa are involved in biological processes, including quorum sensing, regulation of gene expression for virulence factor production, and microbial competition (39). HHQ, HQNO, and PQS and their C9 congeners have been measured from CF samples (29, 56–62). Although PQS is heavily studied, HHQ- and HQNO-type quinolones are the most abundant in CF samples and have been investigated as diagnostic biomarkers of P. aeruginosa infection. Quinolone levels were highest in Romling, Winstanley, and ASMDM ASM samples. These media contain the highest concentrations of aromatic amino acids. The elevated levels of tryptophan in these media increase the available pool of the quinolone precursor anthranilate via the kynurenine pathway, and tyrosine and phenylalanine have been shown to increase quinolone levels by P. aeruginosa in SCFM1 through an undefined mechanism (9, 63, 64). HHQ-type quinolones accounted for ∼90% of the quinolones produced in Soothill ASM cultures, while proportionally lower levels of HHQ-type, but higher levels of PQS-type, quinolones were measured from Romling, Winstanley, and ASMDM ASM samples. HHQ and PQS (and their C9 congeners NHQ and C9-PQS) activate the quinolone receptor PqsR, driving autoinduction of their biosynthesis and subsequent production of virulence factors associated with the quinolone quorum-sensing pathways (65). Although PQS-type quinolones were the lowest proportion of quinolones produced across all ASM, the PQS levels detected may have a significant impact on virulence factor production by P. aeruginosa. PQS regulates the expression of ∼182 genes, acts as an iron trap when associated with the outer membrane, and promotes the formation of membrane vesicles (66, 67). In ASM formulations that promote PQS biosynthesis, virulence factor production may be elevated and influence the interpretation of related phenotypes. Unlike HHQ and PQS, HQNO-type quinolone production by PAO1 across ASM had no distinct pattern. HQNO affects the growth of the CF opportunistic pathogen Staphylococcus aureus through induction of the small colony variant phenotype, increased biofilm formation, and altering sensitivity to antibiotics (68–70). Therefore, interspecies competition between P. aeruginosa and S. aureus in different ASM may yield different results due to variations in HQNO-type quinolone levels.

Unlike phenazine and quinolone production, rhamnolipid levels were similar between ASM cultures except for those in SCFM2 and SCFM3. Rhamnolipids are virulence-associated glycolipids that function as biosurfactants, which aid in the transport of PQS, development of biofilms, microbial competition, and infiltration of airway epithelial cells by P. aeruginosa (71). A recent study correlated higher levels of rhamnolipid production to virulence phenotypes of 35 P. aeruginosa clinical isolates (72). Two di-rhamnolipids have been consistently detected in CF sputum and explant lung tissue, Rha-Rha-C10-C10 and Rha-Rha-C10-C12, although other congeners have been measured (29, 57–60). In ASM, similar di-rhamnolipid production was favored with preferential incorporation of C10-C10 HAAs. The production of rhamnolipids is regulated by the rhl quorum-sensing (QS) system (73). The similar levels of rhamnolipids between most ASM cultures suggest that this QS pathway is not heavily influenced by ASM composition.

The most striking difference in secondary metabolite production by PAO1 among ASM formulations was the production of siderophores. P. aeruginosa produces two siderophores to acquire environmental iron: high-affinity pyoverdine and low-affinity pyochelin. Their production levels are regulated by the ferric uptake regulatory (Fur) system, which senses the amount of bioavailable environmental iron (42, 43). Both siderophores have been measured from CF sputum, although pyoverdine production is strain dependent (58, 59, 74, 75). Pyochelin and pyoverdine were abundant only in SCFM1 cultures. As siderophore production is regulated by the concentration of iron in the medium, the disparate levels of siderophores between SCFM1 and SCFM2 cultures suggest that only SCFM1 was iron limited, although the added iron in the two media was the same.

PGM, the most complex and impure additive of SCFM2, was identified as a source of iron in ASM. Although a variety of CF-relevant nutrients, such as phosphatidylcholines, sphingomyelins, and tryptophan, were detected in PGM, our findings demonstrate that the high concentration of iron (∼30 μM), at least in part, drove changes in levels of PCA, rhamnolipids, and siderophores between SCFM1 and SCFM2 cultures (64, 76, 77). Supporting our conclusion, high concentrations of iron have been shown to repress siderophore production and transcription of rhamnolipid biosynthesis as well as increase PCA levels by P. aeruginosa (43, 54, 78). Importantly, the effect of other components of ASM, including the mucins themselves, may be masked by the impact of the high levels of iron in PGM on P. aeruginosa. We suspect that the concentration of small-molecule components in PGM varies between lots and vendors and may be altered over time by microbial contaminants of the commercial product.

Despite limitations associated with crude commercial mucins, including altered physical and chemical properties due to the purification process and the inclusion of mucin-bound molecules, mucin is a critical medium component for replicating mucosal infection environments in vitro (6, 79–81). Of the mucins available to researchers, PGM is the most economic choice. High-purity commercial bovine submaxillary mucins are cost-prohibitive for high-volume or continuous culture studies. Purification of native mucins is both cost- and time-intensive. Only ∼65 mg of PGM can be purified from a single pig stomach (82). To produce 1 liter of ASM containing 0.5% (wt/vol) mucin, at least 77 pig stomachs would be required. Although several groups have established procedures for purifying mucins from commercial preparations, recently developed nonbiological mucin-like polymers show promise toward reflecting the properties of native mucins and may eventually offer a cost-effective alternative (82–87).

While all ASM led to production of secondary metabolites by PAO1, the variations in metabolite levels likely influence the conclusions that can be drawn between experiments performed in ASM of different compositions. To determine associations of nutrient source or concentration to a specific phenotype, SCFM-series ASM containing dialyzed or purified commercial PGM would likely yield the clearest results while remaining economically viable, despite being arduous to prepare. Future modifications of SCFM-series ASM will need to consider nutrients and macromolecules in the CF lung environment that are not present in the current formulations, including bovine serum albumin (BSA), protein-bound iron sources (e.g., ferritin), and bioactive lipids. Additionally, comprehensive characterization of CF samples from patients at different disease stages would support the development of SCFM formulations that emulate the distinct nutrient profiles associated with disease progression and treatment.

MATERIALS AND METHODS

Identification of commonly used ASM formulations.

The Web of Science Core Collection (Clarivate Analytics) was searched for papers with the “TOPIC” terms “sputum medium,” “sputum media,” “cystic fibrosis medium,” “cystic fibrosis media,” “sputa medium,” “sputa media,” “CF medium,” or “CF medium” in the title abstract or keyword fields on 10 May 2020. All terms were searched as written, including quotation marks. The results of this search were analyzed using the Web of Science tools “Create Citation Report” and “Analyze Results.” The Analyze Results tool was used to download data corresponding to Web of Science Categories. Google Scholar (Alphabet) was used to quantify the number of citations for the 9 most commonly used ASM formulations (see Table S1 in the supplemental material).

Preparation of media.

Miller Luria broth (LB; Millipore Sigma) was prepared according to the manufacturer’s instructions. ASM formulations were prepared using the media components, concentrations, and commercial sources in Data Set S1 at the reported pH. ASM were prepared as reported, with the following modifications: no antibiotics were included, 0.1% (wt/vol) agar was excluded from the Cordwell formulation, and porcine gastric mucin type III (PGM; Millipore Sigma) was used in all ASM where mucin was a component (9–20). Prior to use, PGM was suspended in 1× 3-morpholinopropane-1-sulfonic acid (MOPS) buffer (pH 7.0) to make a 10% (wt/vol) suspension and sterilized using a liquid autoclave cycle (<20 min sterilization time), and an aliquot was checked for sterility. PGM was added at the reported final mucin concentration. All media were stored at 4°C, checked for sterility prior to use, and used within 1 month of preparation. PGM-SCFM1 was created by supplementing SCFM1 with 0.5% (wt/vol) PGM.

Sample preparation of P. aeruginosa growth in ASM and LB.

P. aeruginosa strain PAO1 (MPAO1; University of Washington, Seattle) was inoculated from a streak plate into 5 ml LB and incubated overnight at 37°C with shaking at 220 rpm (88). We inoculated 1.98 ml of LB or ASM with 20 μl of P. aeruginosa LB broth culture (optical density at 600 nm [OD600] of 0.05; ∼106 CFU/ml final concentration) in polystyrene 24-well plates. Each medium was evaluated in its own 24-well plate (12 control wells and 12 growth wells per medium). The media were inoculated, disrupted, and extracted in three batches, comprised of 2 to 4 ASM formulations and a parallel LB plate as a control. Batch 1 consisted of LB, Soothill, SDSU, and Cordwell ASM, batch 2 consisted of LB, Romling, Winstanley, ASMDM, and SCFM1 ASM, and batch 3 consisted of LB, SCFM2, and SCFM3 ASM. All plates were covered and incubated statically at 37°C for 72 h. After incubation, gross phenotypes of P. aeruginosa cultures (Fig. 3 and Fig. S1) were photographed using a top-down view with a Stemi 508 stereo microscope with an Axiocam 105 color camera (Zeiss). Following visualization, 100 μl cellulase (Aspergillus niger; 100 mg/ml in sterile water; Millipore Sigma) was added to each well, and aggregates were mechanically disrupted by pipetting using wide-bore pipette tips (USA Scientific), followed by water bath sonication (Branson) for 10 min and a secondary mechanical disruption by pipetting using wide-bore pipette tips. Disrupted samples were aliquoted for measurements of growth (CFU per milliliter and pH [Fig. 3]) and metabolomics analysis. Three representative biological replicates from each culture condition were serially diluted, spotted onto LB agar, incubated at 37°C overnight, and counted to determine CFU per milliliter.

Sample preparation of P. aeruginosa growth in SCFM1 and PGM-SCFM1.

P. aeruginosa strain PAO1 was inoculated from a streak plate into 5 ml LB and incubated overnight at 37°C with shaking at 220 rpm. We inoculated 1.98 ml of SCFM1 or PGM-SCFM1 (pH 7.0) with 20 μl of P. aeruginosa LB broth culture (OD600 of 0.05; ∼106 CFU/ml final concentration) in a polystyrene 24-well plate (88). Both media were evaluated in one 24-well plate (4 control wells and 4 growth wells per medium). The plate was covered and incubated statically at 37°C for 72 h. After incubation, gross phenotypes of P. aeruginosa growth (Fig. S8) were photographed using a top-down view with a Stemi 508 stereo microscope with an Axiocam 105 color camera (Zeiss). Following visualization, aggregates were mechanically disrupted by pipetting using wide-bore pipette tips (USA Scientific). Disrupted samples were aliquoted for measurements of growth (CFU per milliliter and pH; Fig. S8) and metabolomics analysis. All four biological replicates from each culture condition were serially diluted, spotted onto LB agar, incubated at 37°C overnight, and counted to determine CFU per milliliter.

Metabolomics sample preparation.

HiPerSolv Chromanorm-grade ethyl acetate (EtOAc; VWR) and Optima-grade methanol (MeOH; Fisher Scientific) were used. Each sample aliquot was chemically disrupted with an equal volume of 1:1 EtOAc/MeOH. The samples were dried, resuspended in 100% MeOH containing 1 μM glycocholic acid (Calbiochem; 100.1% pure), diluted as needed in 100% MeOH containing 1 μM glycocholic acid, and centrifuged for 10 min at 4,000 rpm (Thermo; Sorvall ST 40R) to remove nonsoluble particulates. ASM samples were diluted as follows: Soothill, not diluted; Romling, Cordwell, SCFM1, SCFM2, and SCFM3, 5-fold dilution; and Winstanley, ASMDM, SDSU, and LB, 10-fold dilution. SCFM1 and PGM-SCFM1 were diluted 10-fold. Fifty microliters of 10% (wt/vol) PGM in 1× MOPS (morpholinepropanesulfonic acid) was chemically extracted in technical triplicate with 200 μl of 1:1 EtOAc/MeOH following the Bligh-Dyer extraction protocol. The samples were dried and resuspended in 100% MeOH containing 1 μM glycocholic acid, centrifuged for 10 min at 4,000 rpm, and diluted 10-fold.

Metabolomics data collection.

Mass spectrometry data acquisition for all samples was performed using a Bruker Daltonics Maxis II HD quadrupole time of flight (qTOF) mass spectrometer equipped with a standard electrospray ionization (ESI) source. The mass spectrometer was tuned by infusion of tuning mix ESI-TOF (Agilent Technologies) at a 3-μl/min flow rate. For accurate mass measurements, a wick saturated with hexakis (1H,1H,2H-difluoroethoxy)phosphazene ions (Apollo Scientific, m/z 622.1978) located within the source was used as a lock mass internal calibrant. Samples were introduced by an Agilent 1290 ultraperformance liquid chromatography (UPLC) system using a 10-μl injection volume for ASM samples and a 5-μl injection volume for PGM samples. Optima-grade (Fisher Scientific) acetonitrile (ACN), formic acid (FA), and water were used for UPLC separation. Extracts were separated using a Phenomenex Kinetex 2.6-μm C18 column (2.1 mm by 50 mm) using a 9 min, linear water-ACN gradient (from 98:2 to 2:98% water/ACN) containing 0.1% FA at a flow rate of 0.5 ml/min. The mass spectrometer was operated in data-dependent positive ion mode, automatically switching between full-scan MS and MS/MS acquisitions. Full-scan MS spectra (m/z 50 to 1,500) were acquired in the TOF-MS, and the top five most intense ions in a particular scan were fragmented via collision-induced dissociation (CID) using the stepping function in the collision cell. Data for PA mix, a set of P. aeruginosa secondary metabolite external standards, were acquired under identical conditions. PA mix contains 10 μM each of the following: 1-hydroxyphenazine (1-HP; Tokyo Chemical Industry), pyocyanin (PYO; Sigma), phenazine-1-carboxylic acid (PCA; Ark Pharm), phenazine-1-carboxamide (PCN; Ark Pharm), 2-heptylquinolin-4(1H)-1 (HHQ; Ark Pharm), 2-heptyl-4-hydroxyquinoline N-oxide (HQNO, Cayman Chemical), 2-heptyl-3-hydroxyl-4-quinolone (Pseudomonas quinolone signal [PQS]; Chemodex), N-butyryl-l-homoserine lactone (C4-HSL; Cayman Chemical), N-(3-oxododecanoyl)-l-homoserine lactone (3-oxo-C12-HSL; Millipore Sigma), and a rhamnolipid mixture (AGAE Technologies; 90% pure).

FBMN of ASM samples.

Bruker Daltonics CompassXport was used to apply lock mass calibration and convert the liquid chromatography-tandem mass spectrometry (LC-MS/MS) metabolomics data from proprietary format to mzXML format. MZmine (version 2.53) was used to perform feature finding on the mzXML files (89). The output files were a feature quantification table and MS/MS spectral summary (mgf). The feature quantification table of molecular ion abundance across all samples was corrected for sample dilution prior to analysis. The MS/MS spectral summary provided a representative MS/MS spectrum for each molecular ion in the feature quantification table. The Feature Networking workflow (version release 22) was applied to the dilution-corrected feature table, MS/MS spectral summary, and a metadata file using the Global Natural Products Social Molecular Networking (GNPS) analysis platform (35, 36). Briefly, the data were filtered by removing all MS/MS fragment ions ± 17 Da of the precursor mass. MS/MS spectra were window filtered by choosing only the top 6 fragment ions in each ±50-Da window throughout the spectrum. The precursor ion mass tolerance was set to 0.05 Da, and the MS/MS fragment ion tolerance was set to 0.1 Da. A network of nodes connected by edges was then created. Edges were filtered to have a cosine score above 0.75. Further, edges between two nodes were kept in the network only if each of the nodes appeared in each other’s respective top 5 most similar nodes. Finally, the maximum size of a molecular family was set to 50, and the lowest-scoring edges were removed from molecular families until the molecular family was below this threshold. The spectra in the network were then searched against the GNPS spectral libraries. The library spectra were filtered in the same manner as the input data. All matches kept between the network spectra and library spectra were required to have a cosine score above 0.6 and at least 5 matched peaks. Sum normalization and mean aggregation were applied. The molecular network (https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=baba5df5ab1b437980f81b4e64bcffc0) was visualized using Cytoscape (version 3.7.1) (90).

FBMN of SCFM1 and PGM SCFM1 samples.

Data were processed as described for ASM samples with the following changes: MS/MS spectra were window filtered by choosing only the top 4 fragment ions in each ±50-Da window throughout the spectrum, precursor ion mass tolerance was set to 0.02 Da; MS/MS fragment ion tolerance was set to 0.02 Da, edges were filtered to have a cosine score above 0.7, the maximum size of a molecular family was set to 100, and matches kept between the network spectra and library spectra were required to have a cosine score above 0.7 and at least 4 matched peaks. The molecular network (https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=5624bf8b3174483a89cf733d08ce0fff) was visualized using Cytoscape (version 3.7.1) (90).

Annotation of features.

Identification of features in the molecular network as metabolites was accomplished by comparing the experimental data (exact mass, MS/MS fragmentation, and retention time) to (i) the data acquired for the PA Mix standards (level 1 annotation), (ii) reported structural data for features with MS/MS spectral matches to the GNPS libraries (level 2 annotation), and (iii) putative structures for features within in a molecular family (sets of MS/MS spectra from structurally similar molecules). Benzethonium was identified by manually matching the experimental MS/MS data to the data in the Metlin metabolite database (91). Identified features are listed in Table S3 (ASM samples) and Table S5 (SCFM1 and PGM-SCFM1).

Annotation of metabolites, substructure class, and molecular families.

Multiple features (different adduct ions, in-source fragments, and/or in-source dimers) can represent the same metabolite. Therefore, features associated with the same metabolite were identified by coelution and mass defect from predicted exact mass of the compound adduct, fragment, and/or dimer. Abundance values for features from the same metabolite were summed to provide the total abundance for a metabolite. Abundance values for metabolites with characteristic MS/MS fragmentation patterns were summed to provide total abundance of substructure classes (e.g., HQNO-type quinolones or C10-C10-containing rhamnolipids). Abundance values for all identified metabolites that were structurally related were summed to provide the total abundance for a molecular family (e.g., quinolones or rhamnolipids). Identified features are listed in Data Set S3 (ASM samples) and Data Set S4 (SCFM1 and PGM-SCFM1).

Metal quantification from PGM.

We diluted 2.5% (wt/vol) PGM suspended in 1× MOPS 1:1 with nitric acid, and it was then digested overnight and diluted 10-fold. Gallium (Ga), 50 ppb, was added to each sample as an internal standard. Metal concentration was measured by inductively coupled plasma mass spectrometry (ICP-MS) analysis on an Agilent 7500cx at the University of Nebraska-Lincoln Spectroscopy and Biophysics Core facility. Samples were analyzed in technical triplicate with a blank wash run between each sample. Only concentrations above the limit of quantification were retained.

Classical molecular networking of PGM extracts.

Bruker Daltonics CompassXport was used to apply lock mass calibration and convert the metabolomics data from proprietary format to mzXML format. A classical molecular network was created from the mzXML files using the Molecular Networking workflow (version release 27) of the GNPS web platform (36). Briefly, the data were filtered by removing all MS/MS fragment ions ±17 Da of the precursor mass. MS/MS spectra were window filtered by choosing only the top 6 fragment ions in each ±50-Da window throughout the spectrum. The precursor ion mass tolerance was set to 0.05 Da and the MS/MS fragment ion tolerance to 0.5 Da. A network was then created where edges were filtered to have a cosine score above 0.7 and at least 4 matched fragment peaks between corresponding nodes. Further, edges between two nodes were kept in the network only if each of the nodes appeared in each other’s respective top 10 most similar nodes. Finally, the maximum size of a molecular family was set to 100, and the lowest-scoring edges were removed from molecular families until the molecular family was below this threshold. The spectra in the network were then searched against the GNPS spectral libraries. The library spectra were filtered in the same manner as the input data. All matches kept between the network spectra and library spectra were required to have a cosine score above 0.6 and at least 4 matched peaks. The molecular network (https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=fc3ee81541aa4a60af38b15f08789620) was visualized in Cytoscape (version 3.7.1) (90). Spectral matches to the GNPS libraries were manually verified by comparing the experimental data (exact mass, MS/MS fragmentation) to reported structures (Data Set S5).

Statistical analysis.

Statistical comparison of molecular ion abundance across all ASM was conducted on the normalized feature table generated during FBMN analysis in MetaboAnalyst 5.0 using ANOVA with Tukey’s correction for multiple comparisons (92). Statistical comparison of metabolite, molecular family, and substructure class abundance across all ASM was conducted in GraphPad Prism (version 8.4.0) using ANOVA with Tukey’s correction for multiple comparisons. For all analyses, P values of <0.05 were considered statistically significant. All statistical results are summarized in Data Set S2.

Data availability.

All mass spectrometry data, including the raw files, mzXML files, metadata tables, MZmine settings, and Cytoscape files are available via MassIVE as follows: MSV000086721 (LB and ASM formulations media control, P. aeruginosa culture, quality control [QC], and P. aeruginosa [PA] standard mix data), MSV000086723 (SCFM1 and PGM-SCFM1 media control, P. aeruginosa culture, QC, and PA mix data), and MSV000087153 (PGM lipophilic extracts).