# Peptidomics Analysis Discloses That Novel Bioactive Peptides Participate in Necrotizing Enterocolitis in a Rat Model

Authors: Yiwen Liu; Changlin Wang; Renqiang Yu; Jianfeng Fan; Weilai Jin; Yuting Zhu; Yingzuo Shi; Yulei Jing; Xiaolei Wang; Zhengying Li; Jian Zhou; Le Zhang
Journal: BioMed Research International
 (2020)

---

## Abstract

Necrotizing enterocolitis (NEC) is a common devastating gastrointestinal disease in premature infants, the molecular mechanisms of which have not been fully elucidated. Recently, endogenous peptides have garnered much attention owing to their role in diagnosis and treatment. However, changes in the peptide expression of NEC intestinal tissues remain poorly understood. In the present study, a comparative peptidomics profiling analysis was performed between NEC and control intestinal tissues via liquid chromatography-tandem mass spectrometry (LC-MS). In total, 103 upregulated and 73 downregulated peptides were identified in the intestinal tissues (foldchange≥1.5, p<0.05). Bioinformatics analysis revealed that these differentially expressed peptides were significantly associated with NEC pathophysiology, including apoptosis, the TGF-β signaling pathway, the Wnt signaling pathway, and the MAPK signaling pathway. Furthermore, two putative peptides could inhibit apoptosis and promote the migration of intestinal epithelial cells induced by lipopolysaccharide; these peptides were derived from the protein domains MT1 and EZRI, respectively. In conclusion, our study revealed that endogenous peptides are involved in the pathophysiologic mechanism of NEC; nevertheless, further exploration is required in this regard.


---

## Body

## 1. Introduction

Necrotizing enterocolitis (NEC) is one of the most common fatal diseases in preterm infants, with an incidence rate of 5–10% in very low birthweight (<1500 g) cases [1]. Some studies suggest that the administration of steroids, immunoglobulins, or Lactobacillus can reduce the incidence of NEC in neonates. However, the current treatment strategies predominantly focus on medical stabilization and the prevention of disease progression [2]. The mortality rate may rise to 50% when irreversible necrosis or perforation that requires surgery is present [3]. It is known that the disease progresses primarily via the invasion of inflammatory factors and bacteria in the gastrointestinal tract; however, the underlying molecular mechanisms remain to be elucidated completely [4].Peptidomics, a novel branch of proteomics, has garnered considerable attention for its role as biological markers and therapies [5]. Recently, peptidomics analyses have been applied in detecting peptide changes in different diseases, including cancer [6, 7]. Increasing evidence has shown that endogenous peptides play a crucial role in digestive system diseases [8]. For instance, trefoil factor family peptide 3 (TFF3), mainly located in goblet cells, is significantly upregulated in peptic ulcer. TFF3 not only participates in maintaining and repairing mucosal integrity but also promotes the healing of intestinal ischemia-reperfusion injury in weaned rats [9]. Therefore, an analysis of differentially expressed peptides in NEC tissues may provide new insights to explore treatment strategies.To explore the peptide expression profile of NEC, quantitative liquid chromatography-mass spectrometry (LC-MS) was performed, and the peptide content of rat NEC intestinal tissue (NEC group) and normal control intestinal tissue (Ctrl group) was analyzed. Based on preliminary functional verification, two downregulated peptides derived from the EZRI (EDP1) and MT1 (MDP2) domains showed a significant protective effect in a lipopolysaccharide- (LPS-) induced NEC model in vitro and are worthy of further study.


## 2. Materials and Methods

### 2.1. Ethics Statement

The study was conducted with the approval of the Animal Ethics Committee of Nanjing Medical University (No: IACUC-1807004). To minimize the suffering of the animals used in this research, the procedures followed were in strict compliance with the Guide for the Care and Use of Laboratory Animals.


### 2.2. Establishment of the Animal Model

Eighteen neonatal Sprague–Dawley rats provided by the Animal Core Facility of Nanjing Medical University were randomly assigned to the Ctrl and NEC groups. The rats in the Ctrl group were fed by their mothers as needed; the NEC rat model was established as follows [10]: (1) hypoxia-cold stimulation: asphyxia stress via exposure to 95% nitrogen twice a day for 5 min, followed by exposure to a cold temperature (4°C) for 10 min; (2) formula feeding (ionized water (10 mL) and Esbilac (PetAg) canine milk (5 mL) (2 : 1) mixed with Similac Advance infant formula (Abbott Nutrition) 2 g): the animals received perfused hypertonic milk (0.08–0.1 g/kg) immediately after cold stimulation. Four days for establishing the model, the rats were sacrificed, and intestinal tissue was collected. Intestinal tissue (3 cm) from the terminal ileum was harvested and subjected to NEC model determination and peptide extraction.


### 2.3. Peptide Extraction and Purification

Three pairs of intestinal tissue samples from the Ctrl and NEC groups were ground into powder in liquid nitrogen and fully compatible with denaturing buffers including 8 M urea/100 mM TEAB (pH 8.0), 1 mM PMSF, 2 mM EDTA, and 10 mM DTT. The mixture was sonicated on ice for 10 min and centrifuged at 4°C (12,000 rpm) for 30 min; then, the supernatant was transferred to a fresh centrifuge tube. MWCO (10 kDa) filter tubes (Millipore, Billerica, MA, USA) were used to filter the supernatant. The filter tubes were centrifuged at 4°C (10,000 ×g); flow-through containing the peptides was collected and purified using a C18 column. Finally, the flow-through was dried under vacuum and stored at -80°C.


### 2.4. TMT Labeling and LC-MS Analysis

Formic acid (0.1%) was used to dissolve the peptide samples, followed by reducing in 10 mM DTT and alkylating with 50 mM iodoacetamide for 30 min. Following desalting and lyophilization, the peptides were labeled with TMT reagent (TMT 6 Relabeling Reagent, Thermo Fisher Scientific, USA). The labeled peptides were mixed and analyzed via LC-MS/MS [11]; an LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) was used for this mass spectrometry analysis. In total, 2000 single emission spectra were accumulated from 10 random positions of each sample. Using isotope correction factors recommended by the manufacturer, search for MS data using MASCOT database (Matrix Science, Boston, MA, USA) and SwissProt database. Additionally, to reduce the possibility of erroneous peptide identification, only those peptides with a significance score (≥20) within a 99% confidence interval and greater than “identity” were counted via MASCOT probability analysis [12].


### 2.5. Bioinformatics Analysis

