Abstract

Objective. To explore potential biomarkers to accurately diagnose patients with acute pancreatitis (AP) at early stage and to auxiliary clinicians implement the best treatment options. Methods. We selected 3 patients with AP and 3 healthy controls for microarray analysis to obtain differentially expressed circular RNAs (circRNAs). To further verify the results of the microarray analysis, the six differentially expressed circRNAs were confirmed by quantitative polymerase chain reaction (qPCR). The diagnostic accuracy and sensitivity of differentially expressed circRNAs were assessed using the receiver operating characteristic (ROC) curve. A ceRNA network was constructed based on the 6 differentially expressed circRNAs. Results. There were 25 upregulated circRNAs and 26 downregulated circRNAs in the blood of patients with AP. Next, the qPCR verification results further confirmed three downregulated circRNAs, including hsa_circRNA_002532, has_circRNA_059665, and hsa_circRNA_104156, and three upregulated circRNAs including hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470. Among them, hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 increased with the severity of the disease. ROC analysis showed that the three circRNA models show promise to diagnose AP. And a ceRNA network revealed that above six circRNAs could participate in complex regulated network. Conclusions. Elevated hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 could be used as novel biomarkers to diagnose AP patients.

1. Introduction

AP is one of the most common gastrointestinal diseases, and patients need hospitalization [1]. Although most patients are mild, about 20% of patients develop moderate or severe pancreatitis with surrounding pancreas tissue necrosis or multiple organ failure [24]. The overall mortality rate of AP patients is about 2%, but close to 30% in patients with persistent organ system failure. An increased risk of pancreatic cancer in AP patients was observed in a nationwide, population-based matching cohort study [5]. Furthermore, 17%–22% of AP patients have the possibility of recurrence, and 8%–16% may develop chronic pancreatitis [6, 7]. Therefore, early diagnosis is necessary for AP patients.

As a pancreatic inflammatory disease, AP has become one of the leading causes of gastrointestinal disease admission in the United States and many other countries. And the incidence rate demonstrates an increasing trend. There are many factors that cause the disease, such as gallstones, smoking, and alcohol abuse [8, 9]. Of course, about 10% of patients have no cause, and they are idiopathic pancreatitis. The patients have usually a sudden onset with severe persistent abdominal pain, and 80% of patients are accompanied by vomiting. Pain may radiate to the back, usually in the lower chest area. Therefore, it is critical to carry out some treatments within 48 hours, such as fluid resuscitation, analgesia, and nutritional support [10]. Despite the continuous improvement in diagnostic techniques, early diagnosis is still difficult for clinicians.

Therefore, it is critical to explore early diagnostic biomarkers of AP [11]. More and more evidence has confirmed that circRNAs could play an important role in many diseases. CircRNAs are produced by reverse splicing of precursor mRNAs from exons of thousands of genes in eukaryotes [12]. At the same time, circRNAs can be secreted into blood, saliva, and other body fluids as potential biomarkers for disease prediction [13]. However, the function of most circRNAs remains largely unknown. In the ceRNA mechanism, microRNAs (miRNAs) are important posttranscriptional regulators of gene expression that act on target sites in the untranslated region of messenger RNA (mRNA) by direct base pairing. The circRNA can act as a miRNA sponge to affect the activity of miRNAs in the regulation of mRNA expression [14, 15]. At the same time, circRNA abnormalities can lead to a variety of diseases. For example, the circMTO1 could inhibit the progression of hepatocellular carcinoma by promoting p21 expression, acting as a sponge of oncogenic miR-9 [16]. However, to date, no study has explored the expression of circRNAs in the blood of patients with AP. Therefore, our aim was to explore accurate biomarkers and then diagnose patients with AP as soon as possible by detecting biomarkers in the patients’ blood.

2. Materials and Methods

2.1. Patient Information and Diagnostic Criteria

We reviewed 60 patients who were diagnosed with AP from April 2018 to September 2018 at the First Affiliated Hospital of China Medical University (Shenyang, China). In addition, we recruited 30 subjects who underwent routine health checks at the First Affiliated Hospital of China Medical University and showed no signs of disease as a control group. According to the “Guidelines for the diagnosis and treatment of acute pancreatitis (2014),” clinically meets 2 of the following 3 characteristics to diagnose AP: (1) abdominal pain consistent with AP; (2) serum amylase and/or lipase activity is at least 3 times higher than the upper limit of normal; (3) abdominal imaging examination is consistent with AP imaging changes.

