Dysregulated miRNA Associated with Transcription Factors of Insulin Gene Expression in Chronic Pancreatitis

Background/Aim: MicroRNAs with regulatory functions in gene expression are implicated in different diseases. The present study investigated differentially expressed miRNAs that possibly influence transcription factors involved in insulin gene expression in Chronic Pancreatitis (CP) employing bioinformatics approaches. Methods: Pancreatic tissues were collected from CP patients undergoing partial pancreatectomy (n = 16) and controls (n = 15) undergoing resections for non-pancreatic malignancies. MiRNA profiles obtained using microarrays were validated by qRT-PCR. Target search involving miRWalk and TarBase as well as functional annotation employing KEGG (Kyoto encyclopedia of genes and genomes) and DAVID (Database for Annotation) databases were performed. Ingenuity pathway analysis (IPA) was used to construct networks relating miRNAs to their target genes. mRNA and proteins related to insulin gene transcription factors and hormones were evaluated by qRT-PCR and western blotting followed by confirmation upon immunofluorescent staining. Results: Microarray data revealed 10 up-regulated and 15 down-regulated miRNAs in CP as compared to controls (Log2 FC > 2). Bioinformatic analysis showed 8399 target genes and KEGG pathway analysis suggested a role for the dysregulated miRNAs in modulating cytokine signaling, fibrosis, JAK-STAT signaling and insulin synthesis. IPA analysis suggested a simplified network attributing dysregulated miRNAs to NFκBdependent cytokine signaling. Further, associations could be noted between miRNA 200b with Maf A, 138-1 with Neuro D and 27b with FoxO1. Decreases in mRNA levels of Pdx1, Neuro D and increases of Maf A and FoxO1 transcription factors could be noted (P < 0.01) in CP. These results were confirmed by western blotting and immunofluorescence staining. Conclusion: How to cite this paper: Murali Manohar, K., Sasikala, M., Pavan Kumar, P., Rao, G.V. and Nageshwar Reddy, D. (2016) Dysregulated miRNA Associated with Transcription Factors of Insulin Gene Expression in Chronic Pancreatitis. Open Journal of Endocrine and Metabolic Diseases, 6, 205-227. http://dx.doi.org/10.4236/ojemd.2016.610026 Received: August 20, 2016 Accepted: October 18, 2016 Published: October 21, 2016 Copyright © 2016 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access K. Murali Manohar et al. 206 Our results identified dysregulation of miRNAs 138-1, 27b and 200b which were found to be associated with insulin gene transcription factors Neuro D, FoxO1 and Maf A respectively.


