Peptidoglycomics reveals compositional changes in peptidoglycan between biofilm- and planktonic-derived Pseudomonas aeruginosa

Peptidoglycan (PG) is a critical component of the bacterial cell wall and is composed of a repeating β-1,4–linked disaccharide of N-acetylglucosamine and N-acetylmuramic acid appended with a highly conserved stem peptide. In Gram-negative bacteria, PG is assembled in the cytoplasm and exported into the periplasm where it undergoes considerable maturation, modification, or degradation depending on the growth phase or presence of environmental stressors. These modifications serve important functions in diverse processes, including PG turnover, cell elongation/division, and antibiotic resistance. Conventional methods for analyzing PG composition are complex and time-consuming. We present here a streamlined MS-based method that combines differential analysis with statistical 1D annotation approaches to quantitatively compare PGs produced in planktonic- and biofilm-cultured Pseudomonas aeruginosa. We identified a core assembly of PG that is present in high abundance and that does not significantly differ between the two growth states. We also identified an adaptive PG assembly that is present in smaller amounts and fluctuates considerably between growth states in response to physiological changes. Biofilm-derived adaptive PG exhibited significant changes compared with planktonic-derived PG, including amino acid substitutions of the stem peptide and modifications that indicate changes in the activity of amidases, deacetylases, and lytic transglycosylases. The results of this work also provide first evidence of de-N-acetylated muropeptides from P. aeruginosa. The method developed here offers a robust and reproducible workflow for accurately determining PG composition in samples that can be used to assess global PG fluctuations in response to changing growth conditions or external stimuli.

Peptidoglycan (PG) 3 is an important structural component of the bacterial cell wall that is found in almost all bacterial species. In Gram-negative bacteria, PG is composed of a repeating ␤-1,4-linked disaccharide of N-acetylglucosamine (NAG) and N-acetylmuramic acid (NAM) with an appended highly conserved stem peptide consisting of L-alanine, iso-D-glutamate, meso-diaminopimelic acid, followed by two terminating D-alanine residues. During PG synthesis, NAG-NAMpentapeptide monomer subunits are produced in the cytosol, as nucleotide-linked precursors, and transferred to a membranebound undecaprenyl phosphate before being flipped into the periplasmic space through a well-characterized enzymatic pathway (1,2). Within the periplasm, a collection of penicillinbinding proteins create the rigid PG structure by catalyzing linkages between the disaccharides creating a glycan chain, as well as cross-links between adjacent peptide side chains (3)(4)(5)(6).
Once assembled, newly formed PG undergoes considerable maturation, modification, and/or degradation depending on the growth phase or presence of environmental stressors. Numerous lytic enzymes have been identified that cleave at various locations in the PG structure (7,8). These enzymes include peptidases and amidases that cleave the peptide side chain, as well as glucosaminidases, muramidases, and lytic transglycosylases that cleave the glycan chain. They function in diverse cellular processes including PG turnover, cell elongation, cell division, and the membrane insertion of transmembrane protein complexes such as for flagella and the secretion systems (7,9). Other well-characterized PG modifications include acetylation, deacetylation, or N-glycolylation of the PG glycan chain. These modifications are known to provide resistance to the host immune response during infection and during bacterial competition (10 -12). Lastly, the replacement of the terminal D-alanine with D-serine or D-lactate has been shown to impart resistance to some antimicrobials in Gram-positive organisms (13)(14)(15).
Traditional methods for examining PG composition, such as those developed by Glauner et al. in 1988 (16), have remained virtually unchanged for more than 30 years. These methods include purifying PG sacculi from whole bacterial cells, enzymatically digesting them into individual components called muropeptides, and separating muropeptides by reverse-phase HPLC followed by identification using MS. These techniques are complex and time-consuming; specifically, the identification and correlation of individual MS ions with their corresponding muropeptide structure remains a tedious manual step. Because of this complexity, relatively few in vivo PG compositional studies have been undertaken. More recently, several groups have begun to use feature extraction software (17, 18) or isotope-labeling methods (19) to automate the identification of MS ions and to expand the known composition of PG, primarily in Gram-positive organisms.
Here, we further advance these approaches by incorporating feature extraction with a library of known Gram-negative muropeptide structures. We analyzed the resulting spectra with differential analysis approaches that are typically used to detect and analyze large proteomic data sets. This enables the automated assessment of muropeptide identity and the comparison of muropeptide modifications to quantify changes in PG composition between different growth states. In addition, we incorporated a statistical 1D annotation approach (20) that allows for the comparison of global muropeptide compositional changes.
To test this approach, we assessed Pseudomonas aeruginosa PG composition under two distinct physiological growth conditions. P. aeruginosa is an important opportunistic Gram-negative pathogen that causes significant morbidity and mortality in immune-compromised individuals, such as those with cystic fibrosis. A significant decline in the overall prognosis for cystic fibrosis patients is associated with the switch from acute (motile) to chronic (sessile) P. aeruginosa infections. Chronic infections are linked to the ability of P. aeruginosa to grow as an aggregated community of bacteria commonly referred to as a biofilm (21)(22)(23). Once established, biofilms are difficult to eradicate and demonstrate enhanced resistance to antimicrobial treatments compared with their motile "planktonic" counterparts (24). Understanding the mechanisms that drive biofilm formation and the heightened antimicrobial resistance is critical for the treatment of chronic infections. To date, only one study has examined the effects of biofilm growth on PG composition within a Gram-positive Enterococcus bacterium (19).
Our workflow was designed to highlight both large-and small-scale compositional changes that occur during this physiological transition. Based on our results, we present the most detailed compositional analysis of P. aeruginosa PG to date. We identify a core assembly of PG that is represented in high abundance but that does not significantly change between the two growth states and that likely plays a structural role within the cell wall. We also identified an adaptive assembly of PG that is present in smaller amounts but that fluctuates considerably between growth states in response to physiological changes. The PG modifications that showed significant change during biofilm growth included amino acid substitutions of the peptide side chain, as well as modifications that indicate changes in the activity of amidases, deacetylases, and lytic transglycosylases. Our methodology provides a robust and reproducible workflow for accurately determining PG composition between samples that can be used to assess global PG fluctuations in response to changing growth conditions or external stimuli.

