Next Article in Journal
Engineered Mesoporous Silica-Based Core-Shell Nanoarchitectures for Synergistic Chemo-Photodynamic Therapies
Next Article in Special Issue
The Search for Associations of Serum Proteins with the Presence of Unstable Atherosclerotic Plaque in Coronary Atherosclerosis
Previous Article in Journal
Saffron Extract Attenuates Anxiogenic Effect and Improves Cognitive Behavior in an Adult Zebrafish Model of Traumatic Brain Injury
Previous Article in Special Issue
Biomarkers in Neurodegenerative Diseases: Proteomics Spotlight on ALS and Parkinson’s Disease
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 19
10.3390/ijms231911601
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Proteomics Profiling of Stool Samples from Preterm Neonates with SWATH/DIA Mass Spectrometry for Predicting Necrotizing Enterocolitis

by
David Gagné
1,2,3,
Elmira Shajari
1,2,3,
Marie-Pier Thibault
1,2,3,
Jean-François Noël
4,
François-Michel Boisvert
2,3,
Corentin Babakissa
2,5,
Emile Levy
6,
Hugo Gagnon
4,
Marie A. Brunet
2,5,
David Grynspan
7,
Emanuela Ferretti
8,
Valérie Bertelle
2,9 and
Jean-François Beaulieu
1,2,3,*
1
Laboratory of Intestinal Physiopathology, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
2
Centre de Recherche du Centre Hospitalier Universitaire de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
3
Department of Immunology and Cell Biology, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
4
PhenoSwitch Bioscience Inc., 975 Rue Léon-Trépanier, Sherbrooke, QC J1G 5J6, Canada
5
Department of Pediatrics, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
6
Research Center, Centre Hospitalier Universitaire Ste-Justine, Université de Montréal, Montréal, QC H3T 1C5, Canada
7
Department of Pathology and Laboratory Medicine, Faculty of Medicine, University of British Colombia, Vancouver, BC V6T 2B5, Canada
8
Division of Neonatology, Department of Pediatrics, Children’s Hospital of Eastern Ontario (CHEO) and CHEO Research Institute, Ottawa, ON K1H 8L1, Canada
9
Division of Neonatology, Department of Pediatrics, Faculty of Medicine and Health Sciences, Université de Sherbrooke, Sherbrooke, QC J1H 5N4, Canada
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(19), 11601; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms231911601
Submission received: 24 August 2022 / Revised: 15 September 2022 / Accepted: 20 September 2022 / Published: 1 October 2022
(This article belongs to the Special Issue Proteomics and Its Applications in Disease)
Browse Figures
Review Reports Versions Notes

Abstract

:
Necrotizing enterocolitis (NEC) is a life-threatening condition for premature infants in neonatal intensive care units. Finding indicators that can predict NEC development before symptoms appear would provide more time to apply targeted interventions. In this study, stools from 132 very-low-birth-weight (VLBW) infants were collected daily in the context of a multi-center prospective study aimed at investigating the potential of fecal biomarkers for NEC prediction using proteomics technology. Eight of the VLBW infants received a stage-3 NEC diagnosis. Stools collected from the NEC infants up to 10 days before their diagnosis were available for seven of them. Their samples were matched with those from seven pairs of non-NEC controls. The samples were processed for liquid chromatography-tandem mass spectrometry analysis using SWATH/DIA acquisition and cross-compatible proteomic software to perform label-free quantification. ROC curve and principal component analyses were used to explore discriminating information and to evaluate candidate protein markers. A series of 36 proteins showed the most efficient capacity with a signature that predicted all seven NEC infants at least a week in advance. Overall, our study demonstrates that multiplexed proteomic signature detection constitutes a promising approach for the early detection of NEC development in premature infants.

1. Introduction

Necrotizing enterocolitis (NEC) is a gastrointestinal disease affecting preterm newborns. With a 5–16% occurrence in very-low-birth-weight (VLBW) infants [1,2,3,4] and a mortality rate of 20 to 50% [5,6,7], NEC is one of the most life-threatening conditions for premature infants in neonatal intensive care units (NICU). Even in the survivors, NEC is responsible for various long-term clinical sequels [7,8,9,10]. The etiology and pathology of NEC are clearly multifactorial but are still incompletely understood. As addressed in many seminal reviews, the main risk factors of NEC development include functional and immune intestinal immaturity and dysbiosis, formula feeding and low birth weight [4,11,12,13,14,15,16,17].
NEC prediction remains a challenge in the NICU because its diagnosis mostly relies on clinical and basic radiologic characteristics [18] which are associated with features of other neonatal conditions [13]. Considering these limitations, most of the efforts for improving the prediction of NEC development in the NICU have focused on the search for molecular biomarkers. In addition to being able to distinguish NEC from other gastrointestinal conditions prior to clinical presentation [19,20,21,22], the ideal biomarker candidates should be tested from samples obtained under noninvasive conditions [17,19,21,22,23,24]. It is noteworthy that the detection of NEC before its progression to an advanced stage would allow for preventive intervention such as antibiotic treatment and feeding with a human milk diet [11,13,19,25,26,27].
The search and discovery of disease molecular targets for prognostic and diagnostic purposes is a key role in biomarker analysis studies [28]. Only a few studies have performed biomarker research on samples obtained before the onset of NEC clinical manifestations in VLBW infants. In one recent study, a large panel of biomarkers was tested in blood samples obtained from VLBW infants, but none of them identified the newborns at risk of developing NEC [29]. It has also been found that the urinary intestinal-fatty acid binding protein can predict NEC one day before clinical manifestations [30], while a combination of urinary biomarkers was shown to predict disease severity within the first 6 h of NEC suspicion [31]. The use of stool microbiome features in combination with clinical metadata features has been shown to identify NEC-affected infants more than 24 h before disease onset [32]. Specific epigenetic changes (TLR4 gene methylation) [33] and significant fluctuation in total bile acid levels [34] were also observed in stool samples prior to NEC diagnosis. Finally, our group recently investigated a series of biomarkers in the stools of VLBW infants collected prospectively [35]. Among the tested biomarkers—which included fecal calprotectin, lysozyme, haptoglobin and intestinal alkaline phosphatase, which were previously found to be indicative of NEC near diagnosis time [36,37,38,39] and lipocalin-2, a key marker associated with the mucosal samples of NEC infants [40]—only the calprotectin and lipocalin-2 combination was found to provide a predictive signature for NEC development up to a week in advance of the diagnosis.
In the present study, we have explored SWATH mass spectrometry-based proteomics to further investigate the potential of fecal biomarkers for predicting NEC development in stool samples collected from VLBW infants in the context of a multi-center prospective study [35].

2. Results

2.1. Discovery of Up- and Down-Regulated Proteins in Stools of NEC Infants by SWATH-MS Analysis

As summarized in the flow chart (Figure 1), stool samples harvested from the seven VLBW babies over the 10-day period preceding the diagnosis of NEC (identified as day 0), as well as those of the 14 matched non-NEC VLBW infants, were pooled under three sets according to the period preceding the diagnosis. Set A was a pool of the samples obtained 10 to 7 days before diagnosis for each infant, set B was for the samples obtained 6 to 3 days before diagnosis and set C was for the samples obtained between 2 days before diagnosis and the day of the diagnosis.
Two distinct spectral libraries were generated for comparison. The first one was consolidated from the FragPipe analyses of DDA and SWATH data at 5% FDR comprising 8957 peptides (10,029 precursors) and 1044 unique proteins. The second is the reference spectral library built with DIA-NN using the library-free search mode at 1% FDR and the FASTA file of protein targets identified from the FragPipe analyses, comprising 12646 peptides (15,404 precursors) and 1175 protein groups composed of 1046 unique proteins. As expected, most of the unique proteins (Figure 2A) and peptides (Figure 2B) were shared between the two libraries, although DIA-NN identified more peptides at 1% FDR than FragPipe at 5% FDR. To confirm the validity of the DIA-NN library-free approach, the two libraries were tested with DIA-NN for the generation of quantified results.
Using the same parameter for the analysis, both libraries gave comparable results; however, more protein groups were quantified in total (1085 DIA-NN vs. 1022 FragPipe; 957 shared) and with less than 10% missing values (489 DIA-NN vs. 451 FragPipe; 416 shared; Figure 2C) with the reference spectral library built from the DIA-NN library-free search.
Starting with the 489 protein groups with less than 10% missing values from the reference spectral library, potential protein markers were selected with an area under the curve (AUC) of receiving operating characteristic (ROC) curves of 0.7 and above. Among these, 80 (69 protein groups of a single protein gene; 11 of two or more protein genes; 86 human genes in total), which were also filling the other selection criteria, as defined in M&M, were selected, and half of them were found to be increased in the stools of NEC infants, the other half being reduced (Figure 3).
The list of these protein markers along with precursor information is provided in Supplementary Table S1. The principal functions of these selected protein markers were assessed by gene ontology analysis of the biological process (GO_BiologicalProcess-EBI-UniProt-GOA-ACAP-ARAP) using the ClueGO app [41]. In total, 125 genes picked from 69 of the 86 human genes were clustered between eight main groups of related biological process-enriched terms, represented by their most significant term, which were: humoral immune response (32% of the total clustered genes), acute inflammatory response (21%), proteolysis (18%) and antioxidant activity (10%), as well as tight junction assembly (5%) and digestive process (endopeptidase inhibitor activity, carbohydrate process and digestive system process for 6%, 5% and 4%, respectively) (Figure 4 and Supplementary Table S2).

2.2. Optimization of SWATH-MS-Derived Stool Signatures for NEC Infants

