Substantially Altered Expression Profile of Diabetes/Cardiovascular/Cerebrovascular Disease Associated microRNAs in Children Descending from Pregnancy Complicated by Gestational Diabetes Mellitus—One of Several Possible Reasons for an Increased Cardiovascular Risk

Gestational diabetes mellitus (GDM), one of the major pregnancy-related complications, characterized as a transitory form of diabetes induced by insulin resistance accompanied by a low/absent pancreatic beta-cell compensatory adaptation to the increased insulin demand, causes the acute, long-term, and transgenerational health complications. The aim of the study was to assess if alterations in gene expression of microRNAs associated with diabetes/cardiovascular/cerebrovascular diseases are present in whole peripheral blood of children aged 3–11 years descending from GDM complicated pregnancies. A substantially altered microRNA expression profile was found in children descending from GDM complicated pregnancies. Almost all microRNAs with the exception of miR-92a-3p, miR-155-5p, and miR-210-3p were upregulated. The microRNA expression profile also differed between children after normal and GDM complicated pregnancies in relation to the presence of overweight/obesity, prehypertension/hypertension, and/or valve problems and heart defects. Always, screening based on the combination of microRNAs was superior over using individual microRNAs, since at 10.0% false positive rate it was able to identify a large proportion of children with an aberrant microRNA expression profile (88.14% regardless of clinical findings, 75.41% with normal clinical findings, and 96.49% with abnormal clinical findings). In addition, the higher incidence of valve problems and heart defects was found in children with a prior exposure to GDM. The extensive file of predicted targets of all microRNAs aberrantly expressed in children descending from GDM complicated pregnancies indicates that a large group of these genes is involved in ontologies of diabetes/cardiovascular/cerebrovascular diseases. In general, children with a prior exposure to GDM are at higher risk of later development of diabetes mellitus and cardiovascular/cerebrovascular diseases, and would benefit from dispensarisation as well as implementation of primary prevention strategies.

and environmental programming predisposed the affected children to the development of diabetes mellitus, cardiovascular, and cerebrovascular diseases.
The hypothesis of the assessment of potential diabetes/cardiovascular risk in children prenatally exposed to GDM was founded on the fact that particular microRNAs play a key role in the inducement and progress of diabetes mellitus and cardiovascular/cerebrovascular diseases (Table 1)  .
Previously, by searching the Medline database we identified a large scale of microRNAs playing a role in pathogenesis of diabetes mellitus and cardiovascular/cerebrovascular diseases. Finally, we selected for the study a shortlist of 29 microRNAs demonstrated repetitively by numerous scientific teams to be associated with normal stages (development and homeostasis of the cardiovascular system, angiogenesis, and adipogenesis) and pathological conditions and diseases (vascular endothelial dysfunction and inflammation, hypoxia, hypertension and regulation of hypertension-related genes, obesity, dyslipidaemia, atherosclerosis and atherosclerotic plaque formation, insulin resistance, diabetes mellitus and diabetes-related complications, metabolic syndrome, cardiovascular diseases involving the blood vessels (coronary and peripheral artery diseases, carotid artery disease, pulmonary arterial hypertension, cerebrovascular disease, aortic and intracranial aneurysms), cardiovascular diseases involving the heart (congenital heart disease, cardiomyopathies, cardiac dysrhythmias, hypertensive heart disease, myocardial disease, valvular heart disease, inflammatory heart disease, rheumatic heart disease, pulmonary heart disease, and heart failure), chronic kidney disease, ischemia/reperfusion injury, cardiac regeneration and cachexia) ( Table 1)  .
Furthermore, the MiRWalk database (http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/) was utilized to acquire data on the dysregulation of microRNAs, in particular, the abovementioned human disease and phenotype ontologies and OMIM disorders.
To our knowledge, no data on expression profiles of microRNAs associated with diabetes mellitus, cardiovascular, and cerebrovascular diseases in whole peripheral venous blood (white blood cells) of children descending from GDM affected pregnancies are currently available. The only study which was reported in this field was dedicated to profiling of mir-15 family members, affecting insulin signalling pathway, in skeletal muscle biopsies of adult offspring of women with a history of GDM [222].  [55], gestational diabetes mellitus [9,13,20], diabetes mellitus [56][57][58].