Isolation and purification of peptidoglycan from biofilm and planktonic P. aeruginosa cells
To compare muropeptide composition between biofilm and planktonic P. aeruginosa cells, cultures were processed as described in Fig. 1. We grew both types of cells to match previous proteomic studies (25)(26)(27) that showed significant shifts in proteomes at 96 h compared with earlier time points. These proteome shifts included changes in abundance for proteins related to PG synthesis and modification, such as penicillinbinding proteins and lytic transglycosylases (25)(26)(27). Planktonic and biofilm cultures were initially inoculated with 5 ϫ 10 8 cfu of P. aeruginosa PAO1, and after 96 h of growth, four times more cellular biomass was obtained, which is consistent with observed reductions in growth rate in biofilm cultures (25)(26)(27)(28). The quantity of purified sacculi obtained from the planktonic cultures (58.6 Ϯ 5.7 mg of sacculi) when compared with biofilm cultures (13.7 Ϯ 1.6 mg of sacculi) remained constant per weight of cellular biomass, indicating that the reduction in purified sacculi in biofilm cultures was due to reduced growth rate. In addition, biofilm cultures showed considerable bluegreen pigment production compared with their planktonically grown counterparts (Fig. S1, A and B), likely because of Step 1, muropeptide preparation. Planktonic and biofilm P. aeruginosa cultures were grown in liquid medium or on 2% agar-containing medium, respectively. The peptidoglycan sacculi were extracted from all other cellular components and then digested into individual muropeptides using mutanolysin. Scale bar in middle panel, 2 m.
Step 2, muropeptide identification. Unique muropeptide m/z peaks were identified across the total ion chromatograph at all retention times and were compared with a library of predicted muropeptide m/z. MS/MS was then used to compare fragment sizes to predicted muropeptide fragmentation and repeated as necessary to confirm identity of muropeptide.
Step 3, comparison of muropeptide abundances. Bioinformatic approaches were used to compare the differences in abundance of each muropeptide between culture conditions.

Compositional analysis of Pseudomonas peptidoglycan
increased pyocyanin production (25,29,30). Lyophilized sacculi purified from planktonic cells showed a distinct reddish pigment as compared with the overall white color of the biofilm-derived sacculi (Fig. S1, C and D). To ensure that equivalent amounts of sacculi were used for subsequent MS experiments, PG was quantified and normalized based on muramic acid content because one muramic acid represents one monomer subunit (Fig. S1E). Transmission EM was used to assess sacculi purity and integrity (Fig. S1F).

High-resolution muropeptide analysis
Isolated PG sacculi were digested with mutanolysin (31), and muropeptide composition was assessed using LC-Q-TOF-MS. The total ion chromatograph of single replicates of biofilm-and planktonic-derived PG showed differences in the total peak intensity at various retention times ( Fig. 2A). To examine whether these differences represented genuine deviations in PG composition, we analyzed biological triplicate muropeptide samples (in technical triplicate; nine replicates in total per growth condition). Each sample was analyzed in the same sequence to decrease drift and other deviations in retention time between samples. A principal component analysis (PCA) of the data sets was performed. A PCA plot separates the data based on components inherent to the experiment (e.g. biological or technical differences). Here, a PCA of biofilm-and plank-tonic-derived PG muropeptides (Fig. 2B) demonstrates the largest distinction between the data sets to be the state of growth (component 1, 65.04%). Differences in replicate reproducibility is accounted for in component 2 (19.39%), but to a lesser degree, supports the importance of growth conditions to independently segregate the muropeptides.
To identify individual muropeptides, the corresponding m/z was compared with an in-house prepared library comprised of Ͼ6000 potential muropeptide structures (Table S2) and confirmed using MS/MS. In total, 160 muropeptide structures were identified (Table S1), a selected subset of which is presented in Fig. 3. All of the muropeptides identified in this study were detected in both planktonic and biofilm growth conditions. Therefore, we use ion intensity as a measure of the abundance of individual muropeptides between planktonic-and biofilm-derived PG. This is not necessarily a measure of absolute molar abundance but rather an indication of muropeptide compositional changes between the growth states. However, relative abundance can be used within muropeptide groups (monomers, dimers, etc.) to compare molar stoichiometries, as described in previous studies (18). Of the 160 muropeptides identified, 106 have been assigned a structure with 68 confirmed by MS/MS (Table S1, in bold). However, because of the low abundance and uniqueness of the structures, we were unable to definitely assign structures for 54 muropeptides, which when combined represent ϳ6% of the overall PG composition.