Introduction
Chronic pancreatitis (CP) is a progressive inflammatory disorder ultimately culminating in exocrine and endocrine insufficiency contributing to associated maldigestion, abdominal pain and clinical diabetes.Endocrine insufficiency and clinical diabetes are popularly believed to occur as a result of fibrosis and β-cell destruction (apoptosis) in CP [1].Several investigators demonstrated insulin secretory defects in patients with long standing CP [2].It is held that insulin secretory defects arise initially as a result of β-cell dysfunction and overt clinical diabetes manifests from β-cell destruction much later in the course of the disease.Such findings suggested a need to study the mechanism of β-cell dysfunction under conditions of chronic inflammation in CP.Diabetes mellitus, secondary to chronic pancreatitis, is now categorized as Type 3C DM [3].The complex pathophysiology of CP involving secondary diabetes presents distinct clinical and laboratory features, in contrast to autoimmunity and insulin resistance associated respectively with Type 1 and Type 2 DM.Currently it is reported that 5% -10% in western population and 15% -20% of diabetics in South East Asia have Type 3C DM [4].Despite vast knowledge on Type 1 and Type 2 DM, Type 3C diabetes is relatively the less studied form of diabetes.Interestingly, our earlier reports demonstrated reduced islet response to glucose stimulation even in non-diabetic CP patients [5].
MicroRNAs are small, endogenously expressed noncoding RNAs (~21 -25 nucleotides) known to influence basic cellular functions by regulating gene expression.miRNAs bind to the 3'-untranslated region of mRNAs and lead to either transcriptional repression or mRNA degradation [6] [7].Altered miRNA profiles have been implicated in different diseases such as diabetes, inflammatory diseases and various malignancies [8] [9] [10].A few studies have also established basic functional role of miRNAs in islets (hormone secreting cell clusters) including miRNA375 in islet development, miRNA9 in regulating insulin secretory response, miRNA124 in regulating intracellular signaling and miRNA30d in induction of insulin production by targeting MAP4K4 [11] [12] [13] [14].Further, it was shown that miRNA21, 34a and 146 act as negative regulators of insulin signaling [15].Increased expression of miRNA29 family in pre-diabetic NOD mice was shown to be associated with cytokine mediated β-cell dysfunction [16].
Islet response to circulating nutrients involves tight coordination and regula-tion between synthesis and secretion of pancreatic hormones.Expression of insulin gene in response to secretagogues is known to be controlled by various transcription factors such as Pdx1, neuro D, Maf A that bind to enhancer elements in the promoter region of insulin gene [17].Even though several investigators have reported altered miRNA profiles in CP and pancreatic cancer [18], relatively much less is known about the influence of aberrant miRNA(s) on β-cell function in CP; putative interactions between miRNAs and mRNAs of transcription factors of insulin gene expression have also not been elucidated under inflammatory conditions of CP.The established role of miRNAs in inflammation [9] as well as the reported role of miRNAs in β-cell dysfunction in Type 2 DM [19] and in prediabetic NOD mice [16] have led us to hypothesize that aberrant miRNA expression in CP may also have an effect on insulin gene transcription, which may contribute to β-cell dysfunction.
The present study was thus conducted to (a) evaluate altered pattern of miR-NAs in CP and identify their putative targets including genes coding for islet hormones and associated transcription factors; (b) identify association of aberrant miRNAs with signaling pathways involved in insulin gene expression and establish possible network(s) between miRNAs and target genes in CP using bioinformatics approaches.

Patients and Methods
Patients: Pancreatic tissues were obtained either from CP patients (n = 16) who had undergone partial pancreatectomy/lateral pancreaticojejunostomy for pain relief or from patients who had undergone pancreaticoduodenectomy for ampullary/periampullary pathology without any history of CP (controls; n = 15), Histological confirmation of CP involved H&E staining of the resected specimens by a single pathologist experienced in pancreatic pathology who was blinded to the patient groups.Patients with CP having pancreatic malignancies, endocrine tumors, acute on CP were excluded.The study was initiated after the protocols were approved by the Institutional Ethics Committee of Asian Institute of Gastroenterology and all the patients had given informed consent.

MicroRNA profiling and validation of dysregulated miRNAs
Total RNA from the pancreatic tissues was extracted and purified using miRvana RNA isolation kit following the manufacturer's instructions (Applied Biosystems/Ambion, Austin, TX).Qualitative and quantitative assessment of total RNA was obtained using Agilent 2100 bioanalyzer as per the prescribed protocol of the manufacturer (Agilent Technologies, USA).Good quality RNA samples with RNA Integrity Number (RIN) > 7 were subjected to microarray analysis for miRNA profiling and to qRT-PCR for validation of identified miRNA.Micro-RNA profiling, microarray data analysis and validation with qRT-PCR were performed as described earlier [20].The raw data were submitted to the Ar-rayExpress database (accession ID: E -MTAB-3882).

Target gene prediction and pathway analysis of dysregulated miRNAs
Putative target genes for the differentially expressed miRNAs were identified and retrieved from miRWalk 2.0 database [21], which compiles and compares the predictions of other databases.Target genes that were experimentally validated for dysregulated miRNAs were retrieved from Diana Lab Tarbase v.6 [22].
The putative and validated targets of dysregulated miRNAs were subjected to KEGG pathway annotation using DAVID functional annotation tool [23].A two sided Fisher's exact test and Chi-square test were used to classify the enrichment (Re) of pathway and false discovery rate (FDR) was calculated.A p value of <0.01 and an FDR of <0.05 were chosen for identification of enriched pathways.