One objective of this study was to generate stool proteomic signatures that could best predict NEC development. We first chose to use principal component analysis (PCA) [42,43] to explore the discriminating information provided by the 80 potential protein markers and to evaluate sub-selections of this marker set. This multivariate unsupervised dimension reduction method was used to evaluate the capacity of the series of protein markers to cluster the samples adequately based on their outcomes. For the exploratory purpose, only the first two principal components (PC), which summarize the most variability information contained in the original dataset, were selected. As shown in Figure 5, the PC score plot illustrated the relationship between the NEC and non-NEC control samples (Figure 5A), while the protein loading plot (Figure 5B) showed how each protein marker correlated with the score plots. On the loading plot, the clear parting along the PC1 axis of increased (1 to 40) and decreased (41 to 80) markers in NEC indicates that the variation in abundance between these two sets contributes strongly to the divergent distribution of the two types of samples on the same axis as the score plots (Figure 5A). Like the markers (Figure 5B), the sample distribution on the PC2 axis formed more of a continuum (Figure 5A), but, interestingly, more matched pairs (e.g., N2A and C2A) were found on the same side of the axis (13 pairs) than on the opposite side (8 pairs), suggesting that individual characteristics have an impact on marker levels, which could explain at least one part of the variability summarized by the PC2 on both plots. The first two PCs accounted for 33.6% of the total variance of the dataset (Table 1); however, sample separation into discrete clusters was visibly incomplete (Figure 5A), as confirmed by a PC regression of PC score coordinates (PC1 and PC2) [44] that misclassified the membership of two samples (data not shown). Interestingly, parallel analysis [45] generated by Prism to correct for sampling variability and noise selected the first six PCs (Table 1) as those that represent actual variance (eigenvalues greater than the 95th percentile) in the dataset.
To increase the variance explained by the first PCs, we reduced the number of predictor variables using ROC curve information to guide the sub-selection of complementary marker series that can predict the actual outcomes of NEC and non-NEC samples. As expected, higher cumulative proportions of variance, explained by PC1 and PC2, were found for the best 36 (18 up and 18 down; 37.1%), 20 (10 up, 10 down; 41.8%) and 14 (7 up, 7 down; 47.2%) series of protein markers (Table 1), while parallel analysis selected the first 4 PCs of the 36 series and the first PC of the 20 and 14 series (Table 1). Accordingly, the PC score plots of the best 36 (Figure 6A), 20 (Figure 6B) and 14 (Figure 6C) series show a complete separation of NEC and non-NEC sample clusters, as confirmed by PC regression (Table 2).
To further assess the 36, 20 and 14 series’ capacity to discriminate the NEC and non-NEC control samples, we tested two other classical unsupervised clustering methods (hierarchical clustering and K-means clustering) and one supervised classification method (linear discriminant analysis [42,43]) directly on the normalized scores of the expression values. Hierarchical cluster analysis was performed using the farthest neighbor method with the Pearson correlation interval. Discrete clustering was achieved with the 36 and 20 series, whereas the 14 series also formed a two-main-clusters solution, but four samples were misclustered (Table 3). K-means clustering analysis using a two-clusters solution gave similar results. Although none of the series provided complete discrete clustering, the 36 series was the more accurate one, with one sample attributed to the wrong cluster, while the 20 series had two and the 14 series had four misclustered samples (Table 3). Linear discriminant analysis is a multivariate classification method that aims to find an optimal linear function from variables that best explain predetermined classes [42,43]. This method succeeded in correctly classifying all the samples with the 36 series but not with the 20 and the 14 series, which misclassified one (Table 3). The separation of the group centroids (Chi-square test) was significant for the three series, but the lower Wilk’s Lambda value [46] of the 36 series indicated a better ability of the discriminant function to separate the cases into groups compared to the other two series.
Understandably, the counts of samples misclustered or misclassified for the three periods (A, B and C) through these three methods indicate that those from group A are more difficult to differentiate (Table 2 and Table 3). However, overall, the results show that it is feasible to discriminate between all the NEC and non-NEC samples. The signature of the 36 series showed the most efficient discrimination capacity, emphasizing the importance of selecting the right number of complementary markers to optimize the sorting of the samples. The list of these 36 selected protein markers is provided in Table 4.

2.3. SWATH-MS-Derived Stool Signature for Predicting NEC Development

As shown in Figure 7, and as described in the previous section, the discrete clustering of NEC (N1 to N7 group A, B, C: in red) and non-NEC samples (C1 to C7 group A, B, C: in blue) on the dendrogram of the 36 series confirms its capacity to discriminate them according to their respective type. This analysis formed relatively well-defined subclusters founded on individual cases (e.g., N5A, N5B and N5C), which seem more ostensible for the NEC samples, supporting that the sample of the same type also shows variability in association with individual characteristics. The heat map linked to the dendrogram (Figure 7) showed normalized scores of protein levels using the Rankit proportion estimation formula. The x-axis displays the 36 protein markers with the corresponding number from Table 4 (in brackets), and the y-axis displays the 42 samples in the order determined by the analysis. This heat map clearly shows the distinct profiles of the increased and decreased markers and provides a good illustration of how complementarity between and within the group of markers contributes to the discrete clustering.
The relative score of the samples for the 36 stool protein markers was evaluated as a potential tool to predict NEC development at a distinct time before the diagnosis for each subject pair. Relative scores were determined by calculating the means of the standardized values (z-score) of each sample for the 18 increased and 18 decreased protein markers separately, converting the sign of the decreased mean values (multiply by −1) and adding them to the corresponding increased mean values. As shown in Figure 8, the signature from the combination of the 36 stool proteins discriminated all seven NEC from non-NEC infants at least a week in advance. To further test the stool protein signature as a predicting tool for all VLWB infants born between 24 and 29 weeks, the seven NEC and their corresponding non-NEC controls were pooled. As shown in Figure 9, their relative expression scores were found to be significantly different from the control group for both groups of 18 proteins at all stages, as well as in combination, indicating that the 36 protein-based signatures can predict NEC development more than a week in advance for all seven tested cases.

3. Discussion

In this study, we have explored the power of SWATH mass spectrometry-based proteomics to investigate the potential of predicting NEC development from the stools of premature infants. Major efforts toward the discovery of predictive biochemical and clinical markers for this devastating disease have been made over the last decade, but, as mentioned above, they have had limited success in terms of days of predictability prior to the onset of clinical manifestations. This information is a key factor for allowing for an efficient preventive intervention in the NICU for NEC patients [17,47] and for avoiding the overtreatment of non-NEC infants [47,48].
Stool samples present many advantages for investigating NEC-related biomarkers. First, they can be harvested in a completely non-invasive way, even for VLBW, without the significant additional burden for NICU personnel. Second, they are likely to directly reflect the intestinal alterations that precede NEC development, as illustrated by the mucosal release and fecal accumulation of calprotectin and lipocalin2 [35], methylated TLR4 [33], microRNAs [49] and the fluctuation of fecal bile acid levels [34]. Although bowel movements may be delayed in some VLBW infants, most studies were able to successfully deal with this issue for their stool analyses [50,51,52], although, in one study, the low recovery of stool specimens at the time of the NEC acute phase was reported to reduce the clinical value of the test for NEC diagnosis [49]. In agreement with the preceding studies, herein, we have obtained a sufficient number of stool samples for each of the three selected periods to test seven of the eight NEC cases up to 10 days before the diagnosis [35].
There are only a few studies so far that have used mass spectrometry-based proteomics for discovering candidate biomarkers for NEC using samples obtained from blood [53,54], saliva [55], urine [56] or tissues [57], which have led to the identification of proteins displaying accurate diagnostic and prognostic information, but only few of them have conducted investigation before the suspicion of NEC [17]. The only one was a pilot study where proteomics was used to identify predictive biomarkers for NEC from buccal swabs [55]. To our knowledge, this is the first analysis of stool proteins procured from premature infants. Our strategy was to use SWATH-MS in the DIA mode to combine deep proteome coverage with quantitative consistency and accuracy [58]. The strategy for shotgun proteomics was adapted for the complexity of the protein/peptide matrix of the stool. Firstly, DDA and DIA data analysis with FragPipe ensured the completeness of the protein targets that were analyzed with DIA-NN in the library-free search mode. Secondly, the generation of two independent libraries from two software based on different algorithms allowed us to validate the composition of our reference spectral library and of the quantified results. It is worth noting that the latest version of FragPipe integrates DIA quantification with DIA-NN. Indeed, both the FragPipe/msFragger library and DIA-NN library-free mode have been reported to perform well in different aspects of label-free DIA signal extraction and quantification (e.g., the number of identifications, reproducibility) with DIA-NN [59,60]. Of note, the library-free mode that we used gave a higher number of spectral identifications and quantified proteins, even with the FDR estimate considered to be conservative [61,62].
Overall, 80 proteins with an average of 15 peptides (n = 1–100) were found to be consistently differentially detected in the stools of NEC infants relative to those of non-NEC infants. Collectively, the gene ontology analysis of the significantly differentially expressed proteins revealed that more than 60% of them were associated with immune and inflammatory responses and antioxidant defense, an observation consistent with previous observations of NEC intestinal samples analyzed by RNA-Seq, where 60% of the significant pathways were also identified in samples of patients with Crohn’s disease, including gene families related to immunity, infection, antioxidants and antimicrobials [40,63]. That 18.4% of the proteins in the proteolysis categories were observed in the NEC stool samples is also interesting, since fecal proteolytic and elastase activity were reported to be increased before the onset of ulcerative colitis [64]. Interestingly, the increase in fecal proteolytic activity was associated with gut microbiota changes in patients with ulcerative colitis [64], while only human sequences were considered in this work. In this context, it is pertinent to note that a number of markers involved in proteolysis were pancreatic enzymes such as CELA2A, CELA3A, CELA3B and CTRC, which were all found to be reduced in the stools of NEC infants, an observation consistent with the general immaturity of the gut in premature infants who are most susceptible to developing NEC, which may also include limited pancreatic functions [65].
On the other hand, this study also aimed to evaluate the series of markers needed to find the best signature of NEC. While the evaluation of the 36, 20 and 14 series was founded on different unsupervised clustering and supervised classification analyses, the 36 series was able to discriminate more samples on average and showed a lower (better) Wilk’s Lambda value in linear discriminant analysis. Although the variation of abundance between increased and decreased markers seems central to sorting the sample according to their actual outcome, the PCA and hierarchical clustering analysis results indicated that some level of variability within these sets of markers explains the variability between the samples of the same type. In this sense, reducing the number of markers below a certain threshold may limit the primary and complementary information available to sort the samples, which could explain why the 36 series performed better than the 20 and 14 series. However, as the number of samples is limited in this study, it is essential to mention that we do not consider this specific selection of protein markers to be the absolute NEC signature. Indeed, similar work on larger cohorts will be needed to establish more robust conclusions.
It is pertinent to note that, of the 36 selected proteins, 18 were found to be increased in the stools of NEC infants, while the other 18 were reduced. The findings for some of the stool proteins appear consistent with expectations based on the literature, such as the increase in fecal S100A12 in NEC infants [66] and the elevation of mucosal DEFA5 [40,63], Serpin B1 [67] and peroxiredoxin-1 [57] in NEC infants and/or patients with inflammatory bowel disease. Stool accumulation of brush border enzymes such as trehalase and sucrase-isomaltase is also consistent with the extensive destruction of the small intestinal epithelium noted in NEC infants [57]. However, the reduction in the fecal levels for claudin-3 and SLC26A3 observed herein is more puzzling, since claudin-3 was reported to increase in the urine of neonates with NEC [68], and the expression of SLC26A3, an intestinal brush border chloride anion, is drastically reduced in patients with inflammatory bowel disease and has recently been shown to regulate epithelial barrier integrity [69]. Since these two proteins are reduced in the stools of infants more than a week before NEC is diagnosed, one may speculate that the lower fecal levels of these paracellular components could be indicative of a potential mucosal defect in NEC infants.
Another issue that needs to be pointed out is the absence of the calprotectin S100A8 and 100A9 components and lipocalin-2 (LCN2) in the final list of selected proteins, considering that they were previously identified as a robust stool biomarker combination for predicting more than half of the NEC infants [35]. One key aspect of the current study to consider was that each biomarker was selected on the basis of its ability to discriminate NEC vs. non-NEC infants over the 10-day period preceding diagnosis. Indeed, as observed for calprotectin and LCN2 in the immunoassay [35], approximately half of the S100A8, S100A9 and LCN2 in the NEC samples displayed an overlap with the non-NEC samples. S100A8, S100A9 and LCN2 were therefore not retained based on the AUC of ROC analysis. This is consistent with the fact that our approach favored the robustness of the analysis, in line with the current biomarker discovery strategies [70]. Indeed, the accurate identification and quantification of the proteins (ensured here with the hybrid DDA/DIA strategy and by inferring protein libraries from two independent software) is paramount to biomarker discovery and outweighs the identification of as many proteins as possible [71] while reducing the likelihood of false positives [72].
The optimized signatures based on the selection of the 36 proteins found to be increased and decreased in the stools of the NEC infants were able to predict NEC development in all seven VLWB infants at least one week before the diagnosis. This finding is a major achievement in the context of a non-invasive diagnosis method considering that previous published detection methods were predicting either much less in advance [30,31,32,33] or with less sensitivity [34,35].
It is, however, important to consider this analysis as a proof-of-concept study demonstrating that label-free quantification using SWATH/DIA mass spectrometry-based proteomics assisted with DIA-NN software can successfully be applied to stool analyses for the prediction of NEC development. Indeed, despite being derived from a multi-center prospective study where the stools from 132 VLBW infants were collected, NEC is a relatively infrequent condition occurring in the NICU, so only eight of them received a stage-3 diagnosis, including one for which no stools were available for the 7–10 day period before the diagnosis. In this context, considering the relatively small size of the NEC cohort, the set of the 36 protein markers selected for the establishment of an “NEC signature” should not be seen as definitive, since some variations in the levels of proteins may not be related to the disease but rather are a first step toward the usual strategy for mass spectrometry-based proteomics biomarker discovery, where a small number of samples are analyzed in-depth to select the markers to be tested in the next phases [70,73]. The procedure described herein for the selection of the protein markers should therefore be useful for refining this preliminary NEC signature in a future study with a larger cohort required prior to the design of a validation phase for the selected biomarkers by using targeted proteomics analysis. Another aspect that should be considered in a future study would be the inclusion of samples from premature infants affected with other gut-related diseases to optimize an NEC-specific signature.

4. Materials and Methods

4.1. Sample Collection and Preparation for LC-MS/MS Analysis

Stools from 132 very-low-birth-weight infants were collected on a daily basis in the context of a multi-center prospective study aimed at investigating the potential of fecal biomarkers for NEC prediction [35]. Eight infants received a stage-3 NEC diagnosis. The stools collected up to 10 days before the diagnosis were available for seven of them. They were matched with seven pairs of non-NEC controls and grouped according to three test periods (group A: −10 to −7, group B: −6 to −3 and group C: −2 to 0 days before diagnosis), as described before [35], for a total of 42 unique samples. For the preparation, 100 mg of each stool specimen was solubilized in 1 mL of Tris buffer (25 mM Tris pH 7.5), pooled in groups (as detailed above) for NEC and non-NEC infants and centrifuged (16,000× g; 4 °C) for 30 min. The aqueous phase between the pellet and the floating residuals was recovered and stored at −80 °C until preparation for LC-MS/MS analysis.
The concentration of solubilized proteins in the individual samples was measured using a BCA protein assay (ab102536, Abcam, Cambridge, UK). For the 42 individual sample preparations, volumes corresponding to 40 μg of protein were brought up to 100 μL in a 50 mM Tris buffer pH 8.0 for reduction (dithiothreitol (DTT): 10 mM, 10 min at 65 °C), alkylation (15 mM iodoacetamide, 30 min in the dark at room temperature) and quenching (10 mM DTT) steps. The proteins were recovered by precipitation in acetone for 1 h at −80 °C, followed by centrifugation at 16,000× g for 15 min, and the protein pellets were washed with cold methanol. The enzymatic digestion of proteins was carried out in 100 μL of 50 mM Tris buffer pH 8 using 0.7 μg per sample of the proteolytic mix Trypsin/Lys-C, MS Grade (Promega) for 2 h at 37 °C, continued overnight with an additional 0.7 µg of the proteolytic mix (30:1 w/w protein/protease total ratio) and stopped with 2% formic acid. The cleaning and recovery of the peptides were done with a reverse-phase Strata-X polymeric SPE sorbent column (Phenomenex) according to the manufacturer’s instructions. The recovered peptides were dried under nitrogen flow at 37 °C for 45 min and stored at 4 °C until being resuspended in 20 µL of mobile phase solvent A (see the section below) before LC-MS/MS analysis. To generate an ion library, extracted proteins from a representative pool of samples (three NEC and two non-NEC, covering the three test periods) were separated on a 4–20% polyacrylamide gel. The proteins were separated in 12 fractions and then reduced, alkylated and digested in the gel. Peptides were extracted from the gel by successive rounds of dehydration and sonication and purified using reverse-phase SPE.

4.2. LC-MS/MS Acquisition

LC-MS/MS data acquisitions were carried out at the proteomics facilities of PhenoSwitch Bioscience Inc. (Sherbrooke, QC, Canada). The acquisitions were conducted with an ABSciex TripleTOF 6600 (ABSciex, Foster City, CA, USA) equipped with an electrospray interface with a 25 μm iD capillary and coupled to an Eksigent μUHPLC (Eksigent, Redwood City, CA, USA). Analyst TF 1.8 software was used to control the instrument and for data processing and acquisition. The acquisition was performed in the Data-Dependent Acquisition (DDA) mode for the 12 fractions of the ion library. The 42 individual samples (10 µg) were analyzed in the Sequential Window Acquisition of All Theoretical Mass Spectra (SWATH) acquisition mode [74]. The source voltage was set to 5.5 kV and maintained at 325 °C, the curtain gas was set at 35 psi, gas one was set at 27 psi and gas two was set at 10 psi. Separation was performed on a reverse-phase Kinetex XB column 0.3 mm i.d., 2.6 μm particles, 150 mm (Phenomenex), which was maintained at 60 °C. Samples were injected by loop overfilling into a 5 μL loop. For the 60 min LC gradient, the mobile phase consisted of the following: solvent A (0.2% v/v formic acid and 3% DMSO v/v in water) and solvent B (0.2% v/v formic acid and 3% DMSO in EtOH) at a flow rate of 3 μL/min.

4.3. MS Data Preparation and Pre-Analysis