Growth phase-dependent changes in PG composition
A scatterplot of the MS raw intensity values showed that 10 muropeptides were the most abundant and constituted ϳ70% of the overall PG composition within both growth conditions ( Fig. 4A and Table S1, top 10 rows). The most abundant monomers were AEmA (which lacks the fifth alanine), AEm (which lacks the fourth and fifth alanine), and AEmAG (where the fifth alanine is substituted with a glycine) (Fig. 3). The most abundant cross-linked dimer was the muropeptide AEmA-AmEA (3)(4), where (3)(4) represents the formation of a peptide bond between the meso-diaminopimelic acid residue of one monomer and the alanine in the fourth position on the second monomer (Fig. 3). This dimer was also found in high abundance with both an 1,6-anhydro-NAM (anh) modification producing AEmA-AmEA (anh) (3)(4), and ␤-1,4 -linked to a NAG-NAM disaccharide, (NAG-NAM)-AEmA-AmEA (3-4) ( Table S1). The second most abundant dimer in this study was the crosslinked AEmA-mEA (3)(4) which has the fourth alanine on one monomer removed (Fig. 3). The trimer representing three cross-linked AEmA monomers, AEmA-AmEA-AEmA, and the anh-modified form, AEmA-AmEA-AEmA (anh), were also highly abundant in both growth conditions (Table S1). The remaining highest abundant muropeptide was AQm-mEA (ami), where the glutamic acid is changed to glutamine on one monomer, and the NAG-NAM disaccharide on one monomer was removed (Table S1).
The remaining 150 muropeptides identified were present in lower abundance ( Fig. 4A and Table S1). To reduce the influence of the highly abundant muropeptides, the intensity of each muropeptide was normalized to the median intensity of that

Compositional analysis of Pseudomonas peptidoglycan
muropeptide across all replicates of both growth conditions. A scatterplot of the normalized intensity displays the relative abundance of each muropeptide between the two growth conditions (Fig. 4B). Hierarchical clustering analysis of the normalized intensity demonstrates distinct sets of muropeptides that show a change in abundance in biofilm-derived compared with planktonic-derived PG (Fig. 4C). These changes in PG composition in biofilm-cultured P. aeruginosa can be clustered to indicate that ϳ80 of the 160 muropeptides were decreased in abundance, whereas ϳ50 muropeptides were increased in abundance in the biofilm-derived PG (Fig. 4C).
To determine whether the apparent differences in abundance were significant, a Student's t test was performed. To account for the rate of type I errors in null hypothesis testing when conducting multiple comparisons, we used a Benjamini-Hochberg false discovery rate (FDR) correction for the performed t tests (32). The corrected p value was plotted on a volcano plot against the fold change in abundance between the two growth conditions (Fig. 5A). Twenty-seven muropeptides showed no significant change in abundance between the two conditions ( Fig. 5A, gray circles). In contrast, 133 muropeptides showed a significant change in abundance (Fig. 5A, black circles), where 19 muropeptides demonstrated a Ͼ2-fold increase (Fig. 5A, red area), and 43 demonstrated a Ͼ2-fold decrease in abundance (Fig. 5A, blue area) in biofilm-derived PG (Table  S1).
To determine which muropeptides were differentially represented in biofilm cells, muropeptides that showed a significant change in abundance between the two growth conditions were compiled into categories based on their structural characteristics (as described in Fig. 3) and assessed for enrichment using 1D annotation (Table 1 and Fig. 5, B and C). 1D annotation enrichment tests for every annotation term (e.g. muropeptide category) if the corresponding numerical values (e.g. intensity) Figure 3. Graphical representation of selected muropeptides identified in P. aeruginosa PAO1. Typical muropeptide characteristics included monomers, dimers, trimers, and tetramers. Dimers were determined to be either 3-4 or 3-3 cross-linked between muropeptide monomers. Muropeptide modifications included amino acid substitutions in the fourth and/or fifth position of the peptide side chain, as well as the resultant muropeptide products after deacetylase or amidase activity.

Compositional analysis of Pseudomonas peptidoglycan
have a preference to be systematically larger or smaller than the global distribution of the values for all muropeptides (20). The 1D annotation score represents the direction in which the muropeptide tends to deviate from the overall distribution of all muropeptides (i.e. a positive or negative enrichment of the muropeptide category).
Based on these analyses, P. aeruginosa muropeptides could be divided into two distinct groups. The first group consisted of the most common muropeptides identified to date (16,33) and were generally found in high abundance with only minor fluctuations between the growth conditions (Fig. 5, B and C). Therefore, we refer to these muropeptides as being part of the "core PG." The second group consisted of muropeptides with less well-characterized modifications and were present in lower abundance with a higher degree of variability between growth conditions (Fig. 5D). We refer to these muropeptides as being part of the "adaptive PG."