IPA analysis to deduce network of altered miRNAs and target genes
Differentially expressed miRNAs and the respective target genes were uploaded into QIAGEN's Ingenuity ® Pathway Analysis (IPA ® , QIAGEN Redwood City), which uses pooled genes as queries and then retrieves all other gene objects related to query genes/miRNAs stored in knowledge base.Based on the related genes a set of networks were generated with different scores.Networks that are considered to be best fit (implied by the score derived from the p-value) indicate the maximum likelihood of users' list of genes in a network being found.

Evaluation of the network of transcription factors of insulin gene in CP
One μg of total RNA was subjected to cDNA synthesis using oligodT, and M-MuLV Reverse transcriptase enzyme (Invitrogen).qRT-PCR for islet hormones (insulin, glucagon, somatostatin, polypeptide) and transcription factors (Pdx1, NeuroD, MafA, FoxO1) was performed using SYBR Green PCR Master Mix employing Step One Real Time PCR (Applied Bio systems).Primers were designed spanning multiple exons using IDT (Integrated DNA Technologies, Singapore) software available online.All RT-PCRs included no template controls and RT minus controls.For the relative quantification, 10 μL of master mix was prepared with 0.5 μL RT product, 0.2 μL volume each of appropriate forward and reverse primers, 5 μL SYBR green PCR master mix and 4.1 μL RNase free water.Samples were run in duplicates and data were normalized to human GAPDH.mRNA levels were relatively quantified by using ΔΔCT method (List of used Primers is given in Supplementary Table S1).
Western blot analysis and immunofluorescent staining for expression of Pdx1, Neuro D and MafA Proteins were isolated from minced pancreatic tissue using 500 μl of RIPA lysis buffer with protease cocktail inhibitor (Sigma Aldrich, St Louis, USA) and probing the western blots either with rabbit anti-Pdx1 or rabbit anti-NeuroD or rabbit anti-MafA (Invitrogen, Shanghai, China) or anti β-actin (Santa cruz Biotechnology, USA) as per standard protocols [24].
Paraffin embedded pancreatic tissues obtained from CP patients and controls were subjected to immunofluorescent staining.The sections were permeabilized with 0.2% Triton X-100 for 5 minutes at 4˚C and stained with primary antibo-dies for hormones; insulin and glucagon (guinea pig anti-insulin and mouse anti-glucagon at 1:200; Sigma Aldrich), somatostatin, polypeptide: (Invitrogen, Shanghai, China) and transcription factors (pancreatic duodenal homeobox; PDX-1, (rabbit anti PDX-1, rabbit anti MafA and Mouse Neuro D at 1:100 dilution) followed by overnight incubation at 4˚C.After repeated washings with PBS, sections were treated with secondary antibodies tagged with fluorescent dyes guinea pig Alexa 488 for insulin and goat anti mouse Alexa 546 for glucagon at 1:200 dilutions; Invitrogen, Shanghai, China (goat anti rabbit Alexa 546 for PDX-1 and goat anti guinea pig Alexa 488 for insulin; Invitrogen, China, donkey anti goat Alexa 546 for Neuro D).Fluorescence images were captured using a Bioimager (CARV II, BD BioSciences, CA, USA).
Statistical analysis: The data were analyzed using statistical package for social sciences (SPSS version 20.0, Chicago, USA).A two tailed students P value of < 0.05 was considered statistically significant.

Microarray experiments demonstrated differential expression profiles of miRNA in CP
Significance analysis of Microarray (SAM) identified 52 differentially expressed miRNAs in CP tissues as compared to control tissues (Supplementary Table S2).Hierarchical clustering of the normalized expression data showed two distinct clusters of miRNAs (Figure 1(a)).Further refinement with t-test and Bonferroni adjustment revealed a total of 25 differentially expressed miRNAs that includes 10 up-regulated and 15 down-regulated miRNAs with high significance (P < 0.05).The differentially expressed miRNAs detected by highest logarithmic fold changes between CP and adjacent control tissues are shown in

Target gene prediction of aberrant miRNA revealed fibrosis and insulin signaling pathways in CP
In silico analysis was performed to predict target mRNAs of the 25 differentially expressed miRNAs.The number of mRNA targets predicted for each of the differentially expressed miRNAs varied substantially.A three-way intersection of the three algorithms (Targetscan, miR Walk and miRNaDa) revealed 9577 target transcripts corresponding to 8399 unique official gene symbols.Functional annotations of target genes using DAVID identified 12/53 significant pathways of which TGFβ, T cell receptor and cytokine-cytokine receptor interaction pathways were related to inflammation and fibrosis.Similarly, pathways related to insulin synthesis, calcium signaling and Type2DM could be related to β cell specific pathways as shown in the Table 1.S4.

KEGG, ingenuity pathway analysis and putative networks in CP
Glucose responsive MafA expression increased and Pdx1, Neuro D expression decreased in CP mRNA levels of the transcription factors of insulin gene were noted to be al-

Discussion
Micro RNAs are known to play a pivotal role in controlling basic cellular functions by regulating gene expressions.We conducted this study to explore whether there is any dysregulation of pancreatic miRNAs under the inflammatory conditions prevalent in CP and to identify putative links between differentially expressed miRNAs and transcription factors that coordinately regulate insulin gene expression contributing to secondary diabetes occurring in the disease.Our results demonstrate a network of three miRNAs (miR200b, 138 and 27b) regulating insulin gene transcription factors along with cytokines in CP.
The dysregulated profiles of miRNA obtained on microarray and validated on qRT-PCR indicate that inflammation in CP is associated with altered miRNA expression.The observed expression profiles of dysregulated miRNA in the present studies were similar to earlier studies which examined miRNA patterns in the progression of chronic pancreatitis to pancreatic ductal adenocarcinoma [18] which reported up-regulation of miR-339 and down-regulation of miR-34, miR-181, miR-375, miR125b and miR-99.Further, the down regulation patterns (Supplementary Table S3) were akin to those reported in a study comparing the abundance of circulating and tissue miRNAs in CP [25].Subtle differences  between these two studies in miRNA expressions may be ascribed to differences in etiopathogenesis of CP (largely alcoholic in Western countries and idiopathic in tropical countries).Although these studies focused on miRNA expression patterns during progression of CP to pancreatic cancer, no study has been conducted to find an association between dysregulated miRNA and transcription factors involved in regulating insulin gene expression in CP.Such an identification of networks relating aberrant pancreatic miRNA expression to insulin gene transcription is of relevance in unraveling the mechanism of endocrine dysfunction in CP.
Target gene prediction analysis of differentially expressed miRNAs using bioinformatics approaches (miRbase, Target scan and miRwalk) yielded several transcripts (8399) related to various signaling pathways with ≈200 overlapping targets f or each miRNA (Figure S1).It is known that paired miRNA genes present within 10 kb distance (usually named as miRNA clusters) are usually co-expressed and have functionally significant role in repressing the target genes [26].Among the 25 differentially expressed miRNA, a few clusters of miRNA (miR-29a/29b/29c, miR-200a/200b/200c/141, miR-216/217, miR-30a/30b/30c, miR-27a/27b/24/23b) were noted to be associated with target genes of specific pathways involved in fibrosis, inflammation, epithelial-mesenchymal transition (EMT), oxidative stress and dysregulated insulin signaling.Multiple target genes of miR-29 cluster that include collagens (COL7A1), fibrillins (FBN1), and matrix metalloproteins (MMPs) identified in this study were earlier implicated in cardiac, pulmonary and renal fibrosis [27] [28] [29].It is also known that transforming growth factor (TGF)-β-mediated down-regulation of miR-29 would enhance fibrosis and role of TGF-β in fibrogenesis during progression of CP is well established [30].In addition, Wang et al. showed that the miR-200 cluster (mir200 a, b and mir141) contributes directly to fibrosis by regulating of TGF-β2 [31].Reporter gene assays have also confirmed TGF-β2 to be the direct target of miR-141/miR-200a [32].Our results are in agreement with these studies suggesting that down regulation of miR200 and miR29 may be contributing to fibrosis in CP [33].Decreased expression of miR-216/217 cluster which has been shown to be involved in EMT also suggests its role in fibrogenesis [34].Decreased expression of miRNA-27 cluster is known to target PPARγ ligand, which modulates NF-κB dependent pro-inflammatory cytokine production and pancreatic stellate cell activation [35]; this in turn would contribute to increased oxidative stress and inflammation in CP.It is noteworthy that evidences are increasingly indicating that prolonged oxidative stress may be a contributing factor for pathogenesis and progression of CP [36].Results of the present study as well as an earlier report by our group [37] also denote down-regulation of miR-30 family in CP; the miR-30 family is implicated in induction of insulin synthesis by targeting MAP4K4 regulation [38] [39] [40] these observations indicate that miRNA clusters might play a crucial role in inflammation, fibrosis and aberrant β-cell functions associated with CP.
Although comprehensive search was conducted for target gene prediction, networks were generated restricting the analysis to target genes that are associated with oxidative stress, cytokine signaling, Inflammation and insulin signaling since they are highly relevant to the progression of CP and manifestation type 3CDM.An attempt was also made to understand the intricate relations between the dysregulated miRNAs, target genes, and β-cell dysfunction in CP employing Ingenuity Pathway Analysis.A curated network was generated involving miRNAs, target genes of cytokine signaling, oxidative stress, insulin synthesis, β cell function.Interestingly, the resultant network revealed a web consisting of various proinflammatory cytokines (IL6, IFNγ, IL1β, and TNFα), JAK2, NFκB signaling pathways, statins (STAT1 and STAT3), effectors of β-cell function (MafA, NeuroD and Pdx1), β-cell development (MafB), β cell oxidative stress protective transcription factor (FoxO1) and insulin.Earlier results from our laboratory and that of others have shown increased expression of proinflammatory cytokines [41] in CP and a role for NFκB [42] and oxidative stress was implicated in the inflammation in CP [43].The network established in the present study has also identified links between cytokines and miRNAs (IL6 -miR191; IL10-miR27; IFNγ miR125b and 29C; IL1 miR141; TNFα miR 130 and miR19b) involving putative links between miR27-IL10 [44], miR29-IFNγ [45], TNFα-miR130 [46], that were previously established.In the network, NFκB and STAT proteins have network links with IFNγ, IL1A, IL1B TNFα, IL6 and IL10), while binding site in the Pdx1 promoter [47].FoxO1 and Pdx1 have also been reported to show mutually exclusive nuclear localization [48] [49] [50].Increased oxidative stress, arising from inflammation and decreased miR 27b could result in raised FoxO1 and diminished Pdx1 expression in CP as observed in this study could be due to raised FoxO1.Transgenic mice with FoxO1 gain of function in the pancreas exhibited glucose intolerance [51].Conversely, FoxO1 ablation in β-cells resulted in increased insulin secretion [52].Further, in vitro studies showed suppression of miRNA-27 resulted in increased FoxO1 protein levels and a consequent decrease in cell number [53].Hence we believe that increased FoxO1 resulting from increased oxidative stress leads to dysregulation in β cell functions.
Insulin gene transcription is known to be coordinately regulated by transcriptional activators such as Neuro D, Maf A and Pdx1 [54] [55] [56] [57].To our knowledge, this is the first report to show increased expression of Maf A in CP both at mRNA and protein levels.We had earlier shown that reduced expression of Pdx1 is associated with decreased β-cell function [5].It is well established that Pdx1, Neuro D and Maf A bind to the insulin promoter to initiate insulin gene transcription.Even though FoxO1 is known to activate gene expression (Maf A) in response to oxidative stress via acetylation [58], decrease in Pdx1 and Neuro D in CP might result in loss of coordinated transcription of insulin gene resulting in decreased β cell function.Despite acetylated FoxO1 acting as a protective factor to β cells in response to oxidative stress under acute metabolic change, it fails to preserve β cell function in chronic conditions [58].

Conclusion
In summary, present study identified dysregulated miRNAs and indicated the pathways they target in CP.Dysregulated miRNA target genes are significantly     Submit or recommend next manuscript to SCIRP and we will provide best service for you: Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.A wide selection of journals (inclusive of 9 subjects, more than 200 journals) Providing 24-hour high-quality service User-friendly online submission system Fair and swift peer-review system Efficient typesetting and proofreading procedure Display of the result of downloads and visits, as well as the number of cited articles Maximum dissemination of your research work Submit your manuscript at: http://papersubmission.scirp.org/Or contact ojemd@scirp.org

A
total of 16 patients with CP (male n = 12, mean age 34.44 ± 15.85 years), as confirmed by clinical, endoscopic and imaging (CT, MRI, EUS and ERCP) investigations were recruited for the study.All of the CP patients had abdominal pain for 36 ± 21 months with radiological findings revealing the presence of intraductal calculi and main duct dilatation.Histopathological examination confirmed intra and interlobular acinar atrophy and fibrosis in these patients.The control group (n = 15; male = 10, mean age 50.53 ± 9.53 years) were ascertained to be without CP.

Figure 1 .
Figure 1.Differentially expressed miRNAs in CP tissues vs. normal pancreatic tissues.(a) Clustered heat maps illustrate miRNA differential expression profiles across 5 CP samples and 10 controls; (b) 25 miRNAs significantly (P ≤ 0.05) differentially regulated (fold change ± ≥2) in CP patients compared with normal pancreatic tissues are listed in order of their fold change.The accession number is for the mature miRNA sequence; Chr, Chromosome; (c) Quantitative real-time PCR analysis of dysregulated miRNA expression in CP (n = 16) compared with normal pancreatic tissues (n = 15).The miRNA abundance was normalized with U6SnRNA as internal control.Relative expression was analyzed by 2^− ΔΔCT method.

