hsa_circ_0111707 Is Associated With Risk of Stress-Related Type 2 Diabetes via Sponging miR-144-3p

Introduction: Chronic stress plays an important role in the development of type 2 diabetes (T2D). Circular RNAs (circRNAs) play significant roles in regulating the pathogenesis of diseases by regulating gene expression. The aim of the present study was to identify the association between hsa_circ_0111707 and stress-related T2D.

Methods: The present study was performed based on a three-part design. The association between hsa_circ_0111707 in peripheral blood mononuclear cells (PBMCs) and T2D and stress-related variables were assessed in a cross-sectional study. The causal relationship of hsa_circ_0111707 on T2D was further investigated in a nested case-control study. miR-144-3p as the miRNA target of hsa_circ_0111707 was verified by dual-luciferase reporter assay and RNA immunoprecipitation (RIP) assay.

Results: The relative expression of hsa_circ_0111707 was significantly lower in the T2D and impaired fasting glucose (IFG) cases in comparison with controls. The hsa_circ_0111707 expression was significantly negatively correlated with miR-144-3p expression and plasma cortisol concentration and positively correlated with NR3C1 expression. In addition, hsa_circ_0111707 expression was negatively correlated with scores of “demands at work” and “insecurity at work” of Copenhagen Psychosocial Questionnaire (COPSOQ). Decreased hsa_circ_0111707 expression was associated with increased risk of T2D development. Functional analysis demonstrated that hsa_circ_0111707 functions as a sponge for miR-144-3p.

Conclusion: hsa_circ_0111707 is associated with risk of T2D development via sponging miR-144-3p. hsa_circ_0111707 in PBMCs can be considered a potential biomarker of stress-related T2D.

Introduction

Type 2 diabetes (T2D) has become a major public health problem worldwide. In the last decade, the incidence and prevalence of T2D in children and adolescents has increased dramatically (1). In 2015, 415 million people were estimated to have diabetes, more than 90% of whom had type 2 diabetes, with a projected increase to 642 million by 2040 (2). In addition to genetic causes and unhealthful lifestyles, environmental influences, such as psychosocial stress, has been implicated as an important contributor in the development of T2D (3, 4). Psychological stress can affect health through neuroendocrine responses and has been shown to cause impaired glucose homeostasis in earlier experimental and human researches (4, 5). Meta-analytic evidence indicated that job strain, which is defined as high job demands coupled with low control at work, is associated with an increased risk of T2D (6). Association between perceived stress and incident T2D was validated in several cohort studies (7–9). Our previous investigation also showed that chronic stress is associated with insulin resistance (IR), and plasma cortisol is significantly positively correlated with glucose and HOMA-IR (10). Long-term exposure to excess cortisol, such as Cushing syndrome and glucocorticoid treatment, often results in increased susceptibility to hyperglycemia and manifests diabetes (11, 12).

One of the major neuroendocrine systems responding to psychological stress is the hypothalamic-pituitary-adrenal (HPA) axis, with cortisol (glucocorticoid) secretion as its final hormone (13). The HPA axis is highly responsive to the environment, impacts both nervous and other physiological systems, and is critical to health in a wide variety of contexts (14). As a stress hormone, cortisol takes part in the regulation of metabolism, homeostasis, immune reactions, and behavioral responses (15). The effects of cortisol are mediated by intracellular glucocorticoid receptor (GR), coded by the NR3C1 gene (16). The expression levels of NR3C1 influence sensitivity to cortisol in physiological and pathophysiological conditions and modify response to cortisol in relevant target tissues (17). Investigation in an occupational population showed that chronic stress is associated with increased plasma cortisol and decreased NR3C1 expression in circulating monocytes (18). This phenomenon can be explained that reduced NR3C1 expression helps to alleviate excessive stress-related physiological processes (17).

In the past decade, microRNAs (miRNAs) have been recognized as important regulators of gene expression on the posttranscriptional level through sequence-specific binding to mRNA, which results in inhibiting the translation of protein from mRNA or inducing mRNA degradation (19). miRNA expression levels in rodents and human cells have been found to be altered in response to stress conditions, which in turn alters the expression of HPA axis-related target genes (20). Haramati et al. reported that overexpression of miR-34c appears to reduce cell responsiveness to secrete corticotrophin-releasing hormone (CRH) by inhibiting CRHR1 expression and induce an anxiolytic phenotype (21). miR-135a and miR-124 have also been implicated as potential suppressors of NR3C2 protein expression in rodents (22). In our previous study, let-7b, miR-142, miR-144, and miR-29a have been identified as neuroendocrine stress response-related biomarkers for T2D and IR (23). NR3C1 was further identified as a target gene of miR-144 (24). These evidences suggest that stress-related miRNAs may be potential novel targets for treatment of stress-related disorders.

Circular RNAs (circRNAs) were recently discovered as a special novel type of endogenous noncoding RNAs (ncRNAs) that are involved in transcriptional and posttranscriptional gene expression regulation (25). circRNAs differ from linear RNAs in that they are circular molecules with covalently closed loop structures and lack 5′ or 3′ polarity or a polyadenylated tail, which makes them much more stable than most linear RNA (26). As the competitive endogenous RNAs (ceRNAs), circRNAs can compete for miRNA-binding sites and make them function as miRNA sponges, thus suppressing miRNA activity and participating in regulation of miRNA target genes (25–27). Accumulated evidence indicates that circRNAs play crucial roles in the onset and progression of diseases and are expected to become diagnostic biomarkers and treatment targets in human diseases (27–29).

