Abstract Purpose Acute pancreatitis can be classified histologically as interstitial edema pancreatitis (IEP) or as acute necrotizing pancreatitis (ANP). ANP has a higher mortality and long-term or short-term sequelae than IEP. Therefore, this work aims to explore the differences in pathogenesis between ANP and IEP and it has great clinical importance for the treatment and prevention of ANP. Methods In this work, whole blood samples from IEP and ANP patients were analyzed by whole gene sequencing (WGS). Serum samples from IEP and ANP patients were evaluated via enzyme-linked immunosorbent assay (ELISA). Meanwhile, pancreatic tissues of IEP and ANP rat models were subjected to data independent acquisition (DIA) proteomics assays. Then, the WGS analysis and DIA proteomics assay data were analyzed comprehensively. Results Six pathways were found to be significantly different in the ANP/IEP groups through WGS analysis. DIA proteomics found eleven different pathways. In both assays, the complement and coagulation cascades pathway was the most significantly different (p < 0.01) pathway between the two groups. WGS analysis showed base mutations in ten genes in the complement and coagulation cascades pathway. These results were consistent with the ten proteins detected by DIA proteomics analysis, which were significantly upregulated in the ANP/IEP groups. In addition, five of these proteins, complement C3, complement Factor I, alpha-2-macroglobulin, complement C9, and serpin family C member 1, were successfully verified by parallel reaction monitoring analysis and ELISA. Conclusion C3, CFI, A[2]m, C9, and Serpinc1, which belong to complement and coagulation cascades pathway, may promote pancreatic necrosis and aggravate the severity of ANP. Keywords: acute necrotizing pancreatitis, complement and coagulation cascades pathway, data independent acquisition proteomics, interstitial edema pancreatitis, whole gene sequencing genomics Introduction Acute pancreatitis (AP) is a prevalent and possibly fatal gastrointestinal disease. The incidence rate is rising.[42]^1 According to histology, AP is classified into acute necrotizing pancreatitis (ANP) and interstitial edema pancreatitis (IEP). Necrosis of pancreatic parenchyma and/or peripancreatic tissue may occur in about 10–20% of AP cases, which can exacerbate the severity of AP. About 45% of ANP patients will develop organ failure, leading to prolonged intensive care unit admission and a complex course of disease, most often requiring multiple invasive interventions.[43]^2 ANP has a mortality rate of approximately 15%, the necrotic pancreatic (peripancreatic) tissue is susceptible to infection.[44]^3 Mortality can be as high as 30% if infection develops in necrotic collection.[45]^3 Previous studies have also confirmed that the mortality, length of intensive care unit admission, morbidity of complications and need for invasive intervention of ANP are significantly higher than those of IEP under the same conditions.[46]^2^,[47]^4^,[48]^5 Therefore, ANP will lead to a large use of medical resources and an increase in costs. However, due to the incomplete understanding of the pathogenesis of ANP, there is still a lack of specific medicine for the treatment of ANP, and the treatment method and timing of infectious ANP are not clear enough.[49]^6 If we can further understand the pathogenesis of ANP, it may be helpful for clinical intervention in advance to prevent the occurrence of pancreatic necrosis. Therefore, the study of the pathogenesis of ANP has great clinical importance. Previous reports have suggested that AP may be associated with genetic variation.[50]^7–9 It has also been reported that only a portion of patients exposed to AP risk factors develop AP, and only a small number of AP patients develop more severe AP with necrosis or organ failure.[51]^10 This suggests that genetic mutations may regulate the occurrence of necrosis in AP. However, we have not found reports about the difference between IEP and ANP in genomics. As a result, the purpose of this paper is to explore the differences in pathogenesis between ANP and IEP through genomics and proteomics, which may guide the development of specific drugs to treat ANP and clinical treatment of ANP. Materials and Methods Patients This research was authorized by the Ethical Committees of the Affiliated Hospital of North Sichuan Medical College. Written informed consent was obtained from all patients and all procedures involving humans maintained adherence to the Declaration of Helsinki. ANP patients hospitalized in the hospital from November 2016 to August 2017 were consecutively recruited. The diagnosis of AP needs to meet at least two of the following three criteria: (1) meet the abdominal pain characteristics of AP, such as severe, acute and persistent upper abdominal pain radiating into the back; (2) imaging findings were consistent with AP; (3) serum amylase or lipase activity was more than 3 times higher than the normal upper limit. Patients enrolled must undergo abdominal contrast-enhanced magnetic resonance (CEMR)/contrast-enhanced computed tomography (CECT) examination and with relevant laboratory data, medical records within 3 days of admission. The exclusion criteria were as follows: (1) acute attack of chronic pancreatitis; (2) AP due to pancreatic tumor; (3) complicated with cancer, severe chronic wasting diseases; (4) AP patients with a family history of genetic disease (hemophilia, inherited heart disease, etc.) or congenital anatomic abnormalities of the pancreas, such as pancreatic divisum; and (5) unsatisfactory images or incomplete medical records. The flowchart of ANP patients is shown in [52]Figure 1. IEP patients who matched with age and sex of ANP patients during the same period were recruited, and the exclusion and inclusion criteria are the same as ANP patients’ criteria. Diagnosis of IEP and ANP by CECT/CEMRI. The main imaging findings of ANP were local or diffuse hypo-density (hypointensity) non-enhanced areas of pancreas parenchyma and/or heterogeneous density (signal intensity) of peripancreatic fluid collection on CT/MRI. However, IEP showed local or diffuse enlargement of pancreas parenchyma, blurred peripancreatic fat space, and homogeneous enhancement of pancreas parenchyma on CT/MRI.[53]^11 The information of sex, age, etiology, length of hospital stay, C-reactive protein (CRP), serum amylase and lipase of patients was collected. The diagnostic criteria for different etiologies of AP were shown in [54]Supplemental Methods and Materials. Other clinical characteristics, including computed tomography severity index (CTSI)/magnetic resonance severity index (MRSI) score, severity based on the 2012 revised Atlanta Classification (RAC), Acute Physiology and Chronic Health Evaluation (APACHE)-II score, extra-pancreatic inflammation on CT (EPIC) score, systemic inflammatory response syndrome (SIRS) score, bedside index of severity in acute pancreatitis (BISAP) score were calculated within 3 days of admission. Blood samples of the two groups were collected by professionals from the patients who were recruited to participate in the experiment within 3 days after admission and stored at −80°C for whole gene sequencing (WGS). Figure 1. [55]Figure 1 [56]Open in a new tab Flowchart of ANP patients’ recruitment. We also collected serum samples from 84 AP patients during the first 3 days of hospitalization in the Affiliated Hospital of North Sichuan Medical College from June 2020 to May 2021. The enrollment and exclusion criteria were consistent with those of the above patients, and they were classified into ANP and IEP groups by CEMR or CECT, with 42 patients in each group. Serum samples were placed at −80°C for further analysis by enzyme-linked immunosorbent assay (ELISA). WGS Analysis Blood samples of IEP patients (n = 21) and ANP patients (n = 26) were analyzed by WGS. Genomics DNA sequencing was performed using the Illumina NovaSeq 6000 platform (Illumina Inc., San Diego, CA, USA) at Novogene Bioinformatics Technology Co., Ltd. (Beijing, China). For details of sequencing and bioinformatics analysis, see [57]Supplemental Methods and Materials. Establishment of IEP and ANP Rat Models The animal protocols for rat studies were approved by the Animal Care and Use Committee of North Sichuan Medical College. The animal studies were implemented at the Experimental Animal Center of North Sichuan Medical College and were managed in full compliance with China Animal Welfare Act, the China Public Health Service Policy on Humane Care and Use of Laboratory Animals and the National Research Council’s Guide for the Care and Use of Laboratory Animals. Male adult Sprague–Dawley rats (n = 10, 6–8 weeks, 200–250 g) were provided by the Experimental Animal Center of North Sichuan Medical College. They were randomly classified into two groups of five rats each. In one group, IEP was established by retrograde pancreaticobiliary duct injection with 0.9% saline, and another group, ANP was established by similar injections with 5% sodium taurocholate (Solarbio Science Technology, Beijing, China). Then, rats were killed six hours later for the collection of blood samples and pancreatic tissues. Blood was immediately centrifuged at 3000 g for 5 min after collection. The serum was retained and stored at −80°C until detection. The pancreas was removed and divided into two parts for later analyses. One part was stored at −80°C for data independent acquisition (DIA) proteomics assay and validation via parallel reaction monitoring (PRM). One part was fixed in formalin for histological analysis. Serum levels of amylase in AP rats were measured by using an enzymatic colorimetric assay (Roche Diagnostics, Indianapolis, IN). Hematoxylin and eosin (H&E) were used to stain the pancreatic sections to conduct the histopathologic analyses under light microscopy. Schmidt’s scale was used to classify AP degree.[58]^12 This scale scored edema, acinar necrosis, hemorrhage and fat necrosis, and inflammatory and perivascular infiltration extent. DIA Proteomics Analysis and PRM Validation of Target Proteins DIA proteomics analysis and PRM validation were performed on pancreatic tissues of IEP (n = 5) and ANP (n = 5) rats. This process was carried out in Shanghai Applied Protein Technology Co., Ltd., see [59]Supplemental Methods and Materials for details. Validation of Target Proteins in Patient Serum by Enzyme-Linked Immunosorbent Assay (ELISA) The levels of complement C3 (C3), complement Factor I (CFI), alpha-2-macroglobulin (A[2]m), complement C9 (C9), and serpin family C member 1 (Serpinc1) in the serum of AP patients were determined using double antibody ELISA (Shanghai Enzyme-linked Biotechnology Co., Ltd., China) commercial kits (ml058121-J, ml063609-J, ml060470-J, ml060470-J, and ml060054-J, respectively) according to the manufacturer’s instructions. Statistical Analysis The hierarchical and categorical variables were expressed as percentages. The continuous variables were described as medians (range) or mean ± SEM. Etiology and sex were compared through Fisher’s exact test or chi-squared test. On the basis of data distribution, other clinical features of the two groups were compared via Mann–Whitney U-test or independent t-test. Histopathologic investigation and the analysis of enzyme levels were performed using Kruskal–Wallis and analysis of variance (ANOVA) test. All statistical analyses were implemented by applying the statistical analysis software SPSS version 25.0 (SPSS Inc., Chicago, IL). A p value less than 0.05 was regarded as statistically significant. Result Patient Characteristics Statistical analysis results of basic clinical data and various scoring data of the patients in the IEP and ANP groups were shown in [60]Tables 1 and [61]2. Between the ANP and IEP groups, there did not exist statistical differences in etiology, age, sex, BISAP score, or EPIC score. The CTSI/MRSI score, severity based on RAC, CRP level, and SIRS score in the ANP group were much higher than those in the IEP group (p < 0.05). Although there was no significant difference in hospital stay or APACHE II score between the two groups, the average value of the ANP group was higher than that of the IEP group. The levels of amylase and lipase in the IEP group were much higher, compared with ANP group. Table 1. Basic Clinical Characteristics of Patients with IEP and ANP Total (n=47) IEP (n =21) ANP (n = 26) p value Sex 47 (100%) 21 (44.7%) 26 (55.3%) 0.06  Male 30 (63.8%) 12 (57.1%) 18(69.2%)  Female 17 (36.2%) 9 (42.9%) 8 (30.8%) Age (years) 45.96 ± 2.19 46.05 ± 3.61 45.88 ± 2.74 0.97 Etiology 0.97  Biliary 11 (23.4%) 6 (28.6%) 5 (19.2%)  Alcoholic 13 (27.7%) 7 (33.3%) 6 (23.1%)  Hypertriglyceridemia 11 (23.4%) 3 (14.3%) 8 (30.8%)  Idiopathic 12 (25.5%) 5 (23.8%) 7 (26.9%) Hospital stay (days) 13.85 ± 0.86 12.10 ± 0.97 15.27 ± 1.28 0.06 CRP 31.44 (1.75~403.00) 10.06 (1.75~403.00) 37.57 (2.34~307.00) 0.04* Amylase 339.60 (32.50~2256.0) 515.50 (87.70~2256.00) 264.60 (32.50~1321.00) 0.02* Lipase 4431.50 (21.80~6020.00) 960.90 (142.40~6020.00) 358.50 (21.80~2994.30) 0.02* [62]Open in a new tab Note: *p value <0.05. Abbreviations: IEP, interstitial edema pancreatitis; ANP, acute necrotizing pancreatitis; CRP C-reactive protein. Table 2. Severity Scores Comparison Between IEP and ANP Patients IEP (n =21) ANP (n = 26) p value CTSI/MRSI 3 (0~4) 5 (3~10) <0.001* Severity based on RAC <0.01*  Mild 8 (38.10%) 2 (7.70%)  Moderately severe 13 (61.90%) 22 (84.60%)  Severe 0 (0.00%) 2 (7.70%) APACHE II 3.00 ± 0.57 4.46 ± 0.60 0.09 EPIC 5 (0~7) 5 (2~7) 0.25 BISAP 1 (0~3) 1 (0~3) 0.54 SIRS 0 (0~3) 1 (0~3) 0.03* [63]Open in a new tab Note: *p value <0.05. Abbreviations: IEP, interstitial edema pancreatitis; ANP, acute necrotizing pancreatitis; CTSI, computed tomography severity index; MRSI, magnetic resonance severity index; RAC, revised Atlanta classification; APACHE II, acute physiology and chronic health evaluation; EPIC, extra-pancreatic inflammation on computed tomography; BISAP, bedside index of severity in acute pancreatitis; SIRS, systemic inflammatory response syndrome. Genetic Mutations in the Complement and Coagulation Cascades Pathway May Contribute to ANP To study the differences in genomes between ANP and IEP, patients in ANP group (n = 26) or IEP group (n = 21) underwent WGS. Through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, six pathways were significantly different (p < 0.05) between the two groups ([64]Table 3). As seen from the KEGG bubble diagram ([65]Figure 2), the difference in the complement and coagulation cascades pathway was most significant (p < 0.01), with more genes enriched in this pathway. In the complement and coagulation cascades pathway, 22 patients in the ANP group had gene mutations in the exon or splice region, including 39 sites in 29 genes and Nine patients in the IEP group had mutations in their exons, including 16 sites in 14 genes. Among them, 6 mutation sites of ANP groups were located in the splicing region, and there was no annotation information in ANNOVAR. The annotation of one mutation site in each group was unknown in ANNOVAR. Specific mutation sites in the complement and coagulation cascades pathway in the ANP group and IEP group are shown in [66]Table 4 Table 3. Differential Pathways in WGS Analysis Between ANP and IEP Patients Pathway Alt_ANP Ref_ANP Alt_IEP Ref_IEP p value Complement and coagulation cascades 22 4 9 12 <0.01* Alanine aspartate and glutamate metabolism 10 16 16 5 0.02* Apoptosis 18 8 7 14 0.02* Ppar signaling pathway 15 11 19 2 0.02* Allograft rejection 8 18 1 20 0.03* Oxidative phosphorylation 18 8 8 13 0.04* [67]Open in a new tab Notes: Alt_ANP/Alt_IEP the number of patients who occurred mutant bases in this pathway in ANP group/IEP group, Ref_ANP/Ref_IEP the number of patients who consisted with reference bases in this pathway in ANP group/IEP group; *p value <0.05. Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis; WGS, whole genome sequencing. Figure 2. [68]Figure 2 [69]Open in a new tab Differential pathways between IEP and ANP patients analyzed by WGS. Note: The pathway with the most significant difference (p < 0.01) is highlighted by dotted box and arrow. Table 4. Mutations in Complement and Coagulation Cascades Pathway in Patients with ANP and IEP Group Gene ID Region Variation Type Accession Genotype^a ANP A[2]m rs775058537 Splicing graphic file with name JIR_A_351416_O_ILF0001.jpg [70]NC_000012.11 g.9221327_9221330del A[2]m rs200420993 Exonic Missense SNV [71]NC_000012.1 g.9266056G>A C1R rs749274074 Exonic Unknown [72]NC_000012.1 g.7241493C>T C3 rs138900723 Exonic Missense SNV [73]NC_000019.9 g.6690729G>A C4BPA rs140703266 Exonic Missense SNV [74]NC_000001.10 g.207317935T>C C4BPB rs764052227 Exonic Missense SNV [75]NC_000001.10 g.207263768G>T C6 graphic file with name JIR_A_351416_O_ILF0002.jpg Exonic Missense SNV [76]NC_000005.9 g.41159349A>G C7 rs200622924 Exonic Missense SNV [77]NC_000005.9 g.40958328C>T C7 rs375805678 Exonic Missense SNV [78]NC_000005.9 g.40981640A>G C8A rs145228353 Exonic Missense SNV [79]NC_000001.10 g.57341837G>A C8G rs148537245 Exonic Missense SNV [80]NC_000001.10 g.139840422C>T C8G rs200109501 Splicing graphic file with name JIR_A_351416_O_ILF0003.jpg [81]NC_000001.10 g.139841101A>G C9 rs121909592 Exonic Stop gain [82]NC_000005.9 g.39341378G>A CD46 graphic file with name JIR_A_351416_O_ILF0004.jpg Splicing graphic file with name JIR_A_351416_O_ILF0005.jpg [83]NC_000001.10 g.207930549_ 207930550ins AAGTAAACAAA CD46 rs368371683 Exonic Missense SNV [84]NC_000001.10 g.207934692G>C CFH rs367837226 Splicing graphic file with name JIR_A_351416_O_ILF0006.jpg [85]NC_000001.10 g.196646807A>G CFI rs763931500 Exonic Missense SNV [86]NC_000004.11 g.110670410C>A CPB2 rs202020888 Exonic Missense SNV [87]NC_000013.10 g.46629913A>C CR1 rs55906048 Exonic Missense SNV [88]NC_000001.10 g.207680151A>C CR1 graphic file with name JIR_A_351416_O_ILF0007.jpg Splicing graphic file with name JIR_A_351416_O_ILF0008.jpg [89]NC_000001.10 g.207787835G>A CR1 graphic file with name JIR_A_351416_O_ILF0009.jpg Exonic Missense SNV [90]NC_000001.10 g.207791445G>A CR2 rs201760344 Exonic Missense SNV [91]NC_000001.10 g.207642150C>T F11 rs281875279 Exonic Stopgain [92]NC_000004.11 g.187197527G>A F12 rs138580930 Exonic Missense SNV [93]NC_000005.9 g.176829661C>A F12 rs368879882 Exonic Missense SNV [94]NC_000005.9 g.176832982G>A F5 rs774639785 Exonic Missense SNV [95]NC_000001.10 g.169495167C>T F7 rs121964936 Exonic Missense SNV [96]NC_000013.10 g.113773145T>G KLKB1 graphic file with name JIR_A_351416_O_ILF0010.jpg Exonic Missense SNV [97]NC_000004.11 g.187179251G>T KNG1 rs187202560 Splicing graphic file with name JIR_A_351416_O_ILF0011.jpg [98]NC_000003.11 g.186459303C>T KNG1 rs1085307077 Exonic Frameshift deletion [99]NC_000003.11 g.186460048_ 186460051del MASP1 rs755157142 Exonic Missense SNV [100]NC_000003.11 g.186944295G>C PLG rs4252128 Exonic Missense SNV [101]NC_000006.11 g.161152819C>T PROC rs121918145 Exonic Missense SNV [102]NC_000002.11 g.128183754C>T SERPINA1 rs200945035 Exonic Missense SNV [103]NC_000014.8 g.94844968T>C SERPINA5 rs201869477 Exonic Missense SNV [104]NC_000014.8 g.95056500C>G SERPIND1 rs181138395 Exonic Missense SNV [105]NC_000022.10 g.21138310C>T VWF rs773544469 Exonic Missense SNV [106]NC_000012.11 g.6127617G>A VWF graphic file with name JIR_A_351416_O_ILF0012.jpg Exonic Missense SNV [107]NC_000012.11 g.6143942G>A VWF rs189409574 Exonic Missense SNV [108]NC_000012.11 g.6172134T>C IEP C1QC rs560156356 Exonic Missense SNV [109]NC_000001.10 g.22974053C>T C1R rs144141261 Exonic Unknown [110]NC_000012.11 g.7242267G>A C5 rs747276237 Exonic Missense SNV [111]NC_000009.11 g.123752013C>T C5 graphic file with name JIR_A_351416_O_ILF0013.jpg Exonic Missense SNV [112]NC_000009.11 g.123760048C>T C9 rs527488079 Exonic Missense SNV [113]NC_000005.9 g.39308359A>G CR1 rs376980117 Exonic Missense SNV [114]NC_000001.10 g.207785379G>A F10 graphic file with name JIR_A_351416_O_ILF0014.jpg Exonic Missense SNV [115]NC_000013.10 g.113803358A>G F5 rs201292247 Exonic Missense SNV [116]NC_000001.10 g.169512296T>C F9 rs762986355 Exonic Missense SNV [117]NC_000023.10 g.138633223C>T KNG1 rs1085307077 Exonic Frameshift deletion [118]NC_000003.11 g.186460048_ 186460051del MASP1 rs3774266 Exonic Missense SNV [119]NC_000003.11 g.186937924C>T PLAU rs192623012 Exonic Missense SNV [120]NC_000010.10 g.75673753G>T PLG rs139357983 Exonic Missense SNV [121]NC_000006.11 g.161132157C>T Serpinc1 graphic file with name JIR_A_351416_O_ILF0015.jpg Exonic Missense SNV [122]NC_000001.10 g.173876643G>A VWF rs372028373 Exonic Missense SNV [123]NC_000012.11 g.6128632C>T VWF rs139196998 Exonic Missense SNV [124]NC_000012.11 g.6167131G>A [125]Open in a new tab Notes: ^aNomenclature following HGVS recommendations;[126]^42 Inline graphic Indicates no annotation in dbSNP or ANNOVAR. Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis. C3, CFI, A[2]m, C9, and Serpinc1 Were Significantly Upregulated in ANP H&E-stained pancreatic sections of rats from the control group, IEP group, and ANP group were investigated under a light microscope ([127]Figure 3A). The results showed that the control group had normal pancreatic structure without edema and inflammation. The IEP group showed pancreatic edema with fluid accumulation in interlobular septa but no cell necrosis. The ANP group showed pancreatic edema, inflammation, and acinar cell necrosis. Histological scores were 0 in the control group, 2–5 in the IEP group, and 3–6 in the ANP group ([128]Figure 3B). Serum amylase levels of rats in the IEP group (5924 ± 546.30 U/L, n = 5) and ANP group (4369 ± 287.80 U/L, n = 5) were both higher than that in the control group (2386 ± 49.87 U/L, n = 5) ([129]Figure 3C). Figure 3. [130]Figure 3 [131]Open in a new tab Evaluation of IEP and ANP rat models. Notes: (A) Representative H&E staining sections of pancreatic tissues of rats in the three groups were obtained under a 20×10 times light microscope. (B) Histogram of pathological scores of pancreatic tissue sections of rats in the three groups. (C) Histogram of serum amylase levels in three groups of rats. *Statistically significant differences between the IEP and control groups (p < 0.05). ^#Statistically significant differences between the ANP and control groups (p < 0.05). Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis. A total of 367 differentially expressed protiens (DEPs) in pancreas of rats between the ANP and IEP groups were selected by significant difference protein screening. Among them, 254 DEPs were upregulated in the ANP group, whereas 113 DEPs were downregulated when compared with the IEP group. The results are illustrated by a volcano plot ([132]Figure 4A). In addition, the hierarchical clustering heatmap of 367 DEPs was obtained after normalization with z score, and according to the Euclidean distance, the proteins were grouped into four clusters and effectively separated IEP samples from ANP samples ([133]Figure 4B). Figure 4. [134]Figure 4 [135]Open in a new tab DIA proteomics analysis of IEP and ANP rats. Notes: (A) DEPs between the two groups are displayed by a volcano plot. Compared with the IEP group, the upregulated proteins in the ANP group are shown in red, downregulated proteins are shown in blue, and those with no difference are shown in gray. (B) Hierarchical clustering heatmap of DEPs between the IEP and ANP groups. Each column represents a sample, and each row represents a DEP. Red represents significantly upregulated proteins, blue represents significantly downregulated proteins, and gray represents no quantitative protein information. (C) The top 11 pathways of the KEGG pathway enrichment analysis in the IEP and ANP groups, which are marked with a dotted box, were the most significantly different pathways (p < 0.05). (D) The abscissa represents the GO function category, and the ordinate represents the rich factor and number of proteins associated with this GO function. The shade of color indicates the p value, which is how significantly a function is affected. Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology. To further understand the pathological differences between IEP and ANP, KEGG and Gene Ontology (GO) enrichment analysis was used to investigate the DEPs. There were eleven pathways between the IEP and ANP groups with significant differences (p < 0.05), including the complement and coagulation cascades pathway, neutrophil extracellular trap formation, and NF-kappa B signaling pathway ([136]Figure 4C). The difference in the complement and coagulation cascades pathway was most significant, which was consistent with the WGS analysis. GO analysis of DEPs was performed to understand the biological pathways, functions, and localization of proteins in the organism. In [137]Figure 4D, the top twenty GO enrichment analyses are displayed. For a biological process (BP), these proteins principally participate in the inflammatory response and immune response. Molecular functions (MF) mainly included the regulator and inhibitor activities of enzymes. For the cellular components (CC), these DEPs mainly existed in the extracellular space and mitochondria. In combination with KEGG and GO analysis, eighteen proteins in the complement and coagulation cascades pathway were significantly different between the two groups. Ten of the proteins were identical to the gene with base mutations in WGS ([138]Table 5). Compared with the IEP group, these ten proteins were significantly upregulated in the ANP group. Table 5. Differential Proteins of ANP and IEP in Complement and Coagulation Cascades Pathway Protein Gene ANP/IEP p value Alpha-2-macroglobulin A[2]m 7.262 <0.01* Complement C8 gamma chain C8g 2.093 0.02* Coagulation factor XII F12 2.056 <0.01* Complement Factor H CFH 2.038 <0.01* Complement C3 C3 1.992 0.04* Complement component 4 binding protein, alpha C4bpa 1.938 0.03* Complement Factor I CFI 1.714 <0.01* Complement C9 C9 1.613 <0.05* Complement C7 C7 1.578 <0.01* Serpin family C member 1 Serpinc1 1.524 0.04* [139]Open in a new tab Note: *p value <0.05. Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis. Finally, five proteins were successfully verified by PRM. They were C3, CFI, A[2]m, C9, and Serpinc1. These proteins are primarily involved in inflammatory reactions and immune responses. Compared with the IEP group, the expression of C3, CFI, A[2]m, C9, and Serpinc1 in the ANP group was upregulated by 1.71, 1.38, 4.65, 1.48, and 1.06 times, respectively. The upregulation trend of each protein was almost consistent with the results of the DIA proteomics assay. The fold change of the five proteins in the DIA proteomics assay and PRM are shown in [140]Figure 5. Figure 5. [141]Figure 5 [142]Open in a new tab Fold changes of C3, CFI, A[2]m, C9, and Serpinc1 in the ANP/IEP group detected by DIA proteomics and PRM. Abbreviations: PRM, parallel reaction monitoring; DIA, data independent acquisition; ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis. In combination with the results of genomics sequencing, three of the genes consistent with these five proteins were mutated only in ANP group, including C3 ([143]NC_000019.9: g.6690729G>A); CFI ([144]NC_000004.11: g.110670410C>A), and A[2]m ([145]NC_000012.11: g.9266056G>A). One mutation, Serpinc1 ([146]NC_000001.10: g.173876643G>A), occurred only in the IEP group. In addition, C9 was mutated in both groups, but the sites and types of mutations were different. The mutation type of C9 in ANP group was stop gain (C9 ([147]NC_000005.9: g.39341378G>A)), while in IEP group it was missense SNV (C9 ([148]NC_000005.9: g.39308359A>G)). CFI and Serpinc1 Levels in ANP Patients Were Increased, While A[2]m Levels Were Reduced The levels of C3, CFI, A[2]m, C9, and Serpinc1 in the serum of patients with IEP and ANP were detected by ELISA. The results are shown in scatter diagrams in [149]Figure 6. Although there was no significant difference in the level of C3 between the IEP and ANP groups, the proportion of high levels (100–150 µg/mL) in the IEP group (31.03%) was significantly higher than that in the ANP group (18.52%). The levels of CFI, A[2]m, and Serpinc1 in the serum of patients in the two groups were significantly different (p = 0.02, p = 0.04, and p < 0.05, respectively). The average levels of CFI and Serpinc1 in the ANP group were higher than those in the IEP group, while the average level of A[2]m in the ANP group was lower than that in the IEP group. There existed no evident difference in the levels of serum C9 between the ANP and IEP groups, but the ANP group had a higher mean (9.63 ±1.76 µg/mL) than the IEP group (7.89 ±1.76 µg/mL). Figure 6. [150]Figure 6 [151]Open in a new tab The levels of C3, CFI, A[2]m, C9, and Serpinc1 in IEP and ANP patients were detected by ELISA. Abbreviations: ANP, acute necrotizing pancreatitis; IEP, interstitial edema pancreatitis. Discussion This study includes imaging and genomics in patients with AP and proteomics with rat model of AP. Based on imaging classification, the differences in genes and proteins between ANP and IEP were studied to better understand the molecular mechanism of pancreatic necrosis. In our study, the CTSI/MRSI, SIRS scores, severity based on RAC, and CRP of ANP patients were significantly higher than those of IEP patients. Combined with WGS analysis of AP patients and DIA proteomics assays of AP rats and found that the complement and coagulation cascades pathway was significantly different between the ANP and IEP groups. And in this pathway C3 ([152]NC_000019.9: g.6690729G>A); CFI ([153]NC_000004. 11: g.110670410C>A); A[2]m ([154]NC_000012. 11: g. 9266056G>A); C9 ([155]NC_000005.9: g.39341378G>A) mutation occurred in ANP group, and C9 ([156]NC_000005.9: g.39308359A>G); Serpinc1 ([157]NC_000001.10: g.173876643G>A) gene mutation occurred in IEP group. DIA proteomics found that C3, CFI, A[2]m, C9, and Serpinc1 proteins in ANP group was significantly higher than those in IEP group. The accumulated data show that neutrophils can be recruited by the complement system and that the system can also adjust inflammation and tissue injury under various inflammatory conditions.[158]^13–15 Studies have proven that complement system activation in AP is related to the severity of disease, particularly in ANP models.[159]^14 The data show that C3 expression is higher in the pancreatic tissues of ANP rats than in those of IEP rats. Interestingly, although the average levels of C3 in serum of ANP patients were only slightly higher than those in IEP patients, and there was no statistical difference between the two groups, the proportion of high C3 levels in the IEP group was higher than that in the ANP group. This finding is consistent with the findings of Seelig et al, who showed that in a taurocholate-induced rat model of AP, massive C3 deposits could be detected in the vicinity of acinar necroses and necrobioses.[160]^16 Multiple studies have found that serum C3 levels were significantly reduced in ANP patients,[161]^14^,[162]^17 suggesting that C3 consumption may play a role in the development of ANP. A recent study demonstrated that C3 is an effective promoter of systemic inflammation, pancreatic tissue damage, neutrophil recruitment, and neutrophil extracellular trap (NET) generation in the pathogenesis of AP.[163]^15 There is sufficient evidence that neutrophil recruitment together with subsequent NET generation is one of the pivotal components of SAP development.[164]^18 In addition, C3 is involved in all complement activation pathways and can interact with many complement components and factors. For example, CFI can cleave C3 to C3b with assistance from cofactors.[165]^19 C3b is an effective humoral inflammatory mediator. CFI further cleaves C3b to C3bi, which can bind to complement receptors CR3 (and probably CR4) on neutrophils to trigger inflammation.[166]^20 CFI also inactivates C3b, reduces the production of C3 convertase (C3bBb), and prevents excessive consumption of C3 in plasma. Both an increase or decrease in CFI may upset this balance. Our data showed that CFI levels were markedly higher in the ANP group than in the IEP group. This indicates that CFI may be correlated with ANP. However, the exact mechanism is unclear. In addition, trypsin can cleave C3 as well.[167]^21 As a broad-spectrum protease inhibitor in human plasma, A[2]m can inhibit trypsin[168]^22^,[169]^23 and play a central role in defense proteinase-induced release of C3 cleavage products.[170]^24 In AP, trypsin is rapidly inhibited to form a stable compound with A[2]m.[171]^25 Therefore, in patients with AP, concentrations of A[2]m are typically depressed but decrease profoundly and for longer periods in more severe cases.[172]^26–28 This has been interpreted as A[2]m being consumed by activated trypsin and subsequent clearance of the A[2]m protease complexes from the intravascular space by cells of the reticuloendothelial system within approximately 10 minutes.[173]^29 However, clearance of complexes from the peritoneal fluid is probably much slower.[174]^30 This is consistent with our results that A[2]m in the serum of ANP patients was markedly lower than that in IEP patients, while A[2]m was overtly increased in the inflamed pancreatic tissue of ANP rats compared with IEP rats. The experimental data also showed that the content of C9 in the pancreas of ANP rats was observably higher than that of IEP rats. The mean level of serum C9 in ANP patients was higher than that in IEP patients. Many studies have proposed that C9 is associated with necrosis, apoptosis, and the proinflammatory response. As one of the key components of the membrane attack complex (MAC),[175]^31^,[176]^32 the presence of C9 on the cell surface indicates the formation of MAC.[177]^33 MAC is the terminal complement cascade effector.[178]^34 It binds to the target cell surface and leads to cell lysis by generating transmembrane channels or pores on the cell membrane.[179]^35^,[180]^36 These assembled pores are thought to allow lysozymes to degrade the peptidoglycan layer through the outer membrane.[181]^37^,[182]^38 In addition, Triantafilou et al proposed that MAC insertion into the cell membrane could trigger Ca^2+ influx, increase cytosolic Ca^2+ concentration, and then accumulate in the mitochondrial matrix, ultimately resulting in mitochondrial dysfunction, apoptosis, inflammasome activation, and the secretion of IL-1β.[183]^39 MAC-activated cells can express proinflammatory proteins via NF-κB-dependent transcription, assemble inflammasomes, activate processing, and promote IL-1β and IL-18 secretion.[184]^40^,[185]^41 This may explain the higher levels of C9 in the pancreas of the rats in the ANP group in this study. Proteomics analysis of the pancreas of AP rats and ELISA analysis of the serum of AP patients showed that Serpinc1 expression was increased when ANP occurred. This suggests that Serpinc1 is likely to be associated with ANP and may be a biomarker for predicting ANP. However, there have been no reports about the relationship between Serpinc1 and AP. Our study has some limitations. First, our research involved a single-center study, and a few patients underwent WGS analysis, which may lead to selection bias. Second, although we found base mutations in the C3, CFI, A[2]m, C9, and Serpinc1 genes in IEP and ANP patients and verified their protein expression levels in IEP and ANP rat models, the molecular mechanism was not further studied through in vivo and in vitro investigations. Finally, the network of interrelationships among the five genes is not fully understood. Therefore, the molecular mechanism and network of these five genes in the development of ANP should be investigated in detail in future studies. Conclusion C3, CFI, A[2]m, C9, and Serpinc1, which belong to complement and coagulation cascades pathway could play an important role in the pathogenesis of ANP. Moreover, C3, CFI, A[2]m, C9, and Serpinc1 in serum of patients may be used as biomarkers to predict ANP in the future. Acknowledgments