In clinical treatment, mild acute pancreatitis (MAP) patients receive only relatively simple treatment, while severe acute pancreatitis (SAP) patients usually require intensive care. Therefore, we selected 30 MAP and 30 SAP patients to diagnose AP patients with early stage. Currently, common scoring standards are the APPACHE II scoring standard, the MCTSI scoring standard, and BISAP scoring standard. MAP diagnostic criteria were good response to fluid supplementation, without organ failure, and local or systemic complications, recovery within 1-2 weeks. And APACHE-II score <8 points or the MCTSI score <4 points or BISAP <2 are MAP. Diagnostic criteria for SAP were with persistent organ failure (48 h or more). And APACHE-II score ≥8 points or MCTSI score ≥4 points or BISAP ≥2 are SAP.

The exclusion criteria were as follows: <18 years of age, pregnant and lactating women, taking anticoagulant drugs, blood system diseases, tumors, liver disease, and gastrointestinal bleeding patients. The patients’ information is shown in Table 1.

Then we selected 3 patients with MAP and 3 healthy participants for microarray analysis. The study was approved by the Ethics Committee of the First Affiliated Hospital of China Medical University, and informed consent was obtained from all subjects. Figure 1 depicts further experimental design details.

2.2. RNA Extraction and CircRNA Microarray Analysis

We collected whole blood from 3 MAP patients and 3 healthy participants. Total RNA from each sample was quantified using a NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). Sample preparation and microarray hybridization were performed based on the standard protocol of Arraystar. Briefly, total RNA was digested with Rnase R (Epicenter, Inc.) to remove linear RNA and enrich for circular RNA. Then, the enriched circular RNA was amplified by a random priming method (Arraystar Super RNA Labeling Kit; Arraystar) and transcribed into fluorescent cRNA. The labeled cRNA was hybridized to Arraystar Human circRNA Array V2 (8 × 15 K, Arraystar). The labeled cRNA was purified by RNeasy Mini Kit (Qiagen). The concentration and specific activity of the labeled cRNA (pmol Cy3/μgcRNA) were measured by NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). 1 μg of each labeled cRNA was fragmented by adding 5 μl of 10× blocking agent and 1 μl of 25× fragmentation buffer; then the mixture was heated at 60°C for 30 minutes, and finally 25 μl of 2× hybridization buffer was added to dilute the labeled cRNA. 50 μl of the hybridization solution was dispensed into a spacer slide and assembled onto a circRNA expression microarray slide. Slides were incubated for 17 hours at 65°C in Agilent Hybridization Oven. The hybridization array was fixed and scanned using an Agilent Scanner G2505C wash.

Agilent Feature Extraction software (version 11.0.1.1) was utilized to analyze the acquired array images. Quantile normalization and subsequent data processing were performed using the R software limma package (R version 3.1.2).

2.3. Differential Expression Analysis

A scatter plot is a visualization method used to assess circRNA expression variation. Differentially expressed circRNAs with statistical significance (FC ≥ 1.5 and values ≤0.05) were identified utilizing fold change cutoffs and volcano plots, respectively. Among them, value was calculated utilizing the unpaired t-test. Differentially expressed circRNAs between the two samples were identified by fold change filtration. Hierarchical clustering was performed to show the distinguishable circRNA expression patterns in the samples. Shanghai Kangcheng Biological Engineering Co., Ltd. of the People’s Republic of China conducted microarray work.

2.4. qPCR Verification

Total RNA was extracted from the whole blood of 30 healthy participants, 30 MAP patients, and 30 SAP patients according to standard procedures. Total RNA was reverse-transcribed into cDNA kits (Roche, Penzberg, Germany) using random primers and Transcriptor First Strand cDNA Synthesis according to the manufacturer’s instructions. 6 differentially expressed circRNAs were measured by qPCR using a ViiA 7 Real-time PCR System (Applied Biosystems). The reaction conditions were as follows: 95°C for 10 minutes and 40 cycles of 95°C for 10 seconds, 60°C for 60 seconds. RNA levels were normalized to human β-actin. All qPCR reactions were performed in triplicate. Different primers were designed for circRNAs, rather than the more commonly used convergent primers (Figure 2). Table 2 lists all the primers.