Based on differential expression profile of circRNAs in previous microarray analysis (29) and predictive target adsorption relationship between circRNA and miRNA (miRanda), a circulating circRNA, hsa_circ_0111707, was discovered as a potential regulator of NR3C1 by targeting miR-144-3p. In this study, association between hsa_circ_0111707 expression and T2D was investigated, which may unveil new target for the prevention and treatment of stress-related disorders including T2D.

Materials and Methods

Subjects

This study was performed based on a three-part design (Figure 1). Firstly, a cross-sectional study was conducted to compare the relative expression of hsa_circ_0111707 and stress-related variables among T2D, prediabetes cases and health controls. In addition, associations between relative expression of hsa_circ_0111707 and stress-related variables were analyzed. Secondly, a nested case-control study, prospectively collecting blood samples from T2D cases and matched controls before clinical diagnosis, was performed to further assess the predictive efficiency of hsa_circ_0111707 for T2D development. Thirdly, luciferase assay and RNA immunoprecipitation (RIP) assay were used to validate the function of hsa_circ_0111707 as a sponge for miR-144-3p.

FIGURE 1

Figure 1 The flowchart of study design. T2D, type 2 diabetes; IFG, impaired fasting glucose; COPSOQ, Copenhagen Psychosocial Questionnaire; HOMA-IR, homeostasis model assessment of insulin resistance; RR, relative risk; AUC, area under the curve; HEK, human embryonic kidney; HSkMM, human skeletal muscle myoblasts.

In the first part of the study, the participants were composed of 324 individuals, including 115 newly diagnosed T2D cases, 105 health controls, and 104 individuals with impaired fasting glucose (IFG) as prediabetes cases. Controls and IFG cases were frequency matched with T2D cases according to gender and age (± 3 years). All the participants were aged from 30 to 65 and were recruited from January 2016 to January 2017 in a functional community cohort composed of nearly 10,000 employees in Beijing Xicheng district (23). In the second part of the study, cases were newly diagnosed T2D patients in the same cohort during 2020. A total of 95 incident T2D cases were identified. Controls were selected from study participants who did not develop T2D during the follow-up period. For each case, one control was matched with the criteria of gender and age (± 12 months). The relative expression of hsa_circ_0111707 at baseline in 2016 between the two groups was compared. The average time between blood collection and date of diagnosis was 4 years.

The diagnosis of T2D and IFG were in accordance with the 1999 World Health Organization (WHO) criteria (30). T2D: fasting plasma glucose (FPG) levels ≥7.0 mmol/L or 2 h glucose levels ≥11.1 mmol/L in oral glucose tolerance test (OGTT). IFG: 6.1 ≤ FPG <7.0 mmol/L. Patients with a past history of using antidiabetic treatment or any endocrine disease other than T2D, were excluded from the study. Other exclusion criteria were described in previous study (23). This study was approved by the university ethical committee, and informed consent was obtained from each participant.

Data Collection

All subjects attended a standardized examination protocol. A structured questionnaire was used to collect information including demographic data, medical history, current medication, and health-related behaviors (10). The Chinese translation and adaptation of the short standard version of the Copenhagen Psychosocial Questionnaire (COPSOQ) was used to assess job-related psychosocial stress (10, 31, 32). The COPSOQ version used in this study comprises 5 scales including 34 items concerning (1) demands at work, (2) work organization and content, (3) interpersonal relations and leadership, (4) insecurity at work, and (5) job satisfaction. All scales had a 5-point Likert format (except scales measuring insecurity at work), reflecting either intensity or frequency. These scales were transformed into a theoretical range, extending from 0 (do not agree at all) to 100 (fully agree) points. The higher score for scales of “demands at work” and “insecurity at work” and lower score for scales of “work organization and content”, “interpersonal relations and leadership”, and “job satisfaction” indicate higher level of job strain.

Blood Sample Collection and RNA Extraction

Following an at least 8 h overnight fast, venous blood samples were collected in EDTA tubes from each subject between 7:30 and 8:30 a.m. in a calm state; 2 ml of the blood samples were immediately centrifuged to retrieve plasma. Another 2 ml blood samples were used to isolate peripheral blood mononuclear cells (PBMCs) by density gradient centrifugation. Total RNA was extracted from PBMCs using TRIzol reagent (Invitrogen, Waltham, MA, USA), following the manufacturer’s instructions. The purity and concentration of extracted RNA were determined by Nano Drop 2000 spectrophotometer, and the integrity was inspected by Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA).

Anthropometric Measurement and Biochemical Analysis

