1. Microbiology and Infectious Disease

PGFinder, a novel analysis pipeline for the consistent, reproducible, and high-resolution structural analysis of bacterial peptidoglycans

  1. Ankur V Patel  Is a corresponding author
  2. Robert D Turner
  3. Aline Rifflet
  4. Adelina E Acosta-Martin
  5. Andrew Nichols
  6. Milena M Awad
  7. Dena Lyras
  8. Ivo Gomperts Boneca
  9. Marshall Bern
  10. Mark O Collins  Is a corresponding author
  11. Stéphane Mesnage  Is a corresponding author
  1. School of Biosciences, University of Sheffield, United Kingdom
  2. Department of Computer Science, University of Sheffield, United Kingdom
  3. Institut Pasteur, Unité Biologie et Génétique de la Paroi Bactérienne, France
  4. INSERM, Équipe Avenir, France
  5. CNRS, UMR 2001 "Microbiologie intégrative et moléculaire", France
  6. biOMICS Facility, Faculty of Science Mass Spectrometry Centre, University of Sheffield, United Kingdom
  7. Protein Metrics Inc, United States
  8. Infection and Immunity Program, Monash Biomedicine Discovery Institute, Australia
  9. Department of Microbiology, Monash University, Australia
https://doi.org/10.7554/eLife.70597
Share this article
Cite this article
  1. Ankur V Patel
  2. Robert D Turner
  3. Aline Rifflet
  4. Adelina E Acosta-Martin
  5. Andrew Nichols
  6. Milena M Awad
  7. Dena Lyras
  8. Ivo Gomperts Boneca
  9. Marshall Bern
  10. Mark O Collins
  11. Stéphane Mesnage
(2021)
PGFinder, a novel analysis pipeline for the consistent, reproducible, and high-resolution structural analysis of bacterial peptidoglycans
eLife 10:e70597.
https://doi.org/10.7554/eLife.70597
  2. Download BibTeX
  3. Download .RIS

Abstract

Many software solutions are available for proteomics and glycomics studies, but none are ideal for the structural analysis of peptidoglycan (PG), the essential and major component of bacterial cell envelopes. It icomprises glycan chains and peptide stems, both containing unusual amino acids and sugars. This has forced the field to rely on manual analysis approaches, which are time-consuming, labour-intensive, and prone to error. The lack of automated tools has hampered the ability to perform high-throughput analyses and prevented the adoption of a standard methodology. Here, we describe a novel tool called PGFinder for the analysis of PG structure and demonstrate that it represents a powerful tool to quantify PG fragments and discover novel structural features. Our analysis workflow, which relies on open-access tools, is a breakthrough towards a consistent and reproducible analysis of bacterial PGs. It represents a significant advance towards peptidoglycomics as a full-fledged discipline.

Introduction

The characterisation of bacterial cell walls started with the development of electron microscopy techniques (Mudd and Lackman, 1941), and it has ever since been the focus of countless studies. The major and essential component of the bacterial cell envelope is called peptidoglycan (PG). It confers cell shape and resistance to osmotic stress and represents an unmatched target for antibiotics (Mainardi et al., 2008; Vollmer et al., 2008). Some of the most widely used antibiotics to date (beta-lactams and glycopeptides) inhibit the polymerisation of PG.

PG (murein; originally known as mucopeptide) is a giant, insoluble, bag-shaped molecule, and its composition was characterised soon after its discovery (Cummins and Harris, 1956; Rogers and Perkins, 1959; Weidel and Pelzer, 1964). It is composed of glycan chains containing alternating N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) residues linked by β,1–4 bonds. The lactyl group of MurNAc residues is substituted by pentapeptide stems which often has the L-Ala1-γ-D-Glu2-L-DAA3-D-Ala4-D-Ala5 sequence, where DAA is a diamino acid such as meso-diaminopimelic (mDAP) acid or L-lysine (Figure 1a; Vollmer et al., 2008). In some species, a lateral chain (with variable composition and length) can be found attached to the amino acid in position 3. Peptide stem composition and polymerisation can vary amongst bacterial species (Schleifer and Kandler, 1972). Whilst PG building blocks produced in the cytoplasm are always the same, the final structure undergoes constant lysis and modification, a process referred to as ‘remodelling’. Both remodelling and alternative polymerisation modes (Figure 1b) lead to a considerable variation in PG structure during cell growth and division. PG structural plasticity plays a critical role for adaption to environmental conditions during host-pathogen interaction (Boneca et al., 2007; Juan et al., 2018) or to survive exposure to antibiotics (Mainardi et al., 2008).

Figure 1
Download asset Open asset
Diversity of peptidoglycan composition and structure.

(a) Representative peptidoglycan building block made of N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) forming a disaccharide subunit linked to a pentapeptide stem attached to the MurNAc via a lactyl moiety. Peptide stem contains both L and D-amino acids and show a great diversity in composition. Some examples of amino acids found in peptidoglycan are shown for each residue. Modifications of the sugars are also shown. (b) Representation of crosslinking diversity, 4–3 bonds (direct or via peptide crossbridge) and 3–3 bonds are made by D,D- or L,D-transpeptidases, respectively. The enzymes catalysing 1–3 and 4–2 bonds remain unknown. Acceptors stems are shown in blue and donor stems in red. DAA: diamino acid; m-DAP: meso-diaminopimelic acid; D-Lac: D-lactate; X: cell surface polymer (e.g teichoic acid); Z: lateral chain.

PG material is straightforward to purify, but the structural analysis of this molecule is challenging and remains a time-consuming and labour-intensive process. The intact molecule must be broken down into soluble fragments by enzymatic digestion with a glycosyl hydrolase (lysozyme), and individual building blocks (disaccharide peptides, also called muropeptides) are analysed to gain insight into the structure of the intact molecule. A transformative step for the characterisation of disaccharide peptides has been the use of reversed-phase HPLC (rp-HPLC) and mass spectrometry (MS) towards the end of the 1990s (Garcia-Bustos et al., 1988; Glauner, 1988; Glauner et al., 1988; Martin et al., 1987). Combining muropeptides separation by rp-HPLC and MS characterisation has hinted at a more complex structure than previously reported.

Despite tremendous advances in both rp-HPLC-MS instrumentation and software development for the automated analysis of large datasets, ‘peptidoglycomics’ is still in its infancy. The experimental strategy to analyse PG structure has barely changed over the past 30 years. Even though rp-HPLC-MS has been routinely used over the past decade, the analysis of MS data remains a black box. Except for a recent study describing the PG structure of Pseudomonas aeruginosa (Anderson et al., 2020b), no information is available in the literature about the strategy used to identify muropeptides in rp-HPLC-MS datasets. This task relies on searching a subset of expected structures, but the complexity of both the search space and search process is often not described.

We previously provided the proof of concept that shotgun proteomics tools can be used for the automated and unbiased analysis of PG structure (Bern et al., 2017). The analysis of Clostridioides (previously Clostridium) difficile PG led to the identification of many muropeptides never reported before. This work also demonstrated that PG analysis could be carried out with relatively high throughput, opening the possibility to analyse large numbers of samples (such as clinical or environmental isolates) using a minimal amount of material (typically microgram amounts).

Here, we describe a novel software called PGFinder for the analysis of MS data. PGFinder is a versatile and straightforward open-source software tool that allows automated identification of muropeptides based on the creation of dynamic databases. Sharing PGFinder as a Jupyter Notebook provides a robust and consistent pipeline with the potential to accelerate discovery in the field of peptidoglycomics. This workflow described here allows a comprehensive description of the analysis strategy for a consistent and reproducible PG structure analysis by users in the community.