Participants
The study had a prospective design, that ran from August 2016 to October 2019, and included Caucasian children aged 3 to 11 years descending from normal (n = 85) and GDM complicated pregnancies (n = 118). Children of equal age descending from pregnancies with a normal course of gestation were chosen as a control group.
Ninety-eight children were exposed to a GDM pregnancy on diet only, and 20 children descended from a GDM complicated pregnancy on the combination of diet and therapy (in nineteen pregnancies insulin was administrated and in one pregnancy metformin was prescribed). According to the Institutional Review Board (IRB) Guidelines for gestational diabetes mellitus mothers were divided into groups based on preconception BMI (below 18.5, 18.5-24.9, 25.0-29.9, above 30.0), the achievement of optimal total weight gain during pregnancy (12.5-18.0, 11.5-16.0, 7.0-11.5, 5.0-9.0 kg), and the achievement of optimal week weight gain during the second and the third trimesters of gestation (0.5-0.6, 0.4-0.5, 0.2-0.3, 0.2-0.3 kg) were monitored.
The clinical data of children descending from normal and GDM complicated pregnancies are summarized in Table 2. The preconceptional and gestational clinical data of mothers with normal and GDM complicated pregnancies are involved as well.
Children descending from normal pregnancies were healthy infants born after 37 completed weeks of gestation with the weight >2500 g. Normal pregnancies were those ones with the absence of medical, obstetrical, or surgical complications.
Gestational diabetes mellitus was diagnosed if any degree of glucose intolerance appeared for the first-time during gestation [223][224][225]. Concerning the diagnosis and classification of hyperglycemia in pregnancy the guidelines of The International Association of Diabetes and Pregnancy Study Groups (IADPSG) were followed [223]. The first screening was performed during the first trimester of gestation with the aim to detect women with overt diabetes (fasting plasma glucose level is ≥7.0 mmol/L) and women with GDM (fasting plasma glucose level ≥5.1-<7.0 mmol/L). The second screening phase, 2 h 75-g OGTT, was performed during 24-28 weeks of gestation only in those women not previously diagnosed to have overt diabetes or GDM. The second screening phase identified GDM if the fasting plasma glucose level was ≥5.1 mmol/L, or the 1 h plasma glucose was ≥10.0 mmol/L, or the 2-h plasma glucose was ≥8.5 mmol/L [223].
Children with inborn defects or chromosomal disorders, as well as children descending from pregnancies demonstrating other complications were excluded from the inclusion into the study.
Informed consent was obtained from all participants included in the study. Two independent Ethics Committees (one at the Institute for the Care of the Mother and Child, and the second one at the Third Faculty of Medicine, Charles University) approved the study (grant no. AZV 16-27761A, long-term monitoring of complex cardiovascular profile in the mother, foetus, and offspring descending from pregnancy-related complications, dates of approval: 27 March, 2014 and 28 May, 2015). All procedures were also in compliance with the Helsinki Declaration of 1975, as revised in 2000.