2.5. Statistical Analysis

Statistical analyses were performed using SPSS (version 19.0; IBM, Armonk, NY, United States) and GraghPad Prism (version 7.0; GraphPad Software, La Jolla, CA, United States). The relative expression level of each circRNA was expressed by fold change and converted into 2−ΔΔCt. The expression difference in circRNAs among SAP, MAP patients, and healthy individuals and between MAP patients and posttreatment serum samples was assessed using the t-test. To assess the diagnostic value, ROC curve was established. The cut-off value of each circRNA was analyzed by using SPSS software. Area under the ROC curves (AUCs) were calculated to evaluate the ability of the differentially expressed circRNAs. Due to the relative small sample size for ROC analysis, the statistical power was calculated using PASS (version 15.0), under the following conditions: α = 0.05, AUC0 = 0.5, and n = 30. value < 0.05 was considered statistically significant.

2.6. The Construction of ceRNA Network

By combining cotargeted miRNAs, we constructed a ceRNA network by cytoscape package of R language [17, 18]. In addition to measuring the number of common miRNAs, each ceRNA pair was subjected to a hypergeometric test, which was defined by four parameters: (i) N is the total number of miRNAs used to predict the target; (ii) K is the The number of miRNAs interacting with the selected gene; (iii) n is the number of miRNAs that interact with the candidate ceRNA of the selected gene; (iv) the miRNA number common between the two genes [19]. This test used the following formula to calculate the value:

3. Results

3.1. Identification of Differentially Expressed CircRNAs

To investigate the expression profile of circRNAs in AP, we used microarray analysis to perform circRNA expression profiling in the blood of patients with MAP and matched normal humans. The box plot visualized the dataset distribution of circRNAs. After normalization, Figure 3(a) depicts that the log 2 ratios in the three pairs of samples are almost identical. At the same time, we used a scatter plot of the circRNA expression profile to evaluate the change between the two groups (Figure 3(b)), where the values of the X and Y axes are the average normalized signal values (log 2 scaling) of the sample set. And the green line is the fold line, and the circRNA above the top green line and below the bottom green line indicates that the circRNA change between the two samples is more than 1.5 times. Next, volcano maps were used to identify differentially expressed circRNAs that were statistically significant between the two groups and are shown in Figure 3(c). Among them, the vertical lines correspond to 1.5 times up and down, and the horizontal lines represent . The red dot indicates upregulated circRNAs while blue represents downregulated circRNAs. We found 51 differentially expressed circRNAs in patients with MAP (FC ≥ 1.5 and values ≤0.05) compared with normal subjects. Among them, Table 3 shows 25 upregulated circRNAs and 26 downregulated circRNAs are described in Table 4. Next, we clustered all the different circRNAs to characterize the expression pattern of circRNA. The results are exhibited in Figure 4. Among them, red stand for elevated and green represents downregulated circRNAs.

3.2. qPCR Verification

To verify differentially expressed circRNAs in AP, we made qPCR. And we selected three upregulated circRNAs including hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 and three downregulated circRNAs including hsa_circRNA_002532, hsa_circRNA_059665, and hsa_circRNA_104156. We validated the expression of three downregulated circRNAs in MAP and normal humans. And the results are shown in Figure 5. values are 0.008 (hsa_circRNA_002532), <0.0001 (hsa_circRNA_059665), and <0.0001(hsa_circRNA_104156), respectively.

To further explore the relationship between the three upregulated circRNAs and disease severity, we compared the expression levels in healthy individuals, MAP, and SAP patients, respectively. As shown in Figure 6, we found that as the condition worsened, the expression levels of hsa_circRNA_101015 ( value was 0.0003, <0.0001, and <0.0001, respectively), hsa_circRNA_101211 ( value was 0.0014, 0.0478, and <0.0001, respectively), and hsa_circRNA_103470 ( value was 0.001, <0.0001, and <0.0001) also increased significantly. Then we validated the expression levels of three upregulated circRNAs in MAP patients and MAP after treatment. The results are demonstrated in Figure 7. After treatment, the expression levels of hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 in MAP were significantly declined ( value <0.0001).