We applied the PGFinder pipeline to analyse the muropeptides composition of Escherichia coli which has been extensively studied. We demonstrate that PGFinder can capture an unprecedented level of complexity of PG structure, highlighting the limitations of the search strategies reported so far. Finally, we provide evidence that PGFinder can be used in conjunction with freely available MS data deconvolution software, making PG analysis possible using entirely open-access tools. We propose that our approach represents a significant advance towards a consistent and reproducible analysis of PG structure, allowing peptidoglycomics to take the crucial first leap to parity with other omics disciplines.

Results

PGFinder: a dedicated script for bottom-up identification of PG fragments

No pipeline is currently available for the automated analysis of MS PG data. Therefore, we sought to replicate a shotgun proteomics approach to create an analysis pipeline dedicated to PG analysis, referred to as ‘peptidoglycomics’ (Wheeler et al., 2014).

To limit misidentifications due to mass coincidences, we established a search strategy relying on an iterative process (Figure 2). A first search was carried out using a database made of reduced disaccharide peptides (monomers) and their theoretical monoisotopic masses. MS data were deconvoluted using the Protein Metrics Byos software to generate a list of observed monoisotopic masses alongside other parameters including retention times and signal intensity. Individual theoretical masses contained in the monomer database (Figure 2, database 1) were compared with observed masses in the experimental dataset. Any observed mass within 10 ppm tolerance was considered as a match and the corresponding inferred structure and theoretical mass were then added to a list of matched structures (Figure 2, library 1). As a second step, we used the list of matched monomers to build another database in silico (Figure 2, database 2), corresponding to dimers and trimers and their theoretical masses. Two types of polymerisation events are included in the original PGFinder version depending on the type of crosslink either through peptide stems or glycan chains. Individual theoretical masses from the in silico database were compared to observed masses to generate a list of matched dimers and trimers (Figure 2, library 2). As a third step, we combined the lists of matched monomers and multimers to generate a final library of modified muropeptides (Figure 2, library 3). The final library contained only modified muropeptides corresponding to matched monomers, dimers, and trimers. The modifications accounted for include the presence of anhydro groups, deacetylated sugars, amidated amino acids, and modifications resulting from N-acetylglucosaminidase or amidase activities (loss of GlcNAc and lack of peptide stems, respectively). In-source decay products (loss of GlcNAc) and Na+/K+ salt adducts were also added to library 3. All three libraries corresponding to observed monomers, dimers, trimers, and their modified variants were combined to search the MS data for masses matching theoretical values within a 10 ppm mass accuracy window. This search generated results processed by PGFinder to carry out a ‘clean up step’. The intensities of in-source decay products and salt adducts were combined with that from parent ions when found within close retention time (a 0.5 min time window). The output of this final step is a matched table written to a .csv format file. It contained all the inferred structures identified within the specified mass and retention time windows with an extracted-ion chromatogram (XIC) signal intensity for quantification.

Figure 2
Download asset Open asset
Flowchart outlining the algorithm for the matching script.

The identification of muropeptides was carried out using four successive steps, indicated by different colours (orange, green, blue, and red, respectively). As a first step, observed masses in the dataset are compared to a list of theoretical masses corresponding to monomers (database 1). Matched masses within the ppm tolerance set (10 ppm for Orbitrap data) are used to build a list of inferred monomeric structures and their corresponding theoretical masses (library 1). This is then used to generate a list of theoretical multimers (dimers and trimers) and their masses (database 2). A second matching round is carried out to build a list of inferred multimers (library 2). At this stage, matched monomers and multimers are combined to generate a list of modified muropeptides (library 3). Two libraries of matched theoretical masses (monomers and dimers, trimers) and a third library (their modified counterparts) are used to search the dataset. Muropeptide structures are inferred from a match within tolerance between theoretical and observed masses. This data is then ‘cleaned up’ by combining the intensities of ions corresponding to in-source decay and salt adducts to those of parent ions. The final matched mass spectrometry data is then written to a .csv file.

Using PGFinder to investigate PG structure and identify low-abundance muropeptides

The performance of the matching script was tested using the well-characterised PG from E. coli as a proof of concept. UHPLC-MS/MS data were acquired for three independent PG samples (biological replicates; Appendix 1—figure 1). Following MS1 spectral deconvolution (a process calculating masses from observed m/z values), observed masses were matched to theoretical muropeptides masses according to the strategy described above (Figure 2). A first search was carried out using a minimal mass database made of 10 simple PG fragments including three glycan chains (di-, tetra-, and hexasaccharides) and seven monomers (Table 1—source data 1). Due to the .csv format of the database, the diamino acid in position 3 could not be assigned to a symbol or Greek letter and had to be one of the 26 letters already assigned by the IUPAC-IUB Joint Commission on Biochemical Nomenclature. We used the letter J for mDAP for the initial search and replaced it by the letter m in the final table. The output of the automated search is a .csv file per dataset; all files corresponding to biological replicates were collated into one Excel file (Table 1—source data 2). Each search output contained approximately 3000 rows of masses and corresponding parameters. Depending on the dataset analysed, 41–48% of the total ion intensity was assigned to PG structures. As anticipated, inferred structures were frequently found with multiple retention times, reflecting the existence of stereoisomers, with one species accounting for most of the intensity. In some cases, observed masses matched with more than one inferred structure. The output of the automated search was consolidated as described in Supplementary file 1. Retention times were assigned to individual structures based on the elution of the most abundant stereoisomer. For example, >97% of the most abundant monomer GM-AEJA was eluted at an average retention time of 10.04 ± 0.04 min. Data consolidation revealed an unprecedented muropeptide composition complexity compared to recent LC-MS analyses of E. coli (Kühner et al., 2014; Table 1—source data 2). Sixty PG fragments were identified (Table 1): these included glycan chains lacking peptide stems (4.38%), monomers (63.14%), dimers (29.54%), and trimers (2.94%) (Figure 3). Based on the abundance of multimeric PG fragments, we report a crosslinking index of 15.69%, which is slightly lower than the value previously reported of 23.1% (Glauner, 1988).

Figure 3
Download asset Open asset
Distribution of E. coli peptidoglycan fragments identified using automated search workflow.

Breakdown of peptidoglycan is shown by oligomerisation state (left) branching to specific composition (right). Branch size is proportional to percentage. Monomers, dimers, trimers, and glycan chains (left) are broken down into muropeptide composition and structure (right). Individual structures are grouped by colour according to oligomerisation state. Monomers, green; dimers, yellow; trimers, orange. Residues in square brackets are only found in some muropeptides. For example, GM-AEJ[A]-GM-AEJ[A] can represent GM-AEJA-GM-AEJA, GM-AEJA-GM-AEJ, and GM-AEJ-GM-AEJ. G: N-acetylglucosamine; M, N-acetylmuramic acid; A: L- or D- alanine; E: γ-D-glutamic acid; J: meso-diaminopimelic acid; K: D-lysine; R: D-arginine; G: glycine.