BP and Echocardiography Measurements, BMI Assessment
Standardized BP and echocardiography measurements, and BMI assessment were performed as previously described [226]. In brief, the average of the last two systolic blood pressure (SBP) and diastolic blood pressure (DBP) values was taken under consideration for data analyses. A normal BP was defined as SBP and DBP that was below the 90th percentile for gender, age, and height. Hypertension was defined as an average SBP or DBP ≥95th percentile on at least three separate occasions. Prehypertension was defined as an average SBP or DBP within the range of ≥90th percentile and <95th percentile [227].
The age-and sex-specific BMI Percentile Calculator was used to calculate BMI in children (https: //www.cdc.gov/healthyweight/assessing/bmi/childrens_bmi/about_childrens_bmi.html). BMI above the 5th percentile and below the 85th percentile was considered as normal BMI. Children with BMI above the 85th percentile and below the 95th percentile were in the overweight category. Children with BMI equal to or greater than the 95th percentile were in the obese category.
A complete two-dimensional echocardiography was performed using the Philips HD15 ultrasound machine (Philips Ultrasound, Bothell, WA, USA) with a sector array transducer (3-8 MHz) incorporating colour flow, pulse wave, and continuous wave Doppler measurements with an adaptive technology. A complete two-dimensional echocardiography was performed by an investigator experienced with paediatric echocardiography. Children with abnormal findings were referred to a paediatric cardiologist.

Processing of Samples
Samples of unclotted whole peripheral venous blood (200 µL) were processed as previously described [226]. Briefly, homogenized cell lysates were prepared as soon as possible after blood collection using a QIAamp RNA Blood Mini Kit (Qiagen, Hilden, Germany, no. 52304).

Reverse Transcription
Particular microRNAs were transcribed into cDNA using microRNA-specific stem-loop RT primers, components of TaqMan MicroRNA Assays and TaqMan MicroRNA Reverse Transcription Kit (Applied Biosystems, Branchburg, NJ, USA, no. 4366597) as previously described [226]. A total reaction volume of the reaction was 10 µL. The reverse transcriptase reaction was performed following the guidelines for a 7500 Real-Time PCR system (Applied Biosystems, Branchburg, NJ, USA): 30 min at 16 • C, 30 min at 42 • C, 5 min at 85 • C, and then held at 4 • C.

Relative Quantification of microRNAs
Relative quantification of microRNAs by real-time PCR was performed as previously described [226]. cDNA (3 µL) was mixed with specific TaqMan MGB primers and probes (TaqMan MicroRNA Assay, Applied Biosystems, Branchburg, NJ, USA), and the constituents of the TaqMan Universal PCR Master Mix (Applied Biosystems, Branchburg, NJ, USA, no: 4318157) in a total reaction volume of 15 µL. The samples were regarded as positive if Ct (threshold cycle) was below 40 (Ct < 40).
The comparative Ct method was used to determine the expression of each microRNA [228]. The expression of studied microRNAs was normalized to a geometric mean of RNU58A and RNU38B, endogenous controls with the lowest expression variability between studied samples [229]. A reference sample was used throughout the study for relative quantification. As a reference sample we used small RNAs enriched RNA fraction extracted from the fetal part of one randomly selected placenta of a normally ongoing pregnancy.

Statistical Analysis
Logistic regression was used to compare the presence of abnormal clinical findings (BMI in the category overweight or obesity and/or systolic or diastolic BP values in the category prehypertension or hypertension and/or the presence of valve problems and heart defects) among various groups of children.
The Shapiro-Wilk test was used to assess the data normality [230]. Our experimental data did not follow a normal distribution. Therefore, microRNA levels were compared among the groups of children using the Kruskal-Wallis one-way analysis of variance with a post-hoc test (K-W) for the comparison between multiple groups.
Receivers operating characteristic (ROC) curves were used to assess the areas under the curves (AUC), the optimal cut-off points, and the respective sensitivities at 10.0% false positive rate (FPR) for particular microRNAs (MedCalc Software bvba, Ostend, Belgium). The significance level was established at a p-value of p < 0.05.
To identify the optimal combinations of microRNA biomarkers logistic regression combined with the ROC curve analysis was used (MedCalc Software bvba, Ostend, Belgium). In brief, in this setting, the power of the model's predicted values to discriminate between positive and negative cases is quantified by the area under the ROC curve. To perform a full ROC curve analysis the predicted probabilities are first saved and next used as a new variable in the ROC curve analysis. The dependent variable used in logistic regression then acts as the classification variable in the ROC curve analysis dialog box.
Statistica software (version 9.0; StatSoft, Inc., Tulsa, OK, USA) was used to generate box plots of log-normalized gene expression values (RT-qPCR expression, log 10 2 −∆∆Ct ) for particular microRNAs. The box plots display the medians, the 75th and 25th percentiles (the upper and lower limits of the boxes), the maximum and minimum values (the upper and lower whiskers), outliers (circles), and extremes (asterisks). Dot plots, all observations are also displayed in the charts.

Bioinformatics Analysis-microRNA-Target Interactions on Disease Ontologies, Human Phenotype Ontologies, and OMIM Disorders
MiRWalk database (http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/) and the Predicted Target module [231] were utilized to acquire data on predicted targets for microRNAs that have been identified to be dysregulated in whole peripheral venous blood of children descending from GDM complicated pregnancies.
MiRWalk is a comprehensive database that provides in addition to other things also information on human microRNA predicted and/or validated target genes, information on microRNA-target interactions on disease ontologies, human phenotype ontologies and OMIM disorders, and information on possible interactions between microRNAs and genes associated with KEGG, Panther and Wiki pathways.
Only those predicted targets involved in ontologies of human diseases (obesity, hypertension, glucose intolerance, lipid metabolism disease, type 2 diabetes mellitus, heart septal defects and heart valve disease, heart disease, heart failure, venous insufficiency, and pulmonary embolism) reported by previous population-based cohort studies in children and young adults descending from GDM complicated pregnancies [26][27][28][29][30][31][32][33][34][35][36] are presented below. Figure 1 presents the workflow of the study including the particular groups of the participants, the performance of relative quantification of microRNAs and bioinformatics analyses, and the variety of statistical analyses applied.

Higher Incidence of Valve Problems and Heart Defects in a group of Children Descending from GDM Complicated Pregnancies
We distributed children descending from normal and GDM complicated pregnancies into groups based on the results of anamnesis and the results of consecutive clinical examination. The

Higher Incidence of Valve Problems and Heart Defects in a group of Children Descending from GDM Complicated Pregnancies
We distributed children descending from normal and GDM complicated pregnancies into groups based on the results of anamnesis and the results of consecutive clinical examination. The group of children with abnormal clinical findings consisted of those ones who were already dispensarised in the department of paediatric cardiology, or those ones who were overweight/obese, had prehypertension/hypertension, and/or abnormal echocardiogram findings during the visit (in total: FG, n = 35; GDM, n = 57).
In detail, the studied groups of children with abnormal clinical findings consisted of those ones already dispensarised in the department of paediatric cardiology (FG, n = 8/85; GDM, n = 5/118), those ones indicated by the sonographer during the visit to have valve problems and heart defects ( Logistic regression revealed no difference in the incidence of overweight/obesity and/or prehypertension/hypertension between the groups of children descending from normal and GDM complicated pregnancies. Nevertheless, a higher incidence of valve problems and heart defects was observed in a group of children descending from GDM complicated pregnancies (Table 3). Logistic regression was used to compare the presence of abnormal clinical findings between particular groups. The significance level was established at a p-value of p < 0.05. No difference in the incidence of overweight/obesity and/or prehypertension/hypertension was found between the groups of children descending from normal and GDM complicated pregnancies. A higher incidence of valve problems and heart defects was observed in a group of children descending from GDM complicated pregnancies.

Only Just miR-21-5p Indicates a Trend Towards Differentiation between Children with Normal and Abnormal Clinical Findings (BMI, Blood Pressure, and Echocardiogram Findings) Descending from Normal Pregnancies
Furthermore, we compared the microRNA expression profile between individual groups of children with respect to clinical findings. From the ROC curve analysis, it is apparent that only miR-21-5p (a sensitivity of 17.14% at a specificity of 90.0%) trended to differentiate between children descending from normal pregnancies with normal and abnormal clinical findings (BMI, blood pressure, and echocardiogram findings) ( Figure S1). The performance of the ROC curve analysis revealed that miR-29a-3p (12.28%) and miR-92a-3p (12.28%) were able to slightly differentiate at 10.0% FPR within children descending from GDM complicated pregnancies based on the presence of normal and abnormal clinical findings (BMI, blood pressure, and echocardiogram findings) ( Figure S2).

Substantially Altered Expression Profile of Diabetes/Cardiovascular/Cerebrovascular Disease Associated microRNAs in Children Descending from Pregnancy Complicated by Gestational Diabetes Mellitus
With regard to the assessment of the effect of maternal pregnancy complication on postnatal microRNA expression profile, we further compared microRNA gene expression between children descending from normal and GDM complicated pregnancies irrespective of the clinical findings (overweight/obesity, prehypertension/hypertension, and/or valve problems and heart defects). The K-W and subsequent ROC curve analysis revealed a significant dysregulation of multiple microRNAs in children descending from GDM complicated pregnancies (Figure 2, Figure S3).

Discussion
A diabetes/cardiovascular/cerebrovascular disease associated microRNA expression profile was assessed in whole peripheral blood of children at the age of 3 to 11 years with a prior exposure to GDM with the aim to assess to what extent fetal and environmental programming predisposes the affected children to later development of diabetes mellitus, cardiovascular, and cerebrovascular diseases.
The hypothesis of the assessment of potential diabetes/cardiovascular risk in children prenatally exposed to GDM was based on the knowledge that a serious of microRNAs play a role in the pathogenesis of diabetes mellitus and cardiovascular/cerebrovascular diseases (Table 1)  .
Screening based on the combination of microRNAs with a good sensitivity only in the ROC curve analysis was superior over using individual microRNAs for the prediction of potential risk of later development of diabetes mellitus and/or cardiovascular/cerebrovascular diseases, since it was able to identify 88.14% children with an aberrant microRNA expression profile at 10.0% FPR.
Subsequently, we performed analyses to see how the microRNA expression profile differed in relation to the current absence or presence of cardiovascular risk factors and cardiovascular complications (overweight/obesity, prehypertension/hypertension, and/or valve problems and heart defects) and simultaneously to the previous occurrence of maternal pregnancy complication (GDM).
As expected, screening based on the combination of microRNAs with a good sensitivity in the ROC curve analysis only was superior over using individual microRNAs in the assessment of potential risk of later development of diabetes mellitus and cardiovascular/cerebrovascular diseases in groups of children descending from GDM pregnancies with either normal or abnormal clinical findings. At 10.0% FPR, it was able to identify 75.41% and 96.49% children with an aberrant microRNA expression profile.
Subsequently, to assess the impact of cardiovascular risk factors and cardiovascular complications (overweight/obesity, prehypertension/hypertension, and/or valve problems and heart defects) on the microRNA expression profile we further compared the microRNA gene expression profile between equal groups of children descending from pregnancies that had a normal or abnormal course of gestation.
Similarly, as in our previous study [226] dedicated to cardiovascular risk assessment in children descending from pregnancies affected with gestational hypertension, preeclampsia and/or fetal growth restriction, the expression profile of microRNAs was equal between the groups of children descending from uncomplicated pregnancies with normal and abnormal clinical findings, with an exception of miR-21-5p, which showed a trend towards upregulation in a proportion of children with abnormal clinical findings (17.14%).
Parallel, the presence of cardiovascular risk factors (overweight/obesity and/or prehypertension/ hypertension) and cardiovascular complications (valve problems and heart defects) had little impact on the microRNA expression profile of children with a prior exposure to GDM, since only miR-29a-3p (12.28%) and miR-92a-3p (12.28%) differentiated between children with normal and abnormal clinical findings (BMI, blood pressure, and echocardiogram findings) with a poor sensitivity.
Previously, we demonstrated that a proportion of children affected with pregnancy-related complications such as gestational hypertension (GH), preeclampsia (PE), and/or fetal growth restriction (FGR) had alterations in the microRNA expression profile that may predispose these children to later development of cardiovascular/cerebrovascular diseases [226,232,233].
On the other hand, our study revealed that dysregulation of miR-16-5p, miR-92a-3p, miR-100-5p, miR-143-3p, miR-155-5p, miR-221-3p, miR-499a-5p, and miR-574-3p represented a unique feature of an aberrant expression profile of children with a prior exposure to GDM. From these findings, it is obvious, that a large proportion of children prenatally exposed to GDM has a substantially altered microRNA expression profile associated with diabetes mellitus and cardiovascular/cerebrovascular diseases.
Generally, epigenetics refers to DNA sequence independent alterations that are responsible for control of transcription. The classic epigenetic modifications include DNA methylation, post-translational modifications of histone proteins, silencing of the extra copy of the X chromosome in women, and genomic imprinting. In addition, proteins and protein complexes with epigenetic modifying capabilities have been involved under the definition of epigenetics. Furthermore, with the discovery of the RNA interference machinery, several classes of noncoding RNAs (microRNAs, small-interfering RNAs, and long noncoding RNAs) have been added to the definition of epigenetics [234]. Nevertheless, controversy still exists as to whether or not microRNAs should be considered as part of the epigenetic program [234][235][236]. While classical epigenetic mechanisms, such as histone modification and DNA methylation, regulate expression at the transcriptional level, microRNAs putatively function mainly at the posttranscriptional level [234][235][236].
Epigenetics mechanisms appear to be interconnected on multiple levels. DNA methylation may direct histone methylation or vice versa chromatin remodelling drives DNA methylation [235,237]. RNAi also appears to be interconnected with DNA methylation and histone modifications. In some species, the link between microRNAs and epigenetics is strong as shown in vitro by transfection of synthetic siRNAs. However, it remains to be seen whether endogenous microRNAs or other types of noncoding RNAs can be also linked to epigenetic mechanisms in vivo [235].
From our findings, it is obvious, that a large proportion of children prenatally exposed to GDM has substantially altered the microRNA expression profile associated with diabetes mellitus and cardiovascular/cerebrovascular diseases, which may be one of several possible reasons for an increased cardiovascular risk. In general, children with a prior exposure to GDM are at higher risk of later development of diabetes mellitus and cardiovascular/cerebrovascular diseases and would benefit from dispensarisation and implementation of primary prevention strategies.
The identification of valve problems and heart defects within the group of children descending from GDM complicated pregnancies more often than usual may support the suggestion for their dispensarisation and implementation of primary prevention strategies.
Similarly, as other studies [258][259][260] we observed that infertility and an infertility treatment may be associated with an increased risk of GDM onset. The incidence of infertility and infertility treatment was significantly higher in groups of mothers with a GDM complicated course of gestation when compared to the group of mothers with normal pregnancies (Table 2). In addition, GDM occurring after assisted reproductive technology conception had been reported to increase the risk of adverse obstetric and perinatal outcomes [259]. This finding may strengthen the idea of dispensarisation of children descending from GDM complicated pregnancies and implementation of primary prevention strategies in this risky group. With regard to a low number of children descending from GDM complicated pregnancies after an infertility treatment (n = 16), it is misleading to interpret the association between microRNA gene expression in whole peripheral venous blood of children and infertility treatment in mothers. Most of the significant results was achieved when the comparison between children descending from normal pregnancies and children descending from GDM complicated pregnancies regardless of infertility treatment was performed or when the comparison between the groups of children descending from normal and GDM complicated pregnancies without infertility treatment was made (83 NP and 102 GDM).

Conclusions
In conclusion, any of the tissue-specific and circulation-specific (plasma and/or serum) changes in the microRNA expression profile characteristic for patients with diabetes mellitus and cardiovascular/cerebrovascular diseases are also present in whole peripheral blood (white blood cells) of children previously exposed to GDM. This finding indicates that a previous occurrence of GDM may predispose affected children to later development of diabetes mellitus and cardiovascular/cerebrovascular diseases. Consecutive large-scale studies are needed to verify the findings resulting from this pilot study.

Supplementary Materials:
The following are available online at http://www.mdpi.com/2073-4409/9/6/1557/s1, Figure S1: Aberrant microRNA expression profile in children descending from normal pregnancies, Figure S2: Aberrant microRNA expression profile in children descending from GDM complicated pregnancies, Figure S3: Aberrant microRNA expression profile in children descending from GDM complicated pregnancies irrespective of the clinical findings (overweight/obesity, prehypertension/hypertension, and/or valve problems and heart defects), Figure S4: Aberrant microRNA expression profile in children descending from GDM complicated pregnancies with normal clinical findings., Figure S5: Aberrant microRNA expression profile in children descending from GDM complicated pregnancies with abnormal clinical findings.; Table S1: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to obesity using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions), Table S2: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to hypertension using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S3: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to glucose intolerance using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S4: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to lipid metabolism disease using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S5: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to type 2 diabetes using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S6: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to heart septal defect using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S7: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to heart valve disease using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S8: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to heart disease using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S9: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to heart failure using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S10: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to venous insufficiency using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions); Table S11: A list of predicted targets of appropriate microRNAs dysregulated in whole peripheral blood of children descending from GDM complicated pregnancies in relation to pulmonary embolism using miRWalk2.0 database (putative microRNA binding sites predicted by miRWalk algorithm within mRNA selected regions).

Acknowledgments:
We thank the staff of the Institute for the Care of Mother and child (Jana Kumprichtova and Sarka Stranska) for assistance with the collection of patients' biological samples.

Conflicts of Interest:
The authors declare no conflict of interest.