3.3. ROC Analysis of hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 in Patients with AP

The ROC curve was constructed to assess the diagnostic significance of the three elevated circRNAs (Figure 8). RUCs of hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 MAP patients with healthy people were 0.768 (95% CI, 0.651–0.886, , power = 0.97119), 0.731 (95% CI, 0.605–0.857, , power = 0.90168), 0.770 (95% CI, 0.653–0.887, , power = 0.97340), respectively. RUC of the three circRNAs combination was 0.838 (95% CI, 0.738–0.937, , power = 0.99954). The above results show that hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 can be used as biomarkers to diagnose AP at early stage.

3.4. The Construction of ceRNA Network

Recent evidence suggests that circular RNA plays a crucial role in the regulation of miRNA-mediated gene expression regulation by isolating miRNAs. Their interaction with disease-associated miRNAs suggests that circular RNA is important for disease regulation [20]. To assess the potential function of circRNA, we investigated potential miRNAs that bind to circRNA. Then we constructed a ceRNA network for the six differentially expressed circRNAs. In the ceRNA network, there are 241 nodes and 831 edges (Figure 9). Among them, for all nodes, red represents microRNAs, light-blue color represents protein_coding RNAs, and at the same time, brown color represents circular RNAs. Considering all edges, T-shape arrow represents directed relationships, while edges without arrow stand for undirected relationships. Specific details are shown in Supplementary Table 1.

4. Discussion

AP is an inflammatory process of the pancreas and has become an increasingly common clinical disease. In the second or third week of the disease, 40–70% of patients develop infectious necrosis and are the leading cause of late death [21]. In the United States, AP accounts for $2.5 billion in medical expenses, and the number of hospital admissions is about 275,000 per year [22]. Because of the high mortality rate of 20–30% in severe cases of AP, early detection of patients who may need to be transferred to the intensive care unit (ICU) is critical. For patients with AP, accurate diagnosis, appropriate triage, high-quality supportive care, monitoring and treatment of complications, and prevention of recurrence are also critical [23]. The evaluation of biomarkers helps to further improve the identification of high-risk patients. In general, the clinical treatment decisions for AP depend primarily on the severity of the condition. MAP is usually self-limiting and does not cause death. However, SAP can develop rapidly, leading to multiple organ failure and becoming life threatening. Therefore, the treatment of MAP patients is relatively simple and requires only a short hospitalization, while the treatment of SAP patients usually involves intensive care. Whether it is MAP or SAP, quick and accurate diagnosis is still very difficult.

Accurate diagnosis of AP may allow for effective treatment to begin earlier. At present, the accuracy of different scoring systems is not high, and a unique model is needed [24]. There are several laboratory tests such as blood urea nitrogen, creatinine, and hematocrit [25]. However, there is virtually no laboratory test that consistently and accurately predicts the severity of AP earlier [26, 27]. Many predictive systems use CT findings, but CT evidence of SAP lags behind clinical findings, and early CT studies may underestimate the severity of the disease. The scoring system is complex and cumbersome; therefore, these scoring systems are not a substitute for the clinician’s experience assessment.

In recent years, circRNA has received wide attention as a new class of endogenous and regulatory noncoding RNAs. At the same time, with the widespread use of RNA sequencing (RNA-seq) technology and bioinformatics prediction, a large number of circRNAs have been identified. So far, no study has examined the role of circRNA in AP. Considering that circRNA is involved in a variety of diseases, it is necessary to explore differences in the expression of circRNA in patients with AP. Microarrays are an effective tool for analyzing circRNA. Therefore, in our study, we made full use of microarray technology and selected three patients with MAP and three healthy individuals. Our aim was to explore the differential expression of circRNA in the patients’ blood with AP to diagnose AP patients as soon as possible. As a result, we found 25 upregulated circRNAs and 26 downregulated circRNAs. The vast majority of circRNAs have not been studied and require more in-depth exploration.