Muropeptide composition of the core PG
The core PG constitutes the major proportion of the Gramnegative sacculi and consists of muropeptides that have roles in PG synthesis or turnover (16,33). These included modifications indicative of carboxypeptidase reactions producing the AEm and AEmA; D,D-transpeptidase and L,D-transpeptidase reactions producing the 3-4 and 3-3 cross-links, respectively, as well as lytic transglycosylase activity producing the NAM anh that represent terminal ends of the PG glycan chain (34). Based on fold changes and 1D annotation, we did not observe enrichment of these core muropeptide characteristics between the two growth conditions (Fig. 5B). The relatively low variability of these high-abundance muropeptides suggests that they are involved in maintaining the structural viability of the overall PG sacculi.
To delve deeper into the core PG composition, we reassessed the core muropeptides to group them by differences in crosslinking (including 3-4 and 3-3), as well as variations in anh modification (Fig. 5C). When separated by these parameters, higher-fold change differences were seen overall; however, only two categories demonstrated a significant change between the growth conditions. These categories were enriched in biofilmgrown cells and consisted of the AEmA-AmEA (anh) dimer . Each dot represents the intensity of a single muropeptide recorded in both biofilm-derived (x axis) and planktonic-derived (y axis) peptidoglycan. B, scatterplot of the intensity of individual muropeptides normalized to the median intensity value across all samples. Each dot represents the normalized intensity of a single muropeptide recorded in both biofilm-derived (x axis) and planktonic-derived (y axis) peptidoglycan. C, hierarchical clustering analysis of the normalized intensity values. Each bar represents the abundance of one muropeptide in both the biofilm-and planktonic-derived peptidoglycan.

Compositional analysis of Pseudomonas peptidoglycan
and 3-3 linked (anh), which were enriched 1.3-and 3.6-fold, respectively (Table 1 and Fig. 5C). It is possible the anh-modified 3-3-linked AEmA-AmEA was responsible for the enrichment of the AEmA-AmEA (anh) muropeptide category (Table  1). We could not individually test the enrichment of 3-4-and 3-3-linked AEmA-AmEA using 1D annotation because of the limited number of muropeptides with these characteristics. However, the abundance of the single muropeptide, AEmA-AmEA (anh) (3)(4), did not change between the two growth conditions (Table S1). Lytic transglycosylases form the anh product when cleaving the glycan chains to loosen the overall rigid PG structure to carry out numerous biological processes (7,34,35). There are six families of lytic transglycosylases based on sequence similarity (34,36,37), with each family potentially participating in distinct cellular functions (38). These include cellular processes, such as elongation, division septum formation, PG turnover, and the incorporation of macromolecular complexes, such as flagella and secretion systems (34,38). In addition, the resulting anh muropeptides, such as AEm (anh), when removed from the PG during turnover, are well-known signal molecules that regulate the activity of various resistance mechanisms (39,40). Also, the loss of lytic transglycosylase activity has been shown to reduce biofilm formation in a number of bacterial species (41)(42)(43). Therefore, the relevance of the anh formation on the 3-3 cross-link produced by L,D-transpeptidases that are a known bypass resistance mechanism for ␤-lactam antimicrobials (44, 45) may be of interest.