Anthropometric parameters including weight, height, waist circumference (WC), and blood pressure were obtained using standard measurement (26). Body mass index (BMI) was calculated as weight (kg)/height (m)². FPG, glycated hemoglobin (HbAlc), triglyceride (TG), total cholesterol (TCH), high-density lipoprotein cholesterol (HDLC), and low-density lipoprotein cholesterol (LDLC) were tested through standard procedures in the clinical laboratory of the hospital (23). Plasma concentrations of cortisol and insulin were measured by radioimmunoassay (RIA) kit. The homeostasis model assessment index of insulin resistance (HOMA-IR) calculated as “fasting insulin (μIU/ml) × fasting glucose (mmol/L)/22.5.”

Quantitative Real-Time PCR

The expression levels of hsa_circ_0111707, miR-144-3p, and NR3C1 were measured using quantitative real-time PCR (qRT-PCR). RNA was reverse transcribed into cDNA with PrimeScript RT Reagent Kit (Takara, Beijing, China) and miScript Reverse Transcriptase Kit (Qiagen, Hilden, Germany), respectively, according to manufacturer’s instructions. The real-time PCR analyses were performed using SYBR Premix Ex TaqII (Takara) with LightCycler 480 II Real-Time PCR Instrument (Roche, Basel, Switzerland). β-Actin was used as an internal reference for quantification of circRNA and mRNA, while miR-451a for miRNA (miR-451a expression was stable and consistent throughout all the evaluated samples) (23). The specific primers used are listed in Supplementary Table S1. The relative expression of genes was calculated by 2^−ΔΔCT method.

Cell Culture