To further validate the six circRNAs in AP, we performed qPCR analysis. The results of qPCR studies provide novel biomarkers for molecular diagnosis and evaluation in AP. Three circRNAs including hsa_circRNA_002532 hsa_circRNA_059665, and hsa_circRNA_104156 were downregulated in MAP compared with healthy individuals. Three upregulated circRNAs including hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 significantly increased expression levels as the condition worsens. Furthermore, the expression levels of hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 were significantly reduced after treatment. This result suggests that the above three circRNAs may be involved in the development of AP and show promise as potential biomarkers for AP. To evaluate the specificity and sensitivity of the hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470, we performed a ROC curve for further validation. The RUC was 0.768, 0.731, and 0.770, respectively, and the combined RUC values of the three circRNAs were 0.838. Therefore, the three circRNAs can be used as diagnostic biomarkers for AP with high sensitivity.

However, the pathogenesis of AP remains unknown. Therefore, the 6 differentially expressed circRNAs were selected to construct a ceRNA network including three upregulated hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 and three downregulated hsa_circRNA_002532 hsa_circRNA_059665, and hsa_circRNA_104156. In the ceRNA network, we found a complex regulatory network for the above six circRNAs. However, the above six differentially expressed circRNAs excluding hsa_circRNA_104156 have not been reported. Considering that circRNAs can act as sponges of miRNAs, it is of importance to explore miRNAs in AP. More importantly, miRNAs play a critical role in the progression of AP. For example, a previous study found that miR-92b, miR-10a, and miR-7 are downregulated in the blood of patients with AP and can be used to distinguish between patients with AP and healthy cases. In addition, the expression level of miR-551b-5p distinguishes between MAP and SAP [28]. However, these six circRNAs cannot interact with above these miRNAs. In the ceRNA network, hsa_circRNA_101015 could become sponges of several miRNAs such as hsa-miR-135a/b and hsa-miR-543. It has been reported that hsa-miR-135 was upregulated in the serum of mice with acute pancreatitis [29]. Furthermore, it has reported that hsa-miR-135a/b binding polymorphism can reduce the expression of CD133 and reduce the risk of lung cancer and improve the prognosis of patients [30]. Abnormalities in hsa-miR-543 can lead to increased proliferation of multiple cancers, such as esophageal cancer, prostate cancer, and osteosarcoma [3133]. In the ceRNA network, hsa-miR-135 could be a target of hsa_circRNA_101015. It has been found that hsa-miR-216a was downregulated in the serum of mice with acute pancreatitis [34]. Elevated serum hsa-miR-122 has been considered as a noninvasive marker for acute pancreatitis, which could become a target miRNA of hsa_circRNA_059665 as shown in the ceRNA network [35]. Moreover, hsa-miR-148a has been validated to be lowly expressed in acute pancreatitis [36]. Our results showed that hsa-miR-148a could be a potential target miRNA of hsa_circRNA_059665. Furthermore, hsa-miR-148a could inhibit autophagy through IL-6/STAT3 axis in AP in acute pancreatitis. According to the ceRNA network, hsa_circRNA_104156 could be considered as a sponge of hsa-miR-148a.

Therefore, we inferred that the six circRNAs could participate in pathogenic process of AP. However, specific research mechanisms need to be explored. In the future, we would continue to research the molecular mechanisms of hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 in AP.

5. Conclusion

After microarray analysis and qPCR, we found three downregulated circRNAs including hsa_circRNA_002532, hsa_circRNA_059665, and hsa_circRNA_104156 and three upregulated circRNAs including hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 in the blood of AP patients. The ceRNA network revealed that the six circRNAs could play a critical role in AP. Therefore, elevated hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 can be used as biomarkers in the blood to diagnose AP at early stage. The molecular mechanisms of elevated hsa_circRNA_101015, hsa_circRNA_101211, and hsa_circRNA_103470 in the pathological processes of AP would be explored in our further experiments.

Abbreviations

AP:Acute pancreatitis
circRNAs:Circular RNAs
qPCR:Quantitative polymerase chain reaction
ROC:Receiver operating characteristic
SAP:Severe acute pancreatitis
MAP:Mild acute pancreatitis
miRNAs:MicroRNAs
AUCs:Area under the ROC curves.

Data Availability

The data used to support the findings of this study are included within the supplementary information file.

Disclosure

Chang Liu and Xuan Zhu are co-first authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Chang Liu and Xuan Zhu contributed equally to this article.

Supplementary Materials

Supplementary Table 1: the details of the ceRNA network. (Supplementary Materials)