Table 1
Processed match output.
StructureRT (min)Abundance (%)Monoisotopic mass (Da)
Av±SDAv±SDObsTheoΔppm
GM|03.62±0.013.465±0.683498.205498.2062.5
GlycansGM (x2)|010.11±0.030.428±0.349976.384976.3862.2
4.38%±0.35%GM (anhydro) |08.20±1.920.238±0.025478.179478.1802.9
GM (deacetyl) |02.57±0.000.155±0.032456.194456.1963.5
GM (x2) (deacetyl) |06.86±0.020.093±0.012934.372934.3763.2
GM-AEmA|110.04±0.0436.098±2.131941.405941.4082.8
GM-AEm|16.57±0.0114.352±0.397870.368870.3713.0
GM-AEmKR|19.56±0.058.030±0.7741154.5631154.5673.6
GM-AE|19.57±0.041.809±0.231698.284698.2863.1
GM-AEmG|17.85±0.050.689±0.049927.390927.3922.3
GM-AEm (anhydro) |113.98±0.020.668±0.073850.342850.3442.2
MonomersGM-AEmA (anhydro) |116.55±0.010.573±0.100921.380921.3822.0
63.14%±1.13%GM-AEmAG|19.45±0.050.219±0.009998.426998.4293.1
GM-AEmKR (anhydro) |114.83±0.010.160±0.0391134.5371134.5402.9
GM-AEmA (deacetyl) |18.57±0.060.083±0.055899.394899.3973.1
GM-GM-AEmA|113.10±0.020.075±0.0401419.5841419.5882.9
GM-AE (anhydro) |117.44±0.010.069±0.013678.258678.2602.8
M-AEm|14.56±0.010.062±0.064667.289667.2913.8
M-AEmKR|18.16±0.060.061±0.056*951.484951.4873.2
GM-AEmAA|111.38±0.040.059±0.0031012.4421012.4452.4
M-AEmA|18.52±0.050.053±0.015738.325738.3284.0
GM-GM-AEm|111.31±0.040.042±0.0251348.5471348.5512.4
GM-AEm (deacetyl) |14.77±0.010.024±0.014828.358828.3603.0
GM-GM-AEmKR|112.18±0.030.011±0.002*1632.7421632.7473.0
GM-AEmA-GM-AEmA|216.01±0.0217.247±0.7771864.8001864.8052.3
GM-AEmA-GM-AEmKR|214.83±0.024.589±0.5892077.9572077.9643.0
GM-AEmA-GM-AEm|215.09±0.023.207±0.1681793.7631793.7682.6
GM-AEmA-GM-AEmA (anhydro) |220.56±0.010.873±0.0371844.7741844.7782.4
GM-AEm-GM-AEmKR|214.22±0.000.855±0.1012006.9202006.9263.3
GM-AEmA-GM-AEmKR (anhydro) |218.89±0.170.665±0.0792057.9342057.9371.8
GM-AEm-GM-AEm|214.23±0.010.558±0.0621722.7251722.7303.0
GM-AEm-GM-AEmAG|214.68±0.010.416±0.0251850.7851850.7892.4
GM-AEmA-GM-AEm (anhydro) |219.66±0.010.381±0.0281773.7381773.7412.1
DimersGM-AEmA-GM-AEmAG|215.33±0.020.179±0.0051921.8221921.8262.2
29.54%±0.46%GM-AEm-GM-AEmKR (anhydro) |218.07±0.010.170±0.0241986.8961986.9002.1
GM-AEm-GM-AEm (anhydro) |218.77±0.010.141±0.0151702.6971702.7044.5
GM-AEmA-GM-AEmAA|216.54±0.010.075±0.0021935.8381935.8422.1
GM-AEm-GM-AEmG|213.91±0.010.054±0.0031779.7471779.7522.7
GM-GM-AEmA-GM-AEmA|217.51±0.010.046±0.0282342.9762342.9853.6
GM-AEmA-GM-AEmA (deacetyl) |215.17±0.010.029±0.0221822.7891822.7943.0
GM-AEmA-GM-AEmG (anhydro) |219.12±0.010.021±0.0011830.7611830.7630.7
GM-AEmA-GM-AEmAG (anhydro) |219.73±0.010.019±0.0021901.7961901.8002.1
GM-AEmA-GM-AEmAA (anhydro) |221.17±0.020.015±0.0021915.8121915.8161.8
GM-GM-AEmA-GM-AEm|216.85±0.000.003±0.0042271.9432271.9472.1
GM-AEmA-GM-AEmA-GM-AEmA|318.86±0.011.751±0.2212788.1922788.2023.5
GM-AEmA-GM-AEmA-GM-AEm|318.23±0.210.371±0.0312717.1582717.1642.2
GM-AEmA-GM-AEmA-GM-AEmA (anhydro) |322.39±0.020.222±0.0272768.1692768.1752.3
GM-AEmA-GM-AEmA-GM-AEmKR|317.54±0.010.207±0.0283001.3503001.3603.4
GM-AEmA-GM-AEmA-GM-AEm (anhydro) |321.60±0.020.117±0.0032697.1332697.1381.8
TrimersGM-AEmA-GM-AEmA-GM-AEmKR (anhydro) |320.90±0.160.088±0.0262981.3282981.3342.2
2.94%±0.36%GM-AEmA-GM-AEmA-GM-AEmG|317.72±0.010.039±0.0042774.1822774.1861.4
GM-AEmA-GM-AEm-GM-AEm|317.45±0.010.029±0.0052646.1232646.1271.7
GM-AEmA-GM-AEm-GM-AEm (anhydro) |321.16±0.010.025±0.0012626.0962626.1011.9
GM-AEmA-GM-AEm-GM-AEmKR|317.11±0.010.022±0.0022930.3162930.3232.7
GM-AEmA-GM-AEmA-GM-AEmAG|318.24±0.010.021±0.0012845.2172845.2232.0
GM-AEmA-GM-AEmA-GM-AEmAA|319.23±0.010.014±0.002*2859.2352859.2391.3
GM-AEmA-GM-AEm-GM-AEmKR (anhydro) |320.31±0.020.014±0.0052910.2932910.2971.5
GM-AEm-GM-AEmG-GM-AEmAG|317.18±0.000.004±0.0052703.1432703.1492.0
GM-AEmA-GM-AEm-GM-AEmG (anhydro) |321.21±0.020.011±0.003*2754.1572754.1601.1
GM-AEmA-GM-AEmA-GM-AEmAG (anhydro) |321.77±0.010.006±0.0042825.1892825.1972.8

  1. Inferrred dimers and trimers are based on the most abundant monomers and could correspond to alternative structures.

  2. G: GlcNAc; M: MurNAc; m: meso-diaminopimelic acid; the number following the symbol ‘|’ refers to the oligomerisation state (1 for monomers, 2 for dimers, and 3 for trimers).

  3. *

    Calculated from two values.

Table 1—source data 1

E. coli simple mass database.

https://cdn.elifesciences.org/articles/70597/elife-70597-table1-data1-v1.csv
Table 1—source data 2

E. coli matching output and consolidated data.

https://cdn.elifesciences.org/articles/70597/elife-70597-table1-data2-v1.xlsx
Table 1—source data 3

MS/MS analysis of E. coli glycan chains and monomers.

https://cdn.elifesciences.org/articles/70597/elife-70597-table1-data3-v1.pdf
Table 1—source data 4

E. coli complex mass database.

https://cdn.elifesciences.org/articles/70597/elife-70597-table1-data4-v1.csv
Table 1—source data 5

E. coli muropeptide complex table.

https://cdn.elifesciences.org/articles/70597/elife-70597-table1-data5-v1.xlsx

The automated and unbiased search revealed several muropeptides that were expected but never reported to date for E. coli. These included (i) PG fragments resulting from amidase activity (4.55%), found as ‘denuded glycans’ (disaccharides and tetrasaccharides) and modified variants or muropeptide stem with an extra disaccharide residue; (ii) a low-abundance (0.23%) PG fragments containing deacetylated GlcNAc residues; and (iii) PG fragments resulting from glucosaminidase activity (0.12%). Deacetylated muropeptides were not expected since no E. coli PG deacetylase has been identified in this organism to date. All the structures identified for the first time (glycan chains, monomers containing deacetyl groups, and muropeptides lacking a GlcNAc residue) were confirmed by MS/MS analysis (Table 1—source data 3). The proportion of muropeptides with anhydroMurNAc groups identified (4.55%) was in line with previous studies, yielding an average chain length of 36.05. This value is higher than an earlier study that reported a predominant chain length of 5–10 disaccharide units (Harz et al., 1990), but agreed with recent work that reported long glycan chains in E. coli (Turner et al., 2018). Overall, the quantification of muropeptides across biological replicates was very consistent, with Pearson’s correlation coefficients >0.96 (Appendix 1—figure 2). The most pronounced variations in quantification were observed with low-abundance muropeptides accounting for less than 1% of the species identified.

We further explored the structural diversity of E. coli PG, performing a more complex search with a mass database made of glycan chains and all possible monomers containing di-, tri-, tetra-, and pentapeptide stems (Table 1—source data 4; 224 structures in total). Several monomers with tetra- and pentapeptide stems containing unusual amino acids were identified. Only four structures could be confirmed by MS/MS analysis and corresponding multimers were retained (GM-AEJF, -AEJN, AEJK, -AEJAK, and GM-AEJKD). Collectively, muropeptides containing unusual amino acids accounted for ca. 7.5% of the 80 structures identified (Table 1—source data 5) using a complex mass database.

Using PGFinder for the comparative analysis of PG structures

We showed with E. coli data that PGFinder is suitable to characterise the high-resolution structure of PGs using a ‘bottom-up’ approach. However, this requires a careful analysis of the search output to confirm the identity of muropeptides identified and discriminate between multiple structures that can be assigned to a unique observed mass. A more basic application is the use of PGFinder in organisms that have already been studied in detail to either compare PG composition or quantify the abundance of specific structures. This application accounts for most PG analyses described in the literature, comparing PG structures between different isolates, isogenic mutants, or cells grown in different conditions.

We chose the PG of C. difficile as a proof of concept to demonstrate how PGFinder can carry out a straightforward comparative analysis. We prepared PG samples corresponding to biological triplicates from two clinical isolates, R20291 and M7404. To illustrate the versatility of the software, PG samples were digested by mutanolysin, and disaccharide peptides were beta-eliminated to generate lactyl-peptides (Tipper, 2002) and analysed by UHPLC-MS (Figure 4—figure supplement 1). Since the high-resolution PG structure of C. difficile has been described based on the MS/MS analysis of muropeptides (Bern et al., 2017), we used data previously published to build a PG fragment database for a ‘one-off’ matching step. The database included monomers identified by MS/MS, containing unusual amino acids and the corresponding dimers resulting from either D,D- or L,D-transpeptidation (with a AEJA or AEJ peptide as donor stem). To limit the complexity of the search output, we limited the list of trimers, tetramers, and pentamers to those containing the most abundant peptide stems found in dimers (AEJA, AEJ, or AEJG). The database of theoretical masses contained 74 PG structures described in Figure 4—source data 1. MS data were deconvoluted using Byos Feature Finder, and observed monoisotopic masses were matched to theoretical masses using PGFinder. To perform the matching operation in its simplest form, all options offered by the software were deactivated.

All monomers and dimers searched were identified, except for two muropeptides containing a peptide stem with a methionine residue in position four, both present in low abundance in C. difficile strain 630 (ca. 0.10%). 25 out of the 31 possible trimers, tetramers, and pentamers searched were found (Figure 4—source data 2). Comparison of biological replicates revealed a high reproducibility, both in retention times and quantification. A high correlation was found between biological replicates, confirming the robustness of the quantification method (Figure 4a). Both strains contained a similar amount of mono-, tri-, tetra-, and pentamers, but strain M7404 contained a significantly lower proportion of monomers and a higher proportion of dimers than strain R20291 (p=0.022 and p=0.005, respectively; Student’s t-test; Figure 4b). We next performed a Student’s t-test using permutation-based FDR to identify statistically significant differences in the abundance of individual muropeptides between the two strains. The p-value was plotted on a volcano plot against the fold change in abundance between the two samples (Figure 4c). Two muropeptides were significantly less abundant (Lac-AEJ[AG] and Lac-AEJ-Lac-AEJA), and four others were significantly more abundant in strain R20291 (Lac-AEJV- Lac-AEJA, Lac-AEJ[L/I]-Lac-AEJA, Lac-AEJAA-Lac-AEJA and the trimer (Lac-AEJA)3). These differences are likely to reflect different substrate specificities for PBPs and the Ddl ligases in these strains. Therefore, combining the output of PGFinder with statistical analysis of muropeptide abundance offers a robust workflow to identify differences in PG composition.

Figure 4 with 1 supplement see all
Download asset Open asset
Comparative analysis of C. difficile R20291 and M7404 peptidoglycan (PG) composition.

(a) Pearson’s correlation coefficients across biological replicates of R20291 and M7404 C. difficile isolates. Heatmap gradient shows highest value in green to lowest value in red. (b) Muropeptide distribution according to degree of crosslinking. Comparison was carried out using a Student’s t-test; p-value is indicated for each category of muropeptides. (c) Volcano plot, where each dot represents an individual muropeptide, plotted against the significance (Student’s t-test p-value<0.05, FDR < 0.05, S0 = 0.1) and difference (log2). Muropeptides showing a significantly different abundance between strains are highlighted in red. Lac: lactyl group; A: D/L-alanine; E: γ-D-glutamate; J: meso-diaminopimelic acid V: D-valine; L: D-leucine; I: D-isoleucine; G: glycine.

Figure 4—source data 1

C. difficile mass database.

https://cdn.elifesciences.org/articles/70597/elife-70597-fig4-data1-v1.csv
Figure 4—source data 2

C. difficile 20291 versus M7404, list of muropeptides, abundance, RT.

https://cdn.elifesciences.org/articles/70597/elife-70597-fig4-data2-v1.xlsx

Benchmarking the automated PG analysis pipeline using available datasets

The analysis of E. coli PG established a proof of concept, showing that our matching script is suitable for the automated analysis of PG MS data. We next sought to evaluate the robustness of our peptidoglycomics pipeline using available datasets described in the literature. The most suitable publication was a recent study by Anderson et al. describing a PG analysis of P. aeruginosa planktonic cells (Anderson et al., 2020b). Unlike most (if not all) studies published to date, this work provided datasets from biological and technical triplicates. Unlike our E. coli samples, analysed on a Q Exactive Focus Orbitrap (Thermo), P. aeruginosa samples were analysed on an Agilent Q-TOF mass spectrometer. Spectra were deconvoluted using the Byos Feature Finder module, and observed masses were matched using a mass tolerance of 25 ppm as described by Anderson et al. To limit the occurrence of mass coincidences and misidentification at this slightly lower mass accuracy, we carried out a search with PG modifications (anhydroMurNAc residues, deacetylation, lack of peptide stem resulting from amidase activity and amidation) but only using the most frequent combinations of modifications (double anhydroMurNAc and anhydroMurNAc and deacetyl).

PGFinder identified 63 muropeptides out of the 71 reported by Anderson et al., matching our search criteria (Table 2). The eight muropeptides that were not identified were absent from the list of deconvoluted masses, indicating that the problem was not associated with the script, highlighting that the deconvolution step is a source of variability. Interestingly, the observed masses calculated using Byos Feature Finder were closer to the theoretical value (6.5 ppm versus 10.7 ppm on average), reflecting another source of variability associated with data deconvolution. Four muropeptides previously identified containing the AEJAG pentapeptide stem were matched with distinct structures due to a mass coincidence between K and AG. A careful analysis of MS/MS spectra suggested that these muropeptides contained a K residue rather than the AG dipeptide (the y2 ion being 143 ppm away from the expected mass), showing the added value of an unbiased search. This conclusion is supported by the retention times of the corresponding muropeptides since the tetrapeptide AEJK elutes before EAJA whilst the pentapeptide AEJAG elutes later in the chromatography (Bern et al., 2017). It is worth noting that our search also identified a large number of muropeptides that were not reported previously (Table 2—source data 1). These results collectively show that our PG analysis pipeline can identify an unprecedentedly large number of muropeptide structures based on MS1 data, including all those previously reported (Anderson et al., 2020a; Anderson et al., 2020b).