Human embryonic kidney 293T (HEK293T) cells were cultured in Dulbecco’s modified Eagle’s medium (DEME, HyClone, Logan, UT, USA) containing 10% fetal bovine serum (FBS). Human skeletal muscle myoblasts (HSkMM, ScienCell, Carlsbad, CA, USA) were cultured with Skeletal Muscle Cell Medium (Cat. #3501, ScienCell, USA). All these cell lines were maintained at 37°C with 5% CO₂ in a humidified incubator.

Luciferase Reporter Assay

The full length of 346 bp of hsa_circ_0111707 and its corresponding mutant version without miR-144-3p binding site were constructed and inserted downstream of the luciferase reporter gene of the plasmid psiCHECK2 (The fragment sequences were listed in Supplementary Table S2). All these plasmids were validated by sequencing. HEK293T cells (1 × 10⁵ cells/well) were seeded in 96-well plates overnight and transfected with a luciferase reporter plasmid and miR-144-3p mimic or negative control (NC) using the Lipofectamine 2000 transfection reagent (Invitrogen), whereas empty psiCHECK2 was used as control plasmid. After 48 h, the activities of firefly and Renilla luciferase were detected with dual-luciferase reporter assay Kit (Promega, Madison, WI, USA). The relative activity was equal to firefly luciferase activity/Renilla luciferase activity. All experiments were independently in triplicate.

RNA Immunoprecipitation Assay

Skeletal muscle is an important peripheral tissue that utilizes blood glucose. Cortisol impairs glucose uptake in peripheral tissues like muscle and fat, which may facilitate insulin resistance (33). To further support the connection between circRNA and miRNA in skeletal muscle cells, the RIP assay was performed by a Magna RIP Kit (Millipore, Burlington, MA, USA) according to the manufacturer’s methods. HSkMM cells were transfected with miR-144-3p mimics or miR-NC, and then lysed in complete RNA lysis buffer after 48 h. Cell lysates were divided into three fractions. Two of them then incubated with magnetic beads which were conjugated with anti-Argonaute2 (AGO2) or negative control IgG antibody overnight at 4°C. Immunoprecipitated RNAs were extracted to detect the target RNAs including miR-144-3p and hsa_circ_0111707 expressions by RT-qPCR. The other one fraction was used as blank control (input).

Statistical Analysis

The data were represented as the mean ± standard deviation (SD) or percentages when they fit. The normal distribution of the data was checked using the Kolmogorov-Smirnov test. Student’s t-test or one-way analysis of variance (ANOVA) were performed to examine the differences between the study groups if the data obeyed normal distribution, otherwise nonparametric test (Mann-Whitney U test or Kruskal-Wallis H test) was used. Chi-square test was used to compare the categorical data among groups. Spearman’s correlation coefficient was used to evaluate the relationship between relative expression of hsa_circ_0111707 and stress-related variables. The association between circRNA expression and risk of T2D was analyzed using covariate-adjusted logistic regression and relative risk (RR). Receiver operator characteristic (ROC) curve was used to identify the best cut-off value of circRNA expression to predict T2D. Statistical analyses were performed by using SPSS 24.0, GraphPad Prism 5.01 and R 3.5.1. The level of statistical significance was set at two-sided P< 0.05.

Results

Basic Characteristics of the Participants in the First Part of the Study

The demographic and clinical characteristics of the study participants in the first part of the study (the cross-sectional study) are presented in Table 1. There was no significant difference in the distribution of age, gender, smoking, alcohol use, and physical inactivity among the three groups. The average levels of most of the clinical parameters, including BMI, WC, TCH, TG, LDLC, HDLC, and HOMA-IR, were significantly different among the three compared groups (p < 0.05).

TABLE 1

Table 1 Demographic and clinical characteristics of study participants in the first part of the study.

As to the COPSOQ investigation, those with abnormal FPG, including T2D and IFG, showed a higher level of psychosocial stress than those with normal FPG (Table 2). Compared with the control group, the average scores on scales of demands at work and insecurity at work were significantly higher in the T2D and IFG groups (p < 0.05). However, no significant differences of the average scores of these two scales were observed between T2D and IFG groups.

TABLE 2

Table 2 Association between COPSOQ variables and T2D.

Comparison of Relative Expression of hsa_circ_0111707, miR-144-3p, and NR3C1 Among the Study Groups

The relative expression of hsa_circ_0111707, miR-144-3p, and NR3C1 in PBMCs and plasma concentration of cortisol showed significantly differences among the T2D, IFG and control groups (Figures 2A–C; Table 1). Fold changes (FCs) of expression of the compared groups were also analyzed and presented in Supplementary Table S3. The relative expression of hsa_circ_0111707 and NR3C1 were significantly lower in the T2D and IFG cases in comparison with the controls, being lower in the T2D cases in comparison with IFGs (Figures 2A, C). By contrast, as the miRNA target of hsa_circ_0111707 and stress hormone, respectively, relative expression of miR144-3p and plasma cortisol concentration were significantly higher in the T2D and IFG cases in comparison with the controls, being higher in the T2D cases in comparison with IFGs (Figure 2B and Table 1).

FIGURE 2

Figure 2 Relative expression of hsa_circ_0111707, miR-144-3p, and NR3C1 and Spearman’s correlation between hsa_circ_0111707 and miR-144-3p, NR3C1, and cortisol concentration. (A–C) Comparison between T2D, IFG, and control groups. (D–F) Spearman’s correlation between relative expression of hsa_circ_0111707 and miR-144-3p, NR3C1, and cortisol concentration in all the subjects. (G–I) Spearman’s correlation between relative expression of hsa_circ_0111707 and miR-144-3p, NR3C1, and cortisol concentration in subgroups. T2D, type 2 diabetes; IFG, impaired fasting glucose. ^*p < 0.05, ^**p < 0.01, ^***p < 0.001.

The Relationship Between hsa_circ_0111707 Expression and Stress-Related Variables

The Spearman’s correlation analysis revealed that the relative expression of hsa_circ_0111707 was significantly negatively correlated with miR-144-3p and plasma cortisol concentration and positively correlated with the relative expression of NR3C1 (Figures 2D–F). In subgroup analysis, significant correlations were also observed between hsa_circ_0111707 and miR-144-3p, NR3C1, and cortisol in the T2D, IFG, and control groups, respectively (Figures 2G–I). At the same time, significant correlations also existed between miR-144-3p, NR3C1, and cortisol, in all the subjects and subgroups (Supplementary Table S4), which indicated the potential regulation relationship among them. Relative expression of NR3C1 was negatively associated with miR-144-3p and cortisol concentration, while positive association existed between miR-144-3p and cortisol (Supplementary Table S4).

The relationship between relative expression of hsa_circ_0111707 and job-related psychosocial stress assessed by COPSOQ was also analyzed, which showed that relative expression of hsa_circ_0111707 was negatively correlated with scores on scales of demands at work and insecurity at work. Whereas, plasma cortisol concentration was positively correlated with scores of demands at work and insecurity at work (Table 3).

TABLE 3

Table 3 Spearman’s correlation matrix of hsa_circ_ 0111707 expression and COPSOQ variables.

The Relationship Between hsa_circ_0111707 Expression and HOMA-IR

As insulin resistance is the important pathophysiological basis of T2D, the relationship between relative expression of hsa_circ_0111707 and HOMA-IR was assessed. Significantly, negative correlations were observed between hsa_circ_0111707 expression and HOMA-IR in all the subjects (r = −0.447, p < 0.001), T2D (r = −0.278, p = 0.003), and IFG groups (r = −0.203, p = 0.039). Multiple linear regression analysis was further performed to confirm the association between hsa_circ_0111707 expression and HOMA-IR, with potential confounders being adjusted (Supplementary Table S5). The results showed that the relative expression of hsa_circ_0111707 was a significant negative predictor of HOMA-IR (p < 0.001), indicating that decreased hsa_circ_0111707 expression may contribute to insulin resistance. If BMI was replaced by WC among the above covariables in the model, hsa_circ_0111707 expression was still associated with HOMA-IR independent of obesity index.

The Association Between hsa_circ_0111707 Expression and Risk of T2D

In the second part of the study (nested case-control study), the demographic and clinical characteristics of the newly diagnosed T2D cases and their matched controls at 4 years prior to the diagnosis (baseline) were presented in Supplementary Table S6. The relative expression of hsa_circ_0111707 in the T2D cases (4.99 ± 1.69) was significantly lower than the controls (6.21 ± 2.11, p <0.001). The adjusted RR for the per-unit decrease of hsa_circ_0111707 expression was 1.38 (95% CI: 1.17–1.63) in the fully adjusted model (Figure 3). When hsa_circ_0111707 expression was grouped according to the median (5.375), participants with lower expression of hsa_circ_0111707 had a 2.72-fold (95% CI: 1.43–5.18) increased risk for T2D, compared with those with higher expression level (Figure 3). The ROC analyses showed that the area under curve (AUC) was 0.667 (95% CI: 0.590–0.744) for hsa_circ_0111707 (p < 0.001). The relative expression value of 5.845 was identified as the best cutoff point of hsa_circ_0111707 to predict T2D (sensitivity: 0.716, specificity: 0.568).

FIGURE 3

Figure 3 The association between relative expression of hsa_circ_0111707 assessed by univariate and multiple logistic regression analyses. Model 1: adjusted for smoking, drinking, and physical activity. Model 2: further adjusted for BMI based on model 1. Model 3: further adjusted for TC, TG, HDLC, LDLC, SBP, and DBP based on model 2. RR, relative risk; CI, confidence interval; T2D, type 2 diabetes; IFG, impaired fasting glucose.

hsa_circ_0111707 Functions as a Sponge for miR-144-3p

The luciferase assays showed that the average relative luciferase activity of cells cotransfected with miR-144-3p and hsa_circ_0111707-WT was significantly lower than that in cells cotransfected with miR-144-3p and hsa_circ_0111707-MUT (p < 0.001, Figures 4A, B). It has been known that miRNAs regulate target gene expression by binding to AGO2, the key component of RNA-induced silencing complex (RISC). Thus, an anti-AGO2 RIP assay was conducted in HSkMM cells to pull down the RNA transcripts which bind to AGO2 with anti-AGO2 antibody and IgG as a NC. The results showed that both hsa_circ_0111707 and miR-144-3p were significantly immunoprecipitated by anti-AGO2 antibody compared with IgG and enriched by miR-144-3p mimics compared with miR-NC (Figure 4C, D). These data demonstrated that hsa_circ_0111707 acts as a ceRNA to sponge miR-144-3p by directly binding to it.

FIGURE 4

Figure 4 hsa_circ_0111707 functions as a sponge for miR-144-3p. (A) The binding sites of hsa_circ_0111707 and miR-144-3p. (B) Luciferase assays revealed that miR-144-3p is a miRNA target of hsa_circ_0111707 in HEK293T. (C, D) Anti-AGO2 RIP was executed in HSkMM cells after transfection with miR-144-3p mimic or miR-NC, followed by qRT-PCR to detect the enrichment (expression) of hsa_circ_0111707 and miR-144-3p, respectively. HEK, human embryonic kidney; HSkMM, human skeletal muscle myoblasts; WT, wild type; MUT, mutant; Blank: empty plasma; Input: blank control. ^***p < 0.001. ns, not significant.

Discussion

The present study identified the association between relative expression of hsa_circ_0111707 in PBMCs and T2D in an occupational cohort. Functional analysis demonstrated that hsa_circ_0111707 acts as a sponge for miR-144-3p, which is a direct regulator of NR3C1.

In this study, a valid and relevant short version of the COPSOQ instrument, which covers a broad range of aspects of currently leading concepts and theories (① the job characteristics model, ② the Michigan organizational stress model, ③ the demand-control-support model, ④ the sociotechnical approach, ⑤the action-theoretical approach, ⑥the effort-reward-imbalance model, ⑦the vitamin model), was used to assess psychosocial stress at work (31 24). The COPSOQ is a comprehensive instrument for the assessment of psychosocial factors at work. The main advantage of the COPSOQ is its generic usability in all sorts of professions and industries. The user-friendly version of the short COPSOQ, which almost reaches the same level of criterion validity as the long version, is more suitable for large-scale survey (34).

We found that subjects with abnormal increased FPG showed a higher level of work-related psychosocial stress than those with normal FPG. The mean scores of demands at work and insecurity at work” of the COPSOQ were significantly higher in the T2D and IFG groups than those in the control group, suggesting that psychosocial stress is associated with T2D. The significant association observed between COPSOQ and cortisol was consistent with the concept that, as a stimulus, psychological stress activates the HPA axis and results in cortisol rise (13). Decreased relative expression of NR3C1 in PBMCs and increased plasma cortisol concentration in T2D and IFG cases further support the association between psychological stress and T2D. The negative correlation between relative expression of NR3C1 and cortisol concentration can be explained that reduced NR3C1 expression helps to alleviate excessive stress-related physiological processes (17).

In our previous study, NR3C1 was proved to be a direct target gene of miR-144-3p by luciferase assay (24). miR-144 has been identified as a stressor-responsive miRNA by several studies. A human miRNA microarray and subsequent RT-qPCR confirmed that elevated miR-144 in whole blood is a potential biomarker for naturalistic stress in healthy young adults (35). miR-144 has also been reported as a potential signature miRNA for detection of T2D based on the evidence that miR-144 expression was highly upregulated in T2D, and it seemed to exhibit a linear relationship with increasing glycemic status and insulin resistance (23, 36). We found that significantly increased miR144-3p expression and decreased NR3C1 expression in PBMCs among subjects with abnormal increased FPG, and significantly negative correlation between them, which support the relationship that miR-144-3p negatively regulator of NR3C1.

We found that the relative expression of hsa_circ_0111707 significantly decreased in the PBMCs in T2D and IFG cases. At the same time, significant negative relationship existed between hsa_circ_0111707 expression and scores of demands at work and insecurity at work. miR-144-3p was further identified as a miRNA target of hsa_circ_0111707 in this study. We inferred that decreased hsa_circ_0111707 expression might increase the level of target miRNA and miRNA hyper expression promotes degradation of downstream target gene (by binding to 3-untranslated region of mRNA), thereby affecting the susceptibility of T2D development. At population level, the relative expression of hsa_circ_0111707 was significantly negatively correlated with miR-144-3p and cortisol concentration, and positively correlated with expression of NR3C1. These results were in conjunction with our hypothesis that hsa_circ_0111707 can compete for miRNA-binding sites and make it function as miR-144-3p sponge, thus suppressing miR-144-3p activity and participating in regulation of downstream target gene NR3C1 of miR-144-3p.

Based on the cross-sectional study, the predictive efficiency of hsa_circ_0111707 for T2D was further investigated in a nested case-control study, which allows us to estimate the causal relationship of hsa_circ_0111707 on the development of T2D. The measurement of hsa_circ_0111707 expression and other potential risk factors were prior to the diagnosis of T2D. The decreased hsa_circ_0111707 expression 4 years prior to the diagnosis of T2D was associated with increased risk of T2D development. hsa_circ_0111707 expression was an independent predictor of T2D, and insulin resistance with covariates was adjusted. The best cutoff value of hsa_circ_0111707 expression to predict T2D was also identified using ROC analysis.

Over the past few years, accruing evidence has revealed ncRNAs including microRNAs, long noncoding RNAs, and circular RNAs regulate pivotal cellular and molecular processes during all stages of a number of diseases (37). It is known that circulating ncRNAs are involved in cell-to-cell communication and correlated with the expression of target mRNAs in pathways that important to the development of disease (29, 38). The alternation expression of these ncRNAs in circulation has been demonstrated as valuable biomarkers in clinical practices (39, 40). As a window, PBMCs convert psychosocial stress into cellular dysfunction and finally contribute to the pathophysiology of lifestyle-related diseases such as T2D, cardiovascular disease, and atherosclerosis (41). In the present study, by identification of the association between hsa_circ_0111707 and stress-related T2D, it would be helpful for discovering novel targets for the prevention, prediction, and more personalized treatment of neuroendocrine stress-related disorders. Our results are reliable based on the three-part design. However, there are several limitations in our study. Firstly, the present study was performed in an occupational cohort, which affects the extrapolation of the results. Further prospective studies with larger sample sizes are needed to confirm our findings. Secondly, experimental studies are needed to further unveil the molecular mechanisms of the roles of hsa_circ_0111707 in T2D development.

In conclusion, our findings suggest that hsa_circ_0111707 is associated with risk of T2D development and may be involved in downregulation of NR3C1 via sponging miR-144-3p. hsa_circ_0111707 in PBMCs can be considered a potential biomarker of stress-related T2D.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Capital Medical University ethical committee. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

Y-XY designed the study and drafted the manuscript. H-BX, YS, JD, and L-JW collected the data and samples. Y-XY and L-JW conducted statistical analyses. Y-XY, YS, and SW conducted the experiments. All authors interpreted the data and have approved the final manuscript.

Funding

This study was supported by the National Natural Science Foundation (81773511, 81573214) and the Beijing Municipal Natural Science Foundation (7162020).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.790591/full#supplementary-material

Abbreviations

AUC, area under the curve; BMI, body mass index; circRNA, circular RNA; ceRNA: competitive endogenous RNA; COPSOQ, Copenhagen Psychosocial Questionnaire; DBP, diastolic blood pressure; FBS, fetal bovine serum; FPG, fasting plasma glucose; FC, fold change; HbA1c, glycated hemoglobin; HDLC, high-density lipoprotein-cholesterol; HEK, human embryonic kidney; HSkMM, human skeletal muscle myoblasts; HOMA-IR, homeostasis model assessment of insulin resistance; IFG, impaired fasting glucose; IR, insulin resistance; LDLC, low-density lipoprotein-cholesterol; miRNA, microRNA; MUT, mutant; NC, negative control; ncRNA, noncoding RNA; PBMC, peripheral blood mononuclear cell; ROC, receiver operating characteristic curve; RR, relative risk; qRT-PCR, real-time quantitative polymerase chain reaction; RIP: RNA immunoprecipitation; SBP, systolic blood pressure; TCH, total cholesterol; TG, triglyceride; T2D, type 2 diabetes; WC, waist circumference; WT, wide type.

References

1. American Diabetes Association. 2. Classification and Diagnosis of Diabetes: Standards of Medical Care in Diabetes-2020. Diabetes Care (2020) 43(Suppl 1):S14–31. doi: 10.2337/dc20-S002

PubMed Abstract | CrossRef Full Text | Google Scholar

2. International Diabetes Federation. IDF Diabetes Atlas (2015). Available at: http://www.diabetesatlas.org/ (Accessed Oct 6, 2016).

Google Scholar

3. Kolb H, Martin S. Environmental/lifestyle Factors in the Pathogenesis and Prevention of Type 2 Diabetes. BMC Med (2017) 15:131. doi: 10.1186/s12916-017-0901-x

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Morris T, Moore M, Morris F. Stress and Chronic Illness: The Case of Diabetes. J Adult Dev (2011) 18:70–80. doi: 10.1007/s10804-010-9118-3

CrossRef Full Text | Google Scholar

5. Tamashiro KL, Sakai RR, Shively CA, Karatsoreos IN, Reagan LP. Chronic Stress, Metabolism, and Metabolic Syndrome. Stress (2011) 14:468–74. doi: 10.3109/10253890.2011.606341

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Nyberg ST, Fransson EI, Heikkilä K, Ahola K, Alfredsson L, Bjorner JB, et al. Job Strain as a Risk Factor for Type 2 Diabetes: A Pooled Analysis of 124,808 Men and Women. Diabetes Care (2014) 37:2268–75. doi: 10.2337/dc13-2936

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rod NH, Grønbaek M, Schnohr P, Prescott E, Kristensen TS. Perceived Stress as a Risk Factor for Changes in Health Behaviour and Cardiac Risk Profile: A Longitudinal Study. J Intern Med (2009) 266:467–75. doi: 10.1111/j.1365-2796.2009.02124.x

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wiernik E, Nabi H, Thomas F, Pannier B, Hanon O, Simon T, et al. Association Between Current Perceived Stress and Incident Diabetes is Dependent on Occupational Status: Evidence From the IPC Cohort Study. Diabetes Metab (2016) 42(5):328–35. doi: 10.1016/j.diabet.2016.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Novak M, Björck L, Giang KW, Heden-Ståhl C, Wilhelmsen L, Rosengren A. Perceived Stress and Incidence of Type 2 Diabetes: A 35-Year Follow-Up Study of Middle-Aged Swedish Men. Diabetes Med (2013) 30:e8–16. doi: 10.1111/dme.12037

CrossRef Full Text | Google Scholar

10. Yan YX, Xiao HB, Wang SS, Zhao J, He Y, Wang W, et al. Investigation of the Relationship Between Chronic Stress and Insulin Resistance in a Chinese Population. J Epidemiol (2016) 26:355–60. doi: 10.2188/jea.JE20150183

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Lacroix A, Feelders RA, Stratakis CA, Nieman LK. Cushing’s Syndrome. Lancet (2015) 386:913–27. doi: 10.1016/S0140-6736(14)61375-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Clore JN, Thurby-Hay L. Glucocorticoid-Induced Hyperglycemia. Endocr Pract (2009) 15:469–74. doi: 10.4158/EP08331.RAR

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Nicolaides NC, Galata Z, Kino T, Chrousos GP, Charmandari E. The Human Glucocorticoid Receptor: Molecular Basis of Biologic Function. Steroids (2010) 75:1–12. doi: 10.1016/j.steroids.2009.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Matthews SG, McGowan PO. Developmental Programming of the HPA Axis and Related Behaviours: Epigenetic Mechanisms. J Endocrinol (2019) 242:T69–79. doi: 10.1530/JOE-19-0057

PubMed Abstract | CrossRef Full Text | Google Scholar

15. McEwen BS. Central Effects of Stress Hormones in Health and Disease: Understanding the Protective and Damaging Effects of Stress and Stress Mediators. Eur J Pharmacol (2008) 583:174–85. doi: 10.1016/j.ejphar.2007.11.071

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Zhou J, Cidlowski JA. The Human Glucocorticoid Receptor: One Gene, Multiple Proteins and Diverse Responses. Steroids (2005) 70:407–17. doi: 10.1016/j.steroids.2005.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Raison CL, Miller AH. When Not Enough is Too Much: The Role of Insufficient Glucocorticoid Signaling in the Pathophysiology of Stress-Related Disorders. Am J Psychiatry (2003) 160:1554–65. doi: 10.1176/appi.ajp.160.9.1554

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Yan YX, Dong J, Liu YQ, Zhang J, Song MS, He Y, et al. Association of Suboptimal Health Status With Psychosocial Stress, Plasma Cortisol and mRNA Expression of Glucocorticoid Receptor α/β in Lymphocyte. Stress (2015) 18:29–34. doi: 10.3109/10253890.2014.999233

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chen K, Rajewsky N. The Evolution of Gene Regulation by Transcription Factors and microRNAs. Nat Rev Genet (2007) 8:93–103. doi: 10.1038/nrg1990

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Dirven BCJ, Homberg JR, Kozicz T, Henckens MJAG. Epigenetic Programming of the Neuroendocrine Stress Response by Adult Life Stress. J Mol Endocrinol (2017) 59:R11–31. doi: 10.1530/JME-17-0019

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Haramati S, Navon I, Issler O, Ezra-Nevo G, Gil S, Zwang R, et al. MicroRNA as Repressors of Stress-Induced Anxiety: The Case of Amygdalar miR-34. J Neurosci (2011) 31:14191–203. doi: 10.1523/JNEUROSCI.1673-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Mannironi C, Camon J, De Vito F, Biundo A, De Stefano ME, Persiconi I, et al. Acute Stress Alters Amygdala microRNA miR-135a and miR-124 Expression: Inferences for Corticosteroid Dependent Stress Response. PloS One (2013) 8:e73385. doi: 10.1371/journal.pone.0073385

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Liang YZ, Dong J, Zhang J, Wang S, He Y, Yan YX. Identification of Neuroendocrine Stress Response-Related Circulating microRNAs as Biomarkers for Type 2 Diabetes Mellitus and Insulin Resistance. Front Endocrinol (Lausanne) (2018) 9:132. doi: 10.3389/fendo.2018.00132

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Yan YX, Xiao HB, Zhang J, Wang S, Dong J, Wu LJ. Pri-miR-144 Rs9279 is Associated With Type 2 Diabetes and Regulation of Stress Response. J Cell Physiol (2021) 236:561–9. doi: 10.1002/jcp.29883

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, et al. Natural RNA Circles Function as Efficient microRNA Sponges. Nature (2013) 495:384–8. doi: 10.1038/nature11993

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chen LL, Yang L. Regulation of circRNA Biogenesis. RNA Biol (2015) 12(4):381–8. doi: 10.1080/15476286.2015.1020271

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Dong Y, He D, Peng Z, Peng W, Shi W, Wang J, et al. Circular RNAs in Cancer: An Emerging Key Player. J Hematol Oncol (2017) 10:2. doi: 10.1186/s13045-016-0370-2

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Li JJ, Wang W, Wang XQ, He Y, Wang SS, Yan YX. A Novel Strategy of Identifying circRNA Biomarkers in Cardiovascular Disease by Meta-Analysis. J Cell Physiol (2019) 234:21601–12. doi: 10.1002/jcp.28817

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Lu YK, Chu X, Wang S, Sun Y, Zhang J, Dong J, et al. Identification of Circulating Hsa_Circ_0063425 and Hsa_Circ_0056891 as Novel Biomarkers for Detection of Type 2 Diabetes. J Clin Endocrinol Metab (2021) 106:e2688–99. doi: 10.1210/clinem/dgab101

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Department of Noncommunicable Disease Surveillance. Definition, Diagnosis and Classification of Diabetes Mellitus and its Complications: Report of a WHO Consultation. In: Part 1. Diagnosis and Classification of Diabetes Mellitus. Geneva: World Health Organization (1999).

Google Scholar

31. Kristensen TS, Hannerz H, Høgh A, Borg V. The Copenhagen Psychosocial Questionnaire-A Tool for the Assessment and Improvement of the Psychosocial Work Environment. Scand J Work Environ Health (2005) 31:438–49. doi: 10.5271/sjweh.948

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Shang L, Liu P, Fan LB, Gu HK, Li J. Psychometric Properties of the Chinese Version of Copenhagen Psychosocial Questionnaire. J Environ Occup Med (2008) 25:572–6.

Google Scholar

33. Sjöstrand M, Eriksson JW. Neuroendocrine Mechanisms in Insulin Resistance. Mol Cell Endocrinol (2009) 297(1-2):104–11. doi: 10.1016/j.mce.2008.05.010

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Nübling M, Stößel U, Hasselhorn HM, Michaelis M, Hofmann F. Measuring Psychological Stress and Strain at Work - Evaluation of the COPSOQ Questionnaire in Germany. Psychosoc Med (2006) 3:Doc05.

PubMed Abstract | Google Scholar

35. Katsuura S, Kuwano Y, Yamagishi N, Kurokawa K, Kajita K, Akaike Y, et al. MicroRNAs miR-144/144^* and miR-16 in Peripheral Blood are Potential Biomarkers for Naturalistic Stress in Healthy Japanese Medical Students. Neurosci Lett (2012) 516:79–84. doi: 10.1016/j.neulet.2012.03.062

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Karolina DS, Armugam A, Tavintharan S, Wong MT, Lim SC, Sum CF, et al. MicroRNA 144 Impairs Insulin Signaling by Inhibiting the Expression of Insulin Receptor Substrate 1 in Type 2 Diabetes Mellitus. PloS One (2011) 6:e22839. doi: 10.1371/journal.pone.0022839

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Goodall GJ, Wickramasinghe VO. RNA in Cancer. Nat Rev Cancer (2021) 21(1):22–36. doi: 10.1038/s41568-020-00306-0

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Mirra P, Raciti GA, Nigro C, Fiory F, D’Esposito V, Formisano P, et al. Circulating miRNAs as Intercellular Messengers, Potential Biomarkers and Therapeutic Targets for Type 2 Diabetes. Epigenomics (2015) 7:653–67. doi: 10.2217/epi.15.18

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Di Mauro S, Scamporrino A, Fruciano M, Filippello A, Fagone E, Gili E, et al. Circulating Coding and Long non-Coding RNAs as Potential Biomarkers of Idiopathic Pulmonary Fibrosis. Int J Mol Sci (2020) 21:8812. doi: 10.3390/ijms21228812

CrossRef Full Text | Google Scholar

40. Di Mauro S, Scamporrino A, Petta S, Urbano F, Filippello A, Ragusa M, et al. Serum Coding and non-Coding RNAs as Biomarkers of NAFLD and Fibrosis Severity. Liver Int (2019) 39(9):1742–54. doi: 10.1111/liv.14167

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Bierhaus A, Wolf J, Andrassy M, Rohleder N, Humpert PM, Petrov D, et al. A Mechanism Converting Psychosocial Stress Into Mononuclear Cell Activation. Proc Natl Acad Sci USA (2003) 100(4):1920–5. doi: 10.1073/pnas.0438019100

PubMed Abstract | CrossRef Full Text | Google Scholar