Study on inflammation-related genes and microRNAs, with special emphasis on the vascular repair factor HGF and miR-574-3p, in monocytes and serum of patients with T2D
Diabetology & Metabolic Syndrome volume 8, Article number: 6 (2016)
Recently, we reported signs of inflammation (raised IL-8, reduced miR-146a) and signs of vascular repair (raised HGF) in the serum of Ecuadorian patients with type 2 diabetes (T2D). In contrast, we found that the circulating monocytes lacked up-regulation of classical inflammatory genes (IL-1B, IL-6, and TNF) and there was even significant down-regulation of PTGS2. Notably, genes and a microRNA involved in adhesion, cell differentiation and morphology (CD9, DHRS3, PTPN7 and miR-34c-5p) were up-regulated in the T2D monocytes, suggesting a role of the anti-inflammatory cells in adhesion, vascular repair and invasion.
To determine the gene expression of the vascular repair factor HGF in the circulating monocytes of patients with T2D and to investigate the relationship between HGF and the expression of the other previously tested monocyte genes and the contribution to the raised serum level of HGF. In addition, we tested the level of 6 microRNAs, which were previously found abnormal in the circulating monocytes, in the serum of the patients.
A gene and microRNA expression study in monocytes and serum of 64 Ecuadorian patients with T2D (37–85 years) and 44 non-diabetic controls (32–87 years).
The gene expression of HGF was significantly raised in the monocytes of the patients with T2D and associated with the expression of genes involved in adhesion, cell differentiation and morphology. HGF gene expression did not correlate with the serum level of HGF. The monocyte expression of pro-inflammatory cytokine genes was also not associated with the serum levels of these cytokines. The level of miR-574-3p was significantly decreased in the serum of the patients with T2D, and correlated in expression with the decreased well-established inflammation-regulating miR-146a. The level of the microRNAs in serum did not correlate with their expression level in monocytes.
In circulating monocytes of Ecuadorian T2D patients, the microRNA and gene expression of important inflammatory/chemotactic/motility/vascular repair factors differs from the expression in serum. While monocytes show a gene expression profile compatible with an anti-inflammatory state, serum shows a molecular profile compatible with an inflammatory state. Both compartments show molecular signs of vascular repair support, i.e. up-regulated HGF levels.
There is increasing evidence that monocytes, macrophages and related cells are closely involved in the pathogenesis of the metabolic syndrome (MetS) and type 2 diabetes (T2D). Importantly, in obesity the number of macrophages increases from 10–15 % to 50–60 % of total adipose tissue cells [1, 2]. The increase in macrophage number is accompanied by a hyper activation of the cells and leads to a raised secretion of pro-inflammatory cytokines (TNF-α, IL-1β, IL-6, CCL-4, PAI-1) and chemokines (CCL2) causing a state of chronic low-grade inflammation [1, 3–7] and insulin resistance.
Important circulating precursors for the macrophages in the adipose tissue [8–10] are the blood-borne monocytes. There is a relative paucity on the state of inflammatory activation of circulating monocytes in patients with the MetS [11, 12] and T2D , but increases in pattern recognition receptors, oxidative stress and the machinery for the production of pro-inflammatory cytokines have all been described [14–16].
Contrary to this view, we recently reported that monocytes of a group of 64 Ecuadorian patients with T2D were characterized by an anti-inflammatory set point rather than a pro-inflammatory set point when compared to monocytes of a group of 44 non-diabetic controls. We found a decrease in expression of a cluster of 11 mutually correlated core inflammatory cytokine/compound genes (IL-1B, IL-6, TNF, TNFAIP3, PTGS2, CCL20, CCL2, CCL4, PDE4B, DUSP2 and ATF3; reaching significance for PTGS2) in the monocytes of patients with T2D . However, the study on the monocytes of the patients with T2D also showed that there was up-regulated gene expression for genes occurring in a cluster of mutually correlating genes, many of which are involved in adhesion, migration, cell differentiation and cell morphological change . A significant up-regulation as compared to non-diabetic controls was reached for the genes CD9, DHRS3 and PTPN7. Other important genes in this gene cluster were MAPK6, NAB2, STX1A, EMP-1, CDC42, DHRS3, FABP5, BCL2A1, PTX3 and CXCL2. We interpreted these data as indicating that circulating monocytes in our group of patients with T2D were activated, not towards an inflammatory state, but to a state of enhanced adhesion, migration and further differentiation into descendent cell types, most likely into monocyte-derived pro-angiogenic cells, instrumental in vascular repair. This view was further supported by our observation that the expression of miR-146a, a well-known inflammation down-regulating microRNA, was not changed in the T2D monocytes, while a microRNA targeting genes involved in processes of cell morphology and shape change, i.e. miR-34c-5p, was up-regulated as compared to the group of non-diabetic controls .
In the serum of the patients with T2D in whom we performed the monocyte studies, we found clear signs of inflammation . Although there were no increases in the levels of classical cytokines, such as of IL-1β, IL-6 and TNF-α, there was an increase in the level of serum IL-8, and also the level of miR-146a was significantly down-regulated. HGF was increased in the serum of the cases with T2D too . Since HGF is an important vascular repair factor [19–22] and an anti-inflammatory agent [23, 24], and monocyte-derived angiogenic cells are characterized by the expression of HGF , we hypothesized that there was an enhanced monocyte-linked endothelial repair mechanism going on in our patients with T2D.
In the present study, we therefore tested the hypothesis that HGF is increased in the circulating anti-inflammatory monocytes of patients with T2D and we determined the gene expression level of HGF (and the HGF-R, cMET) in the monocytes of our patients with T2D and investigated whether HGF belonged to the cluster of typical inflammatory compound genes or to the clusters of typical adhesion, migration and differentiation genes. In addition, we compared the monocyte gene expression levels of HGF to the serum HGF levels to investigate whether the circulating monocytes could be the main producers of this vascular repair factor in serum.
In addition, we measured the serum level of miR-34c-5p, the microRNA up-regulated in the monocytes of patients with T2D and playing a role in cell shape processes, to see whether this microRNA was also raised in the serum of the patients. Finally, we determined the serum level of the other 5 microRNAs (miR-122, miR-138, miR-410, miR-574-3p and miR-92), which we had previously reported as abnormally expressed in the monocytes of patients with T2D in a finding study .
Patients, materials and methods
A total of 64 subjects diagnosed with diabetes type 2, according to the criteria of The Expert Committee on the diagnosis and classification of Diabetes Mellitus . Patients were recruited in four medical centers of Quito, Ecuador (Eugenio Espejo Hospital, Club de Leones Sur, Fundación Oftalmológica del Valle and Fundación de la Psoriasis) from 2009 until 2012. For demographic and clinical details see Table 1. At the same time, 44 healthy controls with similar ethnical and social background, neither suffering from T2D nor other important medical disorders (including acute infection) served as controls. Controls had to be over 30 years of age, considering the age dependency of T2D , and preferably of the same gender as the patients.
Patients and controls with other immune disorders, other serious medical illnesses, recent infections (last 2 weeks), obvious vascular complications such as diabetic foot and ulcers, fever, pregnancy/postpartum, use of statins and LADA patients (patients positive for GAD-65 Abs) were excluded. The Medical Ethical Review Committee of the Ecuadorian Corporation of Biotechnology Quito, Ecuador and the Ethic Committee of the Central University of Quito approved the study. Written informed consent was obtained of all subjects participating in the study. The Ecuadorian Ministry of Health (MSP) gave the permit to export and process the samples in Erasmus MC, Rotterdam, and The Netherlands.
Blood collection and preparation
In the morning fasting venous blood was collected. Ten mL were collected in a clotting tube and processed within 4 h. Serum was frozen and stored at minus 80 °C for approximately 24 months before testing. Thirty milliliters were collected in tubes containing sodium-heparin for immune cell preparation. From the heparinized blood, peripheral blood mononuclear cell (PBMC) suspensions were prepared in the afternoon by low-density gradient centrifugation, as previously described in detail  within 8 h to avoid activation of the monocytes. PBMCs were frozen in 10 % dimethylsulfoxide and stored in liquid nitrogen. This enabled us to test patient and control serum and immune cells in the same series of experiments later.
Isolation of monocytes
CD14-positive (CD14+) monocytes were isolated from frozen PBMCs by a magnetic cell sorting system (MACS; Miltenyi Biotec, Auburn, CA, USA). The purity of monocytes was >95 % (determined by morphological screening after Trypan Blue staining and flow cytometric analysis). As previously reported; the positive vs. negative selection of immune cells did not influence gene expression profiles .
Real time quantitative PCR (qRT PCR) for monocytes
mRNA expression in monocytes via TaqMan array cards
For the previous report we had determined the expression of 24 mRNAs in the monocytes of the T2D cases and the non-diabetic controls. RNA had been isolated from monocytes using RNeasy columns (Qiagen, Hilden, Germany), and both this method and quantitative RT-PCR has been described in detail elsewhere . One microgram of RNA was reverse-transcribed using the High Capacity cDNA kit (Applied Biosystems, Foster City, CA, USA). qPCR was performed using custom TaqMan Arrays, format 48 (Applied Biosystems), according to the manufacturer’s protocol and validated against the single RT-qPCR method. Per fill port, 400 ng of cDNA (converted from total RNA) was loaded. PCR amplification was performed using an Applied Biosystems Prism 7900HT sequence detection system with TaqMan Array block. Thermal cycler conditions were 2 min at 50 °C, 10 min at 94.5 °C, and then 30 s at 97 °C, and 1 min at 59.7 °C for 40 cycles. The expressions of ATF3, BCL2A1, CCL20, CCL2, CCL7, CD9, CDC42, CXCL2, DHRS3, DUSP2, EMP1, FABP5, HSPA1A/HSPA1B, IL-1B, IL-6, MAPK6, NAB2, PDE4B, PTGS2, PTPN7, PTX3, STX1A, TNF, and TNFAIP3 were determined using this card system. The PCR amplification of the housekeeping gene ABL was performed for each sample to allow normalization between the samples. We chose ABL as the housekeeping gene because it has previously been shown that ABL was the most consistently expressed housekeeping gene in hematopoietic cells. The SDS software (ABI) was used to collect the data and the RQ Manager Program (ABI) was used to assign, check, and standardize CT values. The Data Assist software (ABI) was used to normalize the data against ABL expression. For threshold cycles below 40, the corresponding mRNA was considered detected, higher cycle numbers were not included in calculations. Data were expressed as cycle threshold (CT) values corrected to ABL (ΔCT = CT gene X − CT housekeeping gene ABL), and as fold change values determined via the CT method (User Bulletin 2; Applied Biosystems) (ΔΔCt method, formula 2ΔΔCt) to linear transform the data and to avoid negative values. To correct for inter-assay variance we set the mean of the studied genes found in the healthy control groups in the same assay for each gene to 1, and the fold change (FC) values of the genes in patient monocytes were expressed relative to this set mean of 1 of the healthy controls for the given values.
Individual mRNA qRT–PCR assays for HGF, HGF-R, resistin
For the current report we additionally determined in single assays (not using the card system) the gene expression for HGF, the HGF-R (cMET), and resistin using the same cDNA used in the above described experiments (we measured resistin because resistin had also been found raised in the serum of the T2D cases, though just nor reaching significance, p = 0.07) . TaqMan probes and consensus primers for HGF, HGF-R and resistin were provided by Applied Biosystems. PCR amplification of the housekeeping gene ABL was performed for each sample to allow normalization between the samples. Assays were carried out as described by the manufacturers in 15 µL assays. The fold change values between different groups were determined from the normalized CT values as described above for the outcomes of the card assays.
Individual microRNA qRT–PCR assays for monocytes
Total RNA was isolated from purified monocytes using RNeasy columns (Qiagen, Hilden, Germany) as described by the manufacturer’s manual. Purity and integrity of the RNA were assessed on the Agilent 2100 bioanalyzer with the RNA 6000 Nano LabChip® reagent set (Agilent Technologies, Santa Clara, CA, USA) and the RNA was then stored at −80 °C. Subsequently, specific stem-looped reverse transcription primers were used to obtain cDNA for mature microRNAs. The RNA was reverse transcribed using the TaqMan® MicroRNA Reverse Transcription Kit from Applied Biosystems, The Netherlands (ABI). PCR was performed using pre-designed TaqMan® microRNA assays and TaqMan® Universal Master Mix, NoAmpErase®UNG, with an ABI 7900 HT real-time PCR machine. RNU44 was used as reference microRNA. The PCR conditions were 2 min at 50 °C, 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C, and finally 1 min at 60 °C. The Data Assist software (ABI) was used to normalize the data to RNU44. Calculations of FC were as described above.
The levels of TNFα, IL-1β, IL-6, NGF, HGF, PAI, resistin, CCL2 (MCP-1), adiponectin, leptin, IL-8, and MIP1β (CCL4) were measured by flow cytometry (BD LSR II Biosciences, CA, and EE.UU.) using a commercially available multi-analyte cytometric bead array system (Milliplex® Map, USA). This multi-analyte assay panel was specifically designed for T2D studies by Millipore and is a relative economic assay using minimal quantities of serum (25 μL). The data were analyzed using a 5-parameter logistic method for calculating analyte concentrations in samples (MAGPIX® with xPONENT software, Luminex, Austin, USA). Serum analyte levels below the detection limit of the assay (only found with one analyte, i.e. IL-1β in 45 % of the determinations) were considered as 0 pg/ml and included as such in the statistical analysis (if e.g. half maximal values of the detection limit had been used slightly higher values would have been obtained, but significance levels had not changed for IL-1β). Lipid and hepatic profile was performed according to standard lab procedures in Quito-Ecuador (AMCOR laboratory) and validated in Erasmus MC.
Serum microRNA quantitative real-time PCR (qPCR) in serum
Total RNA was isolated from serum using the same Qiagen miRNeasy kit as used for the monocytes. In order to correct for variations in RNA isolation derived, we spiked-in a non-human (C. elegans) synthetic miRNA cel-miR-39 miRNA Mimic (MSY000010) into the sample before nucleic acid isolation. Subsequently, specific stem-looped reverse transcription primers were used to obtain cDNA for mature microRNAs. The RNA was reverse transcribed using the TaqMan® MicroRNA Reverse Transcription Kit from Applied Biosystems, The Netherlands (ABI). PCR was performed using pre-designed TaqMan® microRNA assays and TaqMan® Universal Master Mix, NoAmpErase®UNG, with an ABI 7900 HT real-time PCR machine. The PCR conditions were 2 min at 50 °C, 10 min at 95 °C, followed by 40 cycles of 15 s at 95 °C, and finally 1 min at 60 °C. The spiked-in syn-cel-miR-39 goes through the entire RNA isolation process and serves as endogenous control for data normalization. The Data Assist software (ABI) was used to normalize the data to syn-cel-miR-39, further calculations of FC were as described above.
Statistical analysis was performed using the SPSS (IBM, Inc.) package for Windows. Data were tested for normal distribution using the Kolmogorov–Smirnov test. The Grubbs’ test for outlier detection was applied (http://graphpad.com/support/faqid/1598/). Depending on the distribution pattern and the total number of subjects, parametric (normal distribution, independent t test) or nonparametric group comparison (Mann–Whitney U test) were applied. Correlations were determined by Spearman’s correlation coefficient. Levels of significance were set at p ≤ 0.05 (two tailed). A dendrogram visualizing associations was constructed in SPSS using hierarchical cluster analysis of the serum cytokines, genes and microRNA expression using the between-groups linkage method. Graphs were designed with Illustrator CS6 for Windows.
HGF is over expressed in monocytes of patients with T2D
Table 1 shows the demographic and clinical data of the Ecuadorian type 2 diabetic (T2D) patients and their controls. It is important to note that we were not completely successful in matching the controls to the patients with regard to age. We therefore corrected all further data for age. It is also worthy to note that controls were equally overweight and obese as our patients and had the same abnormal lipid profiles as the patients. Furthermore patients had a better liver profile than the non-diabetic controls. With regard to medication 31 % were not on anti-diabetic medication (recently discovered patients), while 45 % used metformin and 9 % used insulin (15 % were on both medications). None of the patients and non-diabetic controls used statins, Table 2A and Fig. 1a shows that the HGF expression levels were significantly higher in the monocytes of the patients with T2D as compared to the non-diabetic controls [(fold change T2D vs non-diabetic controls) 1.17 ± SEM 0.62, p = 0.03, n = 59]. Controlling for age, gender, dyslipidemia or liver function via hierarchical clustering showed that these factors did not contribute to the association of monocyte HGF expression with disease. With regard to BMI both BMI and T2D were significantly correlated to monocyte HGF expression. There was a non-significant difference in monocyte HGF expression between patients on metformin (FC 1.21 ± 0.37, n = 23) and un-medicated patients (FC 1.09 ± 0.30, n = 8), although group numbers are too small for a valid statistical evaluation.
The expression levels of the HGF-R (cMET, FC 1.34 ± 0.36 p = 0.40, n = 22) and of resistin (FC 0.47 ± 0.07, p = 0.24, n = 59) were not different between the groups (Table 2A).
Figure 2 shows the heat map and cluster diagram for HGF and resistin with the previously determined genes and the previously determined microRNAs. As can be seen HGF and resistin co-clustered positively with each other and with many genes of the cluster of adhesion/differentiation and shape change genes. Since HGF was significantly over expressed in the T2D monocytes, we focused on this compound. The association of HGF was significant at the p < 0.001 level with DHRS3 (r = .498, p = 0.004, n = 32), CD9 (r = .490, p = 0.004, n = 32), BCL2A1 (r = .503, p = 0.003, n = 32), Resistin (r = .532, p = 0.002, n = 32), HSPA1 (r = .525, p = 0.002, n = 32), but existed also at a lower level for MAPK6 (r = .385, p = 0.03, n = 32) and STX1A (r = .419, p = 0.024, n = 32). It is worthy to note that of these genes HGF, DHRS3 and CD9 were all three significantly higher expressed in the monocytes of cases with T2D as compared to the non-diabetic controls (see ).
It is also worthy to note that HGF expression did significantly negatively correlate with the expression of many genes of the inflammatory cluster in the monocytes (such as CCL4, IL-6, TNF, IL1-β, ATF3, CXCL2 and CCL20), reaching significance for TNFAIP3(r = −.350, p = 0.05, n = 32), supporting the concept that HGF is an anti-inflammatory agent. With regard to clinical parameters HGF expression in monocytes did positively correlate with the BMI (r = .327, p = 0.011, n = 59).
As reported previously, HGF was significantly raised in the serum of the patients with T2D as compared to the non-diabetic controls. Interestingly serum HGF correlated with monocyte DHRS3 gene expression (r = .326, p = 0.008, n = 64), but not with monocyte HGF gene expression (r = .189, p = 0.152, n = 59). Neither was there a positive or negative correlation for resistin expression in the monocytes and in the serum (data not shown), nor for any of the pro-inflammatory compounds, such as IL-1B, TNF, IL-6 and CCL2 (data not shown). This makes it unlikely that monocytes are the prime source of HGF, resistin or pro-inflammatory compounds in serum.
Mir-34c-5p is unaltered, but miR-574-3p is significantly reduced in the serum of patients with T2D
For the previous report on cytokines in T2D serum , we had determined the expression of various cytokines, growth factors, miR-146a and miR-155 in the serum of the cases with T2D and the non-diabetic controls. We found IL-8, HGF and resistin (the latter at a significance level of p = 0.07) raised in the serum of the patients with T2D in comparison to the non-diabetic controls, while miR-146a was down-regulated (see also “Background”).
For the current report, we determined the microRNAs miR-34c-5p, miR-122, miR-138, miR-410, miR-574-3p and miR-92 in the serum of the patients with T2D and the non-diabetic controls, since we had also measured these microRNAs in the monocytes of the patients in the previously reported study on gene and microRNA expression in the monocytes. In that study, we reported that miR-34c-5p was significantly up-regulated in the monocytes of the patients with T2D (see also “Background”).
Table 2B shows that in the current study the serum level of miR-34c-5p was not changed in the patients with T2D as compared to the non-diabetic controls. However, the serum level of microRNA miR-574-3p, was significantly reduced in the T2D serum as compared to the non-diabetic controls (see also Fig. 1b). Controlling for age, gender, BMI and dyslipidemia via hierarchical clustering showed that these factors did not contribute to the association of miR-574-3p with disease. With regard to medication there was a trend that particularly in the metformin-medicated group the miR-574 levels were reduced (FC 0.58 ± 0.35, n = 36 to non-diabetic controls), while in the unmediated group FC were only slightly reduced as compared to non-diabetic controls (FC 0.92 ± 0.77, n = 14).
Figure 3 shows the heat map and cluster diagram of the measured microRNAs and cytokines/growth factors. Since miR-574-3p was found significantly reduced in the serum of the patients with T2D, we focused in particular on this microRNA. It is clear from Fig. 3 that there is a strong clustering and association of miR-574-3p with miR-146a (r = 0.744, p < 0.001, n = 88) and miR-410 (r = 0.324, p = 0.03, n = 80), all microRNAs being decreased in the serum of the patients with T2D. The association of miR-574-3p with the other microRNAs and cytokines/growth factors in serum was not very strong, although there was an association with the serum CCL2 level (r = 0.337, p = 0.001, n = 89).
It is also important to note that correlations between the serum levels of the tested microRNAs and the expression levels of the same microRNAs in monocytes were not present (see Table 3).
Target prediction for miR-574-3p
Since the expression of miR-574-3p was significantly down-regulated in the serum of patients with T2D, we asked if there were in silico indications linking miR-574-3p expression to processes of inflammation or cell adhesion/differentiation/shape change. We used miRecords as a resource for miRNA-target interactions integrating predicted miRNA targets produced by 11 established miRNA target prediction programs (DIANA-microT, MicroInspector, miRanda, MirTarget2, miTarget, NBmiRTar, PicTar, PITA, RNA22, RNAhybrid and TargetScan/TargetScanS, available at http://www.mirecords.bioled.org).
A minimal target gene prediction coverage of three algorithms was used to perform prediction analysis for miR-574-3p. Filtering to a minimum coverage of three algorithms resulted in 934 hits. Ingenuity pathway analysis (Ingenuity® Systems) was used for mapping of the predicted target genes to biological functions.
Interestingly, the top molecular and cellular function of the miR-574-3p predicted target genes was “cell morphology” and “cellular assembly and organization”, while inflammation did not turn up in any of the predicted pathways (see Additional file 1: Ingenuity analysis).
This study showed that the gene expression of the vascular repair factor HGF was significantly raised in the monocytes of patients with T2D as compared to non-diabetic controls. HGF belonged to the cluster of adhesion, differentiation and shape change genes previously described as up-regulated in the T2D monocytes and correlated significantly to the expression of many genes in that cluster. The association of HGF with these differentiation, adhesion and shape change genes is in accordance with a view that the T2D monocytes are differentiating into the elongated vascular support pro-angiogenic cells (CACs) , HGF is a marker of such cells and suggests that the monocytes in our Ecuadorian patients with T2D are instrumental in repairing the vessel walls damaged by T2D related processes.
The up-regulation of HGF in the T2D monocytes might have been instrumental too in the previously reported anti-inflammatory state of the monocytes (see ). Indeed there is ample literature on the anti-inflammatory effects of HGF. It has been shown that monocytes treated with HGF produce high levels of IL-10, a potent immune suppressing cytokine. Mechanistically, HGF modulated IL-10 production in monocytes through the ERK1/2 pathway . With regard to dendritic cells (DC), Molnarfi et al. reported that DC differentiated in the presence of HGF adopt a pro-tolerogenic phenotype with increased ability to generate regulatory T cells , while with regard to endothelial cells Jeong-Ki Min et al. showed that HGF suppresses vascular endothelial growth factor (VEGF)–induced inflammation by inhibiting the nuclear factor kappa B (NFκB) pathway . In support of this anti-inflammatory action of HGF, we found a negative correlation of intra-monocyte HGF expression with the cluster of inflammatory genes, reaching statistical significance for TNFAIP3 (A20, a molecule induced by TNFα signaling, ) expression.
We reported previously that HGF was raised in the serum of the Ecuadorian patients with T2D; however, we did not find a correlation of the serum HGF with the expression of the HGF gene in the monocytes. This suggests that the circulating HGF is not primarily produced by the circulating monocytes, but originates from other sources. This notion also applies to the other tested cytokines and growth factors in the serum of the patients with T2D, for which we could also not find a correlation with intra-monocyte gene expression. Within the serum, the level of HGF correlated with the levels of TNF-α, IL-8 and resistin. Since it is generally thought that these pro-inflammatory compounds and insulin-resistance inducing substances originate from the adipose tissue and the liver [2, 4, 32, 33], we assume that also HGF in the serum of the patients with T2D primarily stems from these sources. There is ample literature on the production of HGF by adipose tissue and the liver [32, 34–37].
Although HGF might originate from sources other than the circulating monocytes, it is nevertheless possible that HGF in the serum could have affected the function of the circulating monocytes in patients with T2D, since the current study shows that monocytes of patients with T2D do express the HGF receptor, though not differently from monocytes of the non-diabetic controls. Interestingly, the level of HGF in serum did not correlate with a reduced inflammatory gene expression in the monocytes, as we found for the intra-monocyte expressed HGF. This suggests that the monocyte-endogenously-produced HGF is more important in the down-regulation of the inflammatory state of the monocytes in the patients with T2D than the serum-borne HGF. The level of circulating HGF did correlate positively to the DHRS3 expression in the monocytes, suggesting that serum-borne HGF might influence the proliferation and differentiation potential of circulating monocytes to pro-angiogenic cells in patients with T2D.
Furthermore, we found that the level of miR-574-3p was significantly reduced in the serum of patients with T2D, similar to miR-146a, of which we reported a down-regulation in the serum of the patients with T2D at an earlier occasion . In the cluster analysis and in correlation studies there was a strong association between the serum level of miR-574-3p and miR-146a, miR-410 and miR-155. This suggests an association of serum miR-574-3p (and also serum miR-410) with inflammatory processes, since miR-146a and miR-155 are important inflammation-regulating microRNAs [38–41]. This notion is further supported by a positive correlation between the serum level of miR-574-3p and the level of CCL2 in serum.
However, when we studied in silico the putative targets of miR-574-3p, ingenuity analysis of the putative targets did not indicate inflammation as an important pathway, whereas cell morphology and cellular assembly and organization were clearly present. The literature on miR-574-3p is in accord with this notion and shows functions of miR-574-3p mainly in the regulation of tumor cell pathology: MiR-574-3p is anti-proliferative, anti-invasive and anti-migratory in gastric and prostate cancer cells [42–44]. In these studies it was found that cullin-1 might be a target of miR-574-3p, and interestingly cullin-1 regulates inflammation via NFκB, thus giving an opening to a relationship with inflammation in connection with cellular assembly and organization .
Our findings suggest a high pro-angiogenic potential of the circulating monocytes and the serum of Ecuadorian patients with T2D. It is important to note that Yang reported that down-regulation of miR-574-3p in pro-angiogenic cells appeared to be a marker of senescence; senescent pro-angiogenic cells have lost their proliferative capacity and are changed in an inflammatory (oxidative radicals) direction . Hence, this might indicate that miR-574-3p is involved in the regulation of the anti-inflammatory and pro-angiogenic state of the circulating monocytes of patients with T2D.
Perhaps the most striking observation in the present study is the absence of correlation between the serum and intracellular monocyte levels of cytokines, growth factors and microRNAs. This suggests that the dynamics of the inflammation-related changes in the monocyte intracellular compartment differ substantially from the dynamics of inflammation-related changes in the serum compartment of patients with T2D. With regard to microRNAs it is possible that these dynamics involve complex processes of micro-vesicle or apoptotic body release from immune, endothelial or other cells  and/or binding to serum lipoproteins , illustrating the complexity of the biological activities of microRNAs in vascular and metabolic disease.
Collectively the data of the previous studies and the current study show that both the monocyte intracellular compartment and the serum compartment of our patients with T2D have undergone inflammation-related changes. However, the monocyte compartment shows in general a reduction in gene expression of typical pro-inflammatory genes, while genes and microRNAs involved in cell adhesion, cell differentiation, growth and vascular repair, such as HGF and miR-34c-5p, are up-regulated. The serum compartment, in contrast to the monocyte compartment, does show signs of high pro-inflammatory activity, e.g. high levels of IL-8 and reduced levels of anti-inflammatory miR-146a and altered levels of miR-574-3p. The serum compartment also shows signs of higher activity of vascular repair and cellular growth induction, yet parameters do not correlate with the monocyte parameters of higher pro-angiogenic cell activity. Most likely different T2D related pathophysiological forces drive the activation and de-activation set points of the circulating monocyte and the serum compartment.
Due to the paucity of material we have not been able to carry out all assays in all patients and thus for some microRNAs (such as e.g. 34c-5p and 576-3p) test numbers might have been too small to detect significant changes. The likelihood that we would have found a significant difference for these microRNAs if we had been able to increase the sample size is possible. However when numbers for miR-574-3p, which we found statistically significant between the groups with n numbers of 53 patients and 37 controls, were recalculated with reduced numbers equal to those for the miRs 34c-5p and 576-3p (with n numbers around half of those for 574-3p), near significant p levels between 0.05 and 0.10 were found. With none of the here tested microRNAs such near significant p levels were found. Nevertheless it is clear that in future studies larger samples need to be tested to come to more solid conclusions.
Although we selected patients and controls on the basis of their diabetic state (glucose levels) to study diabetes as the major determinant for immune differences, we realize that other determinants might have played equally important roles in the outcomes of our study, such as lipid state, liver function, BMI and medication. Correction for lipid state, BMI and liver function did not change our main findings regarding the higher expression of monocyte HGF and reduced serum miR-574-3p. With regard to metformin treatment, which has been shown to alter the immune state , post hoc analyses of our data found non-significant differences between metformin treated and not treated patients, in particular that the reduced serum miR-574-3p values were more clear in the metformin treated patients. However numbers in the test groups were too small for a valid statistical evaluation.
Although focus was on HGF in the here reported studies to evaluate the vessel repair quality of the circulating monocytes, other important analytes could have been studied as well, such as adiponectin and leptin. Levels of these adipokines are important in determining atherosclerosis  and should therefore be studied in future investigations on this subject. This is the more relevant since we found in preliminary unpublished studies that the adiponectin and leptin serum levels were different between Ecuadorian and Dutch T2D patients. Similar observations on an altered leptin/adiponectin ratio have been done by Bribiescas et al. in studying Ache Amerindians versus US individuals .
Despite the limitations of the here described study we conclude that in circulating monocytes of T2D Ecuadorian patients, the microRNA and gene expression of important inflammatory factors, chemotactic/motility factors and a vascular repair factor differs from the expression in serum. While monocytes show a gene expression profile compatible with an anti-inflammatory state, the serum shows a molecular profile of an inflammatory state, suggesting an intricate feedback network. Both compartments show molecular signs of vascular repair support, i.e. up-regulated HGF levels.
tumor necrosis factor a-induced protein (TNFAIP) 3
Abelson murine leukemia viral oncogene homolog 1
cyclic AMP-dependent transcription factor ATF-3
Bcl-2-related protein A1
circulating pro-angiogenic cells
chemokine (C–C motif) ligand 2. Monocyte chemotactic protein 1 (MCP1)
chemokine (C–C motif) ligand 20
chemokine (C–C motif) ligand 4
chemokine (C–C motif) ligand 7
cell division 9
cell division control protein 42 homolog
complementary Deoxyribonucleic acid
hepatocyte growth factor receptor
chemokine (C–X–C motif) ligand 2
short-chain dehydrogenase/reductase 3
dual specificity protein phosphatase 2
epithelial membrane protein 1
fatty acid-binding protein, epidermal
glutamic acid decarboxylase 65
hepatocyte growth factor
hepatocyte growth factor receptor
heat shock 70 kDa protein 1A
heat shock 70 kDa protein 1B
interleukin 1 beta
latent autoimmune diabetes in adults
mitogen-activated protein kinase 6
monocyte chemotactic protein 1 (CCL2)
macrophage inflammatory protein 1 beta (CCL4)
messenger ribonucleic acid
NGFI-A-binding protein 2
nuclear factor kappa-light-chain-enhancer of activated B cells
peripheral blood mononuclear cells
polymerase chain reaction
cAMP-specific 3′,5′-cyclic phosphodiesterase 4B
prostaglandin-endoperoxide synthase 2
protein tyrosine phosphatase non-receptor type 7
pentraxin-related protein 3
quantitative polymerase chain reaction
real-time polymerase chain reaction
tumor necrosis factors
tumor necrosis factor, alpha-induced protein 3
Kamei N, Tobe K, Suzuki R, Ohsugi M, Watanabe T, Kubota N, Ohtsuka-Kowatari N, Kumagai K, Sakamoto K, Kobayashi M, et al. Overexpression of monocyte chemoattractant protein-1 in adipose tissues causes macrophage recruitment and insulin resistance. J Biol Chem. 2006;281(36):26602–14.
Wellen KE, Hotamisligil GS. Obesity-induced inflammatory changes in adipose tissue. J Clin Investig. 2003;112(12):1785–8.
Sierra-Honigmann MR, Nath AK, Murakami C, Garcia-Cardena G, Papapetropoulos A, Sessa WC, Madge LA, Schechner JS, Schwabb MB, Polverini PJ, et al. Biological action of leptin as an angiogenic factor. Science. 1998;281(5383):1683–6.
Hotamisligil GS, Shargill NS, Spiegelman BM. Adipose expression of tumor necrosis factor-alpha: direct role in obesity-linked insulin resistance. Science. 1993;259(5091):87–91.
Hotamisligil GS. Inflammation and endoplasmic reticulum stress in obesity and diabetes. Int J Obes (Lond). 2008;32(Suppl 7):S52–4.
Wajchenberg BL, Nery M, Cunha MR, Silva ME. Adipose tissue at the crossroads in the development of the metabolic syndrome, inflammation and atherosclerosis. Arq Bras de Endocrinol e Metabol. 2009;53(2):145–50.
Heilbronn LK, Campbell LV. Adipose tissue macrophages, low grade inflammation and insulin resistance in human obesity. Curr Pharm Design. 2008;14(12):1225–30.
Swirski FK. The spatial and developmental relationships in the macrophage family. Arterioscler Thromb Vasc Biol. 2011;31(7):1517–22.
Sorisky A, Molgat AS, Gagnon A. Macrophage-induced adipose tissue dysfunction and the preadipocyte: should I stay (and differentiate) or should I go? Adv Nutr. 2013;4(1):67–75.
Dalmas E, Clement K, Guerre-Millo M. Defining macrophage phenotype and function in adipose tissue. Trends Immunol. 2011;32(7):307–14.
Simar D, Jacques A, Caillaud C. Heat shock proteins induction reduces stress kinases activation, potentially improving insulin signalling in monocytes from obese subjects. Cell Stress Chaperones. 2012;17(5):615–21.
Shim WS, Kim HJ, Kang ES, Ahn CW, Lim SK, Lee HC, Cha BS. The association of total and differential white blood cell count with metabolic syndrome in type 2 diabetic patients. Diabetes Res Clin Prac. 2006;73(3):284–91.
Gacka M, Dobosz T, Szymaniec S, Bednarska-Chabowska D, Adamiec R, Sadakierska-Chudy A. Proinflammatory and atherogenic activity of monocytes in type 2 diabetes. J Diabetes Complic. 2010;24(1):1–8.
Dasu MR, Devaraj S, Park S, Jialal I. Increased toll-like receptor (TLR) activation and TLR ligands in recently diagnosed type 2 diabetic subjects. Diabetes Care. 2010;33(4):861–8.
Shiny A, Regin B, Balachandar V, Gokulakrishnan K, Mohan V, Babu S, Balasubramanyam M. Convergence of innate immunity and insulin resistance as evidenced by increased nucleotide oligomerization domain (NOD) expression and signaling in monocytes from patients with type 2 diabetes. Cytokine. 2013;64(2):564–70.
Padmos RC, Schloot NC, Beyan H, Ruwhof C, Staal FJ, de Ridder D, Aanstoot HJ, Lam-Tse WK, de Wit H, de Herder C, et al. Distinct monocyte gene-expression profiles in autoimmune diabetes. Diabetes. 2008;57(10):2768–73.
Baldeón RL, Weigelt K, de Wit H, Ozcan B, van Oudenaren A, Sempértegui F, Sijbrands E, et al. Type 2 diabetes monocyte microRNA and mRNA expression: dyslipidemia associates with increased differentiation-related genes but not inflammatory activation. Plos One. 2015;10(6):e0129421. doi:10.1371/journal.pone.0129421.
Baldeon RL, Weigelt K, de Wit H, Ozcan B, van Oudenaren A, Sempertegui F, Sijbrands E, Grosse L, Freire W, Drexhage HA, et al. Decreased serum level of miR-146a as sign of chronic inflammation in type 2 diabetic patients. PloS One. 2014;9(12):e115209.
Jennische E, Ekberg S, Matejka GL. Expression of hepatocyte growth factor in growing and regenerating rat skeletal muscle. Am J Physiol. 1993;265(1 Pt 1):C122–8.
Ono K, Matsumori A, Shioi T, Furukawa Y, Sasayama S. Enhanced expression of hepatocyte growth factor/c-Met by myocardial ischemia and reperfusion in a rat model. Circulation. 1997;95(11):2552–8.
Aoki M, Morishita R, Taniyama Y, Kida I, Moriguchi A, Matsumoto K, Nakamura T, Kaneda Y, Higaki J, Ogihara T. Angiogenesis induced by hepatocyte growth factor in non-infarcted myocardium and infarcted myocardium: up-regulation of essential transcription factor for angiogenesis, ets. Gene Ther. 2000;7(5):417–27.
Taniyama Y, Morishita R, Aoki M, Nakagami H, Yamamoto K, Yamazaki K, Matsumoto K, Nakamura T, Kaneda Y, Ogihara T. Therapeutic angiogenesis induced by human hepatocyte growth factor gene in rat and rabbit hindlimb ischemia models: preclinical study for treatment of peripheral arterial disease. Gene Ther. 2001;8(3):181–9.
Chen PM, Liu KJ, Hsu PJ, Wei CF, Bai CH, Ho LJ, Sytwu HK, Yen BL. Induction of immunomodulatory monocytes by human mesenchymal stem cell-derived hepatocyte growth factor through ERK1/2. J Leukoc Biol. 2014;96(2):295–303.
Molnarfi N, Benkhoucha M, Juillard C, Bjarnadottir K, Lalive PH. The neurotrophic hepatocyte growth factor induces protolerogenic human dendritic cells. J Neuroimmunol. 2014;267(1–2):105–10.
Bouchentouf M, Paradis P, Forner KA, Cuerquis J, Boivin MN, Zheng J, Boulassel MR, Routy JP, Schiffrin EL, Galipeau J. Monocyte derivatives promote angiogenesis and myocyte survival in a model of myocardial infarction. Cell Transplant. 2010;19(4):369–86.
American Diabetes Association. Diagnosis and classification of diabetes mellitus. Diabetes Care. 2011;34(Suppl 1):S62–9. doi:10.2337/dc11-S062.
Menke A, Casagrande S, Geiss L, Cowie CC. Prevalence of and trends in diabetes among adults in the United States, 1988–2012. JAMA. 2015;314(10):1021–9.
Knijff EM, Breunis MN, van Geest MC, Kupka RW, Ruwhof C, de Wit HJ, Nolen WA, Drexhage HA. A relative resistance of T cells to dexamethasone in bipolar disorder. Bipolar Disord. 2006;8(6):740–50.
Lyons PA, Koukoulaki M, Hatton A, Doggett K, Woffendin HB, Chaudhry AN, Smith KG. Microarray analysis of human leucocyte subsets: the advantages of positive selection and rapid purification. BMC Genom. 2007;8:64.
Min JK, Lee YM, Kim JH, Kim YM, Kim SW, Lee SY, Gho YS, Oh GT, Kwon YG. Hepatocyte growth factor suppresses vascular endothelial growth factor-induced expression of endothelial ICAM-1 and VCAM-1 by inhibiting the nuclear factor-kappaB pathway. Circ Res. 2005;96(3):300–7.
Catrysse L, Vereecke L, Beyaert R, van Loo G. A20 in inflammation and autoimmunity. Trends Immunol. 2014;35(1):22–31.
Bell LN, Ward JL, Degawa-Yamauchi M, Bovenkerk JE, Jones R, Cacucci BM, Gupta CE, Sheridan C, Sheridan K, Shankar SS, et al. Adipose tissue production of hepatocyte growth factor contributes to elevated serum HGF in obesity. Am J Physiol Endocrinol Metab. 2006;291(4):E843–8.
Kammoun HL, Kraakman MJ, Febbraio MA. Adipose tissue inflammation in glucose metabolism. Rev Endocr Metab Disord. 2014;15(1):31–44.
Fain JN. Release of interleukins and other inflammatory cytokines by human adipose tissue is enhanced in obesity and primarily due to the nonfat cells. Vitam Horm. 2006;74:443–77.
Christiaens V, Lijnen HR. Angiogenesis and development of adipose tissue. Mol Cell Endocrinol. 2010;318(1–2):2–9.
Cao Y. Angiogenesis modulates adipogenesis and obesity. J Clin Investig. 2007;117(9):2362–8.
Schirmacher P, Geerts A, Pietrangelo A, Dienes HP, Rogler CE. Hepatocyte growth factor/hepatopoietin A is expressed in fat-storing cells from rat liver but not myofibroblast-like cells derived from fat-storing cells. Hepatology. 1992;15(1):5–11.
Balasubramanyam M, Aravind S, Gokulakrishnan K, Prabu P, Sathishkumar C, Ranjani H, Mohan V. Impaired miR-146a expression links subclinical inflammation and insulin resistance in Type 2 diabetes. Mol Cell Biochem. 2011;351(1–2):197–205.
Taganov KD, Boldin MP, Chang KJ, Baltimore D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci USA. 2006;103(33):12481–6.
Cui JG, Li YY, Zhao Y, Bhattacharjee S, Lukiw WJ. Differential regulation of interleukin-1 receptor-associated kinase-1 (IRAK-1) and IRAK-2 by microRNA-146a and NF-kappaB in stressed human astroglial cells and in Alzheimer disease. J Biol Chem. 2010;285(50):38951–60.
O’Connell RM, Taganov KD, Boldin MP, Cheng G, Baltimore D. MicroRNA-155 is induced during the macrophage inflammatory response. Proc Natl Acad Sci USA. 2007;104(5):1604–9.
Su Y, Ni Z, Wang G, Cui J, Wei C, Wang J, Yang Q, Xu Y, Li F. Aberrant expression of microRNAs in gastric cancer and biological significance of miR-574-3p. Int Immunopharmacol. 2012;13(4):468–75.
Chiyomaru T, Yamamura S, Fukuhara S, Hidaka H, Majid S, Saini S, Arora S, Deng G, Shahryari V, Chang I, et al. Genistein up-regulates tumor suppressor microRNA-574-3p in prostate cancer. PloS One. 2013;8(3):e58929.
Singleton KD, Wischmeyer PE. Glutamine attenuates inflammation and NF-kappaB activation via Cullin-1 deneddylation. Biochem Biophys Res Commun. 2008;373(3):445–9.
Yang DG, Liu L, Zhou SH. MicroRNA alterations in senescent endothelial progenitor cells induced by remnant-like lipoproteins. Chin Med Journal. 2012;125(19):3479–84.
Zampetaki A, Mayr M. MicroRNAs in vascular and metabolic disease. Circ Res. 2012;110(3):508–22.
Norata GD, Sala F, Catapano AL, Fernandez-Hernando C. MicroRNAs and lipoproteins: a connection beyond atherosclerosis? Atherosclerosis. 2013;227(2):209–15.
Lee HM, Kim JJ, Kim HJ, Shong M, Ku BJ, Jo EK. Upregulated NLRP3 inflammasome activation in patients with type 2 diabetes. Diabetes. 2013;62(1):194–204.
Norata GD, Ongari M, Garlaschelli K, Raselli S, Grigore L, Catapano AL. Plasma resistin levels correlate with determinants of the metabolic syndrome. Eur J Endocrinol Eur Fed Endocr Soc. 2007;156(2):279–84.
Bribiescas RG, Hickey MS. Population variation and differences in serum leptin independent of adiposity: a comparison of Ache Amerindian men of Paraguay and lean American male distance runners. Nutr Metab. 2006;3:34.
LBR: study design, data acquisition, carried out the experimental studies, analysis and interpretation of results. Manuscript preparation, and critically evaluate the manuscript. KW: study design, analysis and interpretation of results, drafting of the manuscript and critically evaluated the manuscript. HdW: study design, data acquisition, analysis and interpretation of results, critically evaluated the manuscript, and carried out the experimental studies. BO: contributed to the analysis and interpretation of results, manuscript preparation, and drafting of the manuscript and critically evaluated the manuscript. AvO: data acquisition, analysis and interpretation of results, critically evaluated the manuscript, and carried out the experimental studies. FS: helped in the design of the studies, analysis and interpretation of results, manuscript preparation, and drafting of the manuscript and critically evaluated the manuscript. ES: study design, analysis and interpretation of results, manuscript preparation, and drafting of the manuscript and critically evaluated the manuscript. LG: study design, data acquisition, analysis and interpretation of results, critically evaluated the manuscript, and carried out the experimental studies. AJvZ: study design, manuscript preparation, and drafting of the manuscript and critically evaluated the manuscript. HD: study design, analysis and interpretation of results, manuscript preparation, drafting of the manuscript and critically evaluated the manuscript, supervised the experimental studies. PL: study design, analysis and interpretation of results, manuscript preparation, drafting of the manuscript and critically evaluated the manuscript, supervised the experimental studies. All authors read and approved the final manuscript.
The authors would like to acknowledge Cesar Prócel MD, Internist at Clínica Pasteur, Quito, Ecuador; MDD. Jonhson Viteri, Club de Leones La Villaflora, Quito, Ecuador; MD. Sandra Montalvo, Fundación Oftalmológica del Valle; Corporación Ecuatoriana de Biotecnología, and Fundación Ecuatoriana de la Psoriasis for patients recruitment and the logistic support. This study was supported by Grant 2007.00.035 from the Dutch Diabetes Research Foundation, The Netherlands.
The authors declare that they have no competing interests.
Additional file 1. Ingenuity pathway analysis (Ingenuity® Systems) was used to map the major pathways and processes in which miR-574-3p is involved. The top molecular and cellular function of the miRNA predicted target genes was “cell morphology”; while the second top-associated network was “cell morphology” and “cellular assembly and organization”. We used miRecords as a resource for microRNA-target interactions.
About this article
Cite this article
Baldeón Rojas, L., Weigelt, K., de Wit, H. et al. Study on inflammation-related genes and microRNAs, with special emphasis on the vascular repair factor HGF and miR-574-3p, in monocytes and serum of patients with T2D. Diabetol Metab Syndr 8, 6 (2016). https://doi.org/10.1186/s13098-015-0113-5