Integrative proteomic and metabonomic profiling elucidates amino acid and lipid metabolism disorder in CA-MRSA-infected breast abscesses

Abstract

Objective:

Bacterial culture and drug sensitivity testing have been the gold standard for confirming community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) infection in breast abscess with a long history. However, these tests may delay treatment and increase the risk of nosocomial infections. To handle and improve this critical situation, this study aimed to explore biomarkers that could facilitate the rapid diagnosis of CA-MRSA infection.

Methods:

This study for the first time applied label-free quantitative proteomics and non-targeted metabonomics to identify potential differentially expressed proteins (DEPs) and differentially expressed metabolites (DEMs) in breast abscess infected with CA-MRSA compared to methicillin-susceptible S. aureus (MSSA). The two omics data were integrated and analyzed using bioinformatics, and the results were validated using Parallel Reaction Monitoring (PRM). Receiver operating characteristic (ROC) curves were generated to evaluate the predictive efficiency of the identified biomarkers for diagnosing CA-MRSA infection.

Results:

After using the above-mentioned strategies, 109 DEPs were identified, out of which 86 were upregulated and 23 were downregulated. Additionally, a total of 61 and 26 DEMs were initially screened in the positive and negative ion modes, respectively. A conjoint analysis indicated that the amino acid metabolism, glycosphingolipid biosynthesis, and glycerophospholipid metabolism pathways were co-enriched by the upstream DEPs and downstream DEMs, which may be involved in structuring the related network of CA-MRSA infection. Furthermore, three significant DEMs, namely, indole-3-acetic acid, L-(−)-methionine, and D-sedoheptulose 7-phosphate, displayed good discriminative abilities in early identification of CA-MRSA infection in ROC analysis.

Conclusion:

As there is limited high-quality evidence and multiple omics research in this field, the explored candidate biomarkers and pathways may provide new insights into the early diagnosis and drug resistance mechanisms of CA-MRSA infection in Chinese women.

1 Introduction

Lactational breast abscesses are crucial risk factors for early breastfeeding cessation in postpartum women and poor infant weight gain (Pevzner and Dahan, 2020). Bacterial infection is the main cause of lactational breast abscesses. Through the ruptured nipple, bacteria retrograde to the gland. Coupled with milk stagnation as a good culture medium, bacteria accelerate the reproduction, resulting in mastitis, and inflammation spreads further to form abscesses (Yin et al., 2020). Staphylococcus aureus is the most common pathogenic bacteria causing breast abscesses. Over the past decade, community-acquired methicillin-resistant S. aureus (CA-MRSA), also known as a superbug, has been on the rise, causing breast abscesses and soft tissue infections, accounting for 19.1%–60.0% of S. aureus infections (Boccaccio et al., 2014). CA-MRSA is characterized by high-level drug resistance and complex mechanisms, which can increase infection-caused mortality and prolong hospitalization time (Li Y. et al., 2021). Clinicians may need to intervene more actively for CA-MRSA than methicillin-susceptible S. aureus (MSSA) infection. So far, bacterial culture and drug sensitivity experiments have been the gold standard for diagnosing CA-MRSA infection. However, these methods are time consuming, and the results are easily affected by reagents, experimental operations, sample contamination, and other factors. The application of effective antibiotics is often delayed, which brings severe challenges to clinical treatment and nosocomial infection. Rapid diagnosis and correct management of CA-MRSA infection are essential to avoid complications (Krogerus et al., 2019). Thus, the exploration of more reliable biomarkers for early identification of CA-MRSA-infected breast abscesses is meaningful and urgent.

With the improvement of evidence-based medicine and big multi-omics data, the construction of an accurate inspection system for multiple molecular biomarkers for the early diagnosis and resistance mechanism of CA-MRSA infection would be a strategic breakthrough (Karczewski and Snyder, 2018). There is a general consensus in recent literature that quantitative proteomics could be widely used to identify expression profiles and signal pathways, screen biomarkers, monitor the changes in response to clinical therapeutics, and study protein interaction networks (Cao et al., 2021; Zhang et al., 2021; Dong et al., 2022). Parallel Reaction Monitoring (PRM), based on high-throughput and precision mass spectrometry, has been broadly used for subsequent quantitative verification of target proteins, replacing traditional Western blot technology. This makes the protein validation no longer subject to commercial antibodies, instead of being carried out on large-scale biological samples. Many researchers (Ni et al., 2021; You et al., 2022) combine label-free with PRM to take advantage of wide coverage, converting relative to absolute quantification, to achieve the selection of disease markers and the establishment of diagnostic models. In the process of genetic information transmission from upstream to downstream, metabolites can regulate protein interaction, modify enzyme activity, and stability. Integrating the analysis of proteomics and metabolomics to verify and complement each other, screening differential expressed proteins (DEPs) and differential expressed metabolites (DEMs) with the same change trend and co-participated pathways, are feasible research strategies to fix on potential markers and describe the regulation network of biological organisms (Chen F. et al., 2022; Schumacher-Schuh et al., 2022).

In this study, through the comprehensive proteomic and metabolomic profiling approach, we characterized the potential DEPs and DEMs between CA-MRSA and MSSA and probed into the inflammatory microenvironment features of CA-MRSA-infected breast abscess. The results will provide candidate biomarkers for CA-MRSA early diagnosis and provide clues for multidrug resistance mechanism and the identification of new therapeutic targets.

2 Methods

2.1 Clinical samples and workflow

