Next Article in Journal
Personalised Nutritional Recommendations Based on Individual Post-Prandial Glycaemic Responses Improve Glycaemic Metrics and PROMs in Patients with Type 2 Diabetes: A Real-World Assessment
Next Article in Special Issue
Maternal Dietary Inflammatory Index during Pregnancy Is Associated with Perinatal Outcomes: Results from the IMPACT BCN Trial
Previous Article in Journal
Anticholinergic Load Is Associated with Swallowing Dysfunction in Convalescent Older Patients after a Stroke
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epigenetic Profiling in the Saliva of Obese Pregnant Women

1
Department of Biomedical and Clinical Sciences, Università degli Studi di Milano, 20157 Milan, Italy
2
Department of Dentistry, University Vita-Salute San Raffaele, 20132 Milan, Italy
3
EPIGET LAB, Department of Clinical Sciences and Community Health, Università degli Studi di Milano, 20122 Milan, Italy
4
Department of Pathophysiology and Transplantation, Università degli Studi di Milano, 20054 Segrate, Italy
5
Department of Woman, Mother and Child, Luigi Sacco and Vittore Buzzi Children Hospital, ASST Fatebenefratelli-Sacco, 20154 Milan, Italy
6
Occupational Health Unit, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico, 20122 Milan, Italy
*
Author to whom correspondence should be addressed.
Submission received: 25 April 2022 / Revised: 12 May 2022 / Accepted: 16 May 2022 / Published: 19 May 2022
(This article belongs to the Special Issue The Interplay between Maternal Nutrition and Oxidative Stress)

Abstract

:
Maternal obesity is associated with inflammation and oxidative stress, strongly impacting the intrauterine environment with detrimental consequences for both mother and offspring. The saliva is a non-invasive biofluid reflecting both local and systemic health status. This observational study aimed to profile the epigenetic signature in the saliva of Obese (OB) and Normal-Weight (NW) pregnant women. Sixteen NW and sixteen OB Caucasian women with singleton spontaneous pregnancies were enrolled. microRNAs were quantified by the OpenArray Platform. The promoter region methylation of Suppressor of Cytokine Signaling 3 (SOCS3) and Transforming Growth Factor Beta 1 (TGF-Beta1) was assessed by pyrosequencing. There were 754 microRNAs evaluated: 20 microRNAs resulted in being differentially expressed between OB and NW. microRNA pathway enrichment analysis showed a significant association with the TGF-Beta signaling pathway (miTALOS) and with fatty acids biosynthesis/metabolism, lysine degradation, and ECM–receptor interaction pathways (DIANA–miRPath). Both SOCS3 and TGF-Beta1 were significantly down-methylated in OB vs. NW. These results help to clarify impaired mechanisms involved in obesity and pave the way for the understanding of specific damaged pathways. The characterization of the epigenetic profile in saliva of pregnant women can represent a promising tool for the identification of obesity-related altered mechanisms and of possible biomarkers for early diagnosis and treatment of pregnancy-adverse conditions.

Graphical Abstract

1. Introduction

The obesity epidemic is expected to involve almost 21% of women all over the world by 2025 [1], representing a serious risk for women of reproductive age. Maternal obesity during pregnancy represents a danger for short- and long-term health outcomes for both mothers and children [2,3,4]. Indeed, up to 50% of obese pregnant women develop Gestational Diabetes Mellitus (GDM) [5], and the offspring of obese and diabetic mothers are at higher risk of developing metabolic and cardiovascular diseases in early childhood and later in their adulthood [6], with possible transgenerational effects [7]. Moreover, newborns of obese and diabetic mothers have been recently shown to be hypoxic, acidemic, and to have increased oxidative markers compared to normal-weight pregnancies [8,9].
The metabolic and hormonal balance is dysregulated in obesity and diabetes because of excessive adipose tissue deposition and intracellular fatty acids, leading to alterations of the immune cell profile and inflammation [10]. This metabolic-dependent pro-inflammatory background increases oxidative stress and insulin resistance, which further contribute to metabolic dysregulation, in a deleterious vicious circle. Pregnancy itself is characterized by a physiological generalized low-grade inflammation, and maternal obesity can contribute to additional lipotoxicity, inflammation, and oxidative stress, which may strongly impact the intrauterine environment, with significant effects on fetal health and metabolism [8,11,12]. We previously reported alterations in the placentas of obese mothers, showing oxidative stress and a dysregulated placental metabolome [13,14,15]. These data suggest that obesity-related increased inflammation and oxidative stress may determine a cascade of events, leading to systemic and local dysfunctions and consequent impaired pregnancy outcomes.
Among obesity-related recurring comorbidities, periodontal diseases represent an additional important risk factor for pregnancy outcomes [16]. During pregnancy, both obesity and diabetes can cause oral dysbiosis [17], which is the primary cause of periodontal diseases [18]. In turn, maternal periodontal diseases have been associated with poor obstetric outcomes, such as pre-term birth, premature rupture of membranes, and low birth weight, especially in obese women with pre-gestational obesity and GDM. The systemic inflammatory alterations induced by both obesity/GDM and periodontal diseases can indeed activate adverse responses at the maternal–fetal interface [19,20,21,22]. We recently reported increased total antioxidant capacity and inflammatory levels in both the saliva and blood of obese pregnant women with and without GDM, suggesting a synergic detrimental effect of obesity and periodontal diseases [23].
The saliva is an advantageous, non-invasive, and cost-effective biofluid, reflecting both the local (oral) and systemic health status of the human body [24]. Salivary glands are surrounded by capillaries, resulting in a continuous exchange of molecules between blood flow and saliva secretion from the acinus cells. Therefore, the circulating molecules absorbed by salivary glands can affect the saliva composition, mirroring the systemic status [25]. Among the multitude of molecules contained in the saliva, microRNAs (miRNAs) are proving promising biomarkers, representing important epigenetic regulators. Indeed, about 500 of these short non-coding RNA molecules have been identified in saliva [26]. Interestingly, alterations of several miRNAs’ levels have been reported in the maternal circulation depending on pregestational BMI [27], and a few salivary miRNAs have been identified as predictors of periodontal disease in both non-diabetic and diabetic non-pregnant patients [28]. However, to our knowledge, miRNAs’ profiles in the saliva have never been investigated during pregnancy in the presence of obesity and possible related co-morbidities.
This case–control observational study aimed to assess salivary miRNAs levels in a group of obese pregnant women that were previously deeply characterized for periodontal diseases [23].
Moreover, further epigenetic modifications such as gene methylation are also known to be deregulated in obesity and obesity-related insulin resistance. Therefore, methylation levels of specific genes that are well-known players in this context were assessed in maternal saliva: Suppressor Of Cytokine Signalling-3 (SOCS3), which is upregulated in response to obesity-related stimuli [29], and Transforming Growth Factor-Beta1 (TGF-Beta1), which is a key cytokine in obesity and insulin resistance and has been reported at higher levels in women with a prior history of GDM [30].

2. Materials and Methods

2.1. Study Participants

This was a pilot study performed on a cohort of pregnant women enrolled in the antenatal clinic at the regular checkup in the Obstetric Unit of the L. Sacco Hospital (ASST Fatebenefratelli-Sacco) in Milan. The study was conducted in accordance with the Declaration of Helsinki and in compliance with current Good Clinical Practice guidelines, local laws, regulations, and organizations. The protocol was approved by the hospital ethical committee (Prot. N. 469/2010/52/AP). Participants gave their written informed consent to collect personal data and biological samples. Pregnant women were enrolled during the first trimester and had a regular clinical follow-up during all trimesters. A deep oral health characterization of this population has been previously reported [23].
Only Caucasian women with a singleton spontaneous pregnancy and aged between 18 and 40 years were enrolled. Exclusion criteria were maternal and fetal infections, fetal malformations, chromosomal disorders, maternal alcohol/drugs abuse, and pregestational Body Mass Index (BMI) < 18.5 or > 25. Detailed inclusion and exclusion criteria are reported in Supplementary Table S1.
Thirty-two pregnant women were classified according to their pregestational BMI as Normal-Weight (NW, n = 16: 18.5 ≤ BMI ≤ 24.9 kg/m2) or Obese (OB, n = 16: BMI ≥ 30 kg/m2). Ten out of sixteen OB women were diagnosed with Gestational Diabetes Mellitus (GDM), based on oral glucose tolerance test (75 g), according to our clinical protocol. NW women had uncomplicated pregnancies. Therefore, GDM was an exclusion criterion in this group.
OB patients were followed up in a dedicated antenatal clinic, providing nutritional counseling on diet and lifestyle, particularly focused on pregnancy weight gain. Obese patients with GDM received further advice to control their glycemia according to the clinical protocol.