Fibrosis
, inflammation, insulin signaling pathways and malignancy were found to be enriched on pathway analysis utilizing KEGG database.IPA analysis of 25 dysregulated miRNAs and respective target genes (predicted /validated) generated ten functional networks.Four of these networks that are relevant to CP include those involved in (i) inflammatory disease, fibrosis and cancer (ii) cytokine signaling, insulin synthesis and β cell function (iii) inflammatory response and cell morphology and (iv) immunological disease and cell cycle.Since CP is an inflammatory disease wherein islets have to function in an inflammatory milieu, we conducted further analysis relating miRNA expression to cytokine signaling and insulin gene transcription.Consequently, the final network generated involved cytokines, transcription factors and miRNA.It depicts insulin gene transcription to be coordinately regulated by the three crucial transcription factors (Pdx1, NeuroD, MafA) and FoxO1 which are the confirmed targets of mir19, 130b, 30 family for Neuro D and miR 27 family for FoxO1 respectively.None of the altered miRNA targeted Pdx1 directly.This network also identifies a putative link between miR 200 family and MafA.Inflammatory cytokines (IL1 family, IL10, IFNγ, IL6) were found to be targets of differentially expressed miRNAs 191, 27, 29 and 191 respectively as reflected in the network (Figure 2).Validated targets of networks are shown in Supplementary Table