This study enrolled 60 patients with clinical diagnosis of breast abscesses and confirmed CA-MRSA or MSSA infection by bacterial culture. The case group comprised 24 patients with CA-MRSA infection, while the control group included 36 patients with MSSA infection. The following selection criteria were used: (1) local redness, swelling, fever and pain of the breast, accompanied by a sense of fluctuation; (2) ultrasonographic examination of the lesion showing fluid anechoic areas; and (3) if the abscess is deep, skin redness and swelling are often not obvious; if antipyretics or antibiotics have been applied before hospital visit, the body temperature and white blood cell (WBC) count may not rise. The exclusion criteria were as follows: (1) complicated with other part infections, such as respiratory or urinary system; (2) polymicrobial infection of patients; (3) accompanied by major systemic diseases or tumors; and (4) immunocompromised individuals. The patients’ body mass index (BMI), delivery number, previous history of antibiotic therapy, temperature, WBC, neutrophil, red blood cell (RBC), hemoglobin (HB), platelet (PLT), C-reactive protein (CRP), and ultrasound features of breast abscesses were collected. T-test and chi-square test were used for statistical analysis, and p<0.05 indicated a statistically significant difference.

The CA-MRSA- or MSSA-infected pus samples were collected for the followed dual-omics studies. For label-free analysis, five CA-MRSA and five MSSA samples were used, followed by PRM validation of candidate biomarkers. In the metabolomics analysis, six samples were added to each group, respectively. Finally, the DEPs and DEMs closely related to CA-MRSA infection were screened from the conjoint analysis of proteomics and metabolomics (see Figure 1).

Figure 1

2.2 Label-free quantitative proteomic analysis

A total of 10 samples were minced individually with liquid nitrogen, pretreated with 2mM DL-dithiothreitol (DTT), and lysed in lysis buffer alkylated with indole acetic acid (IAA). Subsequently, the samples were centrifuged and the precipitates were collected, which dissociated in a dissolution buffer containing 0.1 M triethylammonium bicarbonate (TEAB, pH 8.5) and 8 M urea. Protein concentration was quantified by Bradford assay. The supernatant from each sample was digested with Trypsin Gold (Promega) at 1:50 enzyme-to-substrate ratio at 37°C for 16 h. The generated peptides were desalted with C18 cartridge to remove the high concentration urea and dried by vacuum centrifugation.