Table 2
Automated identification of P. aeruginosa peptidoglycan fragments.
Inferred structureMass∆ppmMaxQuant
TheoreticalObservedThis workAnderson et al.
GM (anhydro)478.1799478.17804.0–2.7+
GM498.2061498.20423.9–4.2+
GM (x2) (deacetyl)934.3755934.37065.3–8.6+
GM (x2) (anhydro)956.3598956.35515.06.0+
GM (x2)976.3860976.37946.7–6.1+
GM (x3) (deacetyl)1412.55541412.54904.5–6.2+
GM (x3) (anhydro)1434.53971434.53483.4–7.5+
GM (x3)1454.56591454.55924.6–5.3+
GM (x4)1932.74581932.73525.5–5.1+
GM-AE (anhydro)678.2596678.25674.3–9.1+
GM-AE698.2858698.28303.9–12.9+
GM-AEJ (anhydro)850.3444850.34015.1–10.6+
GM-AEJ870.3706870.36763.5–5.9+
GM-AEJA (anhydro)921.3815921.37655.4–9.9+
GM-AEJG927.3920927.38685.6–8.9+
GM-AEJA941.4077941.40453.4–5.0+
GM-AEJC973.3843973.37638.2–2072.2+
GM-AEJL983.4593983.44989.6–15.5+
GM-AEJK998.4703998.46248.0–10.6+
GM-AEJM1001.41531001.40609.2–13.5+
GM-AEJAA1012.44481012.44133.4–7.8+
GM-AEJY (anhydro)1013.40911013.4242–14.917.8+
GM-AEJF1017.44331017.43478.4–15.0+
GM-AEJY1033.43531033.42787.2–5.3+
GM-AEJAV1040.48081040.47168.8–14.7+
GM-AEJIA1054.49641054.48748.5–11.3+
GM-AEJW1056.43941056.4455–5.84.0+
GM-AEJAM1072.45241072.44605.9–4.3+
GM-AEJKR1154.56671154.56313.1–8.1+
GM-GM-AE1176.48361176.459020.9–24.7+
GM-GM-AEJ1348.56841348.545716.9–24.9+
GM-GM-AEJA1419.60551419.582416.2–23.5+
GM-AEJA-GM-AEJ (amidase product)1313.57211313.56743.5–11.0+
GM-AEJA-GM-AEJA (amidase product)1384.60921384.60374.0–7.4+
GM-AEJ-GM-AEJ (anhydro)1702.70421702.69763.938.3+
GM-AEJ-GM-AEJ1722.73041722.72344.1–8.6+
GM-AEJA-GM-AEJ (double anhydro)1753.71511753.70436.2–7.2+
GM-AEJA-GM-AEJ (anhydro)1773.74131773.73394.2–11.1+
GM-AEJA-GM-AEJ1793.76751793.75964.4–8.8+
GM-AEJA-GM-AEJA (dacetyl)1822.79411822.78087.3–7.4+
GM-AEJA-GM-AEJA (double anhydro)1824.76011824.74478.4–15.6+
GM-AEJA-GM-AEJA (anhydro)1844.77841844.77044.3–8.3+
GM-AEJA-GM-AEJG1850.78891850.8158–14.69.7+
GM-AEJA-GM-AEJA1864.80461864.79624.5–6.6+
GM-AEJA-GM-AEJK (anhydro)1901.84101901.82975.9–14.5+
GM-AEJA-GM-AEJL1906.85621906.84525.8–11.3+
GM-AEJA-GM-AEJK1921.86721921.85864.5–12.0+
GM-AEJA-GM-AEJF1940.84021940.82637.2–8.8+
GM-AEJA-GM-AEJY1956.83221956.82105.7–7.6+
GM-AEJA-GM-AEJAL1977.89331977.88136.0–10.7+
GM-AEJA-GM-AEJKR2077.96362077.95892.2–13.0+
GM-GM-AEJ-GM-AEJ2200.92822200.900012.8–17.7+
GM-GM-AEJA-GM-AEJ2271.96532271.936812.6–18.4+
GM-GM-AEJA-GM-AEJA2343.00242342.973412.4411.4+
GM-AEJA-GM-AEJA-GM-AEJ (double anhydro)2677.11202677.10004.5–10.7+
GM-AEJA-GM-AEJA-GM-AEJ (anhydro)2697.13822697.12594.6–8.6+
GM-AEJA-GM-AEJA-GM-AEJ2717.16442717.15324.1–10.7+
GM-AEJA-GM-AEJA-GM-AEJA (double anhydro)2748.14912748.13634.7–11.0+
GM-AEJA-GM-AEJA-GM-AEJA (anhydro)2768.17532768.16742.9–11.2+
GM-AEJA-GM-AEJA-GM-AEJA2788.20152788.19193.4–9.7+
GM-AEJA-GM-AEJA-GM-AEJK (anhydro)2825.23792825.22056.1–9.3+
GM-GM-AEJA-GM-AEJA-GM-AEJ3195.36223195.326411.2–14.0+
GM-GM-AEJA-GM-AEJA-GM-AEJA3266.39933266.363011.1–12.5+

  1. Alternative structures were matched:

  2. GM-AEJ-GM-AEJK.

  3. GM-AEJ-GM-AEJKA (anhydro).

  4. GM-AEJ-GM-AEJKA.

  5. GM-AEJ-GM-AEJA-GM-AEJKA (anhydro).

Table 2—source data 1

Pseudomonas aeruginosa matched muropeptides not reported previously.

https://cdn.elifesciences.org/articles/70597/elife-70597-table2-data1-v1.xlsx
Table 2—source data 2

Raw output of automated search using MaxQuant and PGFinder.

https://cdn.elifesciences.org/articles/70597/elife-70597-table2-data2-v1.xlsx

A PG analysis workflow using freely available tools

The automated identification of P. aeruginosa muropeptides using PGFinder and the strategy reported previously (Anderson et al., 2020a; Anderson et al., 2020b) relied on commercially available deconvolution software (ProteinMetrics Byos or Agilent MassHunter, respectively). We sought to identify a free alternative software for mass deconvolution to make our PG analysis pipeline accessible to everyone. MaxQuant was selected as a tool of choice since it represents a widely used software package to analyse high-resolution mass-spectrometric data for shotgun proteomics (Cox and Mann, 2008). MaxQuant has native support for Thermo MS data (RAW) and Sciex (WIFF) file formats and supports the open MS data format mzXML. Virtually any proprietary MS data file can be converted to the mzXML open format using freely available tools such as Proteowizard and TOPPAS, making this workflow universally applicable. As proof of concept, we converted P. aeruginosa data to an mzXML file (Appendix 2—figure 1), processed it using MaxQuant for mass deconvolution (Appendix 2—figure 2), and analysed it using PGFinder for PG structure and composition identification. We were able to identify all the expected muropeptides (Table 2). This result confirms that the automated analysis of PG datasets can be carried out using the MaxQuant freeware and our open-source script PGFinder (Table 2 and Table 2—source data 2).

Discussion

This study describes a workflow for the unbiased and automated analysis of bacterial PG using freely available resources. We analysed high-resolution MS datasets corresponding to PG fragments and demonstrated that this approach is a powerful tool to identify muropeptides and carry out comparative analyses based on the MS1 data.

MS analysis of bacterial PG has been carried out since the late 1980s (Garcia-Bustos et al., 1988; Martin et al., 1987). Whilst rp-HPLC-MS is now routinely used to explore PG structure, data analysis remains a manual process. Therefore, this step has become a bottleneck that prevents high-throughput analyses and introduces a series of issues regarding reproducibility. A major issue deals with the definition of the search space used since the list of muropeptides searched is not provided. Our previous work showed that an unbiased approach using shotgun proteomics tools identifies muropeptides containing unusual amino acids in C. difficile (Bern et al., 2017). Our unbiased search identified >106 masses matching E. coli muropeptide structures, representing a number strikingly larger than previously reported (Kühner et al., 2014). This level of complexity was anticipated based on the complement of enzymes involved in E. coli PG synthesis but had never been reported before. Our work therefore highlighted the limitations of search strategies reported so far, especially if we consider the fact that some bacterial PGs like E. coli have been extensively studied over the past 30 years.