Figure 3 .
Figure 3. Islet Hormones and transcription factors in CP.(a) Altered expression patterns of islet hormones and transcription factors in CP tissues identified by qRT-PCR.The mRNA transcript abundance was normalized to that of GAPDH.Relative expression was analyzed by 2^-ΔΔCT method; (b) Protein expression of β-cell specific transcription factors (NeuroD, PDX1 and MafA) in CP by western blot normalized to β-Actin as reference control.The bar graph represents the expression levels as ratio of target protein expression with β-actin; (c) Immunofluorescence staining of islets in CP compared to control tissues confirms altered expressions of islet hormones and transcription factors.Islets from control group and CP patients were stained with fluorescent antibodies to visualize insulin (green), PDX-1, MafA, NeuroD, Glucagon, Somatostatin, Polypeptide (red), and nuclei (blue).All images were obtained at original magnification (400×).

Figure 4 .
Figure 4. Inverse and negative correlation between miRNA and transcription factors.(a) Pearson correlation analysis showed negative correlation of miRNA-target gene pairs (a) miR-138 with NeuroD (r = −0.617,P = 0.012); (b) miR-200b with MafA (r = −0.597,P = 0.015); (c) miR-27b with FoxO1 (r = −0.571P = 0.021).In Graphs aligned towards left vertical axis indicates delta Ct values and are opposite to actual expression values obtained after normalizing with reference internal control genes (U6SnRNA for miRNAs and GAPDH for mRNAs).In Graphs aligned towards right r represents Pearson correlation coefficient; (d) Sequence matching of the predicted binding sites shows miR-138 binds to the 3' UTR region of NeuroD1, miR-200b/200c binds to MafA and miR-27b/27a binds to FoxO1 mRNA.The miRNA seed sequences are highlighted in red and sequences complementary to seed region in 3'UTRs of mRNA are highlighted in green.NeuroD1, Neurogenic differentiation1; MafA, Musculoaponeurotic Fibrosarcoma Oncogene Homolog A; FoxO1, Forkhead box protein O1.
transcription factors and miRNAs (NeuroD-miR138, 19b, 130a, 30c, 148b; Ma-fA-200b; MafB-miR130b, miR29; FoxO1-miR27a; STAT3-miR29b, 130a; STAT1-miR30c) were associated.Interestingly, none of the differentially expressed miRNA targeted Pdx1 (crucial β cell transcription factor).Results of this study established possible links between miRNAs, cytokines, oxidative stress and transcription factors of insulin gene and β-cell functions in CP.Further analysis using target scan, matching the sequences of differentially expressed miRNAs with the target genes identified sequence complementarities between transcription factors of insulin gene and miRNAs.Amongst the 25 differentially expressed miRNA, three miRNAs 200b, 138-1 and 27b were identified to have binding sites in 3'UTRs of MafA, NeuroD and FoxO1 respectively.In addition to this, expression levels of these miRNAs were negatively correlated with Maf A, Neuro D and FoxO1, suggesting a role for these miRNA in regulating the expression of the transcription factors of insulin gene.Decreased mRNA levels of Pdx1, Neuro D, and increased expression of Maf A were confirmed upon probing the proteins on western blot analysis.Although previous studies have demonstrated differential expression of 200b, 27b, 138-1, their role in β cell function were not studied under conditions of inflammation.In vitro assays by earlier investigators demonstrated that FoxO1 inhibits Pdx1 transcription in β cell cultures by competing with FoxA2 (transcription factor) for a common enriched in pathways related to fibrosis and EMT (miR-29, 200, 217 miRNA clusters), inflammation/cytokines (miR-27 family) and insulin signaling pathways (miR-30 family and miR-375).Putative links were also identified upon bioinformatic analysis between miRNAs (miR-138, 200b, 27a/27b) and transcription factors of insulin gene (Neuro D, Maf A, and FoxO1) suggesting that miRNAs, such as miR138-1, 27b and 200b, might be important regulators of crucial transcription factors involved in insulin gene expression in CP.

22045263 ContinuedFigure 2 Figure 2
Figure 2 miR375 decrease JAK2 Not known miR375 -JAK2 miR-375 may function as a tumor suppressor to regulate gastric cancer cell proliferation potentially by targeting the JAK2 oncogene, implicating a role of miR-375 in the pathogenesis of gastric cancer JAK2 mediated pathways were shown to stimulate replication and survival of β-cells 20548334, 22045263

Figure S1 .
Figure S1.Number of mRNA transcripts identified by miRwalk for each differentially expressed miRNA.

Table 1 .
Putative targets and pathways that are associated with dysregulated microRNAs in chronic pancreatitis.

Table S3 .
Comparison of dysregulated miRNA in CP with other studies.

Table S4 .
Interactions of dysregulated miRNAs and target genes: relevance to CP.