Proteomics analyses were performed using EASY-nLCTM 1200 UHPLC system (Thermo Fisher) coupled with a Q Exactive HF-X mass spectrometer (Thermo Fisher). All the resulting spectra were searched from the protein database by the search engine: Proteome Discoverer 2.2 (PD 2.2, Thermo). The protein quantification results were analyzed statistically using the t-test. The proteins that showed significant differences in quantification between the CA-MRSA and MSSA groups were considered as DEPs (p<0.05 and fold change (FC)>1.5 or FC<0.67). The DEPs were further categorized based on their molecular function (MF), cellular component (CC), and biological process (BP) for Gene Ontology (GO) mapping and annotation. Subsequently, the DEPs were searched against the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.geneontology.org/) to obtain information on enriched functional pathways. Possible protein–protein interactions (PPIs) were predicted using the online STRING database (http://string-db.org/) and Cytoscape software.

2.3 Targeted validation by PRM

The DEPs that were significantly upregulated (FC>2.0) in the CA-MRSA group compared to the MSSA group were selected and were annotated as key proteins or pathways in GO and KEGG enrichment analysis. Using the peak area of the internal standard peptide segment DSPSAPVNVTV R labeled with heavy isotopes as the reference, the absolute quantification of target DEPs between CA-MRSA and MSSA groups were calculated in the PRM analysis after correction.

2.4 Untargeted metabolomics analysis (UHPLC-MS/MS)

The 22 pus samples were separately grinded with liquid nitrogen, and the homogenate was resuspended with prechilled 80% methanol and 0.1% formic acid by well vortexing. After centrifugation at 4°C, some of the supernatant was diluted to a final concentration containing 60% methanol by LC-MS grade water. The samples were subsequently transferred to fresh Eppendorf tubes with 0.22 μm filter. Finally, following centrifugation, the filtrate was injected into the LC-MS/MS system analysis.

LC-MS/MS analyses were performed using a Vanquish UHPLC system (Thermo Fisher) coupled with an Orbitrap Q Exactive series mass spectrometer (Thermo Fisher). The raw data files generated by super-high performance liquid chromatography-mass spectrum in series (UHPLC-MS/MS) were processed using the Compound Discoverer 3.1 (CD3.1, Thermo Fisher) to perform peak alignment and peak picking for each metabolite. Then, the peaks of metabolites were matched against the mzCloud (https://www.mzcloud.org/) and ChemSpider (http://www.chemspider.com/) databases to obtain the accurate and relative quantification. We performed multivariate statistical analysis using principal components analysis (PCA) and partial least squares discriminant analysis (PLS‐DA). Univariate analysis (t-test) was used to calculate the statistical significance (p-value). The metabolites with Variable Importance in the Projection (VIP) >1, p<0.05, and FC≥2/FC≤ 0.5 were deemed to be DEMs. Volcano plots were used to filter metabolites of interest, which were based on Log2 (FC) and -log10 (p-value) of metabolites. For clustering heat maps, the data were normalized using z-scores of the intensity areas of DEMs and were plotted by Python-3.5.0 R-3.4.3. We conducted Pearson correlation analysis to describe the relationship between DEMs. These DEMs were annotated using the KEGG database, HMDB database (http://www.hmdb.ca/), and Lipid Maps database (http://www.lipidmaps.org/). The metabolic pathways enriched by DEMs were considered statistically significant when p<0.05.

2.5 Integrated analysis of proteomics and metabolomics

The prominently screened DEPs and DEMs were measured for their degree of association based on the Pearson correlation coefficient, which ranges from -0.99 to +0.99. A negative correlation coefficient indicates negative correlation, while a positive correlation coefficient indicates positive correlation. All the DEPs and DEMs were mapped to the KEGG pathway database to obtain information on their common pathway enrichment.

2.6 Establishment of prediction model of CA-MRSA infection

We selected the DEMs in co-enrichment pathways of proteomics and metabolomics to predict CA-MRSA infection. The area under the curve (AUC) value of the receiver operating characteristic (ROC) was used to calculate the predictive efficiency of individual and combined DEMs. The association between the DEMs quantitation and purulent cavity volume was assessed using the Pearson correlation coefficient. Statistical analyses were performed using GraphPad Prism (version 9.0) and Stata (version 16.0). A p-value of <0.05 was considered statistically significant.

3 Result

3.1 Overview of the participants characteristics

Table 1 presents the clinical characteristics of 24 breast abscesses with CA-MRSA and 36 with MSSA. The average age of the two groups was 29.08 ± 4.596 and 30.22 ± 5.060, respectively, without any confounding bias (p=0.3796). Compared to the MSSA group, the CA-MRSA group had a significantly lower NEU count (6.721 ± 1.847 vs. 5.254 ± 1.839×10⁹/L, p=0.0038), as shown in Supplementary Figure S1A, and a larger purulent cavity volume (15.85 ± 12.26 vs. 30.34 ± 21.70 cm³, p=0.0016) as shown in Supplementary Figure S1B. Additionally, no significant differences were found in the BMI, WBC, RBC, HB, PLT count, CRP, fever, reproductive number, and previous use of antibiotics between the two groups.

From the perspective of ultrasonic images, the abscess cavity of CA-MRSA infection is characterized by thick walls and multiple compartments as shown in Supplementary Figures S1C, D. Concurrently, most patients have enlarged ipsilateral axillary lymph nodes accompanied by pressure pain. The structure of the axillary lymph node in the CA-MRSA group was disordered, with abundant blood flow signals according to Color Doppler Flow Imaging (CDFI). On the contrary, the structure in the MSSA group was relatively clearer, and there was sparser blood flow signal as shown in Supplementary Figures S1E, F.

3.2 Proteomic profiling and enriched pathways

The whole proteomic profiling identified a total of 1,363 proteins. A total of 109 DEPs (86 upregulated and 23 downregulated) between CA-MRSA and MSSA groups were revealed and displayed as a heatmap and volcano plot (Figures 2A, B). In the GO enrichment analysis, the DEPs were mainly involved in the carbohydrate metabolic process, generation of precursor metabolites and energy, calcium ion binding, and GTP binding in the categories of biological process and molecular function (Figure 2C). The top 20 pathways of DEPs enriched in the KEGG analysis are displayed in the scatter plot of Figure 2D. The first 10 pathways (p<0.05), considered significant, were mainly involved in the biosynthesis of amino acids, glycosphingolipid biosynthesis, glycosaminoglycan degradation, AMPK signaling pathway, glutathione metabolism, fatty acid elongation, and glycerophospholipid metabolism.

Figure 2

Each subcellular structure can provide a relatively independent place for a group of specific proteins to exercise their functions, and only proteins in the same or similar subcellular location will have interactions. The subcellular localization of the screened DEPs showed that nuclear proteins accounted for the highest proportion (20.99%), followed by cytoplasmic proteins (19.75%), endoplasmic reticulum proteins (13.58%), and plasma membrane protein (11.11%), respectively, as shown in Figure 2E. Defining the subcellular localization of the DEPs or expression products in the CA-MRSA group can provide the above directions for subsequent in-depth study of the mechanism. A total of 39 DEPs with significant enrichment variance between CA-MRSA and MSSA in GO and KEGG pathway analysis were selected to construct the interaction network. Figure 2F demonstrates the close relationships between four functional clusters with CA-MRSA infection: biosynthesis of amino acids, glycerophospholipid metabolism, AMPK signaling pathway, and tight junction. These clusters were seen to play crucial and interconnected roles in the occurrence of breast abscess infected by CA-MRSA.

3.3 PRM verification of candidate protein biomarkers

PRM was used to determine the abundance variation of 10 candidate DEPs (A0A140VK56, P07686, B2R983, Q5TEC6, Q53EU6, B4DVA7, V9HWE9, P00558, J3KPS3, and Q7L5N7) between CA-MRSA and MSSA groups. As illustrated in Supplementary Figures S2A–J, the expression levels of the above 10 candidate DEPs were increased in CA-MRSA versus MSSA (p<0.05). The relative expression of the target DEPs and peptides determined by label-free and PRM can be displayed in Supplementary Table S1, which indicates consistent upregulation trends (Supplementary Figures S2K–T).

3.4 Metabolomic profiling and related pathways

Because two samples in the CA-MRSA group originating from pus of the right and left breast abscesses in the same case, they were combined in metabolomics to be analyzed as one sample. In positive and negative ion patterns, the relational model between metabolite expression and group category (CA-MRSA and MSSA) was established by partial least squares-discriminate analysis (PLS-DA) (Figures 3A–D). The models provided 61 and 26 DEMs in the positive and negative ion mode (VIP>1.0 and p<0.05), respectively, including 48 upward and 39 downward. The volcano maps and heatmaps visually depicted the overall distribution of DEMs between CA-MRSA and MSSA groups detected in the positive and negative types (Figures 3E–H). As displayed in Figures 3I, J, in the positive ion mode, DEMs such as n-ribosylhistidine, N-(6-amino-5-nitro-4-pyrimidinyl), nopaline, GDP-N-acetyl-alpha-D-perosamine, and L-(-)-methionine had significant positive correlation with each other with the consistency of changing trends. Pathway enrichment analysis of the DEMs was statistically associated with the regulation of actin cytoskeleton, amino sugar and nucleotide sugar metabolism, complement and coagulation cascades, and cGMP-PKG signaling pathway in negative pattern (p<0.05). Other pathways in the positive pattern, such as the sphingolipid signaling pathway and tryptophan metabolism, failed to attain statistical significance, but their p-values were slightly >0.05, suggesting potential contributions as CA-MRSA-infected breast abscesses related metabolic pathways (Figures 3K, L).

Figure 3

3.5 Integrated network analysis of proteomic and metabolomic

The critical DEPs and DEMs screened from proteomic and metabolomic analysis were applied to clarify the correlation degree and inter-reaction (Supplementary Table S2). In Figures 4A, B, the DEPs of A0A140VK56, P00558, J3KPS3, and B4DVA7 were consistent with bradykinin, 3-hydroxypentadecanoic acid, and D-sedoheptulose 7-phosphate, involve in the regulation of actin cytoskeleton, biosynthesis of amino acids, sphingolipid signaling pathway, and pentose phosphate pathway in negative mode. Furthermore, the above-mentioned DEPs were relevant to nopaline, L-(-)-methionine, indole-3-acetic acid, and n-ribosylhistidine in terms of the amino acid metabolism comprised of amino sugar and nucleotide sugar metabolism, arginine and proline metabolism, and tryptophan metabolism in positive mode (Figures 4C, D).

Figure 4

3.6 Identification of CA-MRSA infection and correlation analysis

The ROC curve was plotted for the nine screened DEMs, and the AUC value was used to assess the diagnostic accuracy of individual DEMs for CA-MRSA infection (Figures 5A–I). To better define the optimal diagnostic model, we performed an ROC analysis for a set of biomarkers, including indole-3-acetic acid, L-(-)-methionine, and D-sedoheptulose 7-phosphate, whose AUC values were the top 3. As presented in Figures 5J, K, Model 2 of the mentioned three metabolites showed the best recognition rate of CA-MRSA infection (AUC value, 0.9421; 95% CI, 0.8933–0.9909) compared with Model 1. Furthermore, we selected the three metabolites in Model 2 to continue the association analysis with purulent cavity volume (Figures 5L–Q). The results showed that only the content of L-(-)-methionine in the CA-MRSA group was negatively correlated with purulent cavity volume, but there were no correlations in the MSSA group or other DEMs.

Figure 5

4 Discussion

To the best of our knowledge, this research is the first attempt to link proteomics with metabolomics to investigate the biomarkers and related mechanisms of CA-MRSA infection. The integrated analysis elaborated the pathways of amino acid metabolism and lipid metabolism co-enriched by upstream DEPs and downstream DEMs, which may be involved in structuring the related network of breast abscess infected by CA-MRSA in Chinese women.

4.1 Amino acid metabolism

Amino acids are the basic units of proteins. The amino acid metabolites of microorganisms and host cells are widely implicated in biological processes, such as inflammation and immune response, by regulating the activation of immune cells and the production of antibodies (Costantini et al., 2020; Dhankhar et al., 2020). These have played essential roles in adaptive and innate immunity. Recent literature has documented the important intermediate metabolites in the aromatic amino-acid-related pathways that undergo various changes in sepsis (Chen Q. et al., 2022). The catabolic pathway of converting phenylalanine to phenylpyruvate can impel the escape state of neutrophils in Acinetobacter baumannii (Rodman et al., 2019). In our study, the DEPs between CA-MRSA and MSSA are mostly related to kinases or transferases in the amino acid metabolic pathways. By regulating the synthesis or catabolism of amino acids, the expressions of downstream DEMs are affected correspondingly, leading to the most prominent and widespread alterations in amino acid biosynthesis, phenylalanine metabolism, tryptophan metabolism, arginine, and proline metabolism in the CA-MRSA group (Figure 6). On the other hand, biofilm is an important pathogenic and drug-resistant cause of CA-MRSA. It not only helps bacteria colonize in the host tissue and spreads toxins in vivo but also provides physical and chemical protective barriers for bacterial survival and reinforces the resistance of bacteria to antibiotics and host immune response clearance. Research has announced that transformations in the microenvironment of L-tyrosine residues in phenylalanine metabolism can cause disturbance of cell wall peptidoglycan synthesis and membrane protein conformation in biofilms of CA-MRSA (Zhang et al., 2022).

Figure 6

Tryptophan (Trp) catabolism can lead to Trp depletion and the accumulation of related metabolites, thereby mediating immune escape. Their influences on the host and microorganisms can be achieved through two different catabolic pathways: (1) approximately 90%–95% of Trp is converted into kynurenine (Kyn) and downstream metabolites via the Kyn pathway (Dehhaghi et al., 2019). (2) Approximately 4%–6% of unabsorbed Trp is metabolized by microorganisms into indole and indolic compounds (Wyatt and Greathouse, 2021), consisting of indole-3-acetic acid (IAA), indole ethanol, and indole-3-propionic acid (IPA) (Li J. et al., 2021). Under the catalysis of tryptophan monooxygenase and indole-3-acetamide hydrolase, microorganisms can decompose tryptophan to IAA. An in vivo study substantiated that IAA is beneficial in ameliorating hepatic inflammation that presents in non-alcoholic fatty liver disease (NAFLD) mice, which may be relevant to the lightened inflammatory responses mediated by macrophages (Ji et al., 2019). Mice fed with a high-fat diet (HFD) were infected by MRSA, whose tryptophan metabolite IAA was dramatically reduced, resulting in insensitivity to antibiotic treatment, while combined usage of IAA and ciprofloxacin increased the survival rate for mice (Liu et al., 2021). Similarly, our finding has linked the abnormal metabolism of IAA to CA-MRSA infection and drug resistance, using IAA as one of the potential biomarkers for early diagnosis. More research should decipher the specific mechanisms behind the changes of IAA levels in CA-MRSA infection.

Furthermore, a low expression of L-methionine was found in the CA-MRSA group, and the level was negatively correlated with the volume of breast abscess. Methionine is the most significant methyl donor in the body. By methylating glutathione peroxidase (GPX) and superoxide dismutase (SOD), methionine could increase the synthesis of endogenous phospholipids, thereby stabilizing the lysosomal membrane and reducing the release of acid phosphatase. Thus, a deficiency of methionine can cause local inflammation that is not easily limited and subsequently spreads rapidly. Additionally, L-methionine has the ability to stimulate T-cell activation and maintain T-cell survival and function. L-Methionine intake increases during T-cell activation and contributes to the production of S-adenosine methionine (SAM), which maintains histone and RNA methylation (Sinclair et al., 2019; Bian et al., 2020). However, L-methionine starvation can decrease H3K79me2 in CD8+T cells (Bian et al., 2020) and H3K4me3 in Th17 cells (Roy et al., 2020), thereby limiting the survival and function of T cells. Therefore, it is speculated that the reduced expression of L-methionine in the CA-MRSA group is a potential mechanism for local inflammatory reaction and immune suppression after CA-MRSA infection, which requires further discussion and elucidation.

4.2 Glycosphingolipid biosynthesis and glycerophospholipid metabolism

Glycosphingolipids (GSLs) are binding lipids composed of ceramide and oligosaccharide chains. They are common components of the eukaryotic cell plasma membrane, extending their polar end of hydrophilic sugar chains outward, thus becoming a sign of cell surface bioactivity. The synthesis of GSLs originates from the meshwork of vesicles and tubules in the smooth endoplasmic reticulum (ER) and completes in the Golgi apparatus (Rizzo et al., 2021), where they modulate the activities of membrane proteins, especially signaling receptors. Therefore, GSLs play an important role in molecular signal transduction, immune response, cell adhesion, and recognition (Zhang et al., 2019). Based on the presence of sialic acid, GSLs are divided into neutral and acidic, with the latter known as gangliosides (Chiricozzi et al., 2021). In our data, the pathway of GSLs biosynthesis-ganglio series, in which the DEPs (ID: P07686, B4DVA7) participate, was markedly enriched in the CA-MRSA group (p=0.003). Meanwhile, the DEPs enriched in glycosphingolipids biosynthesis of the CA-MRSA group was localized in subcellular structures such as the endoplasmic reticulum and Golgi apparatus, which, to some extent, supports the subcellular localization finding that larger proportions of DEPs in the CA-MRSA group belonged to cytoplasmic proteins (19.75%), endoplasmic reticulum proteins (13.58%), and plasma membrane proteins (11.11%). The protein annotation of P07686 and B4DVA7 is both β-hexosaminidase B (HexB), which is a group of the most active glycosyl hydrolase isozymes in the lysosomes. HexB catalyzes the conversion from GM2 to GM3 in ganglioside biosynthesis. Substrates for HexB are acidic GSLs, which are implicated in specific infection. Acidic GSLs sugar chains on the host cell surface are binding receptors for many bacterial toxins and determinants of cellular recognition (Aerts et al., 2019). Since an intensified inflammatory state is marked by an increase in the activity of lysosomal exo-glycosidases in tissues and body fluids, especially the most active HexB (Chojnowska et al., 2018), we propose a hypothesis that the combination of CA-MRSA pathogen and GSLs activates HexB with the local inflammatory response to breast abscesses, which is critical for GSLs synthesis and mediates the specific infection process caused by CA-MRSA.

Our findings suggest that lysophosphatidylcholine acyltransferase 2 (LPCAT2, ID: Q7L5N7) is involved in the glycerophospholipid (GPL) metabolism pathway and is significantly enriched in the CA-MRSA breast abscess group (p=0.036). Abnormal GPL metabolism, represented by phosphatidylcholine (PC), has been confirmed to play a considerable role in inflammation, hypertension, and tumors (Sonkar et al., 2019; Sun et al., 2021; Chen Q. et al., 2022). PC is an important component of the biofilm structure. LPCAT2 can acylate lysophosphatidylcholine (LPC) to form PC, while the main source of LPC is the cleavage of PC by phospholipase A2 (PLA2), and the two routes are mutually interconverted (Liu et al., 2020). Recent evidence suggests that LPC possesses advantageous bactericidal clearance and anti-inflammatory effects by enhancing the permeability and perturbation to the cell membrane of CA-MRSA (Miyazaki et al., 2017). LPC reinforces the sensitivity of multidrug-resistant Acinetobacter baumannii to colistin (Miró-Canturri et al., 2021) and strengthens the bactericidal activity of neutrophils by intensifying H₂O₂ production. In Mycobacterium tuberculosis infection, LPC also restrains the inflammatory cytokines production in macrophages via the cAMP-activated PKA-PI3K-P38 MAPK signaling pathway, which can alleviate excessive inflammation (Lee et al., 2018). As the key enzyme in the transition from LPC to PC, high horizontal LPCAT2 in the CA-MRSA group decreased the expression of LPC and weakened its anti-inflammatory action. The lowered LPC/PC ratios are associated with the severity of the inflammatory process (Boldyreva et al., 2021). The above-mentioned may be the clues for CA-MRSA breast abscesses appearing more severe inflammatory response and insensitivity to antibiotics. In addition, LPCAT2 is highly expressed in inflammatory cells and is thought to contribute to the production of lipid mediators in these cells (Tsoupras et al., 2018). LPCAT2 participates in PC remodeling activity in lipid droplet (LD) (Cotte et al., 2018; Jiang et al., 2019) and acts as a key enzyme in the Lands’ Cycle pathway to synthesize PC (Wang and Tontonoz, 2019), which may also affect the plasma membrane composition and the drugs’ pharmacokinetics. Alterations in lipid composition between sensitive and resistant bacteria lead to ineffective drug treatment based on liposome delivery (Saito et al., 2022). In summary, it can be inferred that the resistance of CA-MRSA may be well correlated with the change in plasma membrane structure. This could constitute future directions for untangling the intricate mechanisms of LPCAT2-mediated CA-MRSA inflammation and drug resistance.

4.3 Limitation

Linking proteomics with metabolomics to probe into the biomarkers and related mechanisms of CA-MRSA infection and the biological microenvironment characteristics of CA-MRSA infection in breast abscess were described from the perspectives of amino acid metabolism and lipid metabolism. For all that, several limitations should be acknowledged. First, our sample size was relatively small so as to be considered exploratory until verified in subsequent large cohort studies. From another viewpoint, the small sample size can still provide statistical confidence for our results due to substantial changes in metabolites that could be easily measured. Second, it would be more understandable and useful if we perform follow-up metabolomics profiling of the isolated MRSA and MSSA to validate our findings. Third, the limited number of potential metabolites screened hindered our further analysis using the machine learning methods.

4.4 Conclusion

This research is the first to describe the proteomic and metabolomic profiles of clinically well-characterized breast abscesses caused by CA-MRSA and MSSA in Chinese women. Most of the DEPs and DEMs associated with CA-MRSA infection were found to be enriched in amino acid metabolism, glycosphingolipid biosynthesis, and glycerophospholipid metabolism pathways. More importantly, indole-3-acetic acid, L-(−)-methionine, and D-sedoheptulose 7-phosphate showed good discriminative abilities in identifying early CA-MRSA infection of breast abscesses in a clinical setting. Given the limited high-quality evidence and multi-omics studies in this field, more insights and benefits are required to elucidate the molecular mechanisms of potential biomarkers and pathways for early diagnosis of CA-MRSA infection and drug resistance in the future.

References

• 1

  AertsJ. M. F. G.ArtolaM.van EijkM.FerrazM. J.BootR. G. (2019). Glycosphingolipids and infection. Potential new therapeutic avenues. Front. Cell Dev. Biol.7, 324–339. doi: 10.3389/fcell.2019.00324

  • CrossRef
  • Google Scholar

• 2

  BianY.LiW.KremerD. M.SajjakulnukitP.LiS.CrespoJ.et al. (2020). Cancer SLC43A2 alters T cell methionine metabolism and histone methylation. Nature585 (7824), 277–282. doi: 10.1038/s41586-020-2682-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  BoccaccioC.Verdaguer BabicV.BottoL.CervettoM. M.CetaniS.PaladinoS.et al. (2014). Methicillin-resistant Staphylococcus aureus (MRSA) isolation in breast abscesses in a Public Maternity. Medicina (B Aires)74 (3), 210–215.

  • Google Scholar
  • View reference in article

• 4

  BoldyrevaL. V.MorozovaM. V.SaydakovaS. S.KozhevnikovaE. N. (2021). Fat of the gut: epithelial phospholipids in inflammatory bowel diseases. Int. J. Mol. Sci.22 (21), 11682. doi: 10.3390/ijms222111682

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CaoL.HuangC.Cui ZhouD.HuY.LihT. M.SavageS. R.et al. (2021). Proteogenomic characterization of pancreatic ductal adenocarcinoma. Cell184 (19), 5031–5052. doi: 10.1016/j.cell.2021.08.023

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  ChenF.DaiX.ZhouC. C.LiK. X.ZhangY. J.LouX. Y.et al. (2022). Integrated analysis of the faecal metagenome and serum metabolome reveals the role of gut microbiome-associated metabolites in the detection of colorectal cancer and adenoma. Gut71 (7), 1315–1325. doi: 10.1136/gutjnl-2020-323476

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  ChenQ.LiangX.WuT.JiangJ.JiangY.ZhangS.et al. (2022). Integrative analysis of metabolomics and proteomics reveals amino acid metabolism disorder in sepsis. J. Transl. Med.20 (1), 123. doi: 10.1186/s12967-022-03320-y

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  ChiricozziE.AureliM.MauriL.Di BiaseE.LunghiG.FazzariM.et al. (2021). Glycosphingolipids. Adv. Exp. Med. Biol.1325, 61–102. doi: 10.1007/978-3-030-70115-4_3

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  ChojnowskaS.Cabaj-WiaterI.Mikulska-BaranA.Zalewska-SzajdaB.WaszkiewiczN. (2018). Positive correlation of serum N-Acetyl-β-hexosaminidase with markers of atherosclerosis in diabetes type 2 patients with mild symptoms of depression and anxiety. Dis. Markers2018, 1760592. doi: 10.1155/2018/1760592

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  CostantiniC.BelletM. M.RengaG.StincardiniC.BorghiM.ParianoM.et al. (2020). Tryptophan co-metabolism at the host-pathogen interface. Front. Immunol.11. doi: 10.3389/fimmu.2020.00067

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  CotteA. K.AiresV.FredonM.LimagneE.DerangèreV.ThibaudinM.et al. (2018). Lysophosphatidylcholine acyltransferase 2-mediated lipid droplet production supports colorectal cancer chemoresistance. Nat. Commun.9 (1), 322. doi: 10.1038/s41467-017-02732-5

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  DehhaghiM.Kazemi Shariat PanahiH.GuilleminG. J. (2019). Microorganisms, tryptophan metabolism, and kynurenine pathway: A complex interconnected loop influencing human health status. Int. J. Tryptophan Res.12, 1178646919852996. doi: 10.3390/nu13041189

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  DhankharR.GuptaV.KumarS.KapoorR. K.GulatiP. (2020). Microbial enzymes for deprivation of amino acid metabolism in Malignant cells: Biological strategy for cancer treatment. Appl. Microbiol. Biotechnol.104 (7), 2857–2869. doi: 10.1007/s00253-020-10432-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  DongL.LuD.ChenR.LinY.ZhuH.ZhangZ.et al. (2022). Proteogenomic characterization identifies clinically relevant subgroups of intrahepatic cholangiocarcinoma. Cancer Cell40 (1), 70–87. doi: 10.1016/j.ccell.2021.12.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  JiY.GaoY.ChenH.YinY.ZhangW. (2019). Indole-3-acetic acid alleviates nonalcoholic fatty liver disease in mice via attenuation of hepatic lipogenesis, and oxidative and inflammatory stress. Nutrients11 (9), 2062. doi: 10.3390/nu11092062

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  JiangH.LiZ.HuanC.JiangX. C. (2019). Macrophage lysophosphatidylcholine acyltransfer-ase 3 deficiency-mediated inflammation is not sufficient to induce atherosclerosis in a mouse model. Front. Cardiovasc. Med.5. doi: 10.3389/fcvm.2018.00192

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  KarczewskiK. J.SnyderM. P. (2018). Integrative omics for health and disease. Nat. Rev. Genet.19 (5), 299–310. doi: 10.1038/nrg.2018.4

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  KrogerusC.WernhedenE.HansenL. B. (2019). Mastitis. Ugeskr Laeger181 (47), V07190396.

  • Google Scholar
  • View reference in article

• 19

  LeeH. J.KoH. J.SongD. K.JungY. J. (2018). Lysophosphatidylcholine promotes phagosome maturation and regulates inflammatory mediator production through the protein kinase A-phosphatidylinositol 3 Kinase-p38 mitogen-activated protein kinase signaling pathway during mycobacterium tuberculosis infection in mouse macrophages. Front. Immunol.9. doi: 10.3389/fimmu.2018.00920

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  LiJ.ZhangL.WuT.LiY.ZhouX.RuanZ. (2021). Indole-3-propionic acid improved the intestinal barrier by enhancing epithelial barrier and mucus barrier. J. Agric. Food Chem.69 (5), 1487–1495. doi: 10.1021/acs.jafc.0c05205

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  LiY.MaX. J.HeX. P. (2021). Clinical characteristics of lactational breast abscess caused by methicillin-resistant Staphylococcus aureus: hospital-based study in China. Int. Breastfeed J.16 (1), 80. doi: 10.1186/s13006-021-00429-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  LiuY.YangK.JiaY.ShiJ.TongZ.FangD.et al. (2021). Gut microbiome alterations in high-fat-diet-fed mice are associated with antibiotic tolerance. Nat. Microbiol.6 (7), 874–884. doi: 10.1038/s41564-021-00912-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  LiuP.ZhuW.ChenC.YanB.ZhuL.ChenX.et al. (2020). The mechanisms of lyso-phosphatidylcholine in the development of diseases. Life Sci.247, 117443. doi: 10.1016/j.lfs.2020.117443

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  Miró-CanturriA.Ayerbe-AlgabaR.Jiménez-MejíasM. E.PachónJ.SmaniY. (2021). Efficacy of Lysophosphatidylcholine as Direct Treatment in Combination with Colistin against Acinetobacter baumannii in Murine Severe Infections Models. Antibiotics (Basel)10 (2), 194. doi: 10.3390/antibiotics10020194

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  MiyazakiH.MidorikawaN.FujimotoS.MiyoshiN.YoshidaH.MatsumotoT. (2017). Antimicrobial effects of lysophosphatidylcholine on methicillin-resistant Staphylococcus aureus. Ther. Adv. Infect. Dis.4 (4), 89–94. doi: 10.1177/2049936117714920

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  NiH.PanW.JinQ.XieY.ZhangN.ChenK.et al. (2021). Label-free proteomic analysis of serum exosomes from paroxysmal atrial fibrillation patients. Clin. Proteomics18 (1), 1. doi: 10.1186/s12014-020-09304-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  PevznerM.DahanA. (2020). Mastitis while breastfeeding: prevention, the importance of proper treatment, and potential complications. J. Clin. Med.9 (8), 2328. doi: 10.3390/jcm9082328

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  RizzoR.RussoD.KurokawaK.SahuP.LombardiB.SupinoD.et al. (2021). Golgi maturation-dependent glycoenzyme recycling controls glycosphingolipid biosynthesis and cell growth via GOLPH3. EMBO J.40 (8), e107238. doi: 10.15252/embj.2020107238

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  RodmanN.MartinezJ.FungS.NakanouchiJ.MyersA. L.HarrisC. M.et al. (2019). Human pleural fluid elicits pyruvate and phenylalanine metabolism in Acinetobacter baumannii to enhance cytotoxicity and immune evasion. Front. Microbiol.10. doi: 10.3389/fmicb.2019.01581

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  RoyD. G.ChenJ.MamaneV.MaE. H.MuhireB. M.SheldonR. D.et al. (2020). Methionine metabolism shapes T helper cell responses through regulation of epigenetic reprogramming. Cell Metab.31 (2), 250–266. doi: 10.1016/j.cmet.2020.01.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  SaitoR. F.AndradeL. N. S.BustosS. O.ChammasR. (2022). Phosphatidylcholine-derived lipid mediators: the crosstalk between cancer cells and immune cells. Front. Immunol.13. doi: 10.3389/fimmu.2022.768606

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  Schumacher-SchuhA.BiegerA.BorelliW. V.PortleyM. K.AwadP. S.Bandres-CigaS. (2022). Advances in proteomic and metabolomic profiling of neurodegenerative diseases. Front. Neurol.12. doi: 10.3389/fneur.2021.792227

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  SinclairL. V.HowdenA. J.BrenesA.SpinelliL.HukelmannJ. L.MacintyreA. N.et al. (2019). Antigen receptor control of methionine metabolism in T cells. Elife8, e44210. doi: 10.7554/eLife.44210

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  SonkarK.AyyappanV.TresslerC. M.AdelajaO.CaiR.ChengM.et al. (2019). Focus on the glycerophosphocholine pathway in choline phospholipid metabolism of cancer. NMR Bio Med.32 (10), e4112. doi: 10.1002/nbm.4112

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  SunJ.ZhaoM.YangL.LiuX.PacificoL.ChiesaC.et al. (2021). Identification of potential metabolic markers of hypertension in Chinese children. Int. J. Hypertens.2021, 669173 4. doi: 10.1155/2021/6691734

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  TsouprasA.LordanR.ZabetakisI. (2018). Inflammation, not cholesterol, is a cause of chronic. Disease. Nutrients10 (5), 604–641. doi: 10.3390/nu10050604

  • CrossRef
  • Google Scholar

• 37

  WangB.TontonozP. (2019). Phospholipid remodeling in physiology and disease. Annu. Rev. Physiol.81, 165–188. doi: 10.1146/annurev-physiol-020518-114444

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  WyattM.GreathouseK. L. (2021). Targeting dietary and microbial tryptophan-indole metabolism as therapeutic approaches to colon cancer. Nutrients13 (4), 1189. doi: 10.3390/nu13041189

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  YinY.YuZ.ZhaoM.WangY.GuanX. (2020). Comprehensive evaluation of the risk of lactational mastitis in Chinese women: combined logistic regression analysis with receiver operating characteristic curve. Biosci. Rep.40 (3), BSR20190919. doi: 10.1042/BSR20190919

  • CrossRef
  • Google Scholar

• 40

  YouX.GuoB.WangZ.MaH.ZhangX. (2022). Label-free quantitative proteomic analysis of serum exosomes from patients of renal anemia: The Good and the Bad of Roxadustat. Clin. Proteomics19 (1), 21. doi: 10.1186/s12014-022-09358-w

  • CrossRef
  • Google Scholar

• 41

  ZhangT.de WaardA. A.WuhrerM.SpaapenR. M. (2019). The role of glycosphingolipids in immune cell functions. Front. Immunol.10. doi: 10.3389/fimmu.2019.00090

  • CrossRef
  • Google Scholar

• 42

  ZhangX.MaityT. K.RossK. E.QiY.CultraroC. M.BahtaM.et al. (2021). Alterations in the global proteome and phosphoproteome in third generation EGFR TKI resistance reveal drug targets to circumvent resistance. Cancer Res.81 (11), 3051–3066. doi: 10.1158/0008-5472.CAN-20-2435

  • CrossRef
  • Google Scholar

• 43

  ZhangL.YangW.ChuY.WenB.ChengY.MahmoodT.et al. (2022). The Inhibition Effect of Linezolid with Reyanning Mixture on MRSA and its Biofilm is More Significant than That of Linezolid Alone. Front. Pharmacol.12. doi: 10.3389/fphar.2021.766309

  • CrossRef
  • Google Scholar