The characteristics of the identified peptides were analyzed using ProtParam tool. Presumed functions on peptides and protein–protein interactions were predicted using the UniProt and STRING databases. We used DAVID v6.7 (Database for Annotation, Visualization and Integrated Discovery;https://david-ncifcrf.gov/) for conducting GO and KEGG pathway enrichment analyses to investigate the potential functions of peptides and their precursor proteins.


### 2.6. Hematoxylin-Eosin Staining (H&E Staining)

After sacrifice, intestinal tissue from the rats’ terminal ileum was harvested and fixed with 4% paraformaldehyde for 1 d at room temperature. Following dehydration, embedding, sectioning, defatting, rehydration, and H&E staining, the tissue sections were observed under a light microscope (Zeiss, Oberkochen, Germany).


### 2.7. Cell Culture

IEC-6 enterocytes were purchased from the American Type Culture Collection (ATCC, USA). The growth medium consisted of high glucose Dulbecco’s Modified Eagle’s medium (DMEM, Gibco BRL), 10% fetal bovine serum (FBS, Gibco BRL), 1% penicillin-streptomycin solution, and 0.1 U/mL recombinant human insulin. The cells were incubated at 37°C under 5% CO2 and treated with 100 μg/mL LPS for 6 h to establish the NEC model in vitro.


### 2.8. Wound-Healing Assay

Cell migration was analyzed using the wound-healing assay. IEC-6 cells were seeded in 6-well plates and cultured in DMEM containing 10% FBS. On reaching confluency, the cell monolayer was scraped using sterile 1 mL pipette tips to form a uniform cell-free area. The wells were gently washed with sterile phosphate-buffered saline (PBS) to remove debris from the medium and monitored under an inverted microscope immediately (Eclipse TS100-F, Nikon) (0 h). Next, fresh medium was added to the wells, and the plates were incubated in the presence of EDP1 or MDP2 (20μm/L) at 37°C for 12 h and monitored under an inverted microscope (12 h). Images taken at the 0 h and 12 h timepoints were expressed as the migration area, and the distances were measured using the ImageJ software v4.16.


### 2.9. Detection of Apoptosis

According to the manufacturer’s instructions, an annexin V-fluorescein isothiocyanate (FITC)/propidium iodide (PI) cell apoptosis kit was used to detect the apoptosis level of the IEC-6 cells. In brief, after incubating the NEC in vitro model with EDP1 and MDP2 for 24 h, the IEC-6 cells were collected and washed with PBS, resuspended in 300μL binding buffer at a concentration of 1×106/mL, and mixed with 3 μL annexin V-FITC and 3 μL PI for 15 min. A BD LSRII Flow Cytometer System with the FACS Diva Software was used to analyze the mixtures within 1 h.


### 2.10. Statistical Analysis

All statistical analyses were performed using the SPSS statistical software. Student’st-test was used to assess significant differences between the two groups. Statistical comparisons between multiple groups were analyzed using ANOVA. Statistical significance was determined as p<0.05.


## 2.1. Ethics Statement

The study was conducted with the approval of the Animal Ethics Committee of Nanjing Medical University (No: IACUC-1807004). To minimize the suffering of the animals used in this research, the procedures followed were in strict compliance with the Guide for the Care and Use of Laboratory Animals.


## 2.2. Establishment of the Animal Model

Eighteen neonatal Sprague–Dawley rats provided by the Animal Core Facility of Nanjing Medical University were randomly assigned to the Ctrl and NEC groups. The rats in the Ctrl group were fed by their mothers as needed; the NEC rat model was established as follows [10]: (1) hypoxia-cold stimulation: asphyxia stress via exposure to 95% nitrogen twice a day for 5 min, followed by exposure to a cold temperature (4°C) for 10 min; (2) formula feeding (ionized water (10 mL) and Esbilac (PetAg) canine milk (5 mL) (2 : 1) mixed with Similac Advance infant formula (Abbott Nutrition) 2 g): the animals received perfused hypertonic milk (0.08–0.1 g/kg) immediately after cold stimulation. Four days for establishing the model, the rats were sacrificed, and intestinal tissue was collected. Intestinal tissue (3 cm) from the terminal ileum was harvested and subjected to NEC model determination and peptide extraction.


## 2.3. Peptide Extraction and Purification

Three pairs of intestinal tissue samples from the Ctrl and NEC groups were ground into powder in liquid nitrogen and fully compatible with denaturing buffers including 8 M urea/100 mM TEAB (pH 8.0), 1 mM PMSF, 2 mM EDTA, and 10 mM DTT. The mixture was sonicated on ice for 10 min and centrifuged at 4°C (12,000 rpm) for 30 min; then, the supernatant was transferred to a fresh centrifuge tube. MWCO (10 kDa) filter tubes (Millipore, Billerica, MA, USA) were used to filter the supernatant. The filter tubes were centrifuged at 4°C (10,000 ×g); flow-through containing the peptides was collected and purified using a C18 column. Finally, the flow-through was dried under vacuum and stored at -80°C.


## 2.4. TMT Labeling and LC-MS Analysis

Formic acid (0.1%) was used to dissolve the peptide samples, followed by reducing in 10 mM DTT and alkylating with 50 mM iodoacetamide for 30 min. Following desalting and lyophilization, the peptides were labeled with TMT reagent (TMT 6 Relabeling Reagent, Thermo Fisher Scientific, USA). The labeled peptides were mixed and analyzed via LC-MS/MS [11]; an LTQ-Orbitrap Velos mass spectrometer (Thermo Fisher Scientific) was used for this mass spectrometry analysis. In total, 2000 single emission spectra were accumulated from 10 random positions of each sample. Using isotope correction factors recommended by the manufacturer, search for MS data using MASCOT database (Matrix Science, Boston, MA, USA) and SwissProt database. Additionally, to reduce the possibility of erroneous peptide identification, only those peptides with a significance score (≥20) within a 99% confidence interval and greater than “identity” were counted via MASCOT probability analysis [12].


## 2.5. Bioinformatics Analysis

The characteristics of the identified peptides were analyzed using ProtParam tool. Presumed functions on peptides and protein–protein interactions were predicted using the UniProt and STRING databases. We used DAVID v6.7 (Database for Annotation, Visualization and Integrated Discovery;https://david-ncifcrf.gov/) for conducting GO and KEGG pathway enrichment analyses to investigate the potential functions of peptides and their precursor proteins.


## 2.6. Hematoxylin-Eosin Staining (H&E Staining)

After sacrifice, intestinal tissue from the rats’ terminal ileum was harvested and fixed with 4% paraformaldehyde for 1 d at room temperature. Following dehydration, embedding, sectioning, defatting, rehydration, and H&E staining, the tissue sections were observed under a light microscope (Zeiss, Oberkochen, Germany).


## 2.7. Cell Culture

IEC-6 enterocytes were purchased from the American Type Culture Collection (ATCC, USA). The growth medium consisted of high glucose Dulbecco’s Modified Eagle’s medium (DMEM, Gibco BRL), 10% fetal bovine serum (FBS, Gibco BRL), 1% penicillin-streptomycin solution, and 0.1 U/mL recombinant human insulin. The cells were incubated at 37°C under 5% CO2 and treated with 100 μg/mL LPS for 6 h to establish the NEC model in vitro.


## 2.8. Wound-Healing Assay

Cell migration was analyzed using the wound-healing assay. IEC-6 cells were seeded in 6-well plates and cultured in DMEM containing 10% FBS. On reaching confluency, the cell monolayer was scraped using sterile 1 mL pipette tips to form a uniform cell-free area. The wells were gently washed with sterile phosphate-buffered saline (PBS) to remove debris from the medium and monitored under an inverted microscope immediately (Eclipse TS100-F, Nikon) (0 h). Next, fresh medium was added to the wells, and the plates were incubated in the presence of EDP1 or MDP2 (20μm/L) at 37°C for 12 h and monitored under an inverted microscope (12 h). Images taken at the 0 h and 12 h timepoints were expressed as the migration area, and the distances were measured using the ImageJ software v4.16.


## 2.9. Detection of Apoptosis

According to the manufacturer’s instructions, an annexin V-fluorescein isothiocyanate (FITC)/propidium iodide (PI) cell apoptosis kit was used to detect the apoptosis level of the IEC-6 cells. In brief, after incubating the NEC in vitro model with EDP1 and MDP2 for 24 h, the IEC-6 cells were collected and washed with PBS, resuspended in 300μL binding buffer at a concentration of 1×106/mL, and mixed with 3 μL annexin V-FITC and 3 μL PI for 15 min. A BD LSRII Flow Cytometer System with the FACS Diva Software was used to analyze the mixtures within 1 h.


## 2.10. Statistical Analysis

All statistical analyses were performed using the SPSS statistical software. Student’st-test was used to assess significant differences between the two groups. Statistical comparisons between multiple groups were analyzed using ANOVA. Statistical significance was determined as p<0.05.


## 3. Results

### 3.1. Establishment of the NEC Model

Hypoxia-cold stimulation combined with hypertonic milk gavage was used to establish the NEC rat model. An optical microscope was used to observe stained intestinal tissue sections and evaluate the pathological score. The histopathological results showed shedding of the intestinal villi, submucosal edema, and inflammatory cell infiltration in the NEC group. In contrast, the intestinal mucosa epithelium was intact, the glands were arranged regularly, and edema and inflammatory cell infiltration were not evident in all layers of the intestinal wall in the Ctrl group (Figure1(a)). In the NEC rat pups, the pathological score [13] was significantly higher than that in the Ctrl rat pups (Figure 1(b)). These findings revealed that the NEC model was successfully established.Figure 1
Establishment of the NEC model: (a) H&E staining of representative intestinal tissue sections from the Ctrl and NEC groups (bar=100μm); (b) pathological scores of representative intestinal tissue sections from the Ctrl and NEC groups (∗∗∗p<0.001).
(a)(b)


### 3.2. Differentially Expressed Peptides in Intestinal Tissue Samples

TMT-labeled LC-MS/MS analysis was performed on peptides in the intestinal tissue samples from both the Ctrl and NEC groups. A total of 480 nonredundant peptides were identified, with 176 differentially expressed peptides derived from 91 protein precursors (foldchange≥1.5, p<0.05). Detailed information on these peptides is shown in Table 1. Upregulated or downregulated differentially expressed peptides were visualized in a hierarchical clustering heat map (Figure 2(a)). The number of upregulated or downregulated differentially expressed peptides is shown in Figure 2(b).Table 1
Differential peptides in the intestine tissue of NEC rat model.
AccessionGenePeptideFold changep valueUpregulated peptidesG3V8B3Hist1h2bqAVTKYTSSK57.72612<0.05P69897Tubb5FVFGQSGAGNN31.84287<0.01Q10758Krt8DGKLVSESSDIMSK25.55914<0.001Q10758Krt8LVSESSDIMSK22.05736<0.001P15999Atp5a1VGLKAPGIIPR21.49717<0.001O88989Mdh1VIVVGNPANTNCLTASK12.14385<0.05Q66HT1AldobALQASALAAWGGK11.93121<0.001G3V8B3Hist1h2bqLLLPGELAKH11.86002<0.001P04182OatALQDPNVAAF11.16052<0.05P05197Eef2VFSGVVSTGLKVR10.57569<0.05Q10758Krt8QPGFGSVGGSST10.57101<0.01F1LNF1Hnrnpa2b1GFGFVTF10.55872<0.05K7S2S2Hist2h2aa3AGLQFPVGR10.3368<0.05Q7M0E3DstnFLWAPEQAPLKSK10.01604<0.01Q10758Krt8GGLTSPGFS9.126557<0.001M0R7B4Hist1h1dSGVSLAALKK8.945185<0.001Q10758Krt8FQPGFGSVGGSST8.793551<0.001P69897Tubb5FVFGQSGAGN8.771316<0.05D3ZGY4Gapdh-ps2GAAQNIIPASTGAAKAVGK8.399639<0.001D3ZBN0Hist1h1bSLVSKGTLVQTK8.352781<0.01M0R7B4Hist1h1dGILVQTKGTGASGS8.334644<0.001P02262NAVGAGAPVYL8.301511<0.001Q10758Krt8GMSSFQPGFGSVGGSST8.18755<0.01P10719Atp5bIGLFGGAGVGK7.935978<0.01B5DEM5Rpl14AAIAAAAAAK7.751116<0.01K7S2S2Hist2h2aa3GLQFPVGR7.629358<0.05P62630Eef1a1ESFSDYPPLGR7.607978<0.001O55159EpcamQKDCVCNNYKLTSR7.485879<0.001P15865Hist1h1eSETAPAAPAAPAPAEKTPIKKK7.38236<0.01D3ZGY4Gapdh-ps2VKVGVNGFGR7.332894<0.001F1M577LOC100359671AFIAHPKLGK6.954795<0.05Q6P725DesTFGGAPGFSLGSPLS6.720618<0.05G3V8B3Hist1h2bqLLLPGELAK6.595403<0.01D3ZIE9Aldh18a1VGLGGMEAKVK6.321181<0.001Q925G0Rbm3DYSGSQGGYDRYSGGN6.259407<0.01Q10758Krt8SFQPGFGSVGGSST6.237045<0.01Q10758Krt8SSFQPGFGSVGGSST6.119025<0.01Q9QXQ0Actn4STALYGESDL5.964186<0.001Q10758Krt8GSLGGFGGAGVGGIT5.600793<0.05P05065AldoaALQASALK5.474993<0.01M0RA08Plin3VTGAVDVTCGAVK5.466526<0.001Q6P9V9Tuba1bVAEITNACFEPAN5.315759<0.05A0A0G2K3Z9NASSGNAKIGHPAPSFK4.969645<0.01Q66H80Arcn1VLLAAAVCTK4.865557<0.05D3ZGY4Gapdh-ps2ALNDNFVK4.838244<0.05B0K031Rpl7IALTDNSLVAR4.834916<0.001P62268Rps23KANPFGGASHAK4.801351<0.001D3Z7Y6Krt20LGVAPSVYGGAGGHGT4.677876<0.001Q10758Krt8KLEVDPNIQAVR4.663074<0.001D3ZBN0Hist1h1bSLVSKGTLVQTKGTGASGSF4.551461<0.01M0R7B4Hist1h1dSGVSLAALK4.443164<0.01P05197Eef2VNFTVDQIR4.417121<0.01V5QR27NAFPGQPGGPGA4.341792<0.05P63029Tpt1LDYREDGVTPF4.230072<0.05Q5EB49Eno1SFRNPLAK4.118444<0.05Q80ZE7Cd2apGLPAGGIQPHPQTK3.871932<0.01M0R7B4Hist1h1dASGPPVSELITK3.803902<0.001P18420Psma1LVSLIGSK3.745425<0.05P13437Acaa2GVFIVAAK3.727173<0.01D3ZGY4Gapdh-ps2GAAQNIIPASTGAAK3.718082<0.001Q10758Krt8SLGGFGGAGVGGIT3.670916<0.001D3ZBN0Hist1h1bALAAGGYDVEKNNSR3.654154<0.05Q9WVK7HadhTGEGFYKYK3.578854<0.01K7S2S2Hist2h2aa3AGLQFPVGRVH3.386712<0.05A0A0G2K6S9Myh11STVAALEAK3.278155<0.01P06302PtmaSDAAVDTSSEITTK3.185452<0.05G3V779Lad1TEVLVTPAGVASK3.18017<0.01B2GVB1S100a6LIYNEALK3.138931<0.001K7S2S2Hist2h2aa3AEILELAGNAAR3.121051<0.01P69897Tubb5YNEATGGKYVPR3.110143<0.05D3Z7Y6Krt20SLSSSSQGPALSTSGSL3.061158<0.05P63018Hspa8SKGPAVGIDLGTT3.057092<0.05A0A0G2JWK7TaglnEFTDSQLQEGK3.027585<0.05P62630Eef1a1QTVAVGVIK3.020094<0.05P21913SdhbAQTAAAAAPRIK2.929567<0.01Q5RJK6Inpp1LIQSLGPLKT2.899027<0.01P10860Glud1VVQGFGNVGLH2.760697<0.05P47875Csrp1GFGFGQGAGALVH2.639933<0.001Q63191Mamdc4ADQVTLPESITSNP2.631403<0.01Q07936Anxa2STVHEILCKL2.621438<0.001P06302PtmaAAVDTSSEITTK2.568546<0.01Q3MHS9Cct6aAAVKTLNPKAEVAR2.554117<0.05P38552Lgals4YPSAGYNPPQ2.508311<0.05P05197Eef2VAVEAKNPADLPK2.48379<0.01Q5RKI0Wdr1SVADGYSENNVF2.474195<0.001G3V8B3Hist1h2bqHAVSEGTKAVTKYTSSK2.468693<0.05Q6URK4Hnrnpa3YGGGGNYNDFGN2.394012<0.05Q6VPP3Clca4DTQPVLENFSR2.382552<0.05P12346TfLLEACTFHKS2.333879<0.05P63029Tpt1GAIDDSLIGGN2.312302<0.001Q10758Krt8GGFGGAGVGGIT2.276436<0.05Q07936Anxa2STVHEILCK2.247658<0.01P00731Cpa1ALSTDSFNYATYH2.173098<0.05Q6URK4Hnrnpa3SGSPYGGGYGSGGGSGGYGS2.128543<0.05P38552Lgals4YPSAGYNPPQMN2.11776<0.001Q4FZY0Efhd2ATDELASKLSR2.05165<0.05A0A0G2QC04Pls1EGITAIGGTSSI2.004625<0.01P38552Lgals4PAYPSAGYNPPQMN1.905471<0.001Q9Z144Lgals2SEKFEVTNLNMK1.899939<0.05D3ZVN3Fbxo10HNAEAGVDIR1.77112<0.05G3V779Lad1QQVVGAVQAPGQEKVE1.722176<0.01P46462VcpAVANETGAF1.668905<0.001D4A3Z8Tmco3LRIRPTQSV1.649045<0.05Downregulated peptidesA0A0G2K890EzrFVIKPIDKK-1.50403<0.05Q9JJ19Slc9a3r1EALVEPASESPRPA-1.52587<0.05K7S2S2Hist2h2aa3AQGGVLPNIQAV-1.53578<0.05P15865Hist1h1eAPAAPAAPAPAEKTP-1.57197<0.001F1M4J0RictorSTELLLGV-1.57785<0.05A0A0G2JUA5AhnakDVDVQGPDWH-1.61943<0.01Q63191Mamdc4NADQVTLPESITSNP-1.64595<0.05Q6P9V9Tuba1bDLEPTVIDEVR-1.66594<0.05B5DFA0Vil1EENQVITPRLF-1.69645<0.05A0A0G2K3K2ActbNELRVAPEEHPV-1.70784<0.05Q63191Mamdc4QVTLPESITSNP-1.7079<0.05D3ZDZ1Apbb1ipLPPPPPPP-1.72276<0.01Q5EB49Eno1KLNVVEQEKIDQLMIEMDGTENK-1.7416<0.05P62630Eef1a1GVGEFEAGISKN-1.78783<0.01P63269Actg2FAGDDAPRAVFPS-1.81415<0.001D3ZBN0Hist1h1bGTGASGSFKLN-1.85091<0.001D3ZAU0Muc5bITFCEGSCSGK-1.90446<0.05P06302PtmaSDAAVDTSSEITTKDL-1.93122<0.05D3ZX87LOC100910017KNLQTVNVDEN-1.93203<0.05D3ZUM4Glb1LELCTVEFVDTPVIG-1.97941<0.01Q6P9V9Tuba1bQPPTVVPGGDLAK-2.01585<0.05P63312Tmsb10ADKPDMGEIASFD-2.02372<0.001D4A269NAADEIAKAQVA-2.02826<0.001A0A0G2K890EzrAPPPPPPPVYEPVN-2.10935<0.01F1LQ48HnrnplDEGYGPPPPHYE-2.13291<0.01P08081CltaAPAGAPGGPALGN-2.17074<0.001P06761Hspa5MKETAEAYLGKK-2.21256<0.01Q62803Spam1LEDDLVNTIGEIV-2.21478<0.05P38552Lgals4IPAYPSAGYNPPQ-2.21577<0.01A0A0G2K6S9Myh11AQKGQLSDDEKF-2.22383<0.05P13437Acaa2SACIGGGQGISLIIQNTA-2.23257<0.01P30009MarcksSPEAPPAPVAE-2.24701<0.001P45592Cfl1ASGVAVSDGVIK-2.24734<0.001O55212NADAGAAGGPGGPGGPGLG-2.30987<0.05D4AE06Fkbp15TPSVQPSLQPSHP-2.32989<0.05P62630Eef1a1SGDAAIVDMVPGKP-2.36181<0.001P63269Actg2FAGDDAPRAVFPSIV-2.37939<0.05P62630Eef1a1FAPVNVTTEVK-2.45167<0.01P30009MarcksCSPEAPPAPVAE-2.54825<0.01A0A0G2K3K2ActbDSGDGVTHTVPIYE-2.55819<0.01D3ZFU9MylkAQPVGALKPIGN-2.77191<0.05P15865Hist1h1eSETAPAAPAAPAPAEKTPIK-2.77847<0.001Q53Z83Mt1SCCSCCPVGCSKCAQGCVC-2.79709<0.05D3ZDZ1Apbb1ipNFLPPPPPPP-2.8025<0.001P63312Tmsb10ADKPDMGEIAS-2.84338<0.001Q6URK4Hnrnpa3SPYGGGYGSGGGSGGY-2.93019<0.01Q10758Krt8GFGGAGVGGIT-2.95692<0.01P62630Eef1a1VGEFEAGISKN-2.98559<0.001D3ZFU9MylkAKPAETLKPVGN-2.98675<0.01A0A0G2K890EzrAPPPPPPPVYE-3.14377<0.001Q91W32Casp1KEGGAFPGPSGSLK-3.15783<0.001Q53Z83Mt1SCCSCCPVGCSKCAQGCVCKG-3.27655<0.01P06302PtmaSDAAVDTSSEITTKDLK-3.51799<0.001B6ID08Mt2ASCCSCCPVGCAKCSQGCICK-3.56818<0.05Q6P9V9Tuba1bLEPTVIDEVR-3.6423<0.001P62630Eef1a1AGVGEFEAGISKN-3.79477<0.01D3ZJF8FcgbpCVAESTAVCRAQGD-3.88626<0.001A0A0G2K890EzrAELSSEGILDDRN-3.89919<0.05A0A0G2JUA5AhnakEGPEVDVSLP-4.05973<0.01Q53Z83Mt1CAQGCVCKGASDKCTCCA-4.23021<0.01Q63041A1mAVDQSVLLLKPEAK-4.3155<0.001D3ZBN0Hist1h1bGPPVSEL-4.64319<0.001P13668Stmn1ASSDIQVKELE-4.88393<0.01P97881Muc13GVPSQTSNPYAN-5.02829<0.001Q63041A1mALVQKDTVVKPVIVEPEG-5.23966<0.01Q6P9V9Tuba1bTVVPGGDL-5.71525<0.001P15865Hist1h1eSETAPAAPAAPAPAEKTP-6.05552<0.001Q7M0E3DstnFLWAPEQAPLK-6.43258<0.01P62630Eef1a1ALPGDNVGFN-6.78263<0.001Q5PQK2FusYPTQPGQGYSQQSN-8.52323<0.05Q10758Krt8FGGAGVGGIT-8.58147<0.001P62630Eef1a1APVNVTTEVK-13.4474<0.001P63269Actg2GQKDSYVGDEAQSKRG-15.8726<0.001Figure 2
Differentially expressed intestinal peptides between Ctrl and NEC groups: (a) hierarchical cluster analysis of the differentially expressed peptides; (b) number of upregulated or downregulated peptides. The red color represents upregulation, and the green color represents downregulation; a darker color indicates a higher fold change.
(a)(b)


### 3.3. Characteristics of the Differentially Expressed Peptides

Some specific characteristics of the peptides were detected in this study. Among these characteristics, the amino acid numbers are mainly focused on 9–14 (Figure3(a)), the peptides range largely from 800–1800 (Figure 3(b)), and the pI values are evenly distributed on both sides of acid and basic (Figure 3(c)). In addition, the distribution of MW relative to pI (MW/pI) was investigated based on the amino acid composition and MW distribution specific to a pI (Figure 3(d)). As the generation of peptides mainly depends upon the enzymatic cleavage of precursor proteins [6], the cleavage site patterns of the identified peptides were analyzed (Figures 3(e) and 3(f)). Among the upregulated peptides, lysine (K), arginine (R), alanine (A), and serine (S) were the four dominant amino acids, whereas among the downregulated peptides, alanine (A), lysine (K), asparagine (N), and glutamic acid (E) were the four dominant amino acids. In particular, we found that multiple peptides are derived from same precursor proteins. Among these, the precursor proteins Krt8, Eef1a1, Hist1h1b, and NA generated a total of 31 peptides (Figure 3(g)).Figure 3
Characteristics of the differentially expressed peptides: (a) amino acid number, (b) molecular weights (MWs), and (c) isoelectric points (pI) of the peptides. (d) Scatter plot of MW versus pI. (e) The distribution diagram of the N- and C-terminal cleavage sites of the upregulated peptides. (f) The distribution diagram of the N- and C-terminal cleavage sites of the downregulated peptides. (g) The number of peptides produced by the same precursor proteins.
(a)(b)(c)(d)(e)(f)(g)


### 3.4. GO and KEGG Pathway Analyses of the Differently Expressed Peptides

GO and KEGG pathway analyses were performed to determine the potential roles of these differentially expressed peptides and their corresponding precursor proteins. For biological processes, cellular developmental process, response to oxygen levels, regulation of cell development, NF-κB import into the nucleus, regeneration, and Wnt-activated receptor activity had a high enrichment score (Figure 4(a)). Extracellular exosome, myelin sheath, intracellular ribonucleoprotein complex, cell–cell adherens junction, membrane raft, and brush border as cellular components were found to be closely related to the identified peptides (Figure 4(b)). Regarding molecular function, Wnt-activated receptor activity, RNA polymerase II transcription factor binding, peptide hormone binding, peptide receptor activity, and calcium-dependent protein binding were identified as the most enriched subcategories (Figure 4(c)). Further, KEGG pathway analysis indicated that the potential function of these differentially expressed peptides was related to the TGF-β signaling pathway, riboflavin metabolism, axon guidance, apoptosis, the Wnt signaling pathway, and the MAPK signaling pathway (Figure 4(d)).Figure 4
GO and KEGG pathway analyses of the differentially expressed peptides: (a) the biological process categories; (b) the cellular component categories; (c) the molecular function categories; (d) the KEGG signaling pathways.
(a)(b)(c)(d)


### 3.5. Bioactive Peptides Potentially Related to NEC

The functions of the precursor proteins largely determine the functions of the bioactive peptides [14]. We used the UniProt and STRING databases to analyze the protein interactions of the differentially expressed peptides and the protein–protein interaction of their protein precursors. It is well known that bioactive peptides can be derived from their precursor domain. Hence, we used the UniProt database to analyze peptides derived from their precursor protein functional domain (Table 2). Among the downregulated peptides, the protein precursor Ezrin (Ezr) plays a vital role in the apical integrity of the intestinal epithelium, and its severe defects can lead to incomplete villi morphogenesis and neonatal death [15]. In addition, Metallothio (MT1) has a protective effect on intestinal mucosal inflammation. Therefore, we attempted to detect bioactive peptides with potential NEC protective effects from the downregulated peptides.Table 2
Differential peptides located in functional domains.
GeneAccessionPeptidePeptide domainp valueHist1h2bqG3V8B3AVTKYTSSKHistone H2B<0.05Tubb5P69897FVFGQSGAGNNTubulin/FtsZ family, GTPase domain<0.01Mdh1O88989VIVVGNPANTNCLTASKLdh_1_N<0.05AldobQ66HT1ALQASALAAWGGKGlycolytic<0.001Hist1h2bqG3V8B3LLLPGELAKHHistone H2B<0.001Eef2P05197VFSGVVSTGLKVRGTP_EFTU_D2<0.05Hnrnpa2b1F1LNF1GFGFVTFRNA recognition motif<0.05Hist1h1dM0R7B4SGVSLAALKKDomain in histone families 1 and 5<0.001Tubb5P69897FVFGQSGAGNTubulin/FtsZ family, GTPase domain<0.05Gapdh-ps2D3ZGY4GAAQNIIPASTGAAKAVGKGp_dh_C<0.001Hist1h1bD3ZBN0SLVSKGTLVQTKDomain in histone families 1 and 5<0.01Hist1h1dM0R7B4GILVQTKGTGASGSDomain in histone families 1 and 5<0.001NAP02262VGAGAPVYLHistone 2A<0.001Atp5bP10719IGLFGGAGVGKATPases associated with a variety of cellular activities<0.01Gapdh-ps2D3ZGY4VKVGVNGFGRGlyceraldehyde 3-phosphate dehydrogenase, NAD binding domain<0.001Hist1h2bqG3V8B3LLLPGELAKHistone H2B<0.01Aldh18a1D3ZIE9VGLGGMEAKVKAA_kinase<0.001Actn4Q9QXQ0STALYGESDLCa2+ insensitive EF hand<0.001Krt8Q10758GSLGGFGGAGVGGITKeratin_2_head<0.05AldoaP05065ALQASALKGlycolytic<0.01Plin3M0RA08VTGAVDVTCGAVKPerilipin<0.001NAA0A0G2K3Z9SSGNAKIGHPAPSFK<0.01Gapdh-ps2D3ZGY4ALNDNFVKGp_dh_C<0.05Rps23P62268KANPFGGASHAKRibosom_S12_S23<0.001Krt8Q10758KLEVDPNIQAVRKeratin_2_head<0.001Hist1h1bD3ZBN0SLVSKGTLVQTKGTGASGSFDomain in histone families 1 and 5<0.01Hist1h1dM0R7B4SGVSLAALKDomain in histone families 1 and 5<0.01Tpt1P63029LDYREDGVTPFTCTP<0.05Eno1Q5EB49SFRNPLAKEnolase, C-terminal TIM barrel domain<0.05Psma1P18420LVSLIGSKProteasome<0.05Acaa2P13437GVFIVAAKThiolase_N<0.01Gapdh-ps2D3ZGY4GAAQNIIPASTGAAKGp_dh_C<0.001Krt8Q10758SLGGFGGAGVGGITKeratin_2_head<0.001HadhQ9WVK7TGEGFYKYK3HCDH<0.01Myh11A0A0G2K6S9STVAALEAK<0.01PtmaP06302SDAAVDTSSEITTKProthymosin<0.05Hist2h2aa3K7S2S2AEILELAGNAARHistone 2A<0.01Tubb5P69897YNEATGGKYVPRTubulin/FtsZ family, GTPase domain<0.05Hspa8P63018SKGPAVGIDLGTTHSP70<0.05Eef1a1P62630QTVAVGVIKGTP_EFTU_D3<0.05Glud1P10860VVQGFGNVGLHGlutamate/leucine/phenylalanine/valine dehydrogenase<0.05PtmaP06302AAVDTSSEITTKProthymosin<0.01Eef2P05197VAVEAKNPADLPKEFG_II<0.01Wdr1Q5RKI0SVADGYSENNVFWD40 repeats<0.001Hist1h2bqG3V8B3HAVSEGTKAVTKYTSSKHistone H2B<0.05Tpt1P63029GAIDDSLIGGNTCTP<0.001Krt8Q10758GGFGGAGVGGITKeratin_2_head<0.05Efhd2Q4FZY0ATDELASKLSRMetallothio<0.05Pls1A0A0G2QC04EGITAIGGTSSI<0.01Lgals2Q9Z144SEKFEVTNLNMKGalactoside-binding lectin<0.05Fbxo10D3ZVN3HNAEAGVDIRParallel beta-helix repeats<0.05VcpP46462AVANETGAFATPases associated with a variety of cellular activities<0.001EzrA0A0G2K890FVIKPIDKK<0.05AhnakA0A0G2JUA5DVDVQGPDWH<0.01Tuba1bQ6P9V9DLEPTVIDEVRTubulin/FtsZ family, GTPase domain<0.05Eef1a1P62630GVGEFEAGISKNGTP_EFTU<0.01PtmaP06302SDAAVDTSSEITTKDLProthymosin<0.05Tmsb10P63312ADKPDMGEIASFDThymosin beta actin-binding motif.<0.001Hspa5P06761MKETAEAYLGKKHSP 70<0.01Spam1Q62803LEDDLVNTIGEIVGlyco_hydro_56<0.05Myh11A0A0G2K6S9AQKGQLSDDEKF<0.05Acaa2P13437SACIGGGQGISLIIQNTAThiolase_C<0.01Eef1a1P62630SGDAAIVDMVPGKPGTP_EFTU_D3<0.001Actg2P63269FAGDDAPRAVFPSIVActin<0.05ActbA0A0G2K3K2DSGDGVTHTVPIYE<0.01Tmsb10P63312ADKPDMGEIASThymosin beta actin-binding motif.<0.001Krt8Q10758GFGGAGVGGITKeratin_2_head<0.01Eef1a1P62630VGEFEAGISKNGTP_EFTU<0.001PtmaP06302SDAAVDTSSEITTKDLKProthymosin<0.001Tuba1bQ6P9V9LEPTVIDEVRTubulin/FtsZ family, GTPase domain<0.001Eef1a1P62630AGVGEFEAGISKNGTP_EFTU<0.01FcgbpD3ZJF8CVAESTAVCRAQGDvon Willebrand factor (vWF) type D domain<0.001AhnakA0A0G2JUA5EGPEVDVSLP<0.01Mt1Q53Z83CAQGCVCKGASDKCTCCAMetallothio<0.01A1mQ63041AVDQSVLLLKPEAKAlpha-2-macroglobulin<0.001Stmn1P13668ASSDIQVKELEStathmin<0.01Eef1a1P62630ALPGDNVGFNGTP_EFTU_D2<0.001Krt8Q10758FGGAGVGGITKeratin_2_head<0.001


### 3.6. Protective Effects of EDP1 and MDP2

To further analyze the function of the differentially expressed peptides, we randomly selected downregulated peptide 1 derived from EZRI (EDP1) and downregulated peptide 2 derived from MT1 (MDP2) to investigate its NEC protective ability. Interestingly, EDP1 treatment could partially reverse the inhibition of cell migration caused by LPS (Figures5(a) and 5(b)). In addition, MDP2 treatment not only reversed the inhibition of cell migration in the NEC cell model but also decreased the apoptosis level (Figures 5(a) and 5(b)).Figure 5
Protective effects of EDP1 and MDP2: (a) the effect of EDP1 and MDP2 on the NEC in vitro model wound restitution; (b) the effect of EDP1 and MDP2 on the NEC in vitro model apoptosis ratio.
(a)(b)


## 3.1. Establishment of the NEC Model

Hypoxia-cold stimulation combined with hypertonic milk gavage was used to establish the NEC rat model. An optical microscope was used to observe stained intestinal tissue sections and evaluate the pathological score. The histopathological results showed shedding of the intestinal villi, submucosal edema, and inflammatory cell infiltration in the NEC group. In contrast, the intestinal mucosa epithelium was intact, the glands were arranged regularly, and edema and inflammatory cell infiltration were not evident in all layers of the intestinal wall in the Ctrl group (Figure1(a)). In the NEC rat pups, the pathological score [13] was significantly higher than that in the Ctrl rat pups (Figure 1(b)). These findings revealed that the NEC model was successfully established.Figure 1
Establishment of the NEC model: (a) H&E staining of representative intestinal tissue sections from the Ctrl and NEC groups (bar=100μm); (b) pathological scores of representative intestinal tissue sections from the Ctrl and NEC groups (∗∗∗p<0.001).
(a)(b)


## 3.2. Differentially Expressed Peptides in Intestinal Tissue Samples

TMT-labeled LC-MS/MS analysis was performed on peptides in the intestinal tissue samples from both the Ctrl and NEC groups. A total of 480 nonredundant peptides were identified, with 176 differentially expressed peptides derived from 91 protein precursors (foldchange≥1.5, p<0.05). Detailed information on these peptides is shown in Table 1. Upregulated or downregulated differentially expressed peptides were visualized in a hierarchical clustering heat map (Figure 2(a)). The number of upregulated or downregulated differentially expressed peptides is shown in Figure 2(b).Table 1
Differential peptides in the intestine tissue of NEC rat model.
AccessionGenePeptideFold changep valueUpregulated peptidesG3V8B3Hist1h2bqAVTKYTSSK57.72612<0.05P69897Tubb5FVFGQSGAGNN31.84287<0.01Q10758Krt8DGKLVSESSDIMSK25.55914<0.001Q10758Krt8LVSESSDIMSK22.05736<0.001P15999Atp5a1VGLKAPGIIPR21.49717<0.001O88989Mdh1VIVVGNPANTNCLTASK12.14385<0.05Q66HT1AldobALQASALAAWGGK11.93121<0.001G3V8B3Hist1h2bqLLLPGELAKH11.86002<0.001P04182OatALQDPNVAAF11.16052<0.05P05197Eef2VFSGVVSTGLKVR10.57569<0.05Q10758Krt8QPGFGSVGGSST10.57101<0.01F1LNF1Hnrnpa2b1GFGFVTF10.55872<0.05K7S2S2Hist2h2aa3AGLQFPVGR10.3368<0.05Q7M0E3DstnFLWAPEQAPLKSK10.01604<0.01Q10758Krt8GGLTSPGFS9.126557<0.001M0R7B4Hist1h1dSGVSLAALKK8.945185<0.001Q10758Krt8FQPGFGSVGGSST8.793551<0.001P69897Tubb5FVFGQSGAGN8.771316<0.05D3ZGY4Gapdh-ps2GAAQNIIPASTGAAKAVGK8.399639<0.001D3ZBN0Hist1h1bSLVSKGTLVQTK8.352781<0.01M0R7B4Hist1h1dGILVQTKGTGASGS8.334644<0.001P02262NAVGAGAPVYL8.301511<0.001Q10758Krt8GMSSFQPGFGSVGGSST8.18755<0.01P10719Atp5bIGLFGGAGVGK7.935978<0.01B5DEM5Rpl14AAIAAAAAAK7.751116<0.01K7S2S2Hist2h2aa3GLQFPVGR7.629358<0.05P62630Eef1a1ESFSDYPPLGR7.607978<0.001O55159EpcamQKDCVCNNYKLTSR7.485879<0.001P15865Hist1h1eSETAPAAPAAPAPAEKTPIKKK7.38236<0.01D3ZGY4Gapdh-ps2VKVGVNGFGR7.332894<0.001F1M577LOC100359671AFIAHPKLGK6.954795<0.05Q6P725DesTFGGAPGFSLGSPLS6.720618<0.05G3V8B3Hist1h2bqLLLPGELAK6.595403<0.01D3ZIE9Aldh18a1VGLGGMEAKVK6.321181<0.001Q925G0Rbm3DYSGSQGGYDRYSGGN6.259407<0.01Q10758Krt8SFQPGFGSVGGSST6.237045<0.01Q10758Krt8SSFQPGFGSVGGSST6.119025<0.01Q9QXQ0Actn4STALYGESDL5.964186<0.001Q10758Krt8GSLGGFGGAGVGGIT5.600793<0.05P05065AldoaALQASALK5.474993<0.01M0RA08Plin3VTGAVDVTCGAVK5.466526<0.001Q6P9V9Tuba1bVAEITNACFEPAN5.315759<0.05A0A0G2K3Z9NASSGNAKIGHPAPSFK4.969645<0.01Q66H80Arcn1VLLAAAVCTK4.865557<0.05D3ZGY4Gapdh-ps2ALNDNFVK4.838244<0.05B0K031Rpl7IALTDNSLVAR4.834916<0.001P62268Rps23KANPFGGASHAK4.801351<0.001D3Z7Y6Krt20LGVAPSVYGGAGGHGT4.677876<0.001Q10758Krt8KLEVDPNIQAVR4.663074<0.001D3ZBN0Hist1h1bSLVSKGTLVQTKGTGASGSF4.551461<0.01M0R7B4Hist1h1dSGVSLAALK4.443164<0.01P05197Eef2VNFTVDQIR4.417121<0.01V5QR27NAFPGQPGGPGA4.341792<0.05P63029Tpt1LDYREDGVTPF4.230072<0.05Q5EB49Eno1SFRNPLAK4.118444<0.05Q80ZE7Cd2apGLPAGGIQPHPQTK3.871932<0.01M0R7B4Hist1h1dASGPPVSELITK3.803902<0.001P18420Psma1LVSLIGSK3.745425<0.05P13437Acaa2GVFIVAAK3.727173<0.01D3ZGY4Gapdh-ps2GAAQNIIPASTGAAK3.718082<0.001Q10758Krt8SLGGFGGAGVGGIT3.670916<0.001D3ZBN0Hist1h1bALAAGGYDVEKNNSR3.654154<0.05Q9WVK7HadhTGEGFYKYK3.578854<0.01K7S2S2Hist2h2aa3AGLQFPVGRVH3.386712<0.05A0A0G2K6S9Myh11STVAALEAK3.278155<0.01P06302PtmaSDAAVDTSSEITTK3.185452<0.05G3V779Lad1TEVLVTPAGVASK3.18017<0.01B2GVB1S100a6LIYNEALK3.138931<0.001K7S2S2Hist2h2aa3AEILELAGNAAR3.121051<0.01P69897Tubb5YNEATGGKYVPR3.110143<0.05D3Z7Y6Krt20SLSSSSQGPALSTSGSL3.061158<0.05P63018Hspa8SKGPAVGIDLGTT3.057092<0.05A0A0G2JWK7TaglnEFTDSQLQEGK3.027585<0.05P62630Eef1a1QTVAVGVIK3.020094<0.05P21913SdhbAQTAAAAAPRIK2.929567<0.01Q5RJK6Inpp1LIQSLGPLKT2.899027<0.01P10860Glud1VVQGFGNVGLH2.760697<0.05P47875Csrp1GFGFGQGAGALVH2.639933<0.001Q63191Mamdc4ADQVTLPESITSNP2.631403<0.01Q07936Anxa2STVHEILCKL2.621438<0.001P06302PtmaAAVDTSSEITTK2.568546<0.01Q3MHS9Cct6aAAVKTLNPKAEVAR2.554117<0.05P38552Lgals4YPSAGYNPPQ2.508311<0.05P05197Eef2VAVEAKNPADLPK2.48379<0.01Q5RKI0Wdr1SVADGYSENNVF2.474195<0.001G3V8B3Hist1h2bqHAVSEGTKAVTKYTSSK2.468693<0.05Q6URK4Hnrnpa3YGGGGNYNDFGN2.394012<0.05Q6VPP3Clca4DTQPVLENFSR2.382552<0.05P12346TfLLEACTFHKS2.333879<0.05P63029Tpt1GAIDDSLIGGN2.312302<0.001Q10758Krt8GGFGGAGVGGIT2.276436<0.05Q07936Anxa2STVHEILCK2.247658<0.01P00731Cpa1ALSTDSFNYATYH2.173098<0.05Q6URK4Hnrnpa3SGSPYGGGYGSGGGSGGYGS2.128543<0.05P38552Lgals4YPSAGYNPPQMN2.11776<0.001Q4FZY0Efhd2ATDELASKLSR2.05165<0.05A0A0G2QC04Pls1EGITAIGGTSSI2.004625<0.01P38552Lgals4PAYPSAGYNPPQMN1.905471<0.001Q9Z144Lgals2SEKFEVTNLNMK1.899939<0.05D3ZVN3Fbxo10HNAEAGVDIR1.77112<0.05G3V779Lad1QQVVGAVQAPGQEKVE1.722176<0.01P46462VcpAVANETGAF1.668905<0.001D4A3Z8Tmco3LRIRPTQSV1.649045<0.05Downregulated peptidesA0A0G2K890EzrFVIKPIDKK-1.50403<0.05Q9JJ19Slc9a3r1EALVEPASESPRPA-1.52587<0.05K7S2S2Hist2h2aa3AQGGVLPNIQAV-1.53578<0.05P15865Hist1h1eAPAAPAAPAPAEKTP-1.57197<0.001F1M4J0RictorSTELLLGV-1.57785<0.05A0A0G2JUA5AhnakDVDVQGPDWH-1.61943<0.01Q63191Mamdc4NADQVTLPESITSNP-1.64595<0.05Q6P9V9Tuba1bDLEPTVIDEVR-1.66594<0.05B5DFA0Vil1EENQVITPRLF-1.69645<0.05A0A0G2K3K2ActbNELRVAPEEHPV-1.70784<0.05Q63191Mamdc4QVTLPESITSNP-1.7079<0.05D3ZDZ1Apbb1ipLPPPPPPP-1.72276<0.01Q5EB49Eno1KLNVVEQEKIDQLMIEMDGTENK-1.7416<0.05P62630Eef1a1GVGEFEAGISKN-1.78783<0.01P63269Actg2FAGDDAPRAVFPS-1.81415<0.001D3ZBN0Hist1h1bGTGASGSFKLN-1.85091<0.001D3ZAU0Muc5bITFCEGSCSGK-1.90446<0.05P06302PtmaSDAAVDTSSEITTKDL-1.93122<0.05D3ZX87LOC100910017KNLQTVNVDEN-1.93203<0.05D3ZUM4Glb1LELCTVEFVDTPVIG-1.97941<0.01Q6P9V9Tuba1bQPPTVVPGGDLAK-2.01585<0.05P63312Tmsb10ADKPDMGEIASFD-2.02372<0.001D4A269NAADEIAKAQVA-2.02826<0.001A0A0G2K890EzrAPPPPPPPVYEPVN-2.10935<0.01F1LQ48HnrnplDEGYGPPPPHYE-2.13291<0.01P08081CltaAPAGAPGGPALGN-2.17074<0.001P06761Hspa5MKETAEAYLGKK-2.21256<0.01Q62803Spam1LEDDLVNTIGEIV-2.21478<0.05P38552Lgals4IPAYPSAGYNPPQ-2.21577<0.01A0A0G2K6S9Myh11AQKGQLSDDEKF-2.22383<0.05P13437Acaa2SACIGGGQGISLIIQNTA-2.23257<0.01P30009MarcksSPEAPPAPVAE-2.24701<0.001P45592Cfl1ASGVAVSDGVIK-2.24734<0.001O55212NADAGAAGGPGGPGGPGLG-2.30987<0.05D4AE06Fkbp15TPSVQPSLQPSHP-2.32989<0.05P62630Eef1a1SGDAAIVDMVPGKP-2.36181<0.001P63269Actg2FAGDDAPRAVFPSIV-2.37939<0.05P62630Eef1a1FAPVNVTTEVK-2.45167<0.01P30009MarcksCSPEAPPAPVAE-2.54825<0.01A0A0G2K3K2ActbDSGDGVTHTVPIYE-2.55819<0.01D3ZFU9MylkAQPVGALKPIGN-2.77191<0.05P15865Hist1h1eSETAPAAPAAPAPAEKTPIK-2.77847<0.001Q53Z83Mt1SCCSCCPVGCSKCAQGCVC-2.79709<0.05D3ZDZ1Apbb1ipNFLPPPPPPP-2.8025<0.001P63312Tmsb10ADKPDMGEIAS-2.84338<0.001Q6URK4Hnrnpa3SPYGGGYGSGGGSGGY-2.93019<0.01Q10758Krt8GFGGAGVGGIT-2.95692<0.01P62630Eef1a1VGEFEAGISKN-2.98559<0.001D3ZFU9MylkAKPAETLKPVGN-2.98675<0.01A0A0G2K890EzrAPPPPPPPVYE-3.14377<0.001Q91W32Casp1KEGGAFPGPSGSLK-3.15783<0.001Q53Z83Mt1SCCSCCPVGCSKCAQGCVCKG-3.27655<0.01P06302PtmaSDAAVDTSSEITTKDLK-3.51799<0.001B6ID08Mt2ASCCSCCPVGCAKCSQGCICK-3.56818<0.05Q6P9V9Tuba1bLEPTVIDEVR-3.6423<0.001P62630Eef1a1AGVGEFEAGISKN-3.79477<0.01D3ZJF8FcgbpCVAESTAVCRAQGD-3.88626<0.001A0A0G2K890EzrAELSSEGILDDRN-3.89919<0.05A0A0G2JUA5AhnakEGPEVDVSLP-4.05973<0.01Q53Z83Mt1CAQGCVCKGASDKCTCCA-4.23021<0.01Q63041A1mAVDQSVLLLKPEAK-4.3155<0.001D3ZBN0Hist1h1bGPPVSEL-4.64319<0.001P13668Stmn1ASSDIQVKELE-4.88393<0.01P97881Muc13GVPSQTSNPYAN-5.02829<0.001Q63041A1mALVQKDTVVKPVIVEPEG-5.23966<0.01Q6P9V9Tuba1bTVVPGGDL-5.71525<0.001P15865Hist1h1eSETAPAAPAAPAPAEKTP-6.05552<0.001Q7M0E3DstnFLWAPEQAPLK-6.43258<0.01P62630Eef1a1ALPGDNVGFN-6.78263<0.001Q5PQK2FusYPTQPGQGYSQQSN-8.52323<0.05Q10758Krt8FGGAGVGGIT-8.58147<0.001P62630Eef1a1APVNVTTEVK-13.4474<0.001P63269Actg2GQKDSYVGDEAQSKRG-15.8726<0.001Figure 2
Differentially expressed intestinal peptides between Ctrl and NEC groups: (a) hierarchical cluster analysis of the differentially expressed peptides; (b) number of upregulated or downregulated peptides. The red color represents upregulation, and the green color represents downregulation; a darker color indicates a higher fold change.
(a)(b)


## 3.3. Characteristics of the Differentially Expressed Peptides

Some specific characteristics of the peptides were detected in this study. Among these characteristics, the amino acid numbers are mainly focused on 9–14 (Figure3(a)), the peptides range largely from 800–1800 (Figure 3(b)), and the pI values are evenly distributed on both sides of acid and basic (Figure 3(c)). In addition, the distribution of MW relative to pI (MW/pI) was investigated based on the amino acid composition and MW distribution specific to a pI (Figure 3(d)). As the generation of peptides mainly depends upon the enzymatic cleavage of precursor proteins [6], the cleavage site patterns of the identified peptides were analyzed (Figures 3(e) and 3(f)). Among the upregulated peptides, lysine (K), arginine (R), alanine (A), and serine (S) were the four dominant amino acids, whereas among the downregulated peptides, alanine (A), lysine (K), asparagine (N), and glutamic acid (E) were the four dominant amino acids. In particular, we found that multiple peptides are derived from same precursor proteins. Among these, the precursor proteins Krt8, Eef1a1, Hist1h1b, and NA generated a total of 31 peptides (Figure 3(g)).Figure 3
Characteristics of the differentially expressed peptides: (a) amino acid number, (b) molecular weights (MWs), and (c) isoelectric points (pI) of the peptides. (d) Scatter plot of MW versus pI. (e) The distribution diagram of the N- and C-terminal cleavage sites of the upregulated peptides. (f) The distribution diagram of the N- and C-terminal cleavage sites of the downregulated peptides. (g) The number of peptides produced by the same precursor proteins.
(a)(b)(c)(d)(e)(f)(g)


## 3.4. GO and KEGG Pathway Analyses of the Differently Expressed Peptides

GO and KEGG pathway analyses were performed to determine the potential roles of these differentially expressed peptides and their corresponding precursor proteins. For biological processes, cellular developmental process, response to oxygen levels, regulation of cell development, NF-κB import into the nucleus, regeneration, and Wnt-activated receptor activity had a high enrichment score (Figure 4(a)). Extracellular exosome, myelin sheath, intracellular ribonucleoprotein complex, cell–cell adherens junction, membrane raft, and brush border as cellular components were found to be closely related to the identified peptides (Figure 4(b)). Regarding molecular function, Wnt-activated receptor activity, RNA polymerase II transcription factor binding, peptide hormone binding, peptide receptor activity, and calcium-dependent protein binding were identified as the most enriched subcategories (Figure 4(c)). Further, KEGG pathway analysis indicated that the potential function of these differentially expressed peptides was related to the TGF-β signaling pathway, riboflavin metabolism, axon guidance, apoptosis, the Wnt signaling pathway, and the MAPK signaling pathway (Figure 4(d)).Figure 4
GO and KEGG pathway analyses of the differentially expressed peptides: (a) the biological process categories; (b) the cellular component categories; (c) the molecular function categories; (d) the KEGG signaling pathways.
(a)(b)(c)(d)


## 3.5. Bioactive Peptides Potentially Related to NEC

The functions of the precursor proteins largely determine the functions of the bioactive peptides [14]. We used the UniProt and STRING databases to analyze the protein interactions of the differentially expressed peptides and the protein–protein interaction of their protein precursors. It is well known that bioactive peptides can be derived from their precursor domain. Hence, we used the UniProt database to analyze peptides derived from their precursor protein functional domain (Table 2). Among the downregulated peptides, the protein precursor Ezrin (Ezr) plays a vital role in the apical integrity of the intestinal epithelium, and its severe defects can lead to incomplete villi morphogenesis and neonatal death [15]. In addition, Metallothio (MT1) has a protective effect on intestinal mucosal inflammation. Therefore, we attempted to detect bioactive peptides with potential NEC protective effects from the downregulated peptides.Table 2
Differential peptides located in functional domains.
GeneAccessionPeptidePeptide domainp valueHist1h2bqG3V8B3AVTKYTSSKHistone H2B<0.05Tubb5P69897FVFGQSGAGNNTubulin/FtsZ family, GTPase domain<0.01Mdh1O88989VIVVGNPANTNCLTASKLdh_1_N<0.05AldobQ66HT1ALQASALAAWGGKGlycolytic<0.001Hist1h2bqG3V8B3LLLPGELAKHHistone H2B<0.001Eef2P05197VFSGVVSTGLKVRGTP_EFTU_D2<0.05Hnrnpa2b1F1LNF1GFGFVTFRNA recognition motif<0.05Hist1h1dM0R7B4SGVSLAALKKDomain in histone families 1 and 5<0.001Tubb5P69897FVFGQSGAGNTubulin/FtsZ family, GTPase domain<0.05Gapdh-ps2D3ZGY4GAAQNIIPASTGAAKAVGKGp_dh_C<0.001Hist1h1bD3ZBN0SLVSKGTLVQTKDomain in histone families 1 and 5<0.01Hist1h1dM0R7B4GILVQTKGTGASGSDomain in histone families 1 and 5<0.001NAP02262VGAGAPVYLHistone 2A<0.001Atp5bP10719IGLFGGAGVGKATPases associated with a variety of cellular activities<0.01Gapdh-ps2D3ZGY4VKVGVNGFGRGlyceraldehyde 3-phosphate dehydrogenase, NAD binding domain<0.001Hist1h2bqG3V8B3LLLPGELAKHistone H2B<0.01Aldh18a1D3ZIE9VGLGGMEAKVKAA_kinase<0.001Actn4Q9QXQ0STALYGESDLCa2+ insensitive EF hand<0.001Krt8Q10758GSLGGFGGAGVGGITKeratin_2_head<0.05AldoaP05065ALQASALKGlycolytic<0.01Plin3M0RA08VTGAVDVTCGAVKPerilipin<0.001NAA0A0G2K3Z9SSGNAKIGHPAPSFK<0.01Gapdh-ps2D3ZGY4ALNDNFVKGp_dh_C<0.05Rps23P62268KANPFGGASHAKRibosom_S12_S23<0.001Krt8Q10758KLEVDPNIQAVRKeratin_2_head<0.001Hist1h1bD3ZBN0SLVSKGTLVQTKGTGASGSFDomain in histone families 1 and 5<0.01Hist1h1dM0R7B4SGVSLAALKDomain in histone families 1 and 5<0.01Tpt1P63029LDYREDGVTPFTCTP<0.05Eno1Q5EB49SFRNPLAKEnolase, C-terminal TIM barrel domain<0.05Psma1P18420LVSLIGSKProteasome<0.05Acaa2P13437GVFIVAAKThiolase_N<0.01Gapdh-ps2D3ZGY4GAAQNIIPASTGAAKGp_dh_C<0.001Krt8Q10758SLGGFGGAGVGGITKeratin_2_head<0.001HadhQ9WVK7TGEGFYKYK3HCDH<0.01Myh11A0A0G2K6S9STVAALEAK<0.01PtmaP06302SDAAVDTSSEITTKProthymosin<0.05Hist2h2aa3K7S2S2AEILELAGNAARHistone 2A<0.01Tubb5P69897YNEATGGKYVPRTubulin/FtsZ family, GTPase domain<0.05Hspa8P63018SKGPAVGIDLGTTHSP70<0.05Eef1a1P62630QTVAVGVIKGTP_EFTU_D3<0.05Glud1P10860VVQGFGNVGLHGlutamate/leucine/phenylalanine/valine dehydrogenase<0.05PtmaP06302AAVDTSSEITTKProthymosin<0.01Eef2P05197VAVEAKNPADLPKEFG_II<0.01Wdr1Q5RKI0SVADGYSENNVFWD40 repeats<0.001Hist1h2bqG3V8B3HAVSEGTKAVTKYTSSKHistone H2B<0.05Tpt1P63029GAIDDSLIGGNTCTP<0.001Krt8Q10758GGFGGAGVGGITKeratin_2_head<0.05Efhd2Q4FZY0ATDELASKLSRMetallothio<0.05Pls1A0A0G2QC04EGITAIGGTSSI<0.01Lgals2Q9Z144SEKFEVTNLNMKGalactoside-binding lectin<0.05Fbxo10D3ZVN3HNAEAGVDIRParallel beta-helix repeats<0.05VcpP46462AVANETGAFATPases associated with a variety of cellular activities<0.001EzrA0A0G2K890FVIKPIDKK<0.05AhnakA0A0G2JUA5DVDVQGPDWH<0.01Tuba1bQ6P9V9DLEPTVIDEVRTubulin/FtsZ family, GTPase domain<0.05Eef1a1P62630GVGEFEAGISKNGTP_EFTU<0.01PtmaP06302SDAAVDTSSEITTKDLProthymosin<0.05Tmsb10P63312ADKPDMGEIASFDThymosin beta actin-binding motif.<0.001Hspa5P06761MKETAEAYLGKKHSP 70<0.01Spam1Q62803LEDDLVNTIGEIVGlyco_hydro_56<0.05Myh11A0A0G2K6S9AQKGQLSDDEKF<0.05Acaa2P13437SACIGGGQGISLIIQNTAThiolase_C<0.01Eef1a1P62630SGDAAIVDMVPGKPGTP_EFTU_D3<0.001Actg2P63269FAGDDAPRAVFPSIVActin<0.05ActbA0A0G2K3K2DSGDGVTHTVPIYE<0.01Tmsb10P63312ADKPDMGEIASThymosin beta actin-binding motif.<0.001Krt8Q10758GFGGAGVGGITKeratin_2_head<0.01Eef1a1P62630VGEFEAGISKNGTP_EFTU<0.001PtmaP06302SDAAVDTSSEITTKDLKProthymosin<0.001Tuba1bQ6P9V9LEPTVIDEVRTubulin/FtsZ family, GTPase domain<0.001Eef1a1P62630AGVGEFEAGISKNGTP_EFTU<0.01FcgbpD3ZJF8CVAESTAVCRAQGDvon Willebrand factor (vWF) type D domain<0.001AhnakA0A0G2JUA5EGPEVDVSLP<0.01Mt1Q53Z83CAQGCVCKGASDKCTCCAMetallothio<0.01A1mQ63041AVDQSVLLLKPEAKAlpha-2-macroglobulin<0.001Stmn1P13668ASSDIQVKELEStathmin<0.01Eef1a1P62630ALPGDNVGFNGTP_EFTU_D2<0.001Krt8Q10758FGGAGVGGITKeratin_2_head<0.001


## 3.6. Protective Effects of EDP1 and MDP2

To further analyze the function of the differentially expressed peptides, we randomly selected downregulated peptide 1 derived from EZRI (EDP1) and downregulated peptide 2 derived from MT1 (MDP2) to investigate its NEC protective ability. Interestingly, EDP1 treatment could partially reverse the inhibition of cell migration caused by LPS (Figures5(a) and 5(b)). In addition, MDP2 treatment not only reversed the inhibition of cell migration in the NEC cell model but also decreased the apoptosis level (Figures 5(a) and 5(b)).Figure 5
Protective effects of EDP1 and MDP2: (a) the effect of EDP1 and MDP2 on the NEC in vitro model wound restitution; (b) the effect of EDP1 and MDP2 on the NEC in vitro model apoptosis ratio.
(a)(b)


## 4. Discussion

NEC is the most common gastrointestinal complication in premature infants [16]. Although extensive research has been conducted to investigate this disease, the underlying molecular mechanisms have not yet been fully elucidated. Peptidomics is an innovative branch of proteomics, which has been used to explore the mechanisms of various diseases [17, 18]. In our study, we detected differentially expressed peptides in intestinal tissues that were further analyzed for their potential association with NEC. Moreover, the bioactive peptides EDP1 and MDP2 showed a protective effect in the NEC model in vitro.In total, 176 differentially expressed peptides, including 103 upregulated and 73 downregulated, were identified. GO analysis of precursor proteins revealed that the differentially expressed peptides were mainly involved in the response to oxygen levels, NF-κB import into the nucleus, extracellular exosome, intracellular ribonucleoprotein complex, cell–cell adherens junction, peptide hormone binding, and peptide receptor activity. Among these, accumulation of oxygen free radicals and low antioxidant capacity result in intestinal epithelial cell apoptosis and intestinal inflammation, which are the known causes of NEC [19]. Moreover, the aberrant activation of NF-κB has been shown to be associated with many inflammatory bowel diseases, including NEC, and the inhibition of NF-κB can reduce intestinal damage [20]. Further, exosomes partially assumed the function of the local and long-distance communication of cells. These results revealed that the differentially expressed peptides are involved in the pathological mechanism of NEC.KEGG analysis revealed that the TGF-β signaling pathway, apoptosis, and the Wnt signaling pathway had a high enrichment score. It is well known that the main feature of NEC is hemorrhagic and necrotizing inflammation within all layers of the intestinal wall [21]. As a pleiotropic cytokine, TGF-β can regulate multiple cellular functions and suppress immune responses [22]. In intestinal immunity, TGF-β plays a role in suppressing the strong inflammatory response in the intestinal cavity and promotes the construction of immune tolerance, thereby exerting a protective effect in inflammatory bowel diseases [23]. The destruction of the dynamic barrier between the intestinal epithelial layer and the substances in the intestinal cavity is an important process in NEC pathogenesis [24]. Widespread intestinal epithelial cell death is the key mechanism leading to this damage. Studies have reported high levels of intestinal cell apoptosis in the established NEC animal model; moreover, the effect of caspase inhibitors on apoptosis prevents NEC progression. The renewal of intestinal stem cells is a protective response of the intestine to acute injury. Intestinal renewal as well as stem cell maintenance depends mainly on the Wnt/β-catenin pathway [25]. Studies have found that Wnt/β-catenin signaling decreases during NEC development [26, 27]. This leads to impaired intestinal stem cell activity and poor intestinal regeneration ability. However, the administration of Wnt can maintain intestinal epithelial homeostasis and avoid the intestinal injury observed in NEC [28]. These results indicate that the differentially expressed peptides may be potentially involved in the mechanism underlying NEC.Recently, an increasing number of peptides have been found to exert a protective effect in intestinal diseases. For instance, vasoactive intestinal peptides (as effective anti-inflammatory agents) exhibited the ability to regulate homeostasis of the intestinal epithelial barrier. It serves a therapeutic role by reducing inflammation and destroying tight junctions in NEC [29]. In addition, a peptide derived from Porphyra yezoensis was able to promote IEC-6 cell proliferation through the activation of insulin-like growth factor I receptors [30]. However, there are no reports about the treatment of NEC with endogenous active peptides.In our study, peptides derived from the precursor proteins EZRI and MT1 were significantly downregulated in the NEC group. We initially synthesized two peptides EDP1 and MDP2 from the functional domains of Ezrin and Metallothio, respectively. We found that the peptide EDP1 promotes cell migration after acting on the NEC model in vitro, whereas the peptide MDP2 promotes cell migration and inhibits apoptosis in the NEC model in vitro. These results indicate that bioactive peptides are involved in the progression of NEC; however, the mechanisms underlying the peptide functions need to be further investigated.


## 5. Conclusion

In conclusion, we initially screened the differentially expressed peptides in the intestinal tissues of the NEC rat model based on TMT markers combined with LC-MS/MS analysis. Next, we investigated the pathological mechanism of NEC and explored the potential treatment strategies. Although the peptides EDP1 and MDP2 were confirmed to exhibit protective effects in NEC, we need to explore their therapeutic mechanisms further.



---

*Source: 4705149-2020-12-31.xml*