Several issues and flaws associated with the manual analysis of PG MS data are addressed by our approach, which represents a robust, consistent, and open-access strategy. The PGFinder algorithm has been designed to create dynamic databases that are ultimately combined to perform the final matching process. Optimisation of the search space relies on a preliminary identification of masses matching theoretical monoisotopic masses of monomers, limiting misidentifications based on mass coincidence. Another advantage is that the only information provided by the user is a restricted monomer database rather than a comprehensive one. This avoids a time-consuming operation, prone to human error. PGFinder uses XIC for quantification of PG fragments, providing high resolution, sensitivity, and reproducibility, as indicated by the comparisons across biological replicates (Appendix 1—figure 2 and Figure 4a). It enables the accurate quantification of molecules with overlapping retention times and those present in very low abundance, with a large dynamic range (typically six orders of magnitude). Although a Matlab-based software package (Chromanalysis) has been described to automate the detection and quantification of UV peaks through Gaussian fitting (Desmarais et al., 2015), quantification using XIC is more straightforward. It is worth pointing out that the approach described here does not give absolute quantification of muropeptides. Instead, it allows a relative quantification of muropeptides. Nevertheless, this strategy remains suitable for comparative analyses and overcomes two major limitations associated with UV detection, namely detection threshold and co-elution of molecules.

Our E. coli PG analysis confirmed that PGFinder is a powerful tool that provides a much improved qualitative and quantitative PG analysis (Kühner et al., 2014; Morè et al., 2019). Combining an unbiased search with highly sensitive detection of individual structures is important for two reasons. Firstly, it opens the possibility to identify subtle modifications of PG structure, resulting from either a transient or a localised enzymatic activity such as that taking place at the septum. Secondly, it will permit the identification of previously undetected modifications that may provide new insights into our understanding of PG composition and dynamics. For example, we showed that E. coli PG contains a low abundance of deacetylated sugars. This observation is puzzling because no canonical PG deacetylase genes have been identified in this organism. Although the biological relevance of this property remains to be established, we cannot exclude the possibility that PG deacetylation in E. coli may contribute to PG homeostasis. Another striking outcome resulting from our automated search is the identification of a slightly higher amount of muropeptides containing anhydromuramic acid (4.55%) as compared to 2–3% (Glauner et al., 1988; Liu et al., 2020). To explain the discrepancy between our work and data from the literature, it is tempting to assume that most of the muropeptides containing anhydromuramic acid identified with PGFinder were simply not searched in previous studies. It is worth pointing out that none of the papers describing PG analysis published to date has reported the list of structures searched in the MS data analysed.

One of our objectives was to create an automated PG analysis tool accessible to the broadest audience possible, including people with no prior experience with programming or coding languages. Therefore, we shared PGFinder as a Jupyter Notebook allowing users to customise the search strategy depending on both the question asked and the instrument accuracy. PGFinder is particularly suitable for the characterisation of novel PGs with unknown composition or structural modifications and can be modified by users to add novel functionalities. However, a current limitation of this workflow is that it does not process MS/MS data. Therefore, the fragmentation spectra of individual monomers must be checked using dedicated tools to validate that the inferred structures are correct. We are currently working towards an integrated pipeline that includes MS/MS analysis to our PGFinder pipeline. The ability to disable some PG modifications means that the complexity of the search can be adjusted to focus on specific properties (e.g., the occurrence of acetylation/deacetylation, or amidation) or specific muropeptides resulting from lytic activities (e.g., unsubstituted MurNAc residues resulting from amidase activity). For PG that have already been well characterised (E. coli, P. aeruginosa, or C. difficile), the search parameters are already established, allowing a very straightforward analysis to be performed. Therefore, access to a custom, semi-quantitative sensitive analysis is ideal for comparing PG dynamics or differences in PG structure between a reference strain and isogenic mutants. Both reduced disaccharide peptides or lactyl-peptides (generated by beta-elimination) are identified using PGFinder.

We anticipate that an open access to PGFinder, in conjunction with freely available deconvolution tools, will allow researchers to carry out comparative MS1 analyses. The pipeline defined in this work enables reproducible and consistent data analysis. This represents the first step towards a standardised approach to PG analysis, opening the possibility to reanalyse datasets in repositories. The modular structure of the open-source PGFinder code can be easily integrated into any specific workflow for the automated processing of PG MS data.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background(Escherichia coli)BW25113https://doi.org/10.1073/pnas.120163297RRID:Addgene_72340Model strain for PG analysis
Strain, strain background(Clostridioides difficile)R20291https://doi.org/10.1128/JB.0073107Model strain for PG analysis
Strain, strain background(Clostridioides difficile)M7404https://doi.org/10.1371/journal.ppat.1002317Model strain for PG analysis
Software, algorithmPGFinderv.0.02This workUsed for MS1 analysis of PG structure
Software, algorithmByosv.3.9–32Protein Metrics IncUsed for MS data deconvolution and MS/MS analysis
Software, algorithmMaxQuant v2.0.1.0Cox and Mann, 2008RRID:SCR_014485Used for MS data deconvolution
Software, algorithmPerseusv.1.6.10.53Tyanova et al., 2016RRID:SCR_015753Used statistical analysis of muropeptide abundance

Bacterial strains and culture conditions

Request a detailed protocol

E. coli BW25113 was grown at 37°C in LB under agitation (250 rpm). C. difficile strains were cultured in heart infusion supplemented with yeast extract, L-cysteine, and glucose in an atmosphere of 10% H2, 10% CO2, and 80% N2 at 37°C in a Coy chamber or Don Whitley A300 anaerobic workstation.

PG purification

Request a detailed protocol

PG was purified from exponential (E. coli) or late exponential (C. difficile) phase as described previously (Eckert et al., 2006; Glauner, 1988), freeze-dried, and resuspended in distilled water at a concentration of 5 mg/ml.

Preparation of soluble muropeptides

Request a detailed protocol

PG (1 mg) was digested overnight with 25 µg of mutanolysin at 37°C in 150 µl of 20 mM sodium phosphate buffer (pH 5.5). Soluble disaccharide peptides were recovered in the supernatant following centrifugation (20,000 × g for 20 min at 25°C). To reduce muropeptides, equal volumes (200 µl) of the solution of disaccharide peptides and of borate buffer (250 mM, pH 9.0) were mixed. 2 ml of sodium borohydride was added, and the solution was incubated for 20 min at room temperature. The pH of the solution was adjusted to 4.0 with 20% orthophosphoric acid. Beta-elimination was carried out by mixing 200 µl of muropeptides with 64 µl of 32% (w/v) ammonia. After 5 hr at 37 C, the solution was neutralised with 60 µl of acetic acid glacial, freeze-dried, and resuspended in water (Arbeloa et al., 2004; Eckert et al., 2006).

The reduced muropeptides were desalted by reverse-phase high-performance liquid chromatography (rp-HPLC) on a C18 Hypersil Gold aQ column (3 µm, 2.1 × 200 mm; Thermo Fisher) at a flow rate of 0.4 ml/min. After 1 min in water-0.1% formic acid (v/v) (buffer A), muropeptides were eluted with a 6 min linear gradient to 95% acetonitrile-0.1% formic acid (v/v). Muropeptides were freeze-dried and resuspended in 100 µl. An aliquot of the desalted samples was analysed by rp-HPLC on the same column to measure the UV absorbance of the most abundant monomer (no isocratic step, muropeptides were eluted with a 30 min linear gradient to 15% acetonitrile-0.1% formic acid [v/v]). Samples were diluted to contain 150 mAU/µl of the major monomer and 10 µl were injected. Based on the dry weight of the PG sample, we estimated that this corresponded to approximately 50 µg of material.

UHPLC-MS/MS

Request a detailed protocol

An Ultimate 3000 Ultra High-Performance Chromatography (UHPLC; Dionex/Thermo Fisher Scientific) system coupled with a high-resolution Q Exactive Focus mass spectrometer (Thermo Fisher Scientific) was used for LC/HRMS analysis. Muropeptides were separated using a C18 analytical column (Hypersil Gold aQ, 1.9 µm particles, 150 × 2.1 mm; Thermo Fisher Scientific), column temperature at 50°C. Muropeptides elution was performed by applying a mixture of solvent A (water, 0.1% [v/v] formic acid) and solvent B (acetonitrile, 0.1% [v/v] formic acid). After 10 µl sample injection, MS/MS data were acquired during a 40 min step gradient: 0–12.5% B for 25 min; 12.5–20% B for 5 min; held at 20% B for 5 min, and the column was re-equilibrated for 10 min under the initial conditions.

The Q Exactive Focus was operated under electrospray ionization (H-ESI II)-positive mode. Full scan (m/z 150–2250) used resolution 70,000 (FWHM) at m/z 200, with an automatic gain control (AGC) target of 1 × 106 ions and an automated maximum ion injection time (IT).

Data-dependent MS/MS were acquired on a ‘Top 3’ data-dependent mode using the following parameters: resolution 17,500; AGC 1 × 105 ions, maximum IT 50 ms, NCE 25%, and a dynamic exclusion time 5 s.

MS data deconvolution

Byos search parameters to get .ftrs file Protein Metrics Byos v.3.9–32 was used to identify and compute the XICs. The parameters used for mass deconvolution using MaxQuant v.2.0.1.0 are described in Appendix 2—figure 2.

Data analysis

Request a detailed protocol

The crosslinking index and glycan chain length were calculated as described previously (Glauner, 1988). Label-free relative quantitation of muropeptides from triplicate C. difficile clinical isolates (R20291 and M7404) was performed using Byos 3.11, and statistical analysis of the quantitative data was performed using Perseus v. 1.6.10.53 (Tyanova et al., 2016). Briefly, muropeptide intensities were log2 transformed and normalised by subtraction of the median value. A two-sample Student’s t-test was performed with a permutation-based FDR of 0.05 to determine statistically significant quantitative differences between the strains. Comparisons between R20291 and M7404 muropeptide distribution (mono-, di-, tri-, tetra-pentamer) was evaluated for statistical significance using GraphPad Prism (unpaired t-test).

Runtime environment

Request a detailed protocol

Code is available at https://github.com/Mesnage-Org/PGFinder (Patel, 2021). https://github.com/Mesnage-Org/PGFinder/releases/tag/v0.02 will take you to the archived release used in this paper. We used Python 3 to write the MS1 package and demonstrate its functionality using demo scripts. PGFinder can be run through an interactive Jupyter Notebook hosted on mybinder for ease of use by those less familiar with Python code. A conda environment is provided to ensure reproducible execution. Regression testing has been implemented to ensure changes to code do not cause changes to important results. The GitHub contains an interactive version to run user’s analysis and an end-to-end demo using samples data provided with the script (Interactive PGFinder). The sample data is a MaxQuant deconvolution output from the E. coli MS data analysed in the paper. The current version of the script can handle both .txt (MaxQuant) or .ftrs (Byos) deconvoluted data and offers the possibility for the user to include several modifications in the search. The time window for the ‘clean up step’ (in-source decay and salt adducts) as well as ppm tolerance for matching can also be defined by the user; the default values corresponding to these parameters used in this work are 0.5 min and 10 ppm.

Data availability

Request a detailed protocol

All E. coli and C. difficile MS datasets generated in this study are available through the GlycoPOST repository (GPST000168; Watanabe et al., 2021). P. aeruginosa MS datasets are accessible via Figshare (Anderson et al., 2020b).

Appendix 1

Appendix 1—figure 1
Download asset Open asset
UHPLC-MS chromatogram of E. coli reduced disaccharide peptides.
Appendix 1—figure 2
Download asset Open asset
Consistency of E. coli PG analyses.

(a) Pearson’s correlation coefficients across biological replicates of E. coli BW25113. (b) Muropeptide distribution according to degree of crosslinking. The crosslinking index was calculated as described previously (Glauner, 1988). (c) Pairwise comparisons of intensities corresponding to individual muropeptides identified in biological replicates. WT1, WT2 and WT3 correspond to individual biological replicates; Av., average abundance; SD, standard deviation.

Appendix 2

Appendix 2—figure 1
Download asset Open asset
Workflow for production of MaxQuant compatible MS data files from Agilent QTOF data.

Agilent MS data (data: .d) is converted by Proteowizard to a mzML format (data: XML). Relevant settings for Proteowizard are shown (left). mzML file is then converted by TOPPAS to a mzXML file (data: XML). Relevant settings are shown (right).

Appendix 2—figure 2
Download asset Open asset
Workflow for MS data processing using MaxQuant, before automated analysis.

mzXML (data: XML) is passed to MaxQuant (process) for deconvolution and monoisotopic mass determination. Default values used except where indicated (right). MaxQuant output (data: text file) is then passed to the data parser module (process). This module removes superfluous data and reformats remaining data to be compatible with the matching script as an Excel file (data: xlsx).

Data availability

All raw mass spectrometry data files have are available through the Glycopost repository (ref GPST000168).

The following previously published data sets were used
    1. Anderson EM
    2. Sychantha D
    3. Brewer D
    4. Clarke AJ
    5. Geddes-McAlister J
    6. Khursigara CM
    (2020) figshare
    Peptidoglycomics: Examining compositional changes in peptidoglycan between biofilm- and planktonic-derived Pseudomonas aeruginosa.
    https://doi.org/10.6084/m9.figshare.10277909

References

    1. Glauner B
    2. Höltje JV
    3. Schwarz U
    (1988)
    The composition of the murein of Escherichia coli
    The Journal of Biological Chemistry 263:10088–10095.
    1. Martin SA
    2. Rosenthal RS
    3. Biemann K
    (1987)
    Fast atom bombardment mass spectrometry and tandem mass spectrometry of biologically active peptidoglycan monomers from Neisseria gonorrhoeae
    The Journal of Biological Chemistry 262:7514–7522.

Article and author information

Author details

  1. Ankur V Patel

    School of Biosciences, University of Sheffield, Sheffield, United Kingdom
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Validation, Visualization, Writing – original draft, Writing – review and editing
    For correspondence
    apatel19@sheffield.ac.uk
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8161-3455

  2. Robert D Turner

    Department of Computer Science, University of Sheffield, Sheffield, United Kingdom
    Contribution
    Conceptualization, Data curation, Methodology, Software, Validation, Visualization, Writing – review and editing
    Competing interests
    none

  3. Aline Rifflet

    1. Institut Pasteur, Unité Biologie et Génétique de la Paroi Bactérienne, Paris, France
    2. INSERM, Équipe Avenir, Paris, France
    3. CNRS, UMR 2001 "Microbiologie intégrative et moléculaire", Paris, France
    Contribution
    Investigation, Methodology, Writing – review and editing
    Competing interests
    none

  4. Adelina E Acosta-Martin

    biOMICS Facility, Faculty of Science Mass Spectrometry Centre, University of Sheffield, Sheffield, United Kingdom
    Contribution
    Data curation, Formal analysis, Investigation, Methodology, Resources, Writing – review and editing
    Competing interests
    None

  5. Andrew Nichols

    Protein Metrics Inc, Cupertino, United States
    Contribution
    Methodology, Resources, Funding acquisition
    Competing interests
    Andrew Nichols is affiliated with Protein Metrics Inc. The author has no other competing interests to declare

  6. Milena M Awad

    Infection and Immunity Program, Monash Biomedicine Discovery Institute, Clayton, Australia
    Contribution
    Project administration, Investigation, Funding acquisition, Visualization, Writing – review and editing
    Competing interests
    none

  7. Dena Lyras

    1. Infection and Immunity Program, Monash Biomedicine Discovery Institute, Clayton, Australia
    2. Department of Microbiology, Monash University, Clayton, Australia
    Contribution
    Supervision, Funding acquisition, Visualization, Writing – review and editing
    Competing interests
    none

  8. Ivo Gomperts Boneca

    1. Institut Pasteur, Unité Biologie et Génétique de la Paroi Bactérienne, Paris, France
    2. INSERM, Équipe Avenir, Paris, France
    3. CNRS, UMR 2001 "Microbiologie intégrative et moléculaire", Paris, France
    Contribution
    Project administration, Methodology, Resources, Writing – review and editing
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8122-509X

  9. Marshall Bern

    Protein Metrics Inc, Cupertino, United States
    Contribution
    Conceptualization, Project administration, Methodology, Supervision, Software, Funding acquisition, Validation, Writing – review and editing
    Competing interests
    Marshall Bern is affiliated with Protein Metrics Inc. The author has no other competing interests to declare

  10. Mark O Collins

    1. School of Biosciences, University of Sheffield, Sheffield, United Kingdom
    2. biOMICS Facility, Faculty of Science Mass Spectrometry Centre, University of Sheffield, Sheffield, United Kingdom
    Contribution
    Formal analysis, Project administration, Methodology, Supervision, Resources, Funding acquisition, Validation, Visualization, Writing – review and editing
    For correspondence
    mark.collins@sheffield.ac.uk
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7656-4975

  11. Stéphane Mesnage

    School of Biosciences, University of Sheffield, Sheffield, United Kingdom
    Contribution
    Conceptualization, Data curation, Formal analysis, Project administration, Investigation, Methodology, Supervision, Resources, Funding acquisition, Validation, Visualization, Writing – original draft, Writing – review and editing
    For correspondence
    s.mesnage@sheffield.ac.uk
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1648-4890

Funding

Biotechnology and Biological Sciences Research Council (BB/M011151/1)

  • Ankur V Patel
  • Stephane Mesnage

Medical Research Council (MR/S009272/1)

  • Stéphane Mesnage

Agence Nationale de la Recherche (ANR-10-LABX-62-IBEID)

  • Ivo Gomperts Boneca
  • Aline Rifflet

Agence Nationale de la Recherche (ANR-16-IFEC-0004)

  • Ivo Gomperts Boneca
  • Aline Rifflet

Agence Nationale de la Recherche (ANR-18-CE15-0018)

  • Ivo Gomperts Boneca
  • Aline Rifflet

Australian Research Council (DP210103374)

  • Dena Lyras
  • Milena M Awad

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dominique Mengin-Lecreulx (University Paris XI) for E. coli strain BW25113. Yong Kil and Eric Carslon (Protein Metrics) are acknowledged for their constant support. Motoshi Suzuki (NIH/NIAID) is acknowledged for insightful discussions.

Version history

  1. Received: May 22, 2021
  2. Preprint posted: June 1, 2021 (view preprint)
  3. Accepted: August 8, 2021
  4. Version of Record published: September 28, 2021 (version 1)

Copyright

© 2021, Patel et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,999
    views
  • 226
    downloads
  • 6
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ankur V Patel
  2. Robert D Turner
  3. Aline Rifflet
  4. Adelina E Acosta-Martin
  5. Andrew Nichols
  6. Milena M Awad
  7. Dena Lyras
  8. Ivo Gomperts Boneca
  9. Marshall Bern
  10. Mark O Collins
  11. Stéphane Mesnage
(2021)
PGFinder, a novel analysis pipeline for the consistent, reproducible, and high-resolution structural analysis of bacterial peptidoglycans
eLife 10:e70597.
https://doi.org/10.7554/eLife.70597

Share this article

https://doi.org/10.7554/eLife.70597

Categories and tags

Research organisms

Further reading

    1. Microbiology and Infectious Disease

    Staphylococcus aureus FtsZ and PBP4 bind to the conformationally dynamic N-terminal domain of GpsB

    Michael D Sacco, Lauren R Hammond ... Yu Chen
    Research Article

    In the Firmicutes phylum, GpsB is a membrane associated protein that coordinates peptidoglycan synthesis with cell growth and division. Although GpsB has been studied in several bacteria, the structure, function, and interactome of Staphylococcus aureus GpsB is largely uncharacterized. To address this knowledge gap, we solved the crystal structure of the N-terminal domain of S. aureus GpsB, which adopts an atypical, asymmetric dimer, and demonstrates major conformational flexibility that can be mapped to a hinge region formed by a three-residue insertion exclusive to Staphylococci. When this three-residue insertion is excised, its thermal stability increases, and the mutant no longer produces a previously reported lethal phenotype when overexpressed in Bacillus subtilis. In S. aureus, we show that these hinge mutants are less functional and speculate that the conformational flexibility imparted by the hinge region may serve as a dynamic switch to finetune the function of the GpsB complex and/or to promote interaction with its various partners. Furthermore, we provide the first biochemical, biophysical, and crystallographic evidence that the N-terminal domain of GpsB binds not only PBP4, but also FtsZ, through a conserved recognition motif located on their C-termini, thus coupling peptidoglycan synthesis to cell division. Taken together, the unique structure of S. aureus GpsB and its direct interaction with FtsZ/PBP4 provide deeper insight into the central role of GpsB in S. aureus cell division.

    1. Microbiology and Infectious Disease

    A modified BCG with depletion of enzymes associated with peptidoglycan amidation induces enhanced protection against tuberculosis in mice

    Moagi Tube Shaku, Peter K Um ... Bavesh D Kana
    Research Article

    Mechanisms by which Mycobacterium tuberculosis (Mtb) evades pathogen recognition receptor activation during infection may offer insights for the development of improved tuberculosis (TB) vaccines. Whilst Mtb elicits NOD-2 activation through host recognition of its peptidoglycan-derived muramyl dipeptide (MDP), it masks the endogenous NOD-1 ligand through amidation of glutamate at the second position in peptidoglycan side-chains. As the current BCG vaccine is derived from pathogenic mycobacteria, a similar situation prevails. To alleviate this masking ability and to potentially improve efficacy of the BCG vaccine, we used CRISPRi to inhibit expression of the essential enzyme pair, MurT-GatD, implicated in amidation of peptidoglycan side-chains. We demonstrate that depletion of these enzymes results in reduced growth, cell wall defects, increased susceptibility to antibiotics, altered spatial localization of new peptidoglycan and increased NOD-1 expression in macrophages. In cell culture experiments, training of a human monocyte cell line with this recombinant BCG yielded improved control of Mtb growth. In the murine model of TB infection, we demonstrate that depletion of MurT-GatD in BCG, which is expected to unmask the D-glutamate diaminopimelate (iE-DAP) NOD-1 ligand, yields superior prevention of TB disease compared to the standard BCG vaccine. In vitro and in vivo experiments in this study demonstrate the feasibility of gene regulation platforms such as CRISPRi to alter antigen presentation in BCG in a bespoke manner that tunes immunity towards more effective protection against TB disease.

    1. Microbiology and Infectious Disease

    Tools and methods for high-throughput single-cell imaging with the mother machine

    Ryan Thiermann, Michael Sandler ... Suckjoon Jun
    Tools and Resources

    Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning-based segmentation, ‘what you put is what you get’ (WYPIWYG) – that is, pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother machine-based high-throughput imaging and analysis methods in their research.