For label-free quantification, SWATH/DIA is reported to offer better proteome coverage, reproducibility and quantitative precision compared to DDA [75]. The DIA-NN software [62] can reproduce protein quantification obtained by the semi-automatic but time-consuming curation of precursor spectra and proteins with the software Skyline (MacCoss Lab, University of Washington, Seattle, WA, USA) [76]. However, the computer calculation time from the library-free search mode with DIA-NN increases rapidly with the size of the reference fasta file and the dataset. Trying to build a library from the whole human proteome is not optimal. For this reason, and considering the published workflow on different MS settings and instruments comparing software performance and using cross-compatible workflow [59,60,61], we opted for a hybrid strategy combining the fast and robust FragPipe Graphical User Interface (GUI; v.17.1; https://fragpipe.nesvilab.org/; accessed on 15 March 2022) proteomics platform (integrate MSFragger, Philosopher and EasyPQP) [77,78,79] for pre-analysis of the data and for creating a list of protein targets from preliminary libraries and DIA-NN (GUI; v.1.8) for generating the reference spectral library from this list and for the quantification of the proteins.
To produce an MS data format compatible with FragPipe and DIA-NN, Sciex DDA and SWATH, .WIFF files were converted to .mzML format with MSConvert (GUI) from ProteoWizard (v3.0.22074) [80] using the Peak Picking filter with the “Vendor” algorithm. The .mzML from SWATH data were further converted to pseudo-MS/MS spectra .mzML files using DIA-Umpire [81] in MSConvert. From FragPipe, MSFragger (v.3.4) was used to perform proteomic searches against the human proteome reviewed database (UP000005640; isoforms and contaminants included; www.uniprot.org; accessed on 15 March 2022), with mostly default open search settings (peptide length from 6 to 42, enzyme set to stricttrypsin, missed cleavage set to 1, max fragment charge set to 4, methionine oxidation set as variable modification and carbamidomethylation set as fixed modification) [82,83]. Philosopher was used (v4.1.1) for the proteomic validation (peptide spectrum matches (PSM) validation, protein inference and false discovery rate [FDR] filtering), and EasyPQP (https://github.com/grosenberger/easypqp; accessed on 15 March 2022) was used for spectral library generation from the filtered results.
To maximize the extent of the list of protein targets, three preliminary spectral libraries were built from FragPipe using a permissive 5% FDR to filter the results: one from the DDA data, one from the SWATH pseudo-MS/MS data and one from the DDA and SWATH pseudo-MS/MS data, concurrently. These libraries were combined and filtered to remove the duplicates and counted a total of 1087 protein UniProtKB IDs entries and 1044 unique proteins (isoform IDs of the same protein counted as one protein). Two proteins of interest were added to this list (alkaline phosphatase, germ cell type [ALPG], and bovine beta-lactoglobulin [LBG]) and the “retrieve/ID mapping tool” from Uniprot.org (https://www.uniprot.org/uploadlists/; accessed on 16 March 2022) has been used to create a FASTA file containing the 1089 protein targets and their amino acid sequence.

4.4. Reference Spectral Library Building and Label-Free Quantification Analysis

The reference spectral library was built with DIA-NN using the 1089 protein target database .fasta file as reference. First, an in-silico-predicted spectral library was created from the sequence database with the FASTA digest library-free search/generation and deep learning based-prediction mode enabled. This library was used to analyze the 42 individual .mzML files of the SWATH runs, keeping equivalent settings for precursor length, missed cleavage, charge and modifications, as described above. Match between runs (MBR) was enabled to create a new library and reanalyze the data from this library (precursor FDR and matrix q value FDR set to 1%). The “remove likely interference” and “reannotate” options were activated, and the human proteome reviewed database was used as the reference .fasta file, with protein inference set to genes. The DIA-NN protein inference algorithm (set to genes) creates protein groups according to the gene from which the precursor(s) can be annotated in the reference .fasta file. Precursors from proteotypic peptides are attributed to single protein groups that correspond to a single gene, whereas indistinguishable non-proteotypic precursors are also merged to protein groups that combine the single proteins from which they can originate. The DIA-NN reference spectral library counted 1174 protein groups from 1046 unique protein names, of which 975 protein groups were composed of a single protein, and 199 were composed of two proteins or more. A total of 1085 protein groups were found to have quantified results in the protein groups (pg)_matrix output file, which contains the normalized expression values calculated by DIA-NN. The label-free quantification method applied is based on the MaxLFQ algorithm [84], using only the optimal ion fragments from the top precursors passing the run-specific and global q value threshold (1% FDR) [62].
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium [85] via the PRIDE partner repository [86] with the dataset identifier PXD036178.

4.5. Dataset Composition, Protein Marker Selection and Series Evaluation

From the pg_matrix file, only protein groups with less than 10% missing values were kept for the statistical analysis, for a total of 489 protein groups with ≥38 quantified values on the 42 samples. Since the missing values in the pg_matrix indicate the absence or very low abundance of precursors, zero imputation was applied to missing values. Then, receiver operating characteristic (ROC) curve analysis was performed to assess the binary discriminant value of individual protein groups, and those with an area under the curve (AUC) ≥ 0.7 were kept. To ensure the selection of high-quality and non-biased features, this list of potential protein markers was further refined by eliminating the protein groups with four imputed zero values applied to one particular type of sample (NEC or non-NEC control) and those lacking a precursor with at least 35 quantified values in the precursor (pr)_matrix output file. A maximum of two protein groups containing redundant entries were kept as distinct protein markers. The final list of potential markers thus counted 80 protein groups, 40 with levels increasing in NEC, and 40 with levels decreasing in NEC. Starting with these 80 protein groups, selection for the series of 36, 20 and 14 protein markers was executed successively from the precedent, keeping the protein groups with the best individual AUC, and the others were selected based on their impact on communal ROC curves (increased and decreased taken separately) that collectively show a greater capability to predict the outcomes of NEC and non-NEC samples.
For each series selection, multidimensional analysis was performed using principal component (PC) analysis to observe the behavior of the individual samples and the variables (protein group levels) on the two main dimensions (PC1 and PC2) [87]. PC regressions were based on binary logistic regression of the PC scores from the first two PCs as independent variables and the two a priori group knowledge (NEC and non-NEC control) as dependent variables [44]. Hierarchical cluster analyses were performed with the farthest neighbor (complete linkage based on the similarity of the farthest pair) method and Pearson’s correlation intervals (comparing two vectors of values) [88]. The PC regressions K-means clustering analysis [89] and linear discriminant analysis [90] were used to assess the feasibility of the classification and clustering models from the protein marker series selection. However, no cross-validation methods were tested, since this was only for exploratory purposes, and no validation cohort was available.

4.6. Statistical Analysis and Graph

All calculations and graph plotting were carried out with GraphPad Prism (Prism v.9.3.1; GraphPad Software, Inc) or IBM SPSS Statistics (SPSS v28.0.1.1; IBM), except for the Venn diagrams, which were created with the web-based tool Venny (v.2.1) [91], and the Gene Ontology (GO) analysis, which was conducted with the ClueGO (v.2.5.8) plug-in from the Cytoscape software (v.3.9.1) [41,92].

5. Conclusions

This proof-of-concept study demonstrates that SWATH/DIA mass spectrometry-based proteomics used for analyzing the stools of premature infants represent a promising potential approach for predicting NEC development in NICU. The next steps should include the refinement of this NEC signature in a future study with a larger cohort that should include premature infants affected by non-NEC intestinal diseases and a validation step using targeted proteomics analysis.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms231911601/s1.

Author Contributions

Conceptualization, J.-F.B., E.F., V.B. and D.G. (David Grynspan); methodology, D.G. (David Gagné), J.-F.N., H.G. and M.-P.T.; formal analysis, D.G. (David Gagné) and E.S.; resources, J.-F.B., E.F., V.B. and D.G. (David Grynspan); data curation, D.G. (David Gagné); writing—original draft preparation, D.G. (David Gagné) and J.-F.B.; writing—review and editing, D.G. (David Gagné), E.S., F.-M.B., C.B., E.L., H.G., M.A.B., D.G. (David Grynspan), E.F., V.B. and J.-F.B.; supervision, F.-M.B., E.F. and J.-F.B.; funding acquisition, J.-F.B., D.G. (David Grynspan), C.B., E.F. and V.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Canadian Institutes of Health Research (MOP136991), the Natural Sciences and Engineering Research Council of Canada through a cooperative ENGAGE grant with PhenoSwitch Bioscience Inc. and the Stars Foundation. D.G. (David Gagné) was the recipient of an MITACS postdoctoral fellowship obtained in collaboration with Phenoswitch Bioscience Inc. E.S. was the recipient of a doctoral studentship from the Faculty of Medicine and Health Science of the Université de Sherbrooke. J.F.B. was the recipient of the Canada Research Chair in Intestinal Physiopathology.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board (or Ethics Committee) of the Centre Hospitalier Universitaire de Sherbrooke (protocol code 2007-65, 06-029, last date of approval 4 May 2022).

Informed Consent Statement

Informed consent was obtained from a parent or a guardian for all subjects involved in the study.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository http://www.ebi.ac.uk/pride with the dataset identifier PXD036178.

Acknowledgments

The authors thank the nurses and personnel of the NICUs of the Centre Hospitalier Universitaire de Sherbrooke (CHUS) and The Ottawa Hospital General Campus for cooperation in the daily collection of diapers; Éric Tremblay for the management of the stool sample database; Samuel Lemaire-Paquette from the Biostatistics Service of the Centre de Recherche of the CHUS (CRCHUS) for his assistance in the statistical analysis; and Elizabeth Herring for reviewing the English of the manuscript.

Conflicts of Interest

H.G. and J.F.N. are the CSO and an employee, respectively, of PhenoSwitch Bioscience Inc. The other authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of the data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Neu, J.; Walker, W.A. Necrotizing enterocolitis. N. Engl. J. Med. 2011, 364, 255–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Henry, M.C.; Moss, R.L. Necrotizing enterocolitis. Annu. Rev. Med. 2009, 60, 111–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Heida, F.H.; Stolwijk, L.; Loos, M.H.; van den Ende, S.J.; Onland, W.; van den Dungen, F.A.; Kooi, E.M.; Bos, A.F.; Hulscher, J.B.; Bakx, R. Increased incidence of necrotizing enterocolitis in the Netherlands after implementation of the new Dutch guideline for active treatment in extremely preterm infants: Results from three academic referral centers. J. Pediatr. Surg. 2017, 52, 273–276. [Google Scholar] [CrossRef] [PubMed]
  4. Nino, D.F.; Sodhi, C.P.; Hackam, D.J. Necrotizing enterocolitis: New insights into pathogenesis and mechanisms. Nat. Rev. Gastroenterol. Hepatol. 2016, 13, 590–600. [Google Scholar] [CrossRef] [Green Version]
  5. Fitzgibbons, S.C.; Ching, Y.; Yu, D.; Carpenter, J.; Kenny, M.; Weldon, C.; Lillehei, C.; Valim, C.; Horbar, J.D.; Jaksic, T. Mortality of necrotizing enterocolitis expressed by birth weight categories. J. Pediatr. Surg. 2009, 44, 1072–1075; discussion 1075–1076. [Google Scholar] [CrossRef]
  6. Thyoka, M.; de Coppi, P.; Eaton, S.; Khoo, K.; Hall, N.J.; Curry, J.; Kiely, E.; Drake, D.; Cross, K.; Pierro, A. Advanced necrotizing enterocolitis part 1: Mortality. Eur. J. Pediatr. Surg. 2012, 22, 8–12. [Google Scholar] [CrossRef] [Green Version]
  7. Jones, I.H.; Hall, N.J. Contemporary Outcomes for Infants with Necrotizing Enterocolitis-A Systematic Review. J. Pediatr. 2020, 220, 86–92.e83. [Google Scholar] [CrossRef]
  8. Bazacliu, C.; Neu, J. Necrotizing Enterocolitis: Long Term Complications. Curr. Pediatr. Rev. 2019, 15, 115–124. [Google Scholar] [CrossRef]
  9. Huda, S.; Chaudhery, S.; Ibrahim, H.; Pramanik, A. Neonatal necrotizing enterocolitis: Clinical challenges, pathophysiology and management. Pathophysiology 2014, 21, 3–12. [Google Scholar] [CrossRef]
  10. Pike, K.; Brocklehurst, P.; Jones, D.; Kenyon, S.; Salt, A.; Taylor, D.; Marlow, N. Outcomes at 7 years for babies who developed neonatal necrotising enterocolitis: The ORACLE Children Study. Arch. Dis. Child. Fetal Neonatal Ed. 2012, 97, F318–F322. [Google Scholar] [CrossRef]
  11. Alganabi, M.; Lee, C.; Bindi, E.; Li, B.; Pierro, A. Recent advances in understanding necrotizing enterocolitis. F1000Research 2019, 8, 107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Eaton, S.; Rees, C.M.; Hall, N.J. Current Research on the Epidemiology, Pathogenesis, and Management of Necrotizing Enterocolitis. Neonatology 2017, 111, 423–430. [Google Scholar] [CrossRef] [PubMed]
  13. Bazacliu, C.; Neu, J. Pathophysiology of Necrotizing Enterocolitis: An Update. Curr. Pediatr. Rev. 2019, 15, 68–87. [Google Scholar] [CrossRef] [PubMed]
  14. Samuels, N.; van de Graaf, R.A.; de Jonge, R.C.J.; Reiss, I.K.M.; Vermeulen, M.J. Risk factors for necrotizing enterocolitis in neonates: A systematic review of prognostic studies. BMC Pediatr. 2017, 17, 105. [Google Scholar] [CrossRef]
  15. Berkhout, D.J.C.; Klaassen, P.; Niemarkt, H.J.; de Boode, W.P.; Cossey, V.; van Goudoever, J.B.; Hulzebos, C.V.; Andriessen, P.; van Kaam, A.H.; Kramer, B.W.; et al. Risk Factors for Necrotizing Enterocolitis: A Prospective Multicenter Case-Control Study. Neonatology 2018, 114, 277–284. [Google Scholar] [CrossRef]
  16. Cho, S.X.; Berger, P.J.; Nold-Petry, C.A.; Nold, M.F. The immunological landscape in necrotising enterocolitis. Expert Rev. Mol. Med. 2016, 18, e12. [Google Scholar] [CrossRef] [Green Version]
  17. Agakidou, E.; Agakidis, C.; Gika, H.; Sarafidis, K. Emerging Biomarkers for Prediction and Early Diagnosis of Necrotizing Enterocolitis in the Era of Metabolomics and Proteomics. Front. Pediatr. 2020, 8, 602255. [Google Scholar] [CrossRef]
  18. Bell, M.J.; Ternberg, J.L.; Feigin, R.D.; Keating, J.P.; Marshall, R.; Barton, L.; Brotherton, T. Neonatal necrotizing enterocolitis. Therapeutic decisions based upon clinical staging. Ann. Surg. 1978, 187, 1–7. [Google Scholar] [CrossRef]
  19. Gilfillan, M.; Bhandari, V. Biomarkers for the diagnosis of neonatal sepsis and necrotizing enterocolitis: Clinical practice guidelines. Early Hum. Dev. 2017, 105, 25–33. [Google Scholar] [CrossRef]
  20. Gephart, S.M.; Gordon, P.V.; Penn, A.H.; Gregory, K.E.; Swanson, J.R.; Maheshwari, A.; Sylvester, K. Changing the paradigm of defining, detecting, and diagnosing NEC: Perspectives on Bell’s stages and biomarkers for NEC. Semin. Pediatr. Surg. 2018, 27, 3–10. [Google Scholar] [CrossRef]
  21. Wang, K.; Tao, G.; Sylvester, K.G. Recent Advances in Prevention and Therapies for Clinical or Experimental Necrotizing Enterocolitis. Dig. Dis. Sci. 2019, 64, 3078–3085. [Google Scholar] [CrossRef] [PubMed]
  22. Ng, P.C.; Ma, T.P.; Lam, H.S. The use of laboratory biomarkers for surveillance, diagnosis and prediction of clinical outcomes in neonatal sepsis and necrotising enterocolitis. Arch. Dis. Child. Fetal Neonatal Ed. 2015, 100, F448–F452. [Google Scholar] [CrossRef] [PubMed]
  23. Garg, B.D.; Sharma, D.; Bansal, A. Biomarkers of necrotizing enterocolitis: A review of literature. J. Matern. Fetal Neonatal Med. 2018, 31, 3051–3064. [Google Scholar] [CrossRef] [PubMed]
  24. Beck, L.C.; Granger, C.L.; Masi, A.C.; Stewart, C.J. Use of omic technologies in early life gastrointestinal health and disease: From bench to bedside. Expert Rev. Proteom. 2021, 18, 247–259. [Google Scholar] [CrossRef] [PubMed]
  25. Gephart, S.M.; Newnam, K.M. Closing the Gap between Recommended and Actual Human Milk Use for Fragile Infants: What Will It Take to Overcome Disparities? Clin. Perinatol. 2019, 46, 39–50. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, K.; Tao, G.; Sun, Z.; Sylvester, K.G. Recent Potential Noninvasive Biomarkers in Necrotizing Enterocolitis. Gastroenterol. Res. Pract. 2019, 2019, 8413698. [Google Scholar] [CrossRef]
  27. Floris, I.; Kraft, J.D.; Altosaar, I. Roles of MicroRNA across Prenatal and Postnatal Periods. Int. J. Mol. Sci. 2016, 17, 1994. [Google Scholar] [CrossRef] [Green Version]
  28. Alexovic, M.; Lindner, J.R.; Bober, P.; Longuespee, R.; Sabo, J.; Davalieva, K. Human peripheral blood mononuclear cells: A review of recent proteomic applications. Proteomics 2022, 22, e2200026. [Google Scholar] [CrossRef]
  29. Hoffsten, A.; Markasz, L.; Lilja, H.E.; Olsson, K.W.; Sindelar, R. Early Postnatal Comprehensive Biomarkers Cannot Identify Extremely Preterm Infants at Risk of Developing Necrotizing Enterocolitis. Front. Pediatr. 2021, 9, 755437. [Google Scholar] [CrossRef]
  30. Gollin, G.; Stadie, D.; Mayhew, J.; Slater, L.; Asmerom, Y.; Boskovic, D.; Holden, M.; Angeles, D. Early detection of impending necrotizing enterocolitis with urinary intestinal fatty acid-binding protein. Neonatology 2014, 106, 195–200. [Google Scholar] [CrossRef]
  31. Coufal, S.; Kokesova, A.; Tlaskalova-Hogenova, H.; Frybova, B.; Snajdauf, J.; Rygl, M.; Kverka, M. Urinary I-FABP, L-FABP, TFF-3, and SAA Can Diagnose and Predict the Disease Course in Necrotizing Enterocolitis at the Early Stage of Disease. J. Immunol. Res. 2020, 2020, 3074313. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hooven, T.A.; Lin, A.Y.C.; Salleb-Aouissi, A. Multiple Instance Learning for Predicting Necrotizing Enterocolitis in Premature Infants Using Microbiome Data. In Proceedings of the ACM CHIL ’20: ACM Conference on Health, Inference, and Learning, Toronto, ON, Canada, 2–4 April 2020; pp. 99–109. [Google Scholar] [CrossRef] [Green Version]
  33. Klerk, D.H.; Plosch, T.; Verkaik-Schakel, R.N.; Hulscher, J.B.F.; Kooi, E.M.W.; Bos, A.F. DNA Methylation of TLR4, VEGFA, and DEFA5 Is Associated With Necrotizing Enterocolitis in Preterm Infants. Front. Pediatr. 2021, 9, 630817. [Google Scholar] [CrossRef] [PubMed]
  34. Knapp, S.; Kehring, A.; Stepp, J.; Calton, C.M.; Gephart, S.M.; Bandlamuri, S.; Boyle, K.E.; Dietz, G.I.; Johnson, H.; Romo, R.E.; et al. Elevated Coefficient of Variation in Total Fecal Bile Acids Precedes Diagnosis of Necrotizing Enterocolitis. Sci. Rep. 2020, 10, 249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Thibault, M.P.; Tremblay, E.; Horth, C.; Fournier-Morin, A.; Grynspan, D.; Babakissa, C.; Levy, E.; Ferretti, E.; Bertelle, V.; Beaulieu, J.F. Lipocalin-2 and calprotectin as stool biomarkers for predicting necrotizing enterocolitis in premature neonates. Pediatr. Res. 2022, 91, 129–136. [Google Scholar] [CrossRef]
  36. Pergialiotis, V.; Konstantopoulos, P.; Karampetsou, N.; Koutaki, D.; Gkioka, E.; Perrea, D.N.; Papantoniou, N. Calprotectin levels in necrotizing enterocolitis: A systematic review of the literature. Inflamm. Res. 2016, 65, 847–852. [Google Scholar] [CrossRef]
  37. Salzman, N.H.; Polin, R.A.; Harris, M.C.; Ruchelli, E.; Hebra, A.; Zirin-Butler, S.; Jawad, A.; Martin Porter, E.; Bevins, C.L. Enteric defensin expression in necrotizing enterocolitis. Pediatr. Res. 1998, 44, 20–26. [Google Scholar] [CrossRef] [Green Version]
  38. MohanKumar, K.; Namachivayam, K.; Song, T.; Jake Cha, B.; Slate, A.; Hendrickson, J.E.; Pan, H.; Wickline, S.A.; Oh, J.Y.; Patel, R.P.; et al. A murine neonatal model of necrotizing enterocolitis caused by anemia and red blood cell transfusions. Nat. Commun. 2019, 10, 3494. [Google Scholar] [CrossRef] [Green Version]
  39. Heath, M.; Buckley, R.; Gerber, Z.; Davis, P.; Linneman, L.; Gong, Q.; Barkemeyer, B.; Fang, Z.; Good, M.; Penn, D.; et al. Association of Intestinal Alkaline Phosphatase with Necrotizing Enterocolitis Among Premature Infants. JAMA Netw. Open 2019, 2, e1914996. [Google Scholar] [CrossRef] [Green Version]
  40. Tremblay, E.; Thibault, M.P.; Ferretti, E.; Babakissa, C.; Bertelle, V.; Bettolli, M.; Burghardt, K.M.; Colombani, J.F.; Grynspan, D.; Levy, E.; et al. Gene expression profiling in necrotizing enterocolitis reveals pathways common to those reported in Crohn’s disease. BMC Med. Genom. 2016, 9, 6. [Google Scholar] [CrossRef] [Green Version]
  41. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.H.; Pages, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef] [Green Version]
  42. Lualdi, M.; Fasano, M. Statistical analysis of proteomics data: A review on feature selection. J. Proteom. 2019, 198, 18–26. [Google Scholar] [CrossRef] [PubMed]
  43. Tang, J.; Wang, Y.; Luo, Y.; Fu, J.; Zhang, Y.; Li, Y.; Xiao, Z.; Lou, Y.; Qiu, Y.; Zhu, F. Computational advances of tumor marker selection and sample classification in cancer proteomics. Comput. Struct. Biotechnol. J. 2020, 18, 2012–2025. [Google Scholar] [CrossRef] [PubMed]
  44. Penn, A.M.; Saly, V.; Trivedi, A.; Lesperance, M.L.; Votova, K.; Jackson, A.M.; Croteau, N.S.; Balshaw, R.F.; Bibok, M.B.; Smith, D.S.; et al. Differential Proteomics for Distinguishing Ischemic Stroke from Controls: A Pilot Study of the SpecTRA Project. Transl. Stroke Res. 2018, 9, 590–599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Franklin, S.B.; Gibson, D.J.; Robertson, P.A.; Pohlmann, J.T.; Fralish, J.S. Parallel Analysis: A method for determining significant principal components. J. Veg. Sci. 1995, 6, 99–106. [Google Scholar] [CrossRef] [Green Version]
  46. Chan, Y.H. Biostatistics 303. Discriminant analysis. Singap. Med. J. 2005, 46, 54–62. [Google Scholar]
  47. Neu, J. Necrotizing Enterocolitis: The Future. Neonatology 2020, 117, 240–244. [Google Scholar] [CrossRef]
  48. Kim, J.H.; Sampath, V.; Canvasser, J. Challenges in diagnosing necrotizing enterocolitis. Pediatr. Res. 2020, 88, 16–20. [Google Scholar] [CrossRef]
  49. Ng, P.C.; Chan, K.Y.Y.; Lam, H.S.; Wong, R.P.O.; Ma, T.P.Y.; Sit, T.; Leung, K.T.; Chan, L.C.N.; Pang, Y.L.I.; Cheung, H.M.; et al. A Prospective Cohort Study of Fecal miR-223 and miR-451a as Noninvasive and Specific Biomarkers for Diagnosis of Necrotizing Enterocolitis in Preterm Infants. Neonatology 2020, 117, 555–561. [Google Scholar] [CrossRef]
  50. Lindberg, T.P.; Caimano, M.J.; Hagadorn, J.I.; Bennett, E.M.; Maas, K.; Brownell, E.A.; Matson, A.P. Preterm infant gut microbial patterns related to the development of necrotizing enterocolitis. J. Matern. Fetal Neonatal. Med. 2020, 33, 349–358. [Google Scholar] [CrossRef]
  51. Mai, V.; Young, C.M.; Ukhanova, M.; Wang, X.; Sun, Y.; Casella, G.; Theriaque, D.; Li, N.; Sharma, R.; Hudak, M.; et al. Fecal microbiota in premature infants prior to necrotizing enterocolitis. PLoS ONE 2011, 6, e20647. [Google Scholar] [CrossRef]
  52. Morrow, A.L.; Lagomarcino, A.J.; Schibler, K.R.; Taft, D.H.; Yu, Z.; Wang, B.; Altaye, M.; Wagner, M.; Gevers, D.; Ward, D.V.; et al. Early microbial and metabolomic signatures predict later onset of necrotizing enterocolitis in preterm infants. Microbiome 2013, 1, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ng, P.C.; Ang, I.L.; Chiu, R.W.; Li, K.; Lam, H.S.; Wong, R.P.; Chui, K.M.; Cheung, H.M.; Ng, E.W.; Fok, T.F.; et al. Host-response biomarkers for diagnosis of late-onset septicemia and necrotizing enterocolitis in preterm infants. J. Clin. Investig. 2010, 120, 2989–3000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Chatziioannou, A.C.; Wolters, J.C.; Sarafidis, K.; Thomaidou, A.; Agakidis, C.; Govorukhina, N.; Kuivenhoven, J.A.; Bischoff, R.; Theodoridis, G. Targeted LC-MS/MS for the evaluation of proteomics biomarkers in the blood of neonates with necrotizing enterocolitis and late-onset sepsis. Anal. Bioanal. Chem. 2018, 410, 7163–7175. [Google Scholar] [CrossRef]
  55. Murgas Torrazza, R.; Li, N.; Young, C.; Kobeissy, F.; Chow, M.; Chen, S.; Mai, V.; Sharma, R.; Hudak, M.; Shuster, J.; et al. Pilot study using proteomics to identify predictive biomarkers of necrotizing enterocolitis from buccal swabs in very low birth weight infants. Neonatology 2013, 104, 234–242. [Google Scholar] [CrossRef] [Green Version]
  56. Sylvester, K.G.; Ling, X.B.; Liu, G.Y.; Kastenberg, Z.J.; Ji, J.; Hu, Z.; Wu, S.; Peng, S.; Abdullah, F.; Brandt, M.L.; et al. Urine protein biomarkers for the diagnosis and prognosis of necrotizing enterocolitis in infants. J. Pediatr. 2014, 164, 607–612.e1–7. [Google Scholar] [CrossRef] [Green Version]
  57. Jiang, P.; Smith, B.; Qvist, N.; Nielsen, C.; Wan, J.M.; Sit, W.H.; Jensen, T.K.; Wang, H.; Sangild, P.T. Intestinal proteome changes during infant necrotizing enterocolitis. Pediatr. Res. 2013, 73, 268–276. [Google Scholar] [CrossRef]
  58. Ludwig, C.; Gillet, L.; Rosenberger, G.; Amon, S.; Collins, B.C.; Aebersold, R. Data-independent acquisition-based SWATH-MS for quantitative proteomics: A tutorial. Mol. Syst. Biol. 2018, 14, e8126. [Google Scholar] [CrossRef]
  59. Gotti, C.; Roux-Dalvai, F.; Joly-Beauparlant, C.; Mangnier, L.; Leclercq, M.; Droit, A. Extensive and Accurate Benchmarking of DIA Acquisition Methods and Software Tools Using a Complex Proteomic Standard. J. Proteome Res. 2021, 20, 4801–4814. [Google Scholar] [CrossRef]
  60. Fernández-Vega, A.; Farabegoli, F.; Alonso-Martínez, M.M.; Ortea, I. Comparison of TIMS-PASEF quantitative proteomics data-analysis workflows using FragPipe, DIA-NN, and Spectronaut from a user’s perspective. bioRxiv 2021. [Google Scholar] [CrossRef]
  61. Demichev, V.; Yu, F.; Teo, G.C.; Szyrwiel, L.; Rosenberger, G.A.; Decker, J.; Kaspar-Schoenefeld, S.; Lilley, K.S.; Mülleder, M.; Nesvizhskii, A.I.; et al. High sensitivity dia-PASEF proteomics with DIA-NN and FragPipe. bioRxiv 2021. [Google Scholar] [CrossRef]
  62. Demichev, V.; Messner, C.B.; Vernardis, S.I.; Lilley, K.S.; Ralser, M. DIA-NN: Neural networks and interference correction enable deep proteome coverage in high throughput. Nat. Methods 2020, 17, 41–44. [Google Scholar] [CrossRef] [PubMed]
  63. Tremblay, E.; Ferretti, E.; Babakissa, C.; Burghardt, K.M.; Levy, E.; Beaulieu, J.F. IL-17-related signature genes linked to human necrotizing enterocolitis. BMC Res. Notes 2021, 14, 82. [Google Scholar] [CrossRef] [PubMed]
  64. Galipeau, H.J.; Caminero, A.; Turpin, W.; Bermudez-Brito, M.; Santiago, A.; Libertucci, J.; Constante, M.; Raygoza Garay, J.A.; Rueda, G.; Armstrong, S.; et al. Novel Fecal Biomarkers That Precede Clinical Diagnosis of Ulcerative Colitis. Gastroenterology 2021, 160, 1532–1545. [Google Scholar] [CrossRef]
  65. Lebenthal, E.; Choi, T.S.; Lee, P.C. The development of pancreatic function in premature infants after milk-based and soy-based formulas. Pediatr. Res. 1981, 15, 1240–1244. [Google Scholar] [CrossRef] [Green Version]
  66. Dabritz, J.; Jenke, A.; Wirth, S.; Foell, D. Fecal phagocyte-specific S100A12 for diagnosing necrotizing enterocolitis. J. Pediatr. 2012, 161, 1059–1064. [Google Scholar] [CrossRef]
  67. Uchiyama, K.; Naito, Y.; Takagi, T.; Mizushima, K.; Hirai, Y.; Hayashi, N.; Harusato, A.; Inoue, K.; Fukumoto, K.; Yamada, S.; et al. Serpin B1 protects colonic epithelial cell via blockage of neutrophil elastase activity and its expression is enhanced in patients with ulcerative colitis. Am. J. Physiol. Gastrointest. Liver Physiol. 2012, 302, G1163–G1170. [Google Scholar] [CrossRef]
  68. Thuijls, G.; Derikx, J.P.; van Wijck, K.; Zimmermann, L.J.; Degraeuwe, P.L.; Mulder, T.L.; Van der Zee, D.C.; Brouwers, H.A.; Verhoeven, B.H.; van Heurn, L.W.; et al. Non-invasive markers for early diagnosis and determination of the severity of necrotizing enterocolitis. Ann. Surg. 2010, 251, 1174–1180. [Google Scholar] [CrossRef]
  69. Kumar, A.; Priyamvada, S.; Ge, Y.; Jayawardena, D.; Singhal, M.; Anbazhagan, A.N.; Chatterjee, I.; Dayal, A.; Patel, M.; Zadeh, K.; et al. A Novel Role of SLC26A3 in the Maintenance of Intestinal Epithelial Barrier Integrity. Gastroenterology 2021, 160, 1240–1255.e1243. [Google Scholar] [CrossRef]
  70. Nakayasu, E.S.; Gritsenko, M.; Piehowski, P.D.; Gao, Y.; Orton, D.J.; Schepmoes, A.A.; Fillmore, T.L.; Frohnert, B.I.; Rewers, M.; Krischer, J.P.; et al. Tutorial: Best practices and considerations for mass-spectrometry-based protein biomarker discovery and validation. Nat. Protoc. 2021, 16, 3737–3760. [Google Scholar] [CrossRef]
  71. Mann, M.; Kumar, C.; Zeng, W.F.; Strauss, M.T. Artificial intelligence for proteomics and biomarker discovery. Cell Syst. 2021, 12, 759–770. [Google Scholar] [CrossRef]
  72. Barkovits, K.; Pacharra, S.; Pfeiffer, K.; Steinbach, S.; Eisenacher, M.; Marcus, K.; Uszkoreit, J. Reproducibility, Specificity and Accuracy of Relative Quantification Using Spectral Library-based Data-independent Acquisition. Mol. Cell. Proteom. 2020, 19, 181–197. [Google Scholar] [CrossRef] [Green Version]
  73. Geyer, P.E.; Holdt, L.M.; Teupser, D.; Mann, M. Revisiting biomarker discovery by plasma proteomics. Mol. Syst. Biol. 2017, 13, 942. [Google Scholar] [CrossRef]
  74. Gillet, L.C.; Navarro, P.; Tate, S.; Rost, H.; Selevsek, N.; Reiter, L.; Bonner, R.; Aebersold, R. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: A new concept for consistent and accurate proteome analysis. Mol. Cell. Proteom. 2012, 11, O111.016717. [Google Scholar] [CrossRef] [Green Version]
  75. Bruderer, R.; Bernhardt, O.M.; Gandhi, T.; Xuan, Y.; Sondermann, J.; Schmidt, M.; Gomez-Varela, D.; Reiter, L. Optimization of Experimental Parameters in Data-Independent Mass Spectrometry Significantly Increases Depth and Reproducibility of Results. Mol. Cell. Proteom. MCP 2017, 16, 2296–2309. [Google Scholar] [CrossRef] [Green Version]
  76. Pino, L.K.; Searle, B.C.; Bollinger, J.G.; Nunn, B.; MacLean, B.; MacCoss, M.J. The Skyline ecosystem: Informatics for quantitative mass spectrometry proteomics. Mass Spectrom. Rev. 2020, 39, 229–244. [Google Scholar] [CrossRef]
  77. da Veiga Leprevost, F.; Haynes, S.E.; Avtonomov, D.M.; Chang, H.Y.; Shanmugam, A.K.; Mellacheruvu, D.; Kong, A.T.; Nesvizhskii, A.I. Philosopher: A versatile toolkit for shotgun proteomics data analysis. Nat. Methods 2020, 17, 869–870. [Google Scholar] [CrossRef]
  78. Kong, A.T.; Leprevost, F.V.; Avtonomov, D.M.; Mellacheruvu, D.; Nesvizhskii, A.I. MSFragger: Ultrafast and comprehensive peptide identification in mass spectrometry-based proteomics. Nat. Methods 2017, 14, 513–520. [Google Scholar] [CrossRef] [Green Version]
  79. Teo, G.C.; Polasky, D.A.; Yu, F.; Nesvizhskii, A.I. Fast Deisotoping Algorithm and Its Implementation in the MSFragger Search Engine. J. Proteome Res. 2021, 20, 498–505. [Google Scholar] [CrossRef]
  80. Adusumilli, R.; Mallick, P. Data Conversion with ProteoWizard msConvert. Methods Mol. Biol. 2017, 1550, 339–368. [Google Scholar] [CrossRef]
  81. Tsou, C.C.; Avtonomov, D.; Larsen, B.; Tucholska, M.; Choi, H.; Gingras, A.C.; Nesvizhskii, A.I. DIA-Umpire: Comprehensive computational framework for data-independent acquisition proteomics. Nat. Methods 2015, 12, 258–264. [Google Scholar] [CrossRef] [Green Version]
  82. Chang, H.Y.; Kong, A.T.; da Veiga Leprevost, F.; Avtonomov, D.M.; Haynes, S.E.; Nesvizhskii, A.I. Crystal-C: A Computational Tool for Refinement of Open Search Results. J. Proteome Res. 2020, 19, 2511–2515. [Google Scholar] [CrossRef] [PubMed]
  83. Yu, F.; Teo, G.C.; Kong, A.T.; Haynes, S.E.; Avtonomov, D.M.; Geiszler, D.J.; Nesvizhskii, A.I. Identification of modified peptides using localization-aware open search. Nat. Commun. 2020, 11, 4065. [Google Scholar] [CrossRef] [PubMed]
  84. Cox, J.; Hein, M.Y.; Luber, C.A.; Paron, I.; Nagaraj, N.; Mann, M. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol. Cell. Proteom. MCP 2014, 13, 2513–2526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Deutsch, E.W.; Bandeira, N.; Sharma, V.; Perez-Riverol, Y.; Carver, J.J.; Kundu, D.J.; Garcia-Seisdedos, D.; Jarnuczak, A.F.; Hewapathirana, S.; Pullman, B.S.; et al. The ProteomeXchange consortium in 2020: Enabling ‘big data’ approaches in proteomics. Nucleic Acids Res. 2020, 48, D1145–D1152. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Perez-Riverol, Y.; Bai, J.; Bandla, C.; Garcia-Seisdedos, D.; Hewapathirana, S.; Kamatchinathan, S.; Kundu, D.J.; Prakash, A.; Frericks-Zipper, A.; Eisenacher, M.; et al. The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 2022, 50, D543–D552. [Google Scholar] [CrossRef]
  87. Jolliffe, I.T.; Cadima, J. Principal component analysis: A review and recent developments. Philos. Trans. A Math. Phys. Eng. Sci. 2016, 374, 20150202. [Google Scholar] [CrossRef] [Green Version]
  88. Meunier, B.; Dumas, E.; Piec, I.; Bechet, D.; Hebraud, M.; Hocquette, J.F. Assessment of hierarchical clustering methodologies for proteomic data mining. J. Proteome Res. 2007, 6, 358–366. [Google Scholar] [CrossRef]
  89. Steinley, D. K-means clustering: A half-century synthesis. Br. J. Math. Stat. Psychol. 2006, 59, 1–34. [Google Scholar] [CrossRef] [Green Version]
  90. Huberty, C.J.; Wisenbaker, J.M.; Smith, J.D.; Smith, J.C. Using Categorical Variables in Discriminant Analysis. Multivar. Behav. Res. 1986, 21, 479–496. [Google Scholar] [CrossRef]
  91. Oliveros, J.C. Venny. An Interactive Tool for Comparing Lists with Venn’s Diagrams. Available online: https://bioinfogp.cnb.csic.es/tools/venny/index.html (accessed on 14 March 2022).
  92. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
Ijms 23 11601 g001 550
Figure 1. Flowchart illustrating the main steps of this study.
Figure 1. Flowchart illustrating the main steps of this study.
Ijms 23 11601 g001
Ijms 23 11601 g002 550
Figure 2. Comparison of libraries and quantified results from FragPipe and DIA-NN analysis. (A,B) Venn diagram comparing unique proteins (A) and unique peptides (modified and unmodified; (B)) identified in the two libraries. (C) Venn diagram comparing protein groups with less than 10% missing values quantified with DIA-NN using the combined FragPipe library or DIA-NN library-free mode. FDR: false discovery rate.
Figure 2. Comparison of libraries and quantified results from FragPipe and DIA-NN analysis. (A,B) Venn diagram comparing unique proteins (A) and unique peptides (modified and unmodified; (B)) identified in the two libraries. (C) Venn diagram comparing protein groups with less than 10% missing values quantified with DIA-NN using the combined FragPipe library or DIA-NN library-free mode. FDR: false discovery rate.
Ijms 23 11601 g002
Ijms 23 11601 g003 550
Figure 3. ROC scatterplot summarizing differentially abundant proteins in NEC stools. “Volcano-type” ROC scatter-plot of 489 proteins with ≥38 quantified values. On the plot, the x-axis represents the log2 fold change ratio of the mean protein levels of NEC vs. non-NEC control samples, the y-axis is the area under the curve (AUC) of the receiver operating characteristic (ROC) curve and the orange and green dots show down- and up-regulated protein levels passing the 0.7 AUC threshold (dashed line), respectively. For more clarity, gene symbols are used on the plot (see Supplementary Table S1 for corresponding protein name). Some labels have been skipped.
Figure 3. ROC scatterplot summarizing differentially abundant proteins in NEC stools. “Volcano-type” ROC scatter-plot of 489 proteins with ≥38 quantified values. On the plot, the x-axis represents the log2 fold change ratio of the mean protein levels of NEC vs. non-NEC control samples, the y-axis is the area under the curve (AUC) of the receiver operating characteristic (ROC) curve and the orange and green dots show down- and up-regulated protein levels passing the 0.7 AUC threshold (dashed line), respectively. For more clarity, gene symbols are used on the plot (see Supplementary Table S1 for corresponding protein name). Some labels have been skipped.
Ijms 23 11601 g003
Ijms 23 11601 g004 550
Figure 4. Gene ontology of the biological process. Gene ontology analysis of the 86 human genes associated with the 80 protein markers using the GO_BiologicalProcess-EBI-UniProt-GOA-ACAP-ARAP (22–04-2022) reference term database with the ClueGO (v.2.5.8) plug-in in the Cytoscape software (v.3.9.1). In total, 125 genes picked from 69 of the 86 genes were clustered between the eight main groups of biological process-related terms. On the pie chart, the groups are named based on their leading (most significant) term, and the size of the segment is defined by the % of associated genes per group calculated from the total picked genes. Group-corrected p value calculated by ClueGO with the Bonferroni step-down method: ** p < 0.001.
Figure 4. Gene ontology of the biological process. Gene ontology analysis of the 86 human genes associated with the 80 protein markers using the GO_BiologicalProcess-EBI-UniProt-GOA-ACAP-ARAP (22–04-2022) reference term database with the ClueGO (v.2.5.8) plug-in in the Cytoscape software (v.3.9.1). In total, 125 genes picked from 69 of the 86 genes were clustered between the eight main groups of biological process-related terms. On the pie chart, the groups are named based on their leading (most significant) term, and the size of the segment is defined by the % of associated genes per group calculated from the total picked genes. Group-corrected p value calculated by ClueGO with the Bonferroni step-down method: ** p < 0.001.
Ijms 23 11601 g004
Ijms 23 11601 g005 550
Figure 5. PCA of the 80 protein markers. (A) The PC score plots illustrate the relationship between the NEC (N1 to N7 group A, B, C: in red) and non-NEC control (C1 to C7 group A, B, C: in blue) samples when the dimensionality of the dataset is reduced to its two principal components (PC1: x-axis; PC2: y-axis). (B) Protein loadings show how each protein marker (increased in NEC: from 1 to 40; decreased from 41 to 80; see Supplementary Table S1 for corresponding protein name) is correlated with the two principal components of the PC score plots.
Figure 5. PCA of the 80 protein markers. (A) The PC score plots illustrate the relationship between the NEC (N1 to N7 group A, B, C: in red) and non-NEC control (C1 to C7 group A, B, C: in blue) samples when the dimensionality of the dataset is reduced to its two principal components (PC1: x-axis; PC2: y-axis). (B) Protein loadings show how each protein marker (increased in NEC: from 1 to 40; decreased from 41 to 80; see Supplementary Table S1 for corresponding protein name) is correlated with the two principal components of the PC score plots.
Ijms 23 11601 g005
Ijms 23 11601 g006 550
Figure 6. PCA of the 36, 20 and 14 series of markers. (A–C) PC score plots from the PCA of the (A) 36 (18 up, 18 down), (B) 20 (10 up, 10 down) and (C) 14 (7 up, 7 down) series of protein markers clustering showing the capacity of this series to form discrete clusters of the NEC (N1 to N7 group A, B, C: in red) and non-NEC control (C1 to C7 group A–C: in blue) samples with the two first principal components (PC1: x-axis; PC2: y-axis).
Figure 6. PCA of the 36, 20 and 14 series of markers. (A–C) PC score plots from the PCA of the (A) 36 (18 up, 18 down), (B) 20 (10 up, 10 down) and (C) 14 (7 up, 7 down) series of protein markers clustering showing the capacity of this series to form discrete clusters of the NEC (N1 to N7 group A, B, C: in red) and non-NEC control (C1 to C7 group A–C: in blue) samples with the two first principal components (PC1: x-axis; PC2: y-axis).
Ijms 23 11601 g006
Ijms 23 11601 g007 550
Figure 7. Hierarchical cluster analysis and heatmap of the 36 series. The dendrogram on the left displays the discrete clustering of NEC (N1 to N7 group A, B, C: in red) and non-NEC control samples (C1 to C7 group A, B, C: in blue) between two main clusters. The heat map linked to the dendrogram shows normalized scores of the protein levels for the 36 protein markers (X-axis) with their corresponding number (in brackets; see Table 4 for the protein name) and the 42 individual samples in the order determined by the hierarchical clustering analysis (Y-axis). The colored reference chart on the right indicates the normalized scores (range from −2.26 to 2.26).
Figure 7. Hierarchical cluster analysis and heatmap of the 36 series. The dendrogram on the left displays the discrete clustering of NEC (N1 to N7 group A, B, C: in red) and non-NEC control samples (C1 to C7 group A, B, C: in blue) between two main clusters. The heat map linked to the dendrogram shows normalized scores of the protein levels for the 36 protein markers (X-axis) with their corresponding number (in brackets; see Table 4 for the protein name) and the 42 individual samples in the order determined by the hierarchical clustering analysis (Y-axis). The colored reference chart on the right indicates the normalized scores (range from −2.26 to 2.26).
Ijms 23 11601 g007
Ijms 23 11601 g008 550
Figure 8. Relative score according to paired samples. Plots of the relative score (Y-axis) of the 36 protein marker levels across the three tested periods (X-axis: 10-7, 6-3, 2-0 days before NEC diagnosis) for each of the 21 paired NEC (red dot) and non-NEC control (blue square) samples from the seven paired cases. The two-tailed paired t test on the 36 protein markers (increased and decreased; multiply by −1) gave a significant difference of at least p ≤ 0.05 for all the pairs of the samples.
Figure 8. Relative score according to paired samples. Plots of the relative score (Y-axis) of the 36 protein marker levels across the three tested periods (X-axis: 10-7, 6-3, 2-0 days before NEC diagnosis) for each of the 21 paired NEC (red dot) and non-NEC control (blue square) samples from the seven paired cases. The two-tailed paired t test on the 36 protein markers (increased and decreased; multiply by −1) gave a significant difference of at least p ≤ 0.05 for all the pairs of the samples.
Ijms 23 11601 g008
Ijms 23 11601 g009 550
Figure 9. Relative score by period. The plots show the relative scores (Y-axis) for each of the seven NEC samples (red, yellow and green dots) of the three tested periods (X-axis: 10-7, 6-3, 2-0 days before NEC diagnosis) compared to the average of the three same periods for the non-NEC control samples (Ctrl: 10-0 days; blue dots) to assess the overall relative scores for samples from NEC and non-NEC infants without pairing. The (Increase) and (Decrease) plots were calculated from the 18 increased and the 18 decreased protein markers taken separately, and the (Combined) plot combines the mean values of these two plots for each sample (increased + decreased; multiply by −1). The means ±SD are shown on the plots by horizontal bars. p-values were calculated using one-way ANOVA Dunnett’s multiple comparisons test. ** p < 0.01; *** p < 0.001.
Figure 9. Relative score by period. The plots show the relative scores (Y-axis) for each of the seven NEC samples (red, yellow and green dots) of the three tested periods (X-axis: 10-7, 6-3, 2-0 days before NEC diagnosis) compared to the average of the three same periods for the non-NEC control samples (Ctrl: 10-0 days; blue dots) to assess the overall relative scores for samples from NEC and non-NEC infants without pairing. The (Increase) and (Decrease) plots were calculated from the 18 increased and the 18 decreased protein markers taken separately, and the (Combined) plot combines the mean values of these two plots for each sample (increased + decreased; multiply by −1). The means ±SD are shown on the plots by horizontal bars. p-values were calculated using one-way ANOVA Dunnett’s multiple comparisons test. ** p < 0.01; *** p < 0.001.
Ijms 23 11601 g009
Table
Table 1. PCA information on series’ principal components.
Table 1. PCA information on series’ principal components.
Series80 Markers
PC summaryPC1PC2PC6PC41
Eigenvalue18.238.644.120.04
Proportion of variance22.79%10.80%5.15%0.06%
Cumulative proportion of variance22.79%33.58%62.28%100.00%
Parallel analysis PC selection SelectedSelectedLast selected
Series36 markers
PC summaryPC1PC2PC4PC36
Eigenvalue9.314.052.680.002
Proportion of variance25.86%11.24%7.45%0.00005%
Cumulative proportion of variance25.86%37.10%54.08%100.00%
Parallel analysis PC selection SelectedSelectedLast selected
Series20 markers
PC summaryPC1PC2 PC20
Eigenvalue6.152.21 0.06
Proportion of variance30.77%11.05% 0.28%
Cumulative proportion of variance30.77%41.82% 100.00%
Parallel analysis PC selection Selected
Series14 markers
PC summaryPC1PC2 PC14
Eigenvalue4.721.89 0.08
Proportion of variance33.74%13.49% 0.60%
Cumulative proportion of variance33.74%47.22% 100.00%
Parallel analysis PC selection Selected
PC: principal component
Table
Table 2. PC regression on the first PC score coordinates.
Table 2. PC regression on the first PC score coordinates.
Samples36 Markers20 Markers14 Markers
ClassificationNEC (N)21N21N21N
CTRL (C)21C21C21C
Misclassified
by period		Gr. A000
Gr. B000
Gr. C000
Table
Table 3. Clustering and classification analysis.
Table 3. Clustering and classification analysis.
Hierarchical Clustering36 Markers20 Markers14 Markers
NEC (N)Cluster 121N21N21N/4C
Ctrl (C)Cluster 221C21C17C
Misclustered
by period		Gr. A003
Gr. B001
Gr. C000
K-means clustering36 markers20 markers14 markers
NEC (N)Cluster 121N/1C21N/2C21N/5C
Ctrl (C)Cluster 220C19C16C
Misclustered
by period		Gr. A124
Gr. B001
Gr. C000
Linear discriminant analysis36 markers20 markers14 markers
ClassifiedNEC (N)21N21N/1C21N/1C
Ctrl (C)21C20C20C
Wilk’s Lambda0.0350.1340.169
Significance (chi-square)<0.001<0.001<0.001
Misclassified
by period		Gr. A011
Gr. B000
Gr. C000
Table
Table 4. Selection of the 36 series of increased and decreased protein markers.
Table 4. Selection of the 36 series of increased and decreased protein markers.
No.First Protein Name 1Protein & Protein Groups
(Gene Symbol)		 Quantified Values (42 Samples)Fold Change Ratio 2
NEC/Control		ROC
Curve
AUC		Mann–Whitney Test 3
p-Valueq-Value
More abundant in NEC stools
1Carboxypeptidase A1CPA1423.040.82<0.0010.003
2Peroxiredoxin-1PRDX1;PRDX4392.330.82<0.0010.003
3Glutathione reductase, mitochondrialGSR422.610.81<0.0010.003
4Triosephosphate isomeraseTPI1403.900.800.0010.004
5Peroxiredoxin-4PRDX4412.090.800.0010.004
6Protein S100-A12S100A12394.120.790.0010.004
8Defensin-5DEFA5425.090.760.0030.009
15Leukocyte elastase inhibitorSERPINB1421.620.720.0140.023
17TrehalaseTREH411.630.720.0160.023
22Alpha-2-macroglobulin (1)A2M423.330.83<0.0010.002
23Cathepsin SCTSS421.910.83<0.0010.003
24Alpha-2-macroglobulin (2)A2M;PZP423.330.800.0010.004
25Probable non-functional Ig kappa variable 2D-24IGKV2D-24422.090.780.0020.007
26Angiotensin-converting enzyme 2ACE2422.010.770.0020.007
28Alpha-1-acid glycoprotein 1 (1)ORM1424.820.760.0030.009
32Ig heavy constant gamma 1IGHG1421.960.730.0100.018
35Ig lambda variable 8-61IGLV8-61421.800.720.0160.023
40Sucrase-isomaltase, intestinalSI421.680.700.0520.055
Less abundant in NEC stools
41Pancreatic secretory granule mb. major gp GP2GP2420.360.88<0.001<0.001
42Sulfate transporterSLC26A2410.370.88<0.001<0.001
43Tetraspanin-1TSPAN1420.330.85<0.0010.001
44Chloride anion exchangerSLC26A3400.450.81<0.0010.003
45Chymotrypsin-like elastase family member 3A (1)CELA3A420.430.800.0010.004
46Claudin-3CLDN3,4,6,9420.480.770.0020.007
47Chymotrypsin-CCTRC420.410.770.0020.007
49Intestinal-type alkaline phosphataseALPI420.660.750.0050.012
57Ig lambda variable 3-10IGLV3-10420.440.81<0.0010.003
58Mucin-6MUC6420.530.800.0010.004
59Mucin-1MUC1420.430.800.0010.004
60Ig kappa variable 1D-33IGKV1-33420.640.790.0010.004
61Phospholipase B1, membrane-associatedPLB1380.490.790.0010.004
62Ig heavy variable 2-26IGHV2-26380.500.770.0020.007
64TenascinTNC410.220.760.0030.009
65ColipaseCLPS420.500.750.0060.014
76Ig lambda variable 1-51IGLV1-51420.570.700.0280.032
80Cadherin-related family member 2CDHR2420.680.700.0380.041
Abbreviation: AUC: area under the curve; gp: glycoprotein; Ig: immunoglobulin; MMP: metalloproteinase; NEC: necrotizing enterocolitis; No.: number; ROC: receiver operating characteristic. 1 Duplicated first protein names come from protein groups formed by DIA-NN when indifferentiable non-proteotypic precursors can be attributed to more than one gene. 2 The fold change ratio was calculated from the mean expression value of both groups. 3 p-values and q-values have been calculated on the standardized expression values using the Mann–Whitney test for the p-value and Benjamini–Krieger–Yekutieli’s two-stage step-up method for the q-value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gagné, D.; Shajari, E.; Thibault, M.-P.; Noël, J.-F.; Boisvert, F.-M.; Babakissa, C.; Levy, E.; Gagnon, H.; Brunet, M.A.; Grynspan, D.; et al. Proteomics Profiling of Stool Samples from Preterm Neonates with SWATH/DIA Mass Spectrometry for Predicting Necrotizing Enterocolitis. Int. J. Mol. Sci. 2022, 23, 11601. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms231911601

AMA Style

Gagné D, Shajari E, Thibault M-P, Noël J-F, Boisvert F-M, Babakissa C, Levy E, Gagnon H, Brunet MA, Grynspan D, et al. Proteomics Profiling of Stool Samples from Preterm Neonates with SWATH/DIA Mass Spectrometry for Predicting Necrotizing Enterocolitis. International Journal of Molecular Sciences. 2022; 23(19):11601. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms231911601

Chicago/Turabian Style

Gagné, David, Elmira Shajari, Marie-Pier Thibault, Jean-François Noël, François-Michel Boisvert, Corentin Babakissa, Emile Levy, Hugo Gagnon, Marie A. Brunet, David Grynspan, and et al. 2022. "Proteomics Profiling of Stool Samples from Preterm Neonates with SWATH/DIA Mass Spectrometry for Predicting Necrotizing Enterocolitis" International Journal of Molecular Sciences 23, no. 19: 11601. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms231911601

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop