Prospective insight into the role of benzyl propylene glycoside as a modulator of the cGAS-STING signaling pathway in the management of nonalcoholic fatty pancreas animal model
Biological Research volume 56, Article number: 11 (2023)
Nonalcoholic fatty pancreatitis (NAFP) is one of the metabolic syndrome manifestations that need further studies to determine its molecular determinants and find effective medications. We aimed to investigate the potential effect of benzyl propylene glycoside on NAFP management via targeting the pancreatic cGAS-STING pathway-related genes (DDX58, NFκB1 & CHUK) and their upstream regulator miRNA (miR-1976) that were retrieved from bioinformatics analysis.
The rats were fed either normal chow or a high-fat high-sucrose diet (HFHS), as a nutritional model for NAFP. After 8 weeks, the HFHS-fed rats were subdivided randomly into 4 groups; untreated HFHS group (NAFP model group) and three treated groups which received 3 doses of benzyl propylene glycoside (10, 20, and 30 mg/kg) daily for 4 weeks, parallel with HFHS feeding.
The molecular analysis revealed that benzyl propylene glycoside could modulate the expression of the pancreatic cGAS-STING pathway-related through the downregulation of the expression of DDX58, NFκB1, and CHUK mRNAs and upregulation of miR-1976 expression. Moreover, the applied treatment reversed insulin resistance, inflammation, and fibrosis observed in the untreated NAFP group, as evidenced by improved lipid panel, decreased body weight and the serum level of lipase and amylase, reduced protein levels of NFκB1 and caspase-3 with a significant reduction in area % of collagen fibers in the pancreatic sections of treated animals.
benzyl propylene glycoside showed a potential ability to attenuate NAFP development, inhibit pancreatic inflammation and fibrosis and reduce the pathological and metabolic disturbances monitored in the applied NAFP animal model. The detected effect was correlated with modulation of the expression of pancreatic (DDX58, NFκB1, and CHUK mRNAs and miR-1976) panel.
Although non-alcoholic fatty pancreas (NAFP) was reported early in the 1930s, our knowledge about this disease is still in its infancy and perceived as a relatively new condition . NAFP is described as pancreatic fat cumulation without significant alcohol intake . It was considered a benign incidental finding, and therefore its clinical consequences were ignored. The prevalence of NAFP ranges from 16 to 35% and is increasingly associated with obesity, insulin resistance (IR), deterioration of beta-cell function and metabolic syndrome which might lead to the development of diabetes and pancreatitis . Therefore, its early detection may help to diagnose prediabetic patients to reduce the rising morbidity and mortality due to diabetes mellitus.
Although the shared association between NAFP and non-alcoholic fatty liver disease (NAFLD), the implicating mechanisms still unclear and has led researchers to hypothesize comparable etiologies of NAFP and NAFLD . The metabolic stress in NAFP, including insulin resistance and obesity, can stimulate severe acinar cell injury resulting in progressive acinar cell death, and acute pancreatitis, which acts as a trigger for various signaling mechanisms including the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway . The cGAS-STING pathway was found to be activated in acute pancreatitis and can induce cell injury by activating inflammation and by disturbing glucose and lipid metabolism [5, 6]. STING activation also affects several signaling cascades resulting in the induction of the nuclear factor kappa beta (NF-kB) to produce proinflammatory cytokines and activate fibrogenesis . Therefore, more exploration into this signaling mechanism might help in identifying novel therapies for NAFP disease.
Dysregulation of microRNA (miRNA) may impact the function and status of various tissues, like the pancreas [8,9,10], and liver [11,12,13], contributing to metabolic disorders associated with obesity and insulin resistance-linked diseases including NAFP. miRNAs play a very important role as key regulators of inflammation, insulin signaling, and glucose and lipid metabolism. However, information about the mechanisms of their implication in NAFP progression remains nearly limited, due to the ability of miRNAs to simultaneously affect several gene/pathway networks . This integrated gene (mRNAs)—miRNAs regulatory interaction may provide new early non-invasive diagnostic biomarkers and identification of therapeutic strategies for NAFP. Obviously, bioinformatic analysis facilitates the identification of new candidate RNA species and their interactions as biomarkers for disease screening, diagnosis, and therapy .
Emerging studies have reported that benzyl propylene glycoside (Rosavin), a main constituent of the Rhodiola Rosea plant, possesses several pharmacological effects such as anti-inflammatory, anti-adipogenic and hepato-protective effects on metabolic syndrome and related disorders [16,17,18]. The underlying mechanisms behind these effects may involve inhibition of NF-kB, reducing cell death, inhibition of adipogenesis, and modulation of miRNA expression , and this suggests that the miRNA may be a target for benzyl propylene glycoside treatment. However, its effect on NAFP is not clearly illustrated.
Based on the all previously discussed data, we aimed to investigate the potential therapeutic efficacy of benzyl propylene glycoside on NAFP management via targeting the pancreatic cGAS-STING pathway-related genes (DDX58, NFκB1 & CHUK) and their upstream regulator miRNA (miR-1976) that were retrieved from bioinformatics analysis in NAFP animal model.
Benzyl propylene glycoside—miR-1976 in-silico interaction prediction
The miRNA1976 secondary structure modelling showed MFE for thermodynamic ensemble of -31.85 kcal/mol which was further used to obtain the 3D model (Fig. 1A). Docking scores were recorded as in Table 1. Benzyl propylene glycoside (Rosavin) predicted interactions were sketched as an interaction per nucleotide and type of bond for the top 10 poses (Additional file 1: Fig. S1). HDock calculated confidence score was over 0.5 for one pose with a docking score of − 151.9 and a calculated RMSD of 40.24 (Fig. 1B).
Effect of 8 weeks of HFHS feeding on body weight and blood biochemical parameters
Feeding experimental rats with an HFHS diet for 8 weeks (Table 2) has resulted in a significant increase (p < 0.001) in body weight, HbA1C%, the levels of lipid profile markers, fasting serum glucose, and insulin as well as insulin resistance represented by HOMA-IR, compared to the Sham animals. In addition, a significant elevation in the level of serum amylase and lipase was also recorded in HFHS-fed animals. The results indicated that the HFHS-challenged animals exhibited signs of dyslipidemia and pancreatic damage.
Effect of benzyl propylene glycoside treatment on body weight and blood biochemical parameters
As shown in Table 3, body weight was significantly higher at the end of 12 weeks in untreated HFHS-fed rats than in Sham animals. After 4 weeks from benzyl propylene glycoside treatment, the body weight was significantly attenuated in HFHS rats (R-20 & R-30) compared to the NAFP group. Feeding animals an HFHS diet for the entire 12 weeks (NAFP group) caused a significant upsurge in the levels of serum glucose, insulin, and HbA1C% compared to the normal chow-fed rats (Sham group). Therefore, untreated NAFP rats presented a higher HOMA-IR (p < 0.001) than those of the Sham group. Serum levels of TG, TC, and LDL-C, but not HDL-C, also significantly (p < 0.001) increased in NAFP animals compared to the animals of the Sham group. On the other hand, animals of the three treated groups displayed a significant correction in all previous variables compared to the untreated NAFP group in a dose-dependent manner.
Regarding serum amylase and lipase, a highly significant elevation in the levels of these enzymes was recorded in animals in the untreated NAFP group compared to the Sham group. Daily injection with benzyl propylene glycoside for four weeks, caused a significant reduction in the levels of serum amylase and lipase, in comparison with the NAFP group. This ameliorative effect was more prominent in both the R-20 and R-30 groups.
As compared to the HFHS-8 week group, the benzyl propylene glycoside-treated groups (R-20 and R-30) showed a significant reduction in all previously mentioned biochemical variables that indicated Benzyl Propylene Glycoside has the potential to restore the initial pathological changes induced by HFHS feeding. Moreover, the results demonstrated that the four additional weeks of the HFHS diet resulted in more severe damage in NAFP animals.
The light microscopic examination of H&E sections of the pancreas of the Sham group revealed normal acinar arrangement with basal basophilia and apical acidophlia and the acinar cells have basal open phase nuclei (Fig. 2). The NAFP group showed loss of the normal lobulation of the pancreas. Large areas of pancreatic parenchyma were occupied by fat cells with noticeable areas of fat necrosis. Some intact acini were seen in between fat tissue. Meanwhile, the remaining acini appeared distorted, vacuolated. Rounded structures of variable size were detected which might be regenerating acini. Areas of intense mononuclear cellular infiltration and oedema were also noticed. The interlobular and intralobular connective tissues were relatively thickened. Some blood vessels showed fibrin clots with margination and pavementation of inflammatory cells. While pancreas of R-10 group showed focal structural changes in some lobules as well as interlobular and intralobular connective tissue. The affected acinar cells showed variable structural changes. Some acinar cells were lightly stained with loss of basal basophilia and apical acidophilia. Others showed vacuolated cytoplasm with pyknotic nuclei within oedematous areas. In group R-20, the pancreatic lobules and the pancreatic acini were closely packed, however, noticeable areas of pancreatic affection were still present. There were focal areas of disorganized acini. Localized areas of mononuclear cellular infiltration were still noticed. In R-30 group, pancreas showed normal structure with tightly packed pancreatic acini and thin interlobular septa. Most acini were formed of normal acinar cells with basal basophilia and apical secretory granules and vesicular nuclei. Few acini showed hyalinized cytoplasm.
In Masson’s trichrome stained sections (Fig. 3) showed progressive increase of collagen fibers deposition in all groups to be maximum in NAFP. In the treated group, collagen fibers were still noticed surrounding the blood vessels and thickened interlobular septa in both groups R-10 and R-20, but they were more pronounced in R-10. Despite of that, collagen fibers were apparently less than those in R-30. These results were confirmed by the statistical study. Morphometric and statistical study for area percentage of collagen fibers (Fig. 3F) revealed significant increase in NAFP in relation to other groups. In the treated group R-10 there was a significant increase as compared to the control. However, R-20 showed significant decrease as compared to NAFP group.
The effect of benzyl propylene glycoside on the expression of the pancreatic selected RNA species
Results showed a significant elevation in the expression of pancreatic DDX58, NFκB1 and CHUK mRNAs with a significant reduction in the expression level of miR-1976 in the untreated NAFP group compared to the Sham group (p < 0.001), Fig. 4. Meanwhile, the administration of benzyl propylene glycoside at its two higher dosages 20 & 30 (R-20 and R-30) significantly reduced the significant upregulation in the expression of pancreatic mRNA species manifested in untreated NAFP group animals. Moreover, the data were coupled with a significant increase in the expression of miR-1976 in the treated groups (R-20 and R-30) compared to the untreated NAFP group.
The effect of benzyl propylene glycoside on the pancreatic NFκB1 and Caspase-3
As shown in Fig. 5A–E, caspase-3-stained sections revealed minimal reaction in acinar cells in R-30 group. However moderate reaction was noticed in both the cytoplasm of acinar cells and the rounded structures in R-20. Positive reaction was distinguished in most of acinar cells in NAFP group. Morphometric and statistical study for area percentage of Caspase-3 positive cells reaction revealed significant elevation in NAFP and R-10 in comparison to other groups.
Nuclear factor kappa-stained sections (Fig. 6A–E) revealed minimal reaction in acinar cells in R-30 group. However moderate reaction was noticed in both the acinar cells and the rounded structures in R-20. Maximum Positive reaction was distinguished in most cells in NAFP group. Morphometric and statistical study for area percentage of nuclear factor kappa stained \ cells revealed significant elevation in R-10 group in comparison to other treated groups.
Nonalcoholic fatty pancreatitis (NAFP) considers one of the manifestations of metabolic syndrome that needs further studies to determine molecular determinants of this disorder and find effective medications . Emerging data showed that insulin resistance and dysregulation of the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway are the major driving forces for acute pancreatitis and fibrogenesis in NAFP progression . Thus, herein, we constructed a mRNAs (DDX58, NFκB1& CHUK)—(miR-1976) panel linked to metabolic syndrome and pancreatic cell dysfunction as well as be enrolled in the cGAS-STING pathway via in silico data analysis. Then we evaluated the potential ameliorative effects of benzyl propylene glycoside (Rosavin) treatment, the main constituent of the Rhodiola Rosea plant, on NAFP management and its effects on the constructed RNA panel in the NAFP animal model.
One of the primary mechanisms that explain the incidence of the fatty pancreas is the infiltration of adipocytes into the pancreatic tissue. Obesity and increased body weight are the major contributing factors to this condition. Adipose tissue is an endocrine organ as it emits signals to different organs. During weight gain, the storage of fat in adipose tissue is overridden, resulting in the excess lipid is deposited in visceral and peripheral non-adipose organs including the pancreas . Therefore, fatty infiltration of the pancreas is detected as ectopic adipocytes infiltrating the pancreatic tissue where fats deposit in adipocytes in the pancreatic tissue inducing pathological disorders such as insulin resistance and pancreatic cell injury and ultimately resulting in pancreatitis .
Accordingly, all previously discussed data can illustrate the results we obtained. We have used a high-fat and high-sucrose (HFHS) feeding as a representative experimental animal model of NAFP disease. Accumulating studies have investigated the impacts of HFHS diet on experimental animals, and it has been concluded that consumption of this diet induces obesity and insulin resistance [22,23,24]. In the current study, this nutritional model nearly covered the spectrum of the pathological and metabolic disturbances associated with NAFP. The HFHS diet feeding resulted in increased body weight, hyperglycemia, hyperinsulinemia, insulin resistance, and dyslipidemia in the untreated NAFP group. The animals also showed degrees of pancreatitis manifested by large areas of fat cells with noticeable areas of fat necrosis as well as areas of intense mononuclear cellular infiltration and oedema were also noticed in the pancreatic sections, and elevated serum levels of amylase and lipase.
Recently, it was reported that the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway can be activated by lipotoxicity-induced pancreatic cell injury . cGAS-STING signaling pathway is a crucial regulator of immune responses  and plays an important role in glucolipid metabolic disorders and was found to be activated in the animals fed a high-fat diet, and its gene silencing reversed metabolic dysfunction, insulin resistance, and inflammation . Activated STING can stimulate the phosphorylation of interferons (IRFs). Phosphorylated IRF regulates the expression of target genes, including DEAD Box Protein 58 (DDX58) & Nuclear Factor Kappa B Subunit 1 (NFκB1), to activate diverse downstream signaling pathways and promote the expression of inflammatory and fibrotic genes . Therefore, STING can promote cellular inflammation in several pathological conditions like insulin resistance.
The current results are consistent with the previously discussed data where a significant increase in the expression level of the pancreatic cGAS-STING pathway-related genes (DDX58 and NFκB1) in the untreated HFHS-fed animals (NAFP group) compared to the Sham group (p < 0.001). The results were propped by the results of histological and immunohistochemistry assay which showed a significant prevalent maximum positive immunostaining for NFκB1 marker, coupled with large areas of inflammatory cells infiltration in the pancreatic tissue, in comparison with the Sham group. Parallel studies have also confirmed that NFκB1 and DDX58 are upregulated in pancreatic cell injury [28,29,30]. Moreover, DDX58 was reported to be one of the genes that were differentially expressed in obese patients with type two diabetes mellitus .
The NF-KB activation requires the activity of the upstream serine/threonine protein kinase alpha (IKKα) which is encoded by the CHUK (conserved helix-loop-helix ubiquitous kinase) gene. The phosphorylation of nuclear factor-kB inhibitor (IKB) by IKKα results in its degradation and activation of NFκB1 . The results of the present study revealed a significant upregulation in the expression of the pancreatic CHUK in the NAFP group compared to the Sham group.
Increased expression of caspase-3 has been observed in the absence of apoptosis. Caspase-3 may be implicated in processes other than apoptosis where it can participate in inflammatory responses by cleaving and activating cytokines . Moreover, it was also reported that increased caspase-3 in the high-fat diet-fed animal was associated with a significant elevation in hepatic expression of inflammatory cytokines indicating that increased apoptosis could be an insulting mechanism in hepatic inflammation . Consistently, our results displayed a marked increase in caspase-3 protein expression in the pancreas of NAFP model rats compared to the Sham group.
Herein, all observed disturbances in the untreated NAFP animals were significantly adjusted by treatment of the experimental animals with benzyl propylene glycoside daily for four weeks. The recent studies on benzyl propylene glycoside showed that it exhibits anti-oxidative , anti-cancer , and anti-inflammatory effects . The toxicity of benzyl propylene glycoside has been previously assessed and shown to have a hepatoprotective effect and can alleviate kidney damage [17, 38, 39]. An emerging study that evaluated the toxicity of Rhodiola components showed that LD50 > 5000 mg/kg b.w., considers safe for consumption. This study also showed that long-term administration of Rhodiola doses (100, 250, and 500 mg/kg b.w.) for 28 days didn’t cause any toxic effects in experimental animals. Moreover, all the parameters related to the liver, and kidney were not affected .
It was also reported that benzyl propylene glycoside can attenuate cell injury and fibrosis through inhibition of NF-kB and decreasing the production of pro-inflammatory and fibrotic cytokines . Benzyl propylene glycoside can also improve cellular immunity by inhibiting tissue apoptosis .
The results of the current study were in accordance with the published data. The results revealed that the daily treatment with benzyl propylene glycoside had beneficial actions on the progression of NAFP. It significantly improved the lipid panel, decreased the body weight, lowered the serum insulin and glucose levels, ameliorated the insulin resistance status, and decreased the serum level of lipase and amylase. Surprisingly, as compared to the HFHS-8 week group, the benzyl propylene glycoside-treated groups (R-20 and R-30) revealed significant decreases in all detected biochemical variables that showed benzyl propylene glycoside has the potential to prevent the progression of exocrine pancreatic damage and can recover the initial pathological changes induced by HFHS feeding. Moreover, the applied treatment decreased the expression level of cGAS-STING pathway-related genes, DDX58, NFκB1 & CHUK, coupled with decreased the protein expression of pancreatic inflammatory NFκB1 and caspase-3 as compared with the untreated NAFP animals. Normal structure with tightly packed pancreatic acini and thin interlobular septa with a significant decrease in area percentage of collagen fibers were also detected in the treated groups (R-30) when compared to the NAFP model group. The revealed results indicate that benzyl propylene glycoside could improve pancreatic tissue injury via modulating and inhibiting the cGAS-STING pathway.
The currently identified biomarkers in the early diagnosis of NAFP are insufficient and poorly known. Thus, novel non-invasive biomarkers and precise therapeutic targets are required urgently. MicroRNAs (miRNAs) are a class of small non-coding RNAs that modulate the expression of protein-coding genes . They can be detected in body fluids, like blood and urine, and changes in their levels have been associated with several diseases therefore they can be utilized as diagnostic biomarkers . Accumulating evidence shows that modulation of miRNA expression could be one of the regulatory mechanisms behind the ameliorative activities of benzyl propylene glycoside . miRNAs also play crucial roles in the function and survival of pancreatic cells and have been found to regulate the adaptive responses of pancreatic cells in conditions like obesity and pancreatitis [42, 44].
In the present study, the miRWalk database was utilized to retrieve the upstream regulators, miR-1976, for the selected three mRNAs (CHUK, NFκB1, and DDX58). Regarding the in-silico study of benzyl propylene glycoside-miR-1976 interaction, the confidence scores showed slight significance for successful binding for the first pose with a confidence score > 0.5. However, the high deviation in the RMSD suggests the decreased probability of binding which needs further experimental proof of direct binding . Nevertheless, the alteration of miR-1976 maybe due to indirect effect of benzyl propylene glycoside on miR-1976. Previous emerging studies have demonstrated miR-1976 role as a prognostic indicator and tumor suppressor in non-small lung cancer progression . Moreover, it was reported that the miR-1976 knockdown significantly inhibited cell apoptosis and increased cell proliferation . miR-1976 was also found to be one of the specific downregulated exosomal-miRNA signatures related to pancreatic lesions . Interestingly, the functional enrichment analysis of miR-1976 revealed that it is highly linked to inflammatory cGAS-STING-related and fibrogenic pathways including NF-KB signaling, TGF signaling, and TNF signaling pathways.
Herein, DDX58, NFκB1, and CHUK were screened as target genes of miR-1976 using the mirwalk3 database. miR-1976 can regulate the expression level of these genes via binding to their 3ʹUTR resulting in post-transcriptional inhibition or their degradation . Accordingly, the results showed that there was a significant decrease in the expression level of pancreatic miR-1976 in the untreated NAFP group, in comparison with the Sham control. While benzyl propylene glycoside administration significantly increased its expression, compared to the NAFP group.
Taken all together (Fig. 7), we hypothesized that HFHS-induced lipotoxicity (untreated NAFP) downregulated the expression of miR-1976 which could not exert its inhibitory action on its target genes thereby upregulating the expression of pancreatic DDX58, NFκB1, and CHUK mRNAs. Activating the cGAS-STING signaling pathway stimulated diverse downstream signaling pathways, promoted the expression of inflammatory responses (NFκB1 and Caspase-3), increased the area percentage of collagen fibers (fibrosis), and increased the serum level of lipase and amylase. Consequently, increasing pancreatic cell injury and pancreatitis progression. On treatment, the benzyl propylene glycoside increased the expression of miR-1976 and inhibited the expression of its target genes (DDX58, NFκB1, and CHUK). Inhibiting the cGAS-STING signaling pathway reversed metabolic dysfunction, ameliorated insulin resistance, decreased body weight and obesity, and reduced inflammation, and fibrosis observed in the untreated NAFP group.
The present study may help in better understanding the etiology and pathophysiology of the non-alcoholic fatty pancreas disease (NAFP) and also provides useful information regarding potential molecular targets for NAFP treatment. However, benzyl propylene glycoside may not yet be a suitable fundamental mode of therapy until further preclinical trials are performed.
Benzyl propylene glycoside has demonstrated a potential ability to attenuate NAFP development, inhibit pancreatic cell inflammation and fibrosis and reduce the pathological and metabolic disturbances monitored in the applied NAFP animal model. The detected effect was correlated with upregulation of the expression of pancreatic DDX58, NFκB1, and CHUK mRNAs and downregulation of the expression of pancreatic miR-1976.
Material and methods
Drugs and materials
Sodium pentobarbital was obtained from Sigma Aldrich (St. Louis, Missouri, USA). Rosavin (benzyl propylene glycoside) was supplied from Aktin Chemicals, Inc (Cat. #. APC-380, China).
Animals and treatment
The handling and experimentation protocols were reviewed and approved by the Research Ethics Committee (Number; MoHP0018122017, 1017), Faculty of Medicine, Benha University. The experimental study was performed according to the Declaration of Helsinki guidelines. Male Wistar rats (150–170 g), were housed in cages under standard controlled conditions (12 h light/dark cycles and 21 ± 2 °C) and randomly grouped into normal chow-fed rats (Sham group, n = 8) and high-fat high-sucrose-fed rats (HFHS), Table 4, as a nutritional model for NAFP induction . After 8 weeks of dietary intake, blood samples were drawn to evaluate the effect of the HFHS diet manipulation on the experimental animals. The HFHS-fed animals were then subdivided into 4 groups (n = 8 for each group): untreated HFHS group (NAFP model group) and three benzyl propylene glycoside (rosavin)-treated groups, R-10 group, R-20 group, and R-30 group (Fig. 8). In these treated groups, the rats injected intraperitoneally with10 mg, 20 mg, and 30 mg rosavin/kg body weight, respectively for 4 weeks parallel with HFHS diet . The normal chow-fed rats were given vehicle 0.9% saline intraperitoneally.
Euthanasia and blood and pancreas tissue collection
All the experimental rats were monitored daily for body weight. At the end of the experimental period (12 weeks), the experimental rats were anesthetized with a single dose of sodium pentobarbital (45 mg/kg, intraperitoneally)  and blood samples were rapidly obtained from the retro-orbital vein. Serum was then obtained by centrifugation (1200 g for 10 min) and stored at − 20 °C for the biochemical analyses. The pancreas was carefully removed, weighed, and then rapidly fixed in freshly prepared 10% neutral buffered formaldehyde for analysis by light microscopy.
Serum biochemical analysis
Lipid profile markers and glycated hemoglobin (HbA1C)
Total cholesterol (TC), HDL cholesterol (HDL-C), LDL cholesterol (LDL-C), triglycerides (TG), fasting serum glucose and glycated hemoglobin (HbA1C) were quantitatively determined by the multifunctional biochemistry analyzer (AU680, Beckman Coulter Inc., Brea, CA, USA).
Nonalcoholic fatty pancreas (NAFP)-model markers
Serum insulin was measured using a rat sandwich ELISA kit purchased from Invitrogen (Cat. NO. ERINSX10, Waltham, Massachusetts, USA) according to the manufacturer’s instructions. Serum lipase and amylase were measured using commercial kits obtained from Erba Diagnostics (Miami, Florida, USA) according to the protocol supplied with the respective kits. Homeostasis model assessment-insulin resistance (HOMA-IR) was calculated using the following formula: HOMA-IR = [fasting serum insulin (µU/ml) x fasting serum glucose (mg%)]/405 .
Pancreatic histological and immunohistochemistry assays
The buffered formalin-fixed pancreatic samples were dehydrated using an ascending concentration of alcohol, cleared using methyl benzoate, and mounted in paraffin blocks. Sections were cut at a thickness of 5 μm and stained using hematoxylin and eosin (H&E) and Masson's trichrome stain for the detection of collagen fibers. Other paraffin sections were cut and placed on positively charged slides and were exposed to immune reaction for caspase 3 monoclonal antibody (Cat. No. CPP32 4-1-18, Invitrogen, Waltham, MA, USA) and NFκB1 antibody (Cat. No. BS-3300R, Bioss Antibodies, Woburn, MA, USA). The positive reactions for the caspase 3 and NFκB1 immune-histochemical technique appeared as brown nuclear and cytoplasmic reactions. Negative controls were performed according to the same protocol, but without the usage of the primary antibody. Positive control was performed using a section of tonsils. Finally, the slides were counterstained using Mayer’s hematoxylin. Positive controls were carried out according to the same protocol .
The morphometric study was done using an image analyzer Leica Q win V.3 program installed on a computer which connected to a Leica DM2500 microscope (Wetzlar, Germany). Pancreatic slides from all groups were evaluated by morphometric study. Evaluations were obtained from five different slides taken from each rat. Five non-overlapping fields were selected haphazardly and examined for each slide. The pancreatic slides were used to measure:
I-The mean area percentage (%) of collagen fibers in Masson's trichrome stained sections at objective lens X 20.
II- The mean area percentage (%) of positive reaction of caspase-3 and NFκB1 sections (X20).
Bioinformatics set up
Retrieval of the mRNAs-miRNAs panel
The RNAs species that are related to NAFP development and implicated in obesity and insulin resistance were searched for. Firstly, the differentially expressed genes (mRNAs) associated with pancreatic injury were screened through the Gene Expression Omnibus (GEO) (www.ncbi.nlm.nih.gov/geo/, accessed on 22 Oct 2021) . The screened mRNAs were further filtered according to their significant differential expression (Additional file 1: Fig. S2), their pancreatic tissue-specific expression (Additional file 1: Fig. S3), and their links to the cGAS-STING signaling pathway. From the filtered mRNAs, DEAD Box Protein 58 (DDX58), Nuclear Factor Kappa B Subunit 1 (NFκB1), and Conserved Helix-Loop-Helix Ubiquitous Kinase (CHUK) were selected as they were validated by other microarray databases (Additional file 1: Fig. S4) and by reviews [53,54,55,56,57] to be related to metabolic syndrome and pancreatic cell dysfunction diseases. The selected genes were also mapped and visualized through the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (https://www.genome.jp/kegg/, accessed on 22 Oct 2021) to be enrolled in the cGAS-STING pathway (Additional file 1: Fig. S5). The pathway enrichment analysis using Enrichr (http://amp.pharm.mssm.edu/Enrichr, accessed on 22 Oct 2021)  was primarily enriched in cGAS-STING and NF-kappa B signaling pathways. The top ten terms for pathway enrichment are shown in Fig. 9A. Based on the STRING tool (http://stringdb.org, accessed on 22 Oct 2021) , the protein–protein interaction (PPI) between the three selected genes showed a high confidence level with a combined score > 0.7 (Fig. 9B).
Secondly, miRWalk 3.0 (http://mirwalk.umm.uni-heidelberg.de/, accessed on 22 Oct 2021) was utilized for the retrieval of miRNAs interacting with the three selected mRNAs. miR-1976 (Fig. 9C) was found to target the 3 selected mRNAs with a score ˃ 0.9 (Additional file 1: Fig. S6). DIANA tools mirPath (http://www.microrna.gr/miRPathv3, accessed on 22 Oct 2021) was then used to track pathways of miR-1976. Interestingly, miR-1976 was detected to be related to cGAS-STING-related pathways (Additional file 1: Fig. S7).
All in all, the mRNAs (DDX58, NFκB1& CHUK)—(miR-1976) panel was constructed.
Molecular docking analysis: benzyl propylene glycoside (Rosavin)—miR-1976 in-silico interaction prediction
The molecular docking between the upstream regulator miR-1976 and rosavin was performed. Rosavin ligand was obtained from PubChem with ID: 9,823,887. The miRNA1976 sequence was extracted from the miRbase database with accession number: MI0009986. The secondary structure was computed using RNAFold under ViennaRNA package (Version 2.4.18) . The minimum free energy (MFE) of the secondary structure was computed at 37 °C. The secondary structure was subjected to 3D modelling using RNAComposer web server [61, 62]. The 3D model was used for docking using HDock software which models the protein using two algorithms: template-based and ab initio modelling . The docked forms are ranked upon their docking scores, Root mean standard deviation and confidence score according to the HDock manual. The predicted Rosavin-miRNA1976 interaction for the top 10 poses were calculated using BIOVIA Drug Discovery Studio Visualizer 2021 (version 18.104.22.16898).
Total RNA extraction and quantitative polymerase chain reaction (qPCR)
Total RNA, involving mRNAs and miRNAs, extraction from the 60 mg of frozen pancreas tissue samples was performed using a miRNEasy extraction kit (Qiagen, Hilden, Germany, Cat. No. 217004) according to the protocol supplied with the kit. NanoDrop (Thermo scientific, USA) was utilized to assess the concentration and purity of total RNA and the purity of the isolated RNAs was adjusted to be 1.8–2 (A260/A280). The RNA extracted from the pancreas tissues was then reverse transcribed into complementary DNA using miScript II RT (Cat. No. 218161, Qiagen, Germany).
Relative expression of the selected RNAs species in the pancreatic tissue samples was assessed using a Quantitect SYBR Green Master Mix Kit (Qiagen, Germany, Cat. No. 204143) for DDX58, NFκB1, and CHUK mRNAs and miScript SYBR Green PCR Kit (Qiagen, Germany, Cat no. 218073) for miR-1976 miRNA. Real-time (RT)-qPCR was conducted on 7500 Fast System (Applied Biosystems, Foster City, USA). The GAPDH and SNORD72 were used as housekeeping genes. The primers list used herein was obtained from Qiagen, Germany (Additional file 1: Table S1). The relative quantification of RNA expression was calculated using RQ = 2 –ΔΔCt formula .
GraphPad Prism software, version 8.0 (Inc., CA, USA) was utilized to perform the Statistical analyses. The distribution normality of the data was analyzed using the Kolmogorov–Smirnov test. Data are represented as the mean ± standard deviation (SD). Differences among groups were analyzed by one-way analysis of variance (ANOVA) for statistical significance, followed by Tukeyʼs test.
All data generated during this study are included in this article.
Cyclic GMP-AMP synthase
- CHUK :
Conserved helix-loop-helix ubiquitous kinase
- DDX58 :
DEAD box protein 58
High-fat high-sucrose-fed rats
Non-alcoholic fatty liver disease
Non-alcoholic fatty pancreas
Nuclear factor kappa beta
- NFκB1 :
Nuclear factor kappa b subunit 1
Stimulator of interferon genes
Khoury T, Sbeit W. Fatty pancreas and pancreatic cancer: an overlooked association? J Clin Med. 2022. https://doi.org/10.3390/jcm11030763.
Singh P, Singh K, Dhanju AS, et al. The association of non-alcoholic fatty pancreas disease and glycemic status. Int J Adv Med. 2022. https://doi.org/10.18203/2349-3933.ijam20221350.
Shah N, Rocha JP, Bhutiani N, Endashaw O. Nonalcoholic fatty pancreas disease. Nutr Clin Pract. 2019;112:112.
Sundar V, Dutta A, Ramasamy S, et al. Sting pathway—a futuristic therapeutic target for acute pancreatitis? Gene. 2021. https://doi.org/10.1016/j.gene.2021.145469.
Zhao Q, Wei Y, Pandol SJ, et al. STING signaling promotes inflammation in experimental acute pancreatitis. Gastroenterology. 2018. https://doi.org/10.1053/j.gastro.2018.01.065.
Couillin I, Riteau N. STING Signaling and Sterile Inflammation. Front Immunol. 2021;112:324.
Albadawy R, Agwa SHA, Khairy E, et al. Clinical significance of hspd1/mmp14/itgb1/mir-6881-5p/lnc-sparcl1-1:2 rna panel in nafld/nash diagnosis: Egyptian pilot study. Biomedicines. 2021. https://doi.org/10.3390/biomedicines9091248.
Hu Y, Yu Y. Dysregulation of miR-192-5p in acute pancreatitis patients with nonalcoholic fatty liver and its functional role in acute pancreatitis progression. 2020. Biosci Rep. https://doi.org/10.1042/BSR20194345.
Suzuki R, Asama H, Waragai Y, et al. Fibrosis-related miRNAs as serum biomarkers for pancreatic ductal adenocarcinoma. Oncotarget. 2018. https://doi.org/10.18632/oncotarget.23377.
Eliasson L, Esguerra JLS. MicroRNA networks in pancreatic islet cells: Normal function and type 2 diabetes. Diabetes. 2020. https://doi.org/10.2337/dbi19-0016.
Gadallah SH, Eissa S, Ghanem HM, et al. Probiotic-prebiotic-synbiotic modulation of (YAP1, LATS1 and NF2 mRNAs/miR-1205/lncRNA SRD5A3-AS1) panel in NASH animal model. Biomed Pharmacother. 2021. https://doi.org/10.1016/j.biopha.2021.111781.
Schueller F, Roy S, Vucur M, et al. The role of miRNAs in the pathophysiology of liver diseases and toxicity. Int J Mol Sci. 2018;11:123.
Afshari A, Yaghobi R, Karimi MH, Mowla J. Alterations in MicroRNA gene expression profile in liver transplant patients with hepatocellular carcinoma. BMC Gastroenterol. 2021. https://doi.org/10.1186/s12876-020-01596-2.
Iacomino G, Siani A. Role of microRNAs in obesity and obesity-related diseases. Genes Nutr. 2017;11:231.
Matboli M, Gadallah SH, Rashed WM, et al. Mrna-mirna-lncrna regulatory network in nonalcoholic fatty liver disease. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22136770.
Chen Y, Tang M, Yuan S, et al. Rhodiola rosea: a therapeutic candidate on cardiovascular diseases. Oxid Med Cell Longev. 2022;112:132.
Albadawy R, Hasanin AH, Agwa SHA, et al. Rosavin ameliorates hepatic inflammation and fibrosis in the NASH rat model via targeting hepatic cell death. Int J Mol Sci. 2022. https://doi.org/10.3390/ijms231710148.
Pomari E, Stefanon B, Colitti M. Effects of two different Rhodiola rosea extracts on primary human visceral adipocytes. Molecules. 2015. https://doi.org/10.3390/molecules20058409.
Ling PW, Zhang M, Bai R, et al. Anti-inflammatory effects of Rhodiola rosea L: a review Biomed. Pharmacother. 2020;112:332.
de Silva EL, de Fernandes MS, de Lima EA, et al. Fatty pancreas: disease or finding? Clinics. 2021;112:114.
Guglielmi V, Sbraccia P. Type 2 diabetes: does pancreatic fat really matter? diabetes. Metab Res Rev. 2018;111:234.
Fernandes-Santos C, Evangelista Carneiro R, De Souza ML, et al. Rosiglitazone aggravates nonalcoholic fatty pancreatic disease in C57BL/6 mice fed high-fat and high-sucrose diet. Pancreas. 2009. https://doi.org/10.1097/MPA.0b013e3181987d9d.
Rasool S, Geetha T, Broderick TL, Babu JR. High fat with high sucrose diet leads to obesity and induces myodegeneration. Front Physiol. 2018;112:336.
Pennington KA, van der Walt N, Pollock KE, et al. Effects of acute exposure to a high-fat, high-sucrose diet on gestational glucose tolerance and subsequent maternal health in mice. Biol Reprod. 2017. https://doi.org/10.1095/biolreprod.116.144543.
Hu HQ, Qiao JT, Liu FQ, et al. The STING-IRF3 pathway is involved in lipotoxic injury of pancreatic β cells in type 2 diabetes. Mol Cell Endocrinol. 2020. https://doi.org/10.1016/j.mce.2020.110890.
Rodríguez-García E, Olagüe C, Ríus-Rocabert S, et al. TMEM173 alternative spliced isoforms modulate viral replication through the STING pathway. ImmunoHorizons. 2018. https://doi.org/10.4049/immunohorizons.1800068.
Qiao JT, Cui C, Qing L, et al. Activation of the STING-IRF3 pathway promotes hepatocyte inflammation, apoptosis and induces metabolic disorders in nonalcoholic fatty liver disease. Metabolism. 2018. https://doi.org/10.1016/j.metabol.2017.09.010.
Fawzy H, Fikry E, Fawzy H, Mohammed A. Mito-TEMPO improved L-Arginine- induced acute pancreatitis in rats via TLR-4/ NF-кB/ NLRP3 inflammasome downregulation and antioxidant properties. Azhar Int J Pharm Med Sci. 2021. https://doi.org/10.21608/aijpms.2021.54059.1026.
Farid A, Moussa P, Youssef M, et al. Melatonin relieves diabetic complications and regenerates pancreatic beta cells by the reduction in NF-kB expression in streptozotocin induced diabetic rats. Saudi J Biol Sci. 2022. https://doi.org/10.1016/j.sjbs.2022.103313.
Singh T, Sarmiento L, Luan C, et al. MafA expression preserves immune homeostasis in human and mouse islets. Genes. 2018. https://doi.org/10.3390/genes9120644.
An T, Zhang J, Liu YF, et al. Combined analysis of whole-exon sequencing and lncRNA sequencing in type 2 diabetes mellitus patients with obesity. J Cell Mol Med. 2020. https://doi.org/10.1111/jcmm.14932.
Li X, Hu Y. Attribution of nf-κb activity to chuk/ikkα-involved carcinogenesis. Cancers. 2021. https://doi.org/10.3390/cancers13061411.
Loppnow H, Guzik K, Pryjma J. C HAPTER 8 the role of caspases in modulation of cytokines and other molecules in apoptosis and inflammation. Madame Curie Biosci Database. 2002;1:1–21.
Wang Y, Ausman LM, Russell RM, et al. Increased apoptosis in high-fat diet-induced nonalcoholic steatohepatitis in rats is associated with c-Jun NH2-terminal kinase activation and elevated proapoptotic bax. J Nutr. 2008. https://doi.org/10.1093/jn/138.10.1866.
Huang SC, Lee FT, Kuo TY, et al. Attenuation of long-term Rhodiola rosea supplementation on exhaustive swimming-evoked oxidative stress in the rat. Chin J Physiol. 2009. https://doi.org/10.4077/CJP.2009.AMH029.
Skopińska-Rózewska E, Hartwich M, Siwicki AK, et al. The influence of Rhodiola rosea extracts and rosavin on cutaneous angiogenesis induced in mice after grafting of syngeneic tumor cells. Cent J Immunol. 2008;11:123.
Marchev AS, Dimitrova P, Koycheva IK, Georgiev MI. Altered expression of TRAIL on mouse T cells via ERK phosphorylation by Rhodiola rosea L. and its marker compounds. Food Chem Toxicol. 2017. https://doi.org/10.1016/j.fct.2017.02.009.
Senthilkumar R, Chandran R, Parimelazhagan T. Hepatoprotective effect of Rhodiola imbricata rhizome against paracetamol-induced liver toxicity in rats. Saudi J Biol Sci. 2014. https://doi.org/10.1016/j.sjbs.2014.04.001.
Uyeturk U, Terzi EH, Kemahli E, et al. Alleviation of kidney damage induced by unilateral ureter obstruction in rats by rhodiola rosea. J Endourol. 2013;27:1272–6. https://doi.org/10.1089/end.2013.0319.
Gupta V, Tyagi A, Pandey P, Bhatnagar A. Biological and chemical sciences prophylactic efficacy of rhodiola rosea extract against hypobaric hypoxia- induced pulmonary injuries in rats and its toxicity evaluation. Res J Pharm. 2020;11:34–47.
Xin X, Yao D, Zhang K, et al. Protective effects of Rosavin on bleomycin-induced pulmonary fibrosis via suppressing fibrotic and inflammatory signaling pathways in mice. Biomed Pharmacother. 2019. https://doi.org/10.1016/j.biopha.2019.108870.
Xiang H, Tao X, Xia S, et al. Targeting microRNA function in acute pancreatitis. Front Physiol. 2017;12:326.
Wang J, Ma G, Han X, et al. The low expression of miR-1976 in plasma samples indicating its biological functions in the progression of breast cancer. Clin Transl Oncol. 2020. https://doi.org/10.1007/s12094-020-02361-3.
LaPierre MP, Stoffel M. MicroRNAs as stress regulators in pancreatic beta cells and diabetes. Mol Metab. 2017;11:123.
Yan Y, Tao H, He J, Huang SY. The HDOCK server for integrated protein–protein docking. Nat Protoc. 2020. https://doi.org/10.1038/s41596-020-0312-x.
Chen G, Hu J, Huang Z, et al. MicroRNA-1976 functions as a tumor suppressor and serves as a prognostic indicator in non-small cell lung cancer by directly targeting PLCE1. Biochem Biophys Res Commun. 2016. https://doi.org/10.1016/j.bbrc.2016.04.030.
Wang J, Li M, Han X, et al. MiR-1976 knockdown promotes epithelial–mesenchymal transition and cancer stem cell properties inducing triple-negative breast cancer metastasis. Cell Death Dis. 2020. https://doi.org/10.1038/s41419-020-2711-x.
Vicentini C, Calore F, Nigita G, et al. Exosomal miRNA signatures of pancreatic lesions. BMC Gastroenterol. 2020. https://doi.org/10.1186/s12876-020-01287-y.
Fraulob JC, Ogg-Diamantino R, Fernandes-Santos C, et al. A mouse model of metabolic syndrome: Insulin resistance, fatty liver and non-alcoholic fatty pancreas disease (NAFPD) in C57BL/6 mice fed a high fat diet. J Clin Biochem Nutr. 2010. https://doi.org/10.3164/jcbn.09-83.
Gadallah SH, Ghanem HM, Abdel-Ghaffar A, et al. 4-Phenylbutyric acid and rapamycin improved diabetic status in high fat diet/streptozotocin-induced type 2 diabetes through activation of autophagy. Arch Physiol Biochem. 2021;127:235–44. https://doi.org/10.1080/13813455.2019.1628069.
Suvarna K (2012) Bancroft’s Theory and practice of histological techniques E-Book 7th Edition
Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10. https://doi.org/10.1093/nar/30.1.207.
Chan LK, Tsesmelis M, Gerstenlauer M, et al. Functional IKK/NF-κB signaling in pancreatic stellate cells is essential to prevent autoimmune pancreatitis. Commun Biol. 2022. https://doi.org/10.1038/s42003-022-03371-3.
Meyerovich K, Ortis F, Cardozo AK. The non-canonical NF-κB pathway and its contribution to β-cell failure in diabetes. J Mol Endocrinol. 2018;11:112.
Baker RG, Hayden MS, Ghosh S. NF-κB, inflammation, and metabolic disease. Cell Metab. 2011. https://doi.org/10.1016/j.cmet.2010.12.008.
Yang G, Lee HE, Seok JK, et al. Rig-i deficiency promotes obesity-induced insulin resistance. Pharmaceuticals. 2021. https://doi.org/10.3390/ph14111178.
Hastie E, Cataldi M, Moerdyk-Schauwecker MJ, et al. Novel biomarkers of resistance of pancreatic cancer cells to oncolytic vesicular stomatitis virus. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.11202.
Chen EY, Tan CM, Kou Y, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013. https://doi.org/10.1186/1471-2105-14-128.
Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47:D607–13. https://doi.org/10.1093/nar/gky1131.
Lorenz R, Bernhart SH, Höner C, zu Siederdissen, et al. ViennaRNA package 2.0. Algorithms Mol Biol. 2011. https://doi.org/10.1186/1748-7188-6-26.
Antczak M, Popenda M, Zok T, et al. New functionality of RNAcomposer: An application to shape the axis of miR160 precursor structure. Acta Biochim Pol. 2016. https://doi.org/10.18388/abp.2016_1329.
Popenda M, Szachniuk M, Antczak M, et al. Automated 3D structure composition for large RNAs. Nucleic Acids Res. 2012. https://doi.org/10.1093/nar/gks339.
Yan Y, Zhang D, Zhou P, et al. HDOCK: A web server for protein-protein and protein-DNA/RNA docking based on a hybrid strategy. Nucleic Acids Res. 2017. https://doi.org/10.1093/nar/gkx407.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25:402–8. https://doi.org/10.1006/meth.2001.1262.
Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). This study was supported by Academy of Science and Technology (ASRT), Jesor call 2019, Project ID: jesor_5269.
Ethics approval and consent to participate
The handling and experimentation protocols were approved by the Research Ethics Committee (Number; MoHP0018122017, 1017), Faculty of Medicine, Benha University. The experimental study was performed according to the Declaration of Helsinki guidelines.
Consent for publication
All authors give their consent for publication.
All the authors have no competing interests to disclose.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The docking poses of Rosavin-miRNA1976 interaction. Figures S2. The significant differential expression of the selected candidate genes (DDX58, NFκB1, and CHUK) in pancreatic injury using the Expression Atlas database. Figures S3. Validation of the significant expression of the candidate genes/proteins (DDX58, NFκB1, and CHUK) in the pancreatic tissue. Figure S4. Validation of the implication of DDX58, NFκB1, and CHUK in metabolic syndrome and pancreatic cell dysfunction diseases. Figure S5. The visualization of the selected DDX58, NFκB1, and CHUK genes in the cGAS-STING pathway through KEGG pathway database. Figure S6. Validation of the interaction between the selected m-RNAs and the retrieved miR-1976 from mirwalk3. Figure S7. Validation of the relation of miR-1976 to cGAS-STING-related pathways through DIANA tools mirPath 3. Table S1. List of primer assays.
About this article
Cite this article
Albadawy, R., Hasanin, A.H., Agwa, S.H.A. et al. Prospective insight into the role of benzyl propylene glycoside as a modulator of the cGAS-STING signaling pathway in the management of nonalcoholic fatty pancreas animal model. Biol Res 56, 11 (2023). https://doi.org/10.1186/s40659-023-00423-8
- Nonalcoholic fatty pancreas
- cGAS-STING pathway
- Benzyl propylene glycoside