Muropeptide composition of the adaptive PG
The adaptive PG constitutes a smaller proportion of the overall PG; however, these muropeptides undergo significant fluctuations in abundance between the different culture conditions. It is estimated that Escherichia coli has ϳ3.5 ϫ 10 6 muropeptides/cell (46). Therefore, an individual muropeptide comprising 1.0% of the overall PG (Fig. 4A, light gray box) would amount to ϳ35,000 copies/cell. Even a muropeptide comprising 0.01% of the total PG composition, which represents the lowest abundant muropeptide detected (Table S1), would still equal ϳ350 copies/cell. The low-abundance adaptive muropeptides exhibited the highest fold difference between the biofilm and planktonic cells (Fig. 5D). It is the variability of this dynamic pool of diverse muropeptides that could allow for flexibility in response to growth or stimuli, whereas their low abundance would not affect overall structural viability. In P. aeruginosa, the adaptive muropeptide 1D annotation categories that showed distinct abundance differences between the two growth conditions were the amino acid substitutions of the peptide side chain, amidase activity that separates the NAM from the peptide side chain, and de-N-acetylation of the glycan chain (Fig.  5D).
The majority of Gram-negative muropeptides consist of peptide side chains with D-alanine as the terminal amino acid in positions 4 and 5 (Fig. 3). However, other amino acid substitutions can occur, with glycine as the most well-known (16,(47)(48)(49). In this study, we identified glycine, arginine, cysteine, The red area indicates muropeptides that demonstrate a 2-fold or greater significant increase in biofilm-derived peptidoglycan. The blue area indicates muropeptides that demonstrate a 2-fold or greater significant decrease in biofilm-derived peptidoglycan. B, heat map of fold-change differences between biofilm-and planktonic-derived PG of core muropeptide characteristics. C, heat map of fold-change differences between biofilm-and planktonic-derived PG of detailed core muropeptide characteristics including variations in cross-linkages and variations in anh production. D, heat map of fold-change differences between biofilmand planktonic-derived PG of adaptive muropeptide characteristics. * represents muropeptide characteristics enriched following 1D annotation (Student's t test p value Ͻ 0.05; FDR Ͻ 0.05, s 0 ϭ 1). # represents muropeptide characteristics with only one or two muropeptides within the category and could not be analyzed for significance using 1D annotation.

Compositional analysis of Pseudomonas peptidoglycan
leucine, lysine, methionine, phenylalanine, threonine, tryptophan, tyrosine, and valine incorporated into the peptide side chain at either the fourth and/or fifth position (Table S1). Glycine at position five (AEmAG) was the most abundant amino acid-substituted muropeptide. AEmAG was one of the 10 most abundant muropeptides (Table S1) and decreased in abundance 1.4-fold in biofilm-derived PG. Glycine was also found substituted in position four (AEmG) in very low abundance but demonstrated a much higher 4.0-fold decrease in biofilm-derived PG (Table S1). Overall, amino acid substitutions were reduced 1.5-fold in biofilm-derived PG compared with planktonic-derived PG. Amino acid substitutions on monomers were depleted 1.4-fold (Table 1 and Fig. 5D), whereas dimers containing one amino acid-substituted monomer were depleted 2.4-fold in biofilm-derived PG (Table 1 and Fig. 5D). We did not identify any trimers or tetramers with amino acid substitutions. However, we did identify numerous, as-yet-unknown, m/z peaks that were consistent with potential multimers, some of which may have amino acid substitutions (Table S1). We next determined whether there was a preference for the position on the peptide chain that was substituted. Substitution of the fifth amino acid (AEmAϩ) was approximately three times more prevalent than substitution of the fourth position (AEmϩ) ( Table 1). However, it was the lower abundant AEmϩ that had the greatest variation between the two growth conditions, with 3.4-fold depletion compared with 1.4-fold depletion of the AEmAϩ (Table 1 and Fig. 5D). Only two muropeptides were identified with substitution of both the fourth and fifth amino acid (AEmϩϩ), which was AEmKR, found as both a monomer and as part of a dimer.
The incorporation of alternate amino acids into PG is generally dependent on their concentration within the growth medium (47,50). Biofilms are known to form internal gradients and contain nutrient-limited zones (51)(52)(53). Therefore, the observed reduction of amino acid-substituted muropeptides in the biofilm could be a result of limited availability of these com- Table 1 Abundance differences of core and adaptive muropeptide characteristics identified in P. aeruginosa PAO1 between biofilm-and planktonicderived PG ND indicates categories with only one or two representative muropeptides, 1D annotation could not be performed.

Compositional analysis of Pseudomonas peptidoglycan
pounds within the biofilm. Additionally, exogenous application of D-amino acids can incorporate unique amino acids into muropeptides (50,54,55) and have been shown to reduce biofilm formation or disperse preformed biofilms (54, 56 -58). Exogenous amino acids can also act synergistically with antimicrobials, producing an increased killing effect (56,59). Furthermore, some bacteria (including P. aeruginosa) excrete D-amino acids into the medium in stationary phase (55). The amino acids excreted by P. aeruginosa into the media (55) correlates well with the amino acid-substituted muropeptides that we identified as increased in abundance in the planktonic-grown culture.
Using the MS-based methods presented here, we are unable to determine whether the amino acid-substituted muropeptides contain D-or L-conformations. However, it is possible that the addition of nonstandard amino acids into PG decreases its stability, in turn increasing the susceptibility to antibiotics. This may, in part, contribute to the differences in antibiotic susceptibility observed between planktonic cells and biofilms. Both penicillin-binding proteins and L,D-transpeptidases (responsible for the production of 3-3 cross-links) have been shown to incorporate D-amino acids into PG (55, 60 -62). In Vibrio cholera, L,D-transpeptidases catalyzed the incorporation of D-amino acids preferentially into the fourth position of the muropeptide side chain with other enzymes responsible for incorporation into the fifth position (61). The ortholog of the putative L,Dtranspeptidase, PA2854, in the P. aeruginosa strain PA14 (PA14_27180) was shown to function in the incorporation of D-amino acids into the PG (61). In addition, the abundance of PA2854 increased in the planktonic-cultured P. aeruginosa strain PAO1 proteome at the 96-h time point (25)(26)(27). Therefore, the production of the 3-3 cross-links and the amino acid substitutions may be linked by the activity of a single enzyme, possibly PA2854. Removal of the entire peptide side chain from the NAM is the result of enzymatic cleavage by periplasmic-localized amidase enzymes (7,63). Amidase activity on a dimer results in both a NAG-NAM disaccharide and the remaining cross-linked peptide side chain, e.g. AEmA-AmEA (ami), as residual products (Fig. 3). Previous studies have only detected the amidasecleaved dimer in PAO1 (64), which highlights the sensitivity of our methodology. Within P. aeruginosa PG, muropeptides containing a single NAG-NAM disaccharide did not show a significant change when growing as a biofilm (Fig. 5D). However, muropeptides containing longer NAG-NAM glycan chains of two, three, and four repeats and the amidase residual dimers showed a significant 2.2-and 2.1-fold enrichment, respectively, in biofilm-derived PG (Table 1 and Fig. 5D). Amidases function in combination with other hydrolases to break down the rigid structure of the PG for incorporation of transmembrane complexes, such as the type VI secretion system (9). The type VI secretion system is up-regulated in biofilm and are important for biofilm antimicrobial resistance (65,66). Accordingly, an increase in macromolecular complexes within the biofilm culture could account for the increased PG amidase activity. However, PG amidase activity has been shown to be important for cell division (67), as well as peptidoglycan turnover (68 -70), and amidase abundance is up-regulated by the CpxA-CpxR two component system in response to cell envelope stress in P. aeruginosa (71).
Changes to the acetylation state of the glycan chain are common modifications to Gram-negative PG and generally occurs as either an addition of an acetyl group to the C6 hydroxyl (O-acetylation) or a loss of the acetyl group on the C2 amine (de-N-acetylation) (10,11). P. aeruginosa does not contain the genes necessary for adding the O-acetyl group (72,73). However, PAO1 does possess the pa1517 gene, which is a homolog (26% amino acid identity) of the Helicobacter pylori enzyme responsible for the de-N-acetylation of the NAG (74). In P. aeruginosa, the de-N-acetylation of the glycan chain was depleted 4.0-fold in the biofilm-derived PG (Table 1 and Fig.  5D). NAG de-N-acetylation occurred more frequently, constituting a 4.1-fold decrease in the biofilm-derived PG and was found exclusively on the NAG-NAM disaccharide product of amidase activity (Fig. 3). This could indicate that the P. aeruginosa deacetylase may have an active site that limits substrate binding when the NAM is appended to a peptide side chain. Only a single muropeptide was identified containing a de-Nacetylation modification of the NAM residue, i.e. AEm-mEA (2xNAMdeN) (Fig. 3), which decreased 2.9-fold in biofilm-derived PG. Additional as-yet-unknown muropeptides have m/z values that indicate potential de-N-acetylation modifications, some of which may be NAM de-N-acetylated; however, none were confirmed by MS/MS. Changes to the acetylation of the glycan chain is known to provide protection from hydrolytic enzymes, such as lysozyme, produced by the host immune system during host-pathogen interactions (10,75,76). It is possible that planktonic-grown P. aeruginosa increases the content of the de-N-acetylation in preparation for environmental conditions where host immune stressors and/or bacterial competition may be higher.

Identification of multimeric and unique muropeptides
The high sensitivity of our experimental approach led to the identification of numerous muropeptide trimers and tetramers. The abundance of trimeric muropeptides detected here was ϳ2-fold higher than reported previously for PAO1 (ϳ4%) (3, 77) but was not enriched in either culture condition (Fig. 5B). The fact that we detected not only tetramers but also two anhcontaining forms is a testament to the sensitivity of this methodology. These multimeric muropeptides indicate connection points between several glycan chains and therefore may have significant structural importance. It was also interesting that tetramer formation occurred three times more often at an anh glycan chain termination. Whereas ϳ50% of trimers, ϳ25% of dimers, and only ϳ1% of monomers occurred with an anh modification. The preference for glycan chain termination at a cross-linked muropeptide is also a feature found in H. pylori PG (78) and the PG-derived from stationary phase growth in E. coli (16). Further work into the effects these multimer muropeptides have on the overall ordered structure of PG will be needed.
When examining the abundance of muropeptides grouped by amino acid substitutions, several muropeptides, particularly AEmK, AEmKR, and AQmT-mQA (Fig. 3), were increased in abundance in biofilm-derived PG (Table S1), whereas all other muropeptides containing an amino acid substitution were Compositional analysis of Pseudomonas peptidoglycan depleted as mentioned above (Fig. 5C). AEmK increased in abundance 2.8-fold, AEmKR increased 2.8-fold, and AQmT-mQA increased 2.1-fold in the biofilm-derived PG (Table S1). The uniqueness of these muropeptides could indicate a divergent role for them within the biofilm-grown cultures compared with the remaining amino acid substituted muropeptides. AEmKR is the well-characterized attachment site of PG to the Braun's lipoprotein homolog OprI in P. aeruginosa (16,79,80). The attachment of the PG to the outer membrane through a lipoprotein is thought to increase the stability of the overall cell envelope (81). Therefore, the increased abundance of AEmKR indicates enhanced outer membrane attachment of PG and possibly a more stable cell envelope within biofilm-grown cells, which may increase the overall resistance of the bacteria within a biofilm to antimicrobial treatment.

Conclusions
Coupling high-resolution MS with software capable of handling complex data processing makes high-throughput analysis of PG composition more readily accessible and has been coined "peptidoglycomics" (82,83). However, the use of this technology for the quantitative comparison of PG composition is, as of yet, highly limited (17, 18,84). More specifically, we can now use the semiautomated pipeline described here to compare the detailed makeup of PG under a variety of growth or stress conditions. Using 1D annotation, we can assess the enrichment of muropeptide characteristics to make inferences concerning enzymatic activity. This will potentially provide novel insight as to how physiological and environmental conditions shape the enzymatic pathways involved in the synthesis and modification of a critical cell envelope component.
Gram-negative PG is thought to be relatively consistent in composition between species (33,77,85). We have shown that the Gram-negative bacterium P. aeruginosa PAO1 has a more varied PG composition than previously demonstrated (3,64,77,85). We identified 160 unique muropeptides (Table S1), which represents eight times more muropeptides than previously shown for P. aeruginosa (3,64,77,85) and double the number characterized for E. coli (16). This is the first report of de-Nacetylated muropeptides isolated from P. aeruginosa. In addition, this is the first demonstration of the high diversity and significant quantity of amino acid substituted muropeptides in Gram-negative PG. Therefore, to our knowledge, this approach has produced the most detailed PG compositional analyses of Pseudomonas or any Gram-negative bacterium to date.

Bacterial culture
Starter cultures of P. aeruginosa PAO1 were grown in tryptic soy broth (TSB) at 37°C for 6 h to an approximate optical density at 600 nm of 2.0. The starter culture was used to inoculate 1 liter of TSB for planktonic cells, and 250 ml of TSB ϩ 2% agar (TSA), each with 5 ϫ 10 8 cfu. TSA cultures were allowed to dry briefly after inoculation. Planktonic cultures were grown with shaking at 200 rpm, whereas biofilm cultures were grown static. Both cultures were incubated for 96 h at 37°C. For cell collection, planktonic cultures were centrifuged at 7,000 ϫ g, whereas biofilm cultures were scraped off the surface of the TSA and collected into a 50-ml conical tube. Cell pellets were weighed and stored at Ϫ20°C until further use. All samples were cultured in triplicate.

Isolation of peptidoglycan sacculi
PG sacculi were isolated from planktonic-and biofilmgrown cultures following the method of Glauner (86). Frozen cell pellets were suspended in 250 ml of cold (4°C) 20 mM sodium phosphate buffer, pH 7.5. The suspension was added dropwise to an equal volume of boiling 8% SDS in 20 mM sodium phosphate buffer, pH 7.5, to a final volume of 500 ml in a round-bottomed flask with a water-cooled condenser. The samples were boiled with stirring for 3 h and allowed to cool overnight. The SDS was removed by repeated ultracentrifugation at 70,000 ϫ g for 40 min and washing in 20 mM sodium phosphate buffer, pH 7.5. The concentration of SDS in the supernatant was determined using Stains-All (Sigma-Aldrich), prepared as per manufacturer's instructions; absorbance was read at 510 nm (87) and compared with an SDS standard curve. Washing continued until the concentration of SDS in the supernatant was ϳ0.001%. Washed sacculi were suspended in 20 mM sodium phosphate buffer, pH 7.5, and incubated with 50 g/ml each of amylase, DNase, and RNase for 1 h at 37°C. Next, 100 g/ml Pronase and ϳ0.02% sodium azide were added and incubated for an additional 18 h at 37°C with nutation. Samples were then incubated in 2% SDS at ϳ100°C for 1 h followed by removal of SDS by repeated ultracentrifugation at 70,000 ϫ g for 40 min and washing with H 2 O until SDS concentration was ϳ0.001%. Washed sacculi were suspended in H 2 O and lyophilized. The lyophilized PG sacculi were weighed and stored at room temperature until prepared for quantification and MS.

Quantification of NAM
The quantity of isolated PG was determined for both planktonic and biofilm samples using the method described by Clarke (88). Briefly, 10 mg/ml of purified sacculi were acid hydrolyzed with 6 M HCl for 2 h at 95°C under vacuum in hand-pulled glass ampules. Ampules were opened and dried by heating under vacuum suspended over solid NaOH to neutralize acid. Hydrolyzed PG components were suspended in 10 mM NaOH and injected onto a CarboPac PA-20 column (Dionex Corp., Sunnyvale, CA), equilibrated with 100 mM NaOH for 10 min, and then separated by application of a 0 to 150 mM NaOAc linear gradient over 25 min, with detection using a pulsed-amperometric electrochemical sensor. The concentration of NAM in all samples was determined by peak integration using Chromelean version 7.2.0.3765 (Dionex Corp.) and compared with a muramic acid standard (Sigma-Aldrich).

Preparation of muropeptides for MS
Lyophilized PG sacculi were suspended to 10 mg/ml in H 2 O and sonicated for 20 s at 30% amplitude to disperse clumps. To produce muropeptides, 8 mg/ml of isolated sacculi were digested with 100 g/ml mutanolysin in 100 mM ammonium acetate, pH 5.5, 50 mM MgCl 2 for 18 h at 37°C. Next, an equal volume of 0.5 M borate buffer, pH 9.0, with ϳ10 mg/ml sodium borohydride was added and reacted for 20 min at room temperature. Finally, pH was adjusted to Ͻ4 with diluted (1:5) phos-

Mass spectrometry
LC-MS analyses were performed on an Agilent 1200 HPLC liquid chromatograph interfaced with an Agilent UHD 6540 Q-TOF mass spectrometer at the Mass Spectrometry Facility of the Advanced Analysis Centre, University of Guelph. A C18 column (Agilent AdvanceBio Peptide Map, 50 mm ϫ 2.1 mm ϫ 2.7 m) was used for chromatographic separation with the following solvents; water with 0.1% formic acid for A; and acetonitrile with 0.1% formic acid for B. The mobile phase gradient was as follows: initial conditions 2% B for 2 min increasing to 15% B in 13 min and then to 50% B for an additional 10 min followed by column wash at 98% B and 10 min re-equilibration. The first 2 and last 5 min of the gradient were sent to waste and not the spectrometer. The flow rate was maintained at 0.2 ml/min. The mass spectrometer electrospray capillary voltage was maintained at 4.0 kV, and the drying gas temperature was set at 350°C with a flow rate of 13 liters/min. Nebulizer pressure was 40 p.s.i., and the fragmentor was set to 150 V. Nitrogen was used as both nebulizing and drying gas and collision-induced gas. The mass-to-charge (m/z) ratio was scanned across the m/z range of 200 -2000 m/z in 4GHz (extended dynamic range positive-ion auto MS/MS mode. Three precursor ions per cycle were selected for fragmentation. The instrument was externally calibrated with the ESI TuneMix (Agilent Technologies Inc., Santa Clara, CA). The sample injection volume was 10 l. Triplicate technical replicates were performed for all biological replicates. Selected samples were reinjected for further characterization of ambiguous muropeptides using the same chromatography and general MS conditions as above; however, we incorporated a preferred list of targets for MS/MS at specific m/z values and retention times. These targets were fragmented sequentially with 10, 15, 20, and 25 V collision energy. Raw MS and MS/MS spectra have been deposited into the publicly accessible repository Figshare under number 10277909.

Muropeptide library production
Peptidoglycan muropeptide characteristics selected for incorporation into the library were based on: 1) P. aeruginosa muropeptides published previously (64); 2) known PG compositions in Gram-negative bacteria (11,16,33,64); 3) all possible amino acid substitutions on the stem peptide; and 4) combinations of muropeptide modifications thereof. Structures were produced using ChemDraw Prime version 16.0.1.4 (PerkinElmer Life Science), and the molecular formula for each was manually entered into the MassHunter Personal Compound Database and Library version B.07.00 (Agilent Technologies Inc.). In total Ͼ6000 different muropeptide structures were included (Fig. S2).

Data analysis
Spectra processing was performed using batch untargeted recursive muropeptide feature extraction in MassHunter Profinder version B08.00 (Agilent Technologies Inc.). Profinder recursive feature extraction involves an initial naïve fea-ture finding algorithm, molecular feature extraction, which combines coeluting related ions such as adducts or different charge states into one compound. This list of compounds is then verified in a second round of feature finding with the findby-formula algorithm, which uses the ion m/z values and isotope ratios found by molecular feature extraction to reinterrogate the data. In this way, the Profinder software combines related m/z into single compounds, verifies the extracted compounds, and reduces false positives and negatives. The initial m/z threshold set for feature detection was 300 counts and extraction window of 40 ppm using the Molecular Feature Extraction algorithm. After isotope grouping using the peptide isotope model, the compound threshold was set to 3000 counts in at least two-thirds of samples in one group. For the recursive portion of the feature detection, a list of consensus muropeptides determined from all samples was used to reassess the raw data using the find-by-ion algorithm using a 50 ppm extracted ion chromatograph extraction window. Statistical analysis was then preformed in Mass Profiler Professional (MPP) version 14.9.1 (Agilent Technologies, Inc.). During import of the Profinder prealigned data, coeluting m/z values representing in-source loss of NAG were merged into a single compound. Prior to merging the compounds, MS/MS spectra were assessed to confirm the loss of the NAG. The MPP data were normalized using the NAM quantification and baselined to the median intensity of each individual muropeptide across samples. Significantly different (p value Ͻ 0.05) muropeptides between growth conditions were identified using a moderated t test and corrected using Benjamini-Hochberg FDR. Unsupervised hierarchical clustering of the normalized intensity values was used to group the data. Muropeptides were identified using MPP IDBrowser Identification combined with the muropeptide MassHunter Personal Compound Database and Library and confirmed using the MS/MS spectra (Fig. S2). A list of muropeptides that could not be identified were exported for preferred MS/MS fragmentation and then verified manually. For 1D annotations, the normalized intensity values from MPP of all muropeptides that had a significant FDR Ͻ 0.05 were uploaded into Perseus version 1.6.2.2 (89) and assigned into categories based on the muropeptide characteristics. Each characteristic category was assessed for significance (Student's t test, p value Ͻ 0.05, FDR Ͻ 0.05, s 0 ϭ 1) in Perseus using 1D annotation (20). All graphs were produced in Prism 5.0f (GraphPad Software Inc., San Diego, CA).

Transmission EM of peptidoglycan sacculi
PG sacculi were processed for transmission EM as previously described (90). Sacculi were adhered to 200-mesh carbon-coated copper grids and washed three times with H 2 O, wicking away liquid between washes with #1 Whatman filter paper. The grids were incubated with 0.5% BSA for 5 min before washing again with H 2 O as before. Sacculi were stained with 1% uranyl acetate, air-dried, and viewed with a single-tilt holder on an FEI Tecnai G2 F20 transmission electron microscope operating at an accelerating voltage of 200 kV and equipped with a bottom-mount Gatan 4k chargecoupled device camera.