2.2. Clinical Data and Biological Samples Collection

Clinical data: Maternal, neonatal, and placental data were recorded at recruitment and after delivery.
Oral and periodontal health: A complete clinical oral health examination was performed in the third trimester, as previously reported [23,31]. Pregnant patients were classified as healthy with an overall good oral care (no dental soft debris—dental plaque or materia alba or calculus—or nor gingival Bleeding On Probing (BOP) or Probing Pocket Depth (PPD > 3 mm), or periodontally diseased with an overall neglected oral care, which grouped gingivitis (soft debris and calculus in ≥ 6 teeth and gingival bleeding ≤ 5 teeth) and periodontitis (probing pocket depth ≥ 4 mm, soft debris, calculus, and bleeding ≥ 6 teeth).
Saliva collection: Saliva sampling was performed during the third trimester without any stimulation (passive drool technique), thus avoiding the use of solutions that could alter physico-chemical features and preventing the employment of tools requiring a specific site of withdrawal. Choosing the same timing of withdrawal for all subjects eludes any fluctuations in the saliva composition related to circadian rhythms. Enrolled women were asked to refrain from eating, drinking sugar or alcoholic beverages, smoking, and performing invasive oral care procedures for at least 1h before collection. They were then asked to rinse their mouth (1 min) with a physiological solution to remove any food residues, which could alter pH or promote bacterial growth, and afterward, not to swallow for 3 min and spit in a sterile tube. Collected samples were temporarily stored at −80 °C.

2.3. Isolation of Extracellular Vesicles and miRNA

Saliva samples (1 mL) were thawed on ice to avoid Extracellular Vesicles’ (EVs) thermal damage. Then, they were centrifuged at 4000× g for 30 min at 4 °C to remove any cell debris and aggregates. Supernatants were ultracentrifuged at 110,000× g for 75 min at 4 °C in order to pellet EVs, then stored at −20 °C.
miRNAs’ isolation from the obtained EVs was performed with the combination of the miRNeasy Kit and RNeasy Cleanup Kit (Qiagen, Hilden, Germany), according to the manufacturer’s protocol. They were eluted in 20 µL of nuclease-free water and stored at −80 °C until used. EV miRNAs’ quality and integrity were assessed through the “2100 Bioanalyzer RNA system” with the Pico Kit (Agilent Technologies, Santa Clara, CA, USA), but the concentration (ng/µL) was assessed by a Quantus Fluorometer (Promega, Milan, Italy).

2.4. Screening of miRNA Expression

There were 15 ng miRNAs Reverse Transcribed (RT), preamplified (16 cycles), and analyzed by real-time PCR with the QuantStudio™ 12K Flex OpenArray® Platform (Applied Biosystems, Waltham, MA, USA), as previously described in detail [32]. Gene Expression Suite Software v.1.0.3 (Applied Biosystems, Waltham, MA, USA) was used to process miRNA profiling expression data from the “TaqMan™ OpenArray™ human MicroRNA panel” (ThermoFisher Scientific, Waltham, MA, USA). As described in Figure 1, we obtained 758 Relative Threshold Cycle (Crt) values for each subject, which included 754 unique miRNAs and four internal controls (ath-miR159a, RNU48, RNU44, and U6). For each amplification curve, we obtained an AmpScore value, a quality measurement that indicates the low signal in the amplification curve linear phase (range: 0–2). miRNAs with Crt value > 28, or AmpScore < 1.24, or missing were considered not amplified, and thus, the corresponding Crt value was set to 29. miRNAs that were not amplified in all subjects (n = 420) were excluded, resulting in 334 miRNAs being included in the statistical analysis. miRNA expression was determined using the relative quantification 2-ΔΔCrt [33]. Among the endogenous controls analyzed (i.e., ath-miR159a, U6-rRNA, RNU44, RNU48), U6-rRNA and RNU48 were chosen for normalization by the NormFinder algorithm [34], due to their stability in the comparison between samples.

2.5. DNA Isolation

Total DNA was extracted using the QIAamp® DNA Mini and Blood Mini Kit (Qiagen; Valencia, CA, USA) following the supplementary protocol (“Isolation of genomic DNA from saliva and mouthwash using the QIAamp® DNA Blood Mini Kit”). Briefly, saliva (160 μL) was thawed on ice and added with 640 μL of Phosphate Buffered Saline (PBS). Samples were then centrifuged at 1800× g for 10 min and their supernatant was discarded. The pellets were lysed and DNA extracted with a spin-column-based method. DNA was eluted in 44 μL of Buffer AE (10 mM Tris·Cl; 0.5 mM EDTA; pH 9.0) for a final concentration of 25 ng/μL.

2.6. Bisulfite Conversion and DNA Methylation

500 ng di genomic DNA was bisulfite converted using the EZ DNA Methylation Direct Kit (Zymo Research, Orange, CA, USA), in accordance with the manufacturer’s protocol. Converted DNA was eluted in 30 μl of Elution Buffer.
Analysis of DNA methylation was performed by PyroMark MD Pyrosequencing System (Qiagen, Milan, Italy) as previously described [35,36].
The bisulfite-treated genomic DNA samples were amplified with PCR; CpG sites were analyzed within the promoter regions of the following genes: Suppressor of Cytokine Signalling 3 (SOCS3) and Transforming Growth Factor Beta 1 (TGF-Beta1).
Detailed information concerning primer sequences and genomic regions is listed in Table 1. The percentage of 5- methylcytosine (% 5mC) was reported as the percentage of methylated cytosine divided by the sum of methylated and unmethylated cytosines. Every sample was tested twice for each marker to confirm reproducibility and to increase the precision of the findings.

2.7. Statistical Analysis and Prediction Tools

Maternal, neonatal, and placental characteristics and periodontal health parameters were compared between groups by the Student t-test or Mann–Whitney U test for independence samples according to the data distribution (assessed by the Kolmogorov–Smirnov test). Differences between groups in the frequencies of periodontal disease, mode of delivery, neonatal weight centile, and fetal sex were evaluated by the chi-squared test for independence or Fisher’s exact test.
Multivariable linear regression models were applied to verify the association between BMI groups and miRNA expression in the screening phases. miRNA expression values were log2-transformed to achieve a normal distribution. Regression models were adjusted for age, smoking habit, GDM, and gestational age at saliva withdrawal. Due to the high number of comparisons, we applied a multiple comparison correction method based on the Benjamini–Hochberg False Discovery Rate (FDR) to calculate the FDR p-value. The criteria used to identify the top miRNAs were a p-value < 0.05, an FDR p-value < 0.20, and a Fold Change (FC) < 0.5 or FC > 2. A Volcano plot was produced to select miRNAs characterized by more than ±2-fold case–control differences (FC 2 or < 0.5) with a p-value < 0.05 from linear regression models.
To evaluate the associations between methylation levels (SOCS3, TGF-Beta1) and BMI groups, we used multiple linear regression models adjusted for maternal age, smoking habits, and GDM. We report marginal means and 95% CI.
Statistical analyses were performed using the statistical package SPSS, v.27 (IBM, Armonk, NY, USA) and SAS 9.4 statistical software (SAS Institute Inc., Cary, NC, USA).
miRNA pathway enrichment analysis: miTALOS v2 tools (http://mips.helmholtz-muenchen.de/mitalos (accessed on 15 March 2022)) were used to gain insights into tissue-specific miRNA regulation of biological pathways. Indeed, miTALOS integrates five different miRNA target prediction tools and two different resources (KEGG and NCI), being able to identify a specific association between miRNAs and signaling pathways [37].
We integrated our analysis by using DIANA–miRPath v3 (https://dianalab.e-ce.uth.gr/html/mirpathv3/index.php?r=mirpath (accessed on 15 March 2022)), which uses data from the KEGG and TarBase resources, to perform miRNA pathway enrichment analysis. We thus investigated how these pathways are combined in a posteriori bioinformatic analysis (union and intersection of pathways search). The results are shown as miRNA/pathways’ interaction heat maps [38].
For both miTALOS v2 and DIANA–miRPath v3, the significance of the association (p-value) between miRNAs and the resulting pathways was calculated with Fisher’s exact test; then, the results for multiple pathways were corrected using the EASE score and false discovery rate (Benjamini–Hochberg) [39,40,41].

3. Results

3.1. Maternal and Oral Health Characteristics and Delivery Data

Table 2 resumes the maternal and oral health characteristics.
Pregestational BMI was significantly higher in the OB group compared to NW (p < 0.001), according to the inclusion criteria. Gestational Weight Gain (GWG) was significantly lower in OB vs. NW (p = 0.035), in compliance with IOM indications (upper limits: +9 kg for OB, +16 kg for NW). However, OB women gained more overall excessive weight than recommended for OB, as shown by the percentages reported in Table 2 (maternal GWG to IOM upper limit). Maternal basal glycemia was higher in the OB group (p = 0.034).
Unstimulated saliva was collected at a similar gestational age, and its flow rate had a normal range in both groups [42].
Concerning the overall oral status, a higher percentage of mothers with periodontal disease was observed in the OB group (68.8%) compared to NW (56.2%), though not significantly. Among the recorded measures of dental and periodontal health, the number of teeth was significantly lower (p = 0.027) in the OB vs. NW group. The dental plaque index percentage was significantly higher in OB women vs. NW (p = 0.024).
All the enrolled women delivered at term, without obstetrical complications. No difference was recorded in the mode of delivery, although a higher percentage of cesarean section was observed in the OB group (NW: 12.5%; OB: 43.8%), in agreement with other epidemiological data from obese pregnancies.
Table 3 reports neonatal and placental data at delivery.
Neonates were similar in terms of gestational age and biometric parameters (weight and centiles); as expected, a higher number, though not significant, of Large for Gestational Age (LGA) babies was observed from OB mothers. Fetal sex frequencies did not differ between the two study groups.
Placentas were significantly heavier (p = 0.035) and thicker (p = 0.050) in the OB than in the NW group.

3.2. miRNAs’ Profile in Maternal Saliva

In maternal saliva, the presence of 754 miRNAs was evaluated by a simultaneous quantitation with the TaqMan OpenArray Human MicroRNA Panel, which contains the most biologically relevant miRNAs in miRBase (https://www.mirbase.org/, accessed on 1 March 2011). Data were adjusted for gestational age at saliva withdrawal, maternal age, and smoking habits, as well as GDM association as covariates.
The expression of 334 miRNAs was detectable in the saliva of at least one subject of our study population. Among these miRNAs, 31 were differentially expressed (FC < 0.5 or > 2.0) with statistical significance (p < 0.05) between OB and NW mothers, as shown by the volcano plot (Figure 2).
Following FDR correction, 20 miRNAs were accepted for the downstream analysis (FDR < 0.20). In particular, 2 miRNAs out of 20 (hsa-miR-505 and hsa-miR616) had FC < 0.5, thus resulting in being downregulated in the OB group. FC > 2.0 indicates upregulated levels for 18 miRNAs in OB women. FC, p-value, and FDR are summarized in Table 4.
To gain insights into the miRNA profiling results, miTALOS and DIANA–miRPath were used for miRNA pathway enrichment analysis.
First, all 20 miRNAs were inserted together in miTALOS. Interestingly, the union analysis (i.e., combinations of pathways related to the inserted miRNAs) showed a significant association between the 20 miRNAs and a specific pathway involving TGF-Beta signaling (p = 0.040).
An additional analysis was performed by DIANA–miRPath v3. A heat map depicting the force of association was obtained after inserting the 20 miRNAs differentially expressed in OB vs. NW maternal saliva. Four pathways matches were pointed out by the DIANA–miRPath pathway enrichment analysis (Figure 3 and Table 5).

3.3. DNA Methylation in Maternal Saliva: TGF-Beta1 and SOCS3

Methylation data were corrected for maternal age, smoking habits, and GDM association. Both TGF-Beta1 and SOCS3 were differentially methylated, with significantly decreased methylation levels in OB mothers compared to NW (p = 0.019 and p = 0.025, respectively) (Table 6).

4. Discussion

Maternal obesity and Gestational Diabetes (GDM) have been associated with lipotoxicity, inflammation, and increased oxidative stress both at systemic and intrauterine levels [8,10,44]. We previously reported placental modifications in obese (OB) pregnancies, indicating an excessive oxidative intrauterine/systemic environment [13,15,45]. Increased antioxidant capacity and C-reactive protein levels in both maternal plasma and the saliva of obese women with and without GDM have also been reported, confirming saliva as an effective fluid mirroring the systemic status, representing an interesting non-invasive tool for the evaluation of metabolic signatures during pregnancy [23].
The obesogenic environment can influence pregnancy outcomes via epigenetic mechanisms, including microRNAs (miRNAs) and DNA methylation [46,47]. Alterations occurring in maternal obesity and GDM can affect epigenetic modifications of different maternal and fetal tissues, leading to alterations in several cellular pathways, which can affect the future development in the offspring of obesity, diabetes, and other metabolic and cardiovascular diseases [3,47,48,49]. Studies on maternal circulating miRNAs are currently in progress aiming at identifying specific miRNA profiles associated with pregnancy pathologies [50,51,52,53], some of them proposing specific miRNAs as early biomarkers for gestational diseases such as preeclampsia [50], obesity [51,52], or GDM [53]. Nevertheless, to our knowledge, this is the first study characterizing the miRNAs’ profile and DNA methylation in the saliva of obese compared to normal-weight pregnant women.
Mothers included in this study population had previously been carefully classified for periodontal comorbidities (gingivitis and periodontitis) [23]. Clinical conditions of the enrolled patients were well characterized during regular prenatal checks, with obese women receiving specific nutritional and lifestyle counseling in a dedicated clinic. Indeed, OB gained significantly less weight during pregnancy compared to normal-weight mothers. Nevertheless, obese women exceeded gestational weight gain IOM recommendations, possibly worsening the obesity-related adverse systemic and local environment. In fact, the poorest overall periodontal health status was observed in the OB group. Moreover, both placental weight and thickness were increased in OB, suggesting a trend to a lower placental efficiency, giving rise to a higher, though not significant, mean neonatal weight centile. This obesogenic context also resulted in alterations of the saliva epigenetic profile. In particular, we reported a set of 20 up regulated or down expressed miRNAs and decreased methylation levels of the genes SOCS3 and TGF-Beta1 in the saliva of obese vs. normal-weight mothers.

4.1. miRNAs’ Profile

The saliva shows a very high number of detectable miRNAs, having a unique spectrum of these non-coding RNAs [54].
We performed a pathway enrichment analysis by integrating different miRNA target prediction tools, to identify specific associations with signaling pathways. In particular, saliva miRNAs showing significant differences between OB and NW mothers, matched pathways involving fatty acids biosynthesis and metabolism, Extracellular Matrix (ECM)–receptor interaction, and lysine degradation.

4.1.1. Fatty Acids Biosynthesis and Metabolism

We found significant associations between miR-1254 and fatty acids biosynthesis and metabolism. Interestingly, our group recently reported lower levels of Long Chain-Polyunsaturated Fatty Acid (LC-PUFA) derivatives, arachidonic acid, and DHA, opposite significantly increased saturated palmitic acid levels, in obese placentas [15]. These data can indicate a disruption of the physiologic LC-PUFA biomagnification linked to maternal obesity, similar to alterations occurring in pregnancy pathologies characterized by placental dysfunction, such as intrauterine growth restriction. Moreover, changes in fatty acids biosynthesis and metabolism can damage cells’ and organelles’ membranes, leading to excessive ROS production and oxidative stress, characterizing the systemic and intrauterine environment of obese and GDM pregnancies [3,13,44]. Although this evidence is derived from different tissues, the present results may be suggestive of a disarranged fatty acid metabolism in obese pregnancies mediated by miRNA epigenetic alterations, which can also be detected in the oral non-invasive fluid during the third trimester of pregnancy.

4.1.2. Extracellular Matrix

Pathway enrichment analysis showed a significant match between hsa-miR-184, hsa-miR-206, hsa-miR-302c-5p, hsa-miR-376a-3p, hsa-miR-505-3p, hsa-miR-512-3p hsa-miR-635, and the ECM–receptor interaction pathway.
Intriguingly, human endothelial cell culture exposed to high levels of saturated fatty acids showed alterations in the structure and functionality of the components of the ECM [55]. Similarly, hyperglycemia induces excessive production of ECM collagens and altered proteoglycans by endothelial cells. Indeed, also in patients with type II diabetes, the ECM of blood vessels presented an altered ratio of proteoglycans and a thicker basement membrane vs. normoglycemic controls. These changes could represent pathogenic mechanisms for vascular diabetic complications, which have a strong inflammatory basis [56]. Physiologically, ECM remodeling consists of structural and compositional changes and is critical for the differentiation and expansion of adipocytes. In obesity-related adipose tissue hypertrophy, ECM synthesis and degradation can be altered, causing inflammation and mechanical stress [57]. Both ECM degradation and pro-inflammatory milieu establishment are promoted by the activity of Matrix Metalloproteinases (MMPs), which are well-known targets of several miRNAs, suggesting their epigenetic regulation [58]. In turn, pro-inflammatory cytokines and pathogens associated with periodontal disease raise MMPs’ secretion, thus favoring disintegration of the ECM [59]. Therefore, the association between ECM pathways and OB-related salivary miRNAs’ alterations reported in the present study might be speculated to be part of the dysregulation of the ECM remodeling, due to the obesogenic environment, the linked inflammatory setting, and the oral disease condition.

4.1.3. Lysine Degradation

hsa-miR-302d-3p, hsa-miR-302c-3p, hsa-miR-376a-5p, hsa-miR-505-3p, hsa-miR-505-5p, hsa-miR-512-3p, and hsa-miR-616-3p resulted in being associated with the lysine degradation pathway.
Amino Acids (AAs) are necessary building blocks for protein biosynthesis. Lysine is crucial for cell growth and plays a fundamental role in the production of carnitine, which shuttles long-chain fatty acids into mitochondria for energy production and assists in lowering cholesterol levels [60].
Coherent with the results reported in this manuscript, decreased lysine levels in the adipose tissue of obese non-pregnant subjects have been shown [61], suggesting an altered metabolic pattern typical of the obese condition. Desert and colleagues reported increased lysine levels in the urinary excretion of obese mothers vs. normal-weight, hypothesizing an enhancement in lysine renal secretion [62]. In a metabolomics analysis of OB vs. NW placentas, we previously showed amino acid changes, with marked lower amounts of lysine [15]. Similar amino acid alterations were previously reported in the serum of hyperglycemic mothers, showing a metabolic profile consistent with insulin resistance [63,64].
Despite the consistency of this evidence, the underlying mechanism explaining the relation between obesity and lysine catabolism is still unknown. It is tempting to hypothesize that miRNAs could be involved in these mechanisms. Indeed, the lysine degradation pathway has been shown as a possible target of specific miRNAs related to a diabetic microenvironment in a beta-human cell line [65]. Targeted studies are needed to confirm this hypothesis.

4.2. DNA Methylation

4.2.1. TGF-Beta1

Transforming Growth Factor-Beta 1 (TGF-Beta1) is an anti-inflammatory cytokine with widespread and multi-organ effects. It plays an active role in glucose, lipid, amino acid, and redox metabolism and represents a key cytokine in obesity and insulin resistance [30,66,67,68]. Noteworthy evidence showed that TGF-Beta1 signaling inhibition ameliorated the metabolic profile in mice by increasing glucose and insulin tolerance. In turn, elevated TGF-Beta1 levels were associated with glucose intolerance and increased adiposity [69] and have been reported in women with a prior history of GDM [70]. Accordingly, in the present study, TGF-Beta1 methylation levels were significantly decreased in the saliva of obese vs. lean mothers, supporting its upregulation, even when adjusted for maternal age, smoking habits, and the presence of GDM. It might be hypothesized that OB’s TGF-Beta1 higher systemic levels are mirrored in the saliva, confirming this biofluid as a striking non-invasive tool for the investigation of local and systemic alterations. Moreover, the pathway enrichment analysis by miTALOS showed a noticeable significant association between the 20 miRNAs differentially expressed in OB and the TGF-Beta signaling pathway. This evidence supports previous results showing a miRNA-based regulation of TGF-Beta signaling [71] and represents an interesting starting point to unravel the alterations of the molecular regulatory mechanisms of this pathway in obesity, leading to the possible impairment of nutrients and redox metabolism, with the consequent promotion of oxidative stress.

4.2.2. SOCS3

Over the past few years, SOCS3 has emerged as an interesting target to treat metabolic disorders. Indeed, SOCS3 association with inflammation, cumulative stress, and insulin resistance in obesity and diabetes is well established, and its inhibition is considered a promising strategy for the treatment of metabolic disorders in different conditions [29]. In pregnancy, the possible participation of SOCS3 in the attenuation of leptin and insulin signaling has been recently reported, with increased resistance associated with SOCS3 upregulation during physiological pregnancy-induced metabolic changes [72].
Here, we report a significant decrease in SOCS3 methylation levels in the saliva of OB women compared to NW, likely accounting for exacerbated SOCS3 expression in pregnant women with higher BMI. In non-pregnant individuals, SOCS3 DNA methylation has been previously shown to be inversely related to body weight and to modulate the impact of overall stress on obesity [73]. Furthermore, obesity-related inflammation leads to the upregulation of SOCS3 proteins in several tissues, leading to the inhibition of insulin signaling and, therefore, representing an interesting mechanism connecting the immune and metabolic system in a delicate and balanced cross-talk [29,74]. Results reported in the present study suggest a disruption of this balance, and the altered SOCS3 DNA methylation may represent an intriguing mechanism undergoing OB metabolic alterations, deserving further future investigations both as a biomarker of the dysmetabolic context in pregnancy and as a possible future target for its treatment.

4.3. Strengths and Limitations

This work is a pilot study, therefore presenting a small sample size. Nevertheless, the analyzed population was very well defined in terms of inclusion and exclusion criteria, clinical characteristics, and associated conditions, which strongly reduced any additional clinical bias impacting the results. The strict statistical analysis that was applied included several covariates, allowing us to focus on the effect of obesity on the analyzed epigenetic modifications.
Molecular investigations were performed on the saliva. Indeed, this is an advantageous, non-invasive, and cost-effective biofluid. Saliva composition reflects perturbations occurring in the oral cavity, which might vary depending on the local conditions. However, it can also mirror the systemic status, representing an intriguing tool for biomarker non-invasive examination. Moreover, a careful procedure was employed to remove any contaminants and to avoid chemical alterations or the release of inflammatory mediators. Indeed, the passive drool technique for saliva collection is recommended when the target concentration is low, and the unstimulated collection allows obtaining high purity and great sample quality [54].

5. Conclusions

To our knowledge, this is the first study exploring the epigenetic status in the saliva of obese pregnant women. This pilot study can pave the way for future studies on pregnancy-related epigenetics. Among the mechanisms undergoing obese-related injuries, epigenetic modifications have gained increasing attention over the past years for their high potential as both disease biomarkers and treatment targets. We characterized the epigenetic signatures of a group of clinically well-characterized obese and normal-weight pregnant women by analyzing their saliva, which represents a promising non-invasive biofluid, which can mirror the systemic status of the individual. We found alterations in the expression of 20 different miRNAs and methylation levels of two genes involved in obesity-related inflammation and oxidative stress. These results help to clarify some of the mechanisms involved in the impairment occurring in obesity, which could contribute to the increase in the risk of negative outcomes in OB pregnancies. Furthermore, the underlying association between miRNAs and their related molecular pathways will help understanding specific damaged pathways. Since miRNAs could be promising biomarkers in the early diagnosis and treatment of pregnancy-adverse conditions, these results place the analyzed salivary miRNAs in a clinical context for future studies and propose possible molecular mechanisms in which they could be involved, representing future promising tools for the identification of obesity-related alterations. Further studies are needed for investigating maternal miRNA profiling related to post-natal follow-up data.

Supplementary Materials

The following Supporting Information can be downloaded at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/nu14102122/s1, Table S1: Inclusion and Exclusion Criteria of the Study Population.

Author Contributions

Conceptualization, I.C., S.A., C.M. and V.B.; investigation, G.M.A., B.A. and L.D.; formal analysis, C.F., G.M.A. and C.M.; resources, M.Z.; writing—original draft, C.M., G.M.A. and A.S.; writing—review and editing, I.C., S.A., C.M. and V.B.; supervision, I.C., C.M. and V.B.; funding acquisition, I.C. All authors have read and agreed to the published version of the manuscript.

Funding

PRIN 2010-2011 Prot. 20102chst5_005 “Parto pre-termine: markers molecolari, biochimici e biofisici dell’ unità feto-placentare”, grant from the Italian Ministry of University and Research to Irene Cetin. Progetto di Ricerca Finalizzata RF-2016-02362165, “Epigenetic impact of maternal obesity and nutritional status. Nutritional/lifestyle intervention for the improvement of pregnancy outcomes” (EPI-MOM), grant from the Italian Ministry of Health to Irene Cetin. The authors are thankful to Associazione Studio Malformazioni (ASM) for an unconditioned grant to the Laboratory of Maternal-Fetal Translational Research “Giorgio Pardi.”

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the “Medical Ethical and Institutional Review Board” (Prot. N. 469/2010/52/AP).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All data that support the findings of this study are available from the corresponding author (Chiara Mandò) upon reasonable request.

Acknowledgments

We thank all the midwives and nurses of the Obstetric Unit of Luigi Sacco Hospital (ASST Fatebenefratelli-Sacco Milano) for their expertise and cooperation. We are particularly grateful to all the pregnant women that contributed to the study with their clinical and biological data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NCD Risk Factor Collaboration (NCD-RisC). Trends in adult body-mass index in 200 countries from 1975 to 2014: A pooled analysis of 1698 population-based measurement studies with 19·2 million participants. Lancet 2016, 387, 1377–1396. [Google Scholar] [CrossRef] [Green Version]
  2. Schaefer-Graf, U.; Napoli, A.; Nolan, C.J. Diabetic Pregnancy Study Group. Diabetes in pregnancy: A new decade of challenges ahead. Diabetologia 2018, 61, 1012–1021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Catalano, P.M.; Shankar, K. Obesity and pregnancy: Mechanisms of short term and long term adverse consequences for mother and child. BMJ 2017, 356, j1. [Google Scholar] [CrossRef]
  4. Poston, L.; Caleyachetty, R.; Cnattingius, S.; Corvalán, C.; Uauy, R.; Herring, S.; Gillman, M.W. Preconceptional and maternal obesity: Epidemiology and health consequences. Lancet Diabetes Endocrinol. 2016, 4, 1025–1036. [Google Scholar] [CrossRef]
  5. Bordin, P.; Dotto, L.; Battistella, L.; Rosso, E.; Pecci, L.; Valent, F.; Collarile, P.; Vanin, M. Gestational diabetes mellitus yesterday, today and tomorrow: A 13 year italian cohort study. Diabetes Res. Clin. Pract. 2020, 167, 108360. [Google Scholar] [CrossRef]
  6. Monteiro, L.J.; Norman, J.E.; Rice, G.E.; Illanes, S.E. Fetal programming and gestational diabetes mellitus. Placenta 2016, 48 (Suppl. 1), S54–S60. [Google Scholar] [CrossRef] [Green Version]
  7. Barker, D.J.; Thornburg, K.L.; Osmond, C.; Kajantie, E.; Eriksson, J.G. Beyond birthweight: The maternal and placental origins of chronic disease. J. Dev. Orig. Health Dis. 2010, 1, 360–364. [Google Scholar] [CrossRef] [PubMed]
  8. Parisi, F.; Milazzo, R.; Savasi, V.M.; Cetin, I. Maternal Low-Grade Chronic Inflammation and Intrauterine Programming of Health and Disease. Int. J. Mol. Sci. 2021, 22, 1732. [Google Scholar] [CrossRef]
  9. Bianchi, C.; Taricco, E.; Cardellicchio, M.; Mandò, C.; Massari, M.; Savasi, V.; Cetin, I. The role of obesity and gestational diabetes on placental size and fetal oxygenation. Placenta 2021, 103, 59–63. [Google Scholar] [CrossRef]
  10. Zavatta, A.; Parisi, F.; Mandò, C.; Scaccabarozzi, C.; Savasi, V.M.; Cetin, I. Role of Inflammaging on the Reproductive Function and Pregnancy. Clin. Rev. Allergy Immunol. 2022, 1–16. [Google Scholar] [CrossRef]
  11. Kelly, A.C.; Powell, T.L.; Jansson, T. Placental function in maternal obesity. Clin. Sci. 2020, 134, 961–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Martino, J.; Sebert, S.; Segura, M.T.; García-Valdés, L.; Florido, J.; Padilla, M.C.; Marcos, A.; Rueda, R.; McArdle, H.J.; Budge, H.; et al. Maternal Body Weight and Gestational Diabetes Differentially Influence Placental and Pregnancy Outcomes. J. Clin. Endocrinol. Metab. 2016, 101, 59–68. [Google Scholar] [CrossRef] [Green Version]
  13. Mandò, C.; Anelli, G.M.; Novielli, C.; Panina-Bordignon, P.; Massari, M.; Mazzocco, M.I.; Cetin, I. Impact of Obesity and Hyperglycemia on Placental Mitochondria. Oxid. Med. Cell. Longev. 2018, 2018, 2378189. [Google Scholar] [CrossRef]
  14. Diceglie, C.; Anelli, G.M.; Martelli, C.; Serati, A.; Lo Dico, A.; Lisso, F.; Parisi, F.; Novielli, C.; Paleari, R.; Cetin, I.; et al. Placental Antioxidant Defenses and Autophagy-Related Genes in Maternal Obesity and Gestational Diabetes Mellitus. Nutrients 2021, 13, 1303. [Google Scholar] [CrossRef] [PubMed]
  15. Fattuoni, C.; Mandò, C.; Palmas, F.; Anelli, G.M.; Novielli, C.; Parejo Laudicina, E.; Savasi, V.M.; Barberini, L.; Dessì, A.; Pintus, R.; et al. Preliminary metabolomics analysis of placenta in maternal obesity. Placenta 2018, 61, 89–95. [Google Scholar] [CrossRef] [PubMed]
  16. Foratori-Junior, G.A.; da Silva, B.M.; da Silva Pinto, A.C.; Honório, H.M.; Groppo, F.C.; de Carvalho Sales-Peres, S.H. Systemic and periodontal conditions of overweight/obese patients during pregnancy and after delivery: A prospective cohort. Clin. Oral Investig. 2020, 24, 157–165. [Google Scholar] [CrossRef]
  17. Ye, C.; Kapila, Y. Oral microbiome shifts during pregnancy and adverse pregnancy outcomes: Hormonal and Immunologic changes at play. Periodontol. 2000 2021, 87, 276–281. [Google Scholar] [CrossRef]
  18. Kinane, D.F.; Stathopoulou, P.G.; Papapanou, P.N. Periodontal diseases. Nat. Rev. Dis. Primers 2017, 3, 17038. [Google Scholar] [CrossRef]
  19. Byun, J.S.; Lee, H.Y.; Tian, J.; Moon, J.S.; Choi, J.; Lee, S.H.; Kim, Y.G.; Yi, H.S. Effect of Salivary Exosomal miR-25-3p on Periodontitis with Insulin Resistance. Front. Immunol. 2022, 12, 775046. [Google Scholar] [CrossRef]
  20. Ren, H.; Du, M. Role of Maternal Periodontitis in Preterm Birth. Front. Immunol. 2017, 8, 139. [Google Scholar] [CrossRef] [Green Version]
  21. Cetin, I.; Pileri, P.; Villa, A.; Calabrese, S.; Ottolenghi, L.; Abati, S. Pathogenic mechanisms linking periodontal diseases with adverse pregnancy outcomes. Reprod. Sci. 2012, 19, 633–641. [Google Scholar] [CrossRef] [PubMed]
  22. Srinivas, S.K.; Parry, S. Periodontal disease and pregnancy outcomes: Time to move on? J. Womens Health 2012, 21, 121–125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zambon, M.; Mandò, C.; Lissoni, A.; Anelli, G.M.; Novielli, C.; Cardellicchio, M.; Leone, R.; Monari, M.N.; Massari, M.; Cetin, I.; et al. Inflammatory and Oxidative Responses in Pregnancies with Obesity and Periodontal Disease. Reprod. Sci. 2018, 25, 1474–1484. [Google Scholar] [CrossRef] [PubMed]
  24. Kaczor-Urbanowicz, K.E.; Martin Carreras-Presas, C.; Aro, K.; Tu, M.; Garcia-Godoy, F.; Wong, D.T. Saliva diagnostics—Current views and directions. Exp. Biol. Med. 2017, 242, 459–472. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Yoshizawa, J.M.; Schafer, C.A.; Schafer, J.J.; Farrell, J.J.; Paster, B.J.; Wong, D.T. Salivary biomarkers: Toward future clinical and diagnostic utilities. Clin. Microbiol. Rev. 2013, 26, 781–791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Weber, J.A.; Baxter, D.H.; Zhang, S.; Huang, D.Y.; Huang, K.H.; Lee, M.J.; Galas, D.J.; Wang, K. The microRNA spectrum in 12 body fluids. Clin. Chem. 2010, 56, 1733–1741. [Google Scholar] [CrossRef]
  27. Tsamou, M.; Martens, D.S.; Winckelmans, E.; Madhloum, N.; Cox, B.; Gyselaers, W.; Nawrot, T.S.; Vrijens, K. Mother’s Pre-pregnancy BMI and Placental Candidate miRNAs: Findings from the ENVIRONAGE Birth Cohort. Sci. Rep. 2017, 7, 5548, Erratum in Sci. Rep. 2018, 8, 6063. [Google Scholar] [CrossRef]
  28. Al-Rawi, N.H.; Al-Marzooq, F.; Al-Nuaimi, A.S.; Hachim, M.Y.; Hamoudi, R. Salivary microRNA 155; 146a/b and 203: A pilot study for potentially non-invasive diagnostic biomarkers of periodontitis and diabetes mellitus. PLoS ONE 2020, 15, e0237004. [Google Scholar] [CrossRef]
  29. Pedroso, J.A.B.; Ramos-Lobo, A.M.; Donato, J., Jr. SOCS3 as a future target to treat metabolic disorders. Hormones 2019, 18, 127–136. [Google Scholar] [CrossRef]
  30. Yener, S.; Demir, T.; Akinci, B.; Bayraktar, F.; Kebapcilar, L.; Ozcan, M.A.; Biberoglu, S.; Yesil, S. Transforming growth factor-beta 1 levels in women with prior history of gestational diabetes mellitus. Diabetes Res. Clin. Pract. 2007, 76, 193–198. [Google Scholar] [CrossRef]
  31. Page, R.C.; Eke, P.I. Case definitions for use in population-based surveillance of periodontitis. J. Periodontol. 2007, 78 (Suppl. 7), 1387–1399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Pergoli, L.; Cantone, L.; Favero, C.; Angelici, L.; Iodice, S.; Pinatel, E.; Hoxha, M.; Dioni, L.; Letizia, M.; Albetti, B.; et al. Extracellular vesicle-packaged miRNA release after short-term exposure to particulate matter is associated with increased coagulation. Part. Fibre Toxicol. 2017, 14, 32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  34. Andersen, C.L.; Jensen, J.L.; Ørntoft, T.F. Normalization of real-time quantitative reverse transcription-PCR data: A model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64, 5245–5250. [Google Scholar] [CrossRef] [Green Version]
  35. Bollati, V.; Baccarelli, A.; Hou, L.; Bonzini, M.; Fustinoni, S.; Cavallo, D.; Byun, H.M.; Jiang, J.; Marinelli, B.; Pesatori, A.C.; et al. Changes in DNA methylation patterns in subjects exposed to low-dose benzene. Cancer Res. 2007, 67, 876–880. [Google Scholar] [CrossRef] [Green Version]
  36. Tarantini, L.; Bonzini, M.; Tripodi, A.; Angelici, L.; Nordio, F.; Cantone, L.; Apostoli, P.; Bertazzi, P.A.; Baccarelli, A.A. Blood hypomethylation of inflammatory genes mediates the effects of metal-rich airborne pollutants on blood coagulation. Occup. Environ. Med. 2013, 70, 418–425. [Google Scholar] [CrossRef] [Green Version]
  37. Preusse, M.; Theis, F.J.; Mueller, N.S. miTALOS v2: Analyzing Tissue Specific microRNA Function. PLoS ONE 2016, 11, e0151771. [Google Scholar] [CrossRef]
  38. Vlachos, I.S.; Kostoulas, N.; Vergoulis, T.; Georgakilas, G.; Reczko, M.; Maragkakis, M.; Paraskevopoulou, M.D.; Prionidis, K.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA miRPath v.2.0: Investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012, 40, W498–W504. [Google Scholar] [CrossRef]
  39. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef]
  40. Hosack, D.A.; Dennis, G., Jr.; Sherman, B.T.; Lane, H.C.; Lempicki, R.A. Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4, R70. [Google Scholar] [CrossRef]
  41. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Soc. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  42. Iorgulescu, G. Saliva between normal and pathological. Important factors in determining systemic and oral health. J. Med. Life 2009, 2, 303–307. [Google Scholar] [PubMed]
  43. Bertino, E.; Spada, E.; Occhi, L.; Coscia, A.; Giuliani, F.; Gagliardi, L.; Gilli, G.; Bona, G.; Fabris, C.; De Curtis, M.; et al. Neonatal anthropometric charts: The Italian neonatal study compared with other European studies. J. Pediatr. Gastroenterol. Nutr. 2010, 51, 353–361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Saben, J.; Lindsey, F.; Zhong, Y.; Thakali, K.; Badger, T.M.; Andres, A.; Gomez-Acevedo, H.; Shankar, K. Maternal obesity is associated with a lipotoxic placental environment. Placenta 2014, 35, 171–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Anelli, G.M.; Cardellicchio, M.; Novielli, C.; Antonazzo, P.; Mazzocco, M.I.; Cetin, I.; Mandò, C. Mitochondrial content and hepcidin are increased in obese pregnant mothers. J. Matern. Fetal Neonatal Med. 2018, 31, 2388–2395. [Google Scholar] [CrossRef] [PubMed]
  46. Franzago, M.; Fraticelli, F.; Stuppia, L.; Vitacolonna, E. Nutrigenetics, epigenetics and gestational diabetes: Consequences in mother and child. Epigenetics 2019, 14, 215–235. [Google Scholar] [CrossRef] [Green Version]
  47. Kerr, B.; Leiva, A.; Farías, M.; Contreras-Duarte, S.; Toledo, F.; Stolzenbach, F.; Silva, L.; Sobrevia, L. Foetoplacental epigenetic changes associated with maternal metabolic dysfunction. Placenta 2018, 69, 146–152. [Google Scholar] [CrossRef]
  48. Guarino, E.; Delli Poggi, C.; Grieco, G.E.; Cenci, V.; Ceccarelli, E.; Crisci, I.; Sebastiani, G.; Dotta, F. Circulating MicroRNAs as Biomarkers of Gestational Diabetes Mellitus: Updates and Perspectives. Int. J. Endocrinol. 2018, 2018, 6380463. [Google Scholar] [CrossRef] [Green Version]
  49. Jing, J.; Wang, Y.; Quan, Y.; Wang, Z.; Liu, Y.; Ding, Z. Maternal obesity alters C19MC microRNAs expression profile in fetal umbilical cord blood. Nutr. Metab. 2020, 17, 52. [Google Scholar] [CrossRef]
  50. Luque, A.; Farwati, A.; Crovetto, F.; Crispi, F.; Figueras, F.; Gratacós, E.; Aran, J.M. Usefulness of circulating microRNAs for the prediction of early preeclampsia at first-trimester of pregnancy. Sci. Rep. 2014, 4, 4882. [Google Scholar] [CrossRef] [Green Version]
  51. Sørensen, A.E.; van Poppel, M.N.M.; Desoye, G.; Simmons, D.; Damm, P.; Jensen, D.M.; Dalgaard, L.T.; The Dali Core Investigator Group. The Temporal Profile of Circulating miRNAs during Gestation in Overweight and Obese Women with or without Gestational Diabetes Mellitus. Biomedicines 2022, 10, 482. [Google Scholar] [CrossRef] [PubMed]
  52. Carreras-Badosa, G.; Bonmatí, A.; Ortega, F.J.; Mercader, J.M.; Guindo-Martínez, M.; Torrents, D.; Prats-Puig, A.; Martinez-Calcerrada, J.M.; Platero-Gutierrez, E.; De Zegher, F.; et al. Altered Circulating miRNA Expression Profile in Pregestational and Gestational Obesity. J. Clin. Endocrinol. Metab. 2015, 100, E1446–E1456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Poirier, C.; Desgagné, V.; Guérin, R.; Bouchard, L. MicroRNAs in Pregnancy and Gestational Diabetes Mellitus: Emerging Role in Maternal Metabolic Regulation. Curr. Diabetes Rep. 2017, 17, 35. [Google Scholar] [CrossRef] [PubMed]
  54. Wren, M.E.; Shirtcliff, E.A.; Drury, S.S. Not all biofluids are created equal: Chewing over salivary diagnostics and the epigenome. Clin. Ther. 2015, 37, 529–539. [Google Scholar] [CrossRef] [Green Version]
  55. Camejo, G.; Hurt-Camejo, E.; Olsson, U.; Bondjers, G. Lipid mediators that modulate the extracellular matrix structure and function in vascular cells. Curr. Diabetes Rep. 1999, 1, 142–149. [Google Scholar] [CrossRef]
  56. Olsson, U.; Bondjers, G.; Camejo, G. Fatty acids modulate the composition of extracellular matrix in cultured human arterial smooth muscle cells by altering the expression of genes for proteoglycan core proteins. Diabetes 1999, 48, 616–622. [Google Scholar] [CrossRef]
  57. Pessentheiner, A.R.; Ducasa, G.M.; Gordts, P.L.S.M. Proteoglycans in Obesity-Associated Metabolic Dysfunction and Meta-Inflammation. Front. Immunol. 2020, 11, 769. [Google Scholar] [CrossRef]
  58. Brovkina, O.; Nikitin, A.; Khodyrev, D.; Shestakova, E.; Sklyanik, I.; Panevina, A.; Stafeev, I.; Menshikov, M.; Kobelyatskaya, A.; Yurasov, A.; et al. Role of MicroRNAs in the Regulation of Subcutaneous White Adipose Tissue in Individuals with Obesity and Without Type 2 Diabetes. Front. Endocrinol. 2019, 10, 840. [Google Scholar] [CrossRef]
  59. Chang, Y.C.; Yang, S.F.; Lai, C.C.; Liu, J.Y.; Hsieh, Y.S. Regulation of matrix metalloproteinase production by cytokines; pharmacological agents and periodontal pathogens in human periodontal ligament fibroblast cultures. J. Periodontal Res. 2002, 37, 196–203. [Google Scholar] [CrossRef]
  60. Tomé, D.; Bos, C. Lysine requirement through the human life cycle. J. Nutr. 2007, 137 (Suppl. 6), 1642S–1645S. [Google Scholar] [CrossRef] [Green Version]
  61. Hanzu, F.A.; Vinaixa, M.; Papageorgiou, A.; Párrizas, M.; Correig, X.; Delgado, S.; Carmona, F.; Samino, S.; Vidal, J.; Gomis, R. Obesity rather than regional fat depots marks the metabolomic pattern of adipose tissue: An untargeted metabolomic approach. Obesity 2014, 22, 698–704. [Google Scholar] [CrossRef] [PubMed]
  62. Desert, R.; Canlet, C.; Costet, N.; Cordier, S.; Bonvallot, N. Impact of maternal obesity on the metabolic profiles of pregnant women and their offspring at birth. Metabolomics 2015, 11, 1896–1907. [Google Scholar] [CrossRef] [Green Version]
  63. Scholtens, D.M.; Muehlbauer, M.J.; Daya, N.R.; Stevens, R.D.; Dyer, A.R.; Lowe, L.P.; Metzger, B.E.; Newgard, C.B.; Bain, J.R.; Lowe, W.L., Jr.; et al. Metabolomics reveals broad-scale metabolic perturbations in hyperglycemic mothers during pregnancy. Diabetes Care 2014, 37, 158–166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Jansson, N.; Rosario, F.J.; Gaccioli, F.; Lager, S.; Jones, H.N.; Roos, S.; Jansson, T.; Powell, T.L. Activation of placental mTOR signaling and amino acid transporters in obese women giving birth to large babies. J. Clin. Endocrinol. Metab. 2013, 98, 105–113. [Google Scholar] [CrossRef] [Green Version]
  65. Jeffery, N.; Harries, L.W. miRNAs responsive to the diabetic microenvironment in the human beta cell line EndoC-βH1 may target genes in the FOXO, HIPPO and Lysine degradation pathways. Exp. Cell Res. 2019, 384, 111559, Erratum in Exp. Cell Res. 2019, 385, 111654. [Google Scholar] [CrossRef]
  66. Liu, H.; Chen, Y.G. The Interplay Between TGF-β Signaling and Cell Metabolism. Front. Cell Dev. Biol. 2022, 10, 846723. [Google Scholar] [CrossRef]
  67. Lee, J.H.; Mellado-Gil, J.M.; Bahn, Y.J.; Pathy, S.M.; Zhang, Y.E.; Rane, S.G. Protection from β-cell apoptosis by inhibition of TGF-β/Smad3 signaling. Cell Death Dis. 2020, 11, 184. [Google Scholar] [CrossRef] [Green Version]
  68. Toyoda, S.; Shin, J.; Fukuhara, A.; Otsuki, M.; Shimomura, I. Transforming growth factor β1 signaling links extracellular matrix remodeling to intracellular lipogenesis upon physiological feeding events. J. Biol. Chem. 2022, 298, 101748. [Google Scholar] [CrossRef]
  69. Yadav, H.; Devalaraja, S.; Chung, S.T.; Rane, S.G. TGF-β1/Smad3 Pathway Targets PP2A-AMPK-FoxO1 Signaling to Regulate Hepatic Gluconeogenesis. J. Biol. Chem. 2017, 292, 3420–3432. [Google Scholar] [CrossRef] [Green Version]
  70. Zhu, Y.; Tian, F.; Li, H.; Zhou, Y.; Lu, J.; Ge, Q. Profiling maternal plasma microRNA expression in early pregnancy to predict gestational diabetes mellitus. Int. J. Gynecol. Obstet. 2015, 130, 49–53. [Google Scholar] [CrossRef]
  71. Kim, Y.J.; Hwang, S.J.; Bae, Y.C.; Jung, J.S. MiR-21 regulates adipogenic differentiation through the modulation of TGF-beta signaling in mesenchymal stem cells derived from human adipose tissue. Stem Cells 2009, 27, 3093–3102. [Google Scholar] [CrossRef] [PubMed]
  72. Zampieri, T.T.; Ramos-Lobo, A.M.; Furigo, I.C.; Pedroso, J.A.; Buonfiglio, D.C.; Donato, J., Jr. SOCS3 deficiency in leptin receptor-expressing cells mitigates the development of pregnancy-induced metabolic changes. Mol. Metab. 2014, 4, 237–245. [Google Scholar] [CrossRef] [PubMed]
  73. Xu, K.; Zhang, X.; Wang, Z.; Hu, Y.; Sinha, R. Epigenome-wide association analysis revealed that SOCS3 methylation influences the effect of cumulative stress on obesity. Biol. Psychol. 2018, 131, 63–71. [Google Scholar] [CrossRef] [PubMed]
  74. Lebrun, P.; Van Obberghen, E. SOCS proteins causing trouble in insulin action. Acta Physiol. 2008, 192, 29–36. [Google Scholar] [CrossRef]
Figure 1. Data processing workflow. OB: Obese, NW: Normal-Weight, Crt: Relative Threshold Cycle, FDR: False Discovery Rate; FC: Fold Change; Bold: highlights the most relevant terms.
Figure 1. Data processing workflow. OB: Obese, NW: Normal-Weight, Crt: Relative Threshold Cycle, FDR: False Discovery Rate; FC: Fold Change; Bold: highlights the most relevant terms.
Nutrients 14 02122 g001
Figure 2. Volcano plot showing differentially expressed miRNAs between OB and NW women. The graph reports the 334 detectable miRNAs (all black and red dots). The horizontal line states the –log10 of the p-value 0.05; the two vertical lines indicate the log2 of the FC values 0.5 and 2.0. As a result, miRNAs with a significantly different expression vs. NW are displayed in red (n = 31). In particular, the upper left quadrant reports the downregulated miRNAs, while the upper right one all the upregulated miRNAs.
Figure 2. Volcano plot showing differentially expressed miRNAs between OB and NW women. The graph reports the 334 detectable miRNAs (all black and red dots). The horizontal line states the –log10 of the p-value 0.05; the two vertical lines indicate the log2 of the FC values 0.5 and 2.0. As a result, miRNAs with a significantly different expression vs. NW are displayed in red (n = 31). In particular, the upper left quadrant reports the downregulated miRNAs, while the upper right one all the upregulated miRNAs.
Nutrients 14 02122 g002
Figure 3. Heat map by DIANA–miRPath v.3. The analyzed miRNAs are listed on the right. The 4 specific pathways, resulting in being significant from miRNA pathway enrichment analysis, are listed at the bottom. The force of association (Log(p-value)) between the analyzed miRNAs and the 4 pathways is indicated by the color key reported at the upper left corner. ECM: ExtraCellular Matrix, Log: log10 p-value.
Figure 3. Heat map by DIANA–miRPath v.3. The analyzed miRNAs are listed on the right. The 4 specific pathways, resulting in being significant from miRNA pathway enrichment analysis, are listed at the bottom. The force of association (Log(p-value)) between the analyzed miRNAs and the 4 pathways is indicated by the color key reported at the upper left corner. ECM: ExtraCellular Matrix, Log: log10 p-value.
Nutrients 14 02122 g003
Table 1. Pyrosequencing assays information.
Table 1. Pyrosequencing assays information.
GeneChromosome
Position
CpG SitesPrimers:
Forward (F)
Reverse (R)
Sequencing (S)
Sequencing
Length (BP)
Annealing Temperature (°C)
TGF-Beta1chr17:6354421-763548212F: TTGGGTGATTTTTTTATAGGAGTT
R: bio-TCCCCCCAAAAAAACCTATT
S: GAGATGTTGAAGAGTGGTTA
2552
SOCS3chr19:41353024-413534582F: GGTTTGTTTTTTGAGTTTT
R:bio-CTACAAAAACTAAAAATCTCCC
S: TATTTATTTTTTGGTATTAG
2354
BP: Base Pairs, °C: Celsius degrees.
Table 2. Maternal and periodontal health characteristics.
Table 2. Maternal and periodontal health characteristics.
NW
n = 16
All OB
n = 16
p-Value
Age, years B30.94 ± 3.8933.25 ± 4.61ns
Pregestational BMI, kg/m2 B21.02 ± 2.2536.47 ± 5.53<0.001
GWG, kg B13.38 ± 4.419.13 ± 5.630.035
GWG to IOM upper limit, % A83.61 ± 27.54101.37 ± 62.57ns
Basal glycemia, mg/dL B81.50 ± 5.4892.69 ± 15.180.034
Gestational diabetes, n (%)0 (0)10 (62.5)-
GA at saliva withdrawal, mL A33.24 ± 1.7332.99 ± 2.52ns
Saliva flow rate, mL/min A0.48 ± 0.220.47 ± 0.17ns
Periodontal health C
Healthy, n (%)7 (43.8)5 (31.2)ns
Periodontal disease, n (%)9 (56.2)11 (68.8)
Number of teeth B27.60 ± 0.7425.31 ± 3.010.027
BOP, % sites B17.98 ± 19.7639.57 ± 38.10ns
PPD, mean A2.38 ± 0.452.63 ± 0.88ns
Plaque index, % A25.61 ± 18.9752.90 ± 40.320.024
Calculus, % A27.79 ± 23.6544.53 ± 35.81ns
Values expressed as the mean ± standard deviation were analyzed according to their distribution with independent samples A: Student t-test, or B: Mann–Whitney U test, or C: chi-squared test for independence or Fisher’s exact test; statistical significance vs. normal-weight women. BMI: Body Mass Index; GWG: Gestational Weight Gain at delivery; IOM: Institute of Medicine; maternal basal glycemia: maternal fasting glycemia refers to the first value of the Oral Glucose Tolerance Test (OGTT; physiological value ≤ 92 mg/dL) performed between 24 and 28 weeks of gestation; GA: Gestational Age; flow rate: ratio between mL of saliva and minutes of withdrawal; oral disease: gingivitis and/or periodontitis; PPD: Probing Pocket Depth; BOP: gingival Bleeding On Probing; ns: not statistically significant.
Table 3. Neonatal and placental data at delivery.
Table 3. Neonatal and placental data at delivery.
NW
n = 16
All OB
n = 16
p-Value
GA at delivery, wks A39.76 ± 0.8739.33 ± 1.18ns
Neonatal weight, gr A3331.25 ± 294.603488.75 ± 356.62ns
Neonatal weight centile A47.56 ± 24.1961.00 ± 28.72ns
AGA, n (%) C15 (93.8)11 (68.8)ns
LGA, n (%) C1 (6.3)5 (31.3)
Neonatal sex, n C ns
Males, n (%)9 (56.3)7 (43.8)
Females, n (%)7 (43.8)9 (56.3)
Placental weight, gr A439.00 ± 79.05509.29 ± 80.810.035
Neonatal/placental weight B7.78 ± 1.666.97 ± 1.63ns
Placental area, cm2 B281.94 ± 64.06266.16 ± 45.33ns
Placental thickness, cm A1.61 ± 0.371.97 ± 0.500.050
Values are expressed as the mean ± standard deviation. Data were analyzed according to their distribution with independent samples A: Student t-test, or B: Mann–Whitney U test, or C: chi-squared test for Independence or Fisher’s exact test; statistical significance vs. normal-weight women. GA: Gestational Age; Neonatal Weight Centile: calculated using INeS Charts (http://www.inescharts.com/index.aspx (accessed on 15 March 2022)); refer to [43]; AGA: Appropriate for Gestational Age; LGA: Large for Gestational Age.
Table 4. Upregulated and downregulated miRNAs in maternal saliva.
Table 4. Upregulated and downregulated miRNAs in maternal saliva.
miRNA NameFold Change
(FC)
p-ValueFalse Discovery Rate
p-Value (FDR)
hsa-miR-505 ↓0.3350.00020.0396
hsa-miR-616 ↓0.4340.01230.1873
hsa-miR-618 ↑7.1410.00020.0396
hsa-miR-206 ↑3.9850.00050.0514
hsa-miR-376a ↑5.8700.00080.0671
hsa-miR-517c ↑2.7880.00150.0989
hsa-miR-133a ↑3.7030.00410.1599
hsa-miR-512 ↑23.2950.00440.1599
hsa-miR-302d ↑14.4270.00460.1599
hsa-miR-520b ↑2.7550.00510.1599
hsa-miR-1254 ↑5.4430.00540.1599
hsa-miR-133b ↑5.6340.00570.1599
hsa-miR-1285 ↑40.9560.00950.1859
hsa-miR-635 ↑2.1780.01020.1859
hsa-miR-551b ↑2.6140.01060.1859
hsa-miR-548b-5p ↑5.6920.01100.1859
hsa-miR-1256 ↑2.1690.01100.1859
hsa-miR-302c ↑2.6220.01110.1859
hsa-miR-184 ↑4.8170.01180.1873
hsa-miR-548c-5p ↑4.9120.01310.1909
FC (Fold Change), p-value, and FDR (False Discovery Rate) of the 20 miRNAs with significantly different expression levels in OB vs. NW. Data were analyzed with a multiple linear regression analysis adjusted for gestational age at saliva withdrawal, age, smoking habits, and GDM association. Results are presented in descending order of statistical significance. (↓: downregulated; ↑: upregulated miRNA in OB vs. NW).
Table 5. Union analysis performed with DIANA–miRPath v.3.
Table 5. Union analysis performed with DIANA–miRPath v.3.
Pathways Informationp-ValuemiRNAs
Fatty acids
biosynthesis
Synthesis of fatty acids.
Fatty acids are generally excessive in the Western diet.
2.1444 × 10−11hsa-miR-1254
Fatty acids
metabolism
Anabolic and catabolic processes involving fatty acids or related molecules.
Fatty acids are generally excessive in the Western diet.
0.0002hsa-miR-1254
Lysine
degradation
Catabolism of lysine (from dietary up-taken or intracellular proteins).
Lysine is generally excessive in the Western diet.
1.7153 × 10−7hsa-miR-505-5p;
hsa-miR-302d-3p;
hsa-miR-376a-5p;
hsa-miR-302c-3p;
hsa-miR-512-3p;
hsa-miR-505-3p;
hsa-miR-616-3p
ECM–receptor
interaction
Tissue and organ morphogenesis; maintenance of cell and tissue structure and function; adhesion, migration, differentiation, proliferation, and apoptosis; force-transmitting physical link with the cytoskeleton.<1 × 10−325hsa-miR-302c-5p;
hsa-miR-184;
hsa-miR-206;
hsa-miR-635; hsa-miR-512-3p;
hsa-miR-505-3p;
hsa-miR-376a-3p
ECM: ExtraCellular Matrix.
Table 6. TGF-Beta1 and SOCS3 methylation levels.
Table 6. TGF-Beta1 and SOCS3 methylation levels.
Mean OB
(95% CI)
Mean NW
(95% CI)
p-Value
TGF-Beta1
(% 5mC)
0.44
(0.24–0.63)
0.86
(0.61–1.11)
0.019
SOCS3
(% 5mC)
60.95
(57.20–64.70)
68.50
(63.75–73.25)
0.025
Results from a multiple logistic regression analysis adjusted for age, smoking habits, and GDM association. Data are shown as the mean with the 95% Confidence Interval (CI) of methylated cytosines’ percentage (% 5mC); p < 0.05 vs. NW. TGF-Beta1: Transforming Growth Factor-Beta 1; SOCS3: Suppressor of Cytokine Signalling 3.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mandò, C.; Abati, S.; Anelli, G.M.; Favero, C.; Serati, A.; Dioni, L.; Zambon, M.; Albetti, B.; Bollati, V.; Cetin, I. Epigenetic Profiling in the Saliva of Obese Pregnant Women. Nutrients 2022, 14, 2122. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14102122

AMA Style

Mandò C, Abati S, Anelli GM, Favero C, Serati A, Dioni L, Zambon M, Albetti B, Bollati V, Cetin I. Epigenetic Profiling in the Saliva of Obese Pregnant Women. Nutrients. 2022; 14(10):2122. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14102122

Chicago/Turabian Style

Mandò, Chiara, Silvio Abati, Gaia Maria Anelli, Chiara Favero, Anaïs Serati, Laura Dioni, Marta Zambon, Benedetta Albetti, Valentina Bollati, and Irene Cetin. 2022. "Epigenetic Profiling in the Saliva of Obese Pregnant Women" Nutrients 14, no. 10: 2122. https://0-doi-org.brum.beds.ac.uk/10.3390/nu14102122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop