Construction of a novel miRNA regulatory network and identification of target genes in gestational diabetes mellitus by integrated analysis

Abstract

Backgrounds: Given the roles of microRNA (miRNA) in human diseases and the high incidence of gestational diabetes mellitus (GDM), the aim of the study was to examine miRNA signatures and crucial pathways, as well as possible biomarkers for GDM diagnosis.

Methods: We conducted a two-stage study to explore functional miRNA and those target genes. Twelve participants (6 GDM and 6 non-GDM) were first enrolled and performed RNA sequencing analysis. The overlapped candidate genes were further screened in combination with differentially expressed genes (DEGs) of GEO datasets (GSE87295, GSE49524 and GSE19649) and potential target genes of DEMs. Candidate genes, critical pathways, small molecular compounds and regulatory networks were identified using bioinformatic analysis. The potential candidate genes were then investigated using the GEO dataset (GSE103552) of 19 participants in the validation stage (11 GDM and 8 non-GDM women).

Results: Briefly, blood samples were sequenced interrogating 50 miRNAs, including 20 upregulated and 30 downregulated differentially expressed microRNAs(DEMs) in our internal screening dataset. After screening GEO databases, 123 upregulated and 70 downregulated genes were overlapped through DEGs of GEO datasets and miRNA-target genes. MiR-29b-1-5p-TGFB2, miR-142-3p-TGFB2, miR-9-5p-FBN2, miR-212-5p-FBN2, miR-542-3p-FBN1, miR-9-5p-FBN1, miR-508-3p-FBN1, miR-493-5p-THBS1, miR-29b-3p-COL4A1, miR-432-5p-COL5A2, miR-9-5p-TGFBI, miR-486-3p-SLC7A5 and miR-6515-5p-SLC1A5 were revealed as thirteen possible regulating pathways by integrative analysis.

Conclusion: Overall, thirteen candidate miRNA-target gene regulatory pathways representing potentially novel biomarkers of GDM diseases were revealed. Ten chemicals were identified as putative therapeutic agents for GDM. This study examined a series of DEGs that are associated with epigenetic alternations of miRNA through an integrated approach and gained insight into biological pathways in GDM. Precise diagnosis and therapeutic targets of GDM would be further explored through putative genes in the future.

Introduction

Gestational diabetes mellitus (GDM) is a comprehensive form of pregnancy-specific glucose intolerance or hyperglycemia that manifests certain degree of glucose intolerance (Plows et al., 2018). Hyperglycemia occurs in approximately 16.7% in pregnancies worldwide, 75%–90% of which is caused by GDM, implying that GDM has become a significant public health concern (IDA, 2021). Although blood glucose level in GDM usually returns to normal after delivery, women are at high risk of acquiring type 2 diabetes later in life, and the risk of metabolic syndrome and insulin resistance in offspring have also increased (Damm et al., 2016), which presents a vicious intergenerational cycle. Compelling evidence suggest that advanced maternal age, family history of diabetes, diet, physical activity or emerging environmental factors are likely to have an impact on the risk of developing GDM (Zhang et al., 2016). However, knowledge concerning the detailed processes governing the initiation and progression of GDM remains unknown.

Although Genome-wide association studies (GWAS) have revealed several genetic loci correlated with the complexity of the disease, the underlying mechanism remain unclear (Kwak et al., 2012). Non-coding RNAs (ncRNAs) are important players in metabolic processes and their deregulated expressions have been observed in several metabolic diseases, including GDM. MicroRNA (miRNA), as a type of generally ubiquitous and multifunctional short non-coding RNAs with 19–22 nucleotides that regulate post-transcriptional gene expression (Brennecke et al., 2005), and participate in a range of biological functions (Du and Zamore, 2005).

Currently, in spite of the presence of discordant data, numerous researches indicate that circulating miRNAs are involved in mediating the key pathophysiological features of GDM, including glucose homeostasis, inflammation, insulin resistance, metabolic adaptations and β-cell dysfunction (Sliwinska et al., 2017; Filardi et al., 2020; Abu Samra et al., 2022). Aberrant levels of miRNA and gene expression in GDM have laid a favorable foundation for personalized target therapy and potential drugs. Compelling studies show that miR-195-5p overexpression in GDM women play an important role in insulin insensitivity regulation (Tagoma et al., 2018; Wang et al., 2020). A study of Filardi et al. identified that miR-222-3p and miR-409-3p signatures were significantly up-regulated in GDM, which are correlated with fasting plasma glucose (Filardi et al., 2022). Concerted efforts are made to gain knowledge about the mechanisms by which metabolic pathways are coordinated by acquired and genetic factors to explore novel insights into GDM treatment.

Although common biomarkers for GDM have been identified previously, the specific molecular mechanisms remain unclear. Genetic susceptibility is represented by gene variants conferring individual differences in response to metabolic-related chronic diseases. Most importantly, few studies have used transcriptome sequencing to explore GDM biomarkers, especially in South-East Chinese population. To the best of our knowledge, systematic regulatory network and pathways have not been well profiled about how alternations of microRNAs and mRNA are related to GDM development integrating with chemical compounds. The understanding of their functions might help unravel the complex pathophysiological mechanisms and identify novel clinical treatment. Early prevention and diagnosis are of great important to avoid adverse effects.

In this present study, a complete transcriptome RNA based on GDM in peripheral blood leukocytes (PBL) was sequenced, and putative hub genes were identified using a stepwise screening approach. The target genes of differentially expressed miRNA (DEMs) were synchronously predicted via online miRNA-target databases. The transcriptome were utilized to reveal molecular pathways and protein-protein interaction (PPI) network of candidate genes and confirmed those hub genes, followed by an external validation on the expression levels through the GSE103552 dataset (11 GDM and 8 normal controls). Simultaneously, by querying through CMap, registered chemicals could be screened and the likelihood of drugs could be obtained, showing gene expression profile that discover potential small molecules targeting diseases. This study aimed to identify pivotal pathways for complex pathophysiological mechanisms for GDM and contextually explore novel potential diagnostic biomarkers for the disease.

Materials and methods

Sample collection

In the first RNA-sequencing screening stage, blood samples were collected form 6 GDM from the Nantong Maternal and Child Health Hospital (NMCHH) in October 2020, and in the meantime 6 GDM-free women were randomly selected from a pool of more than 100 individuals who participated in routine healthcare examination in the same hospital. Non-GDM women were matched with GDM cases according to the age.

The methods for diagnosing GDM were mentioned in our study previously (Shen et al., 2020). Those participants conform to the diagnosis criteria of GDM [International Association of Diabetes Pregnancy Study Group (IADPSG)] plus complete demographic information recruited (Metzger et al., 2010). GDM patients with complications such as diabetes mellitus, chronic hypertension, pre-eclampsia and inflammatory diseases were excluded.

The study was reviewed and approved by the Ethics Committee of Shanghai University of Medicine & Health Sciences. All participants signed written informed consent.

Library preparation and sequencing

White blood cells were extracted by centrifugation at 1,500 g for 20 min with 2 ml whole blood. We utilized TRIzol (Invitrogen, Carlsbad, CA, United States ) to extract total RNAs following the manufacture’s protocol.

YueDa Biotechnology Co., Ltd. (Shanghai, China) performed RNA-seq using 150 ng of total RNA as input, and the results were analyzed by an Illumina HiSeq 2,500 sequencing platform with 10 M reads (Illumina, San Diego, CA, United States ). Bowtie were used to compare DEMs (Langmead, 2010), and feature Counts were adopted to annotate and quantify miRNAs (Liao et al., 2014). Counts were assessed and DEMs were filtered using the DESeq2 package in R (http://bioconductor.org/packages/release/bioc/html/DESeq2.html). Herein, miRNAs with fold changes (FC) > 1.5 or <0.667 and p < 0.05 were considered as significance.

Prediction of downstream target genes of DEMs

Target genes were concurrently predicted using Targetscan 7.2 (http://www.targetscan.org/vert_72/), miRDB (http://mirdb.org/), and miRwalk (http://mirwalk.umm.uni-heidelberg.de/) based on the aforementioned main DEM analysis. The anticipated target genes that better fit among three databases were regarded as those target genes of DEMs for a really steady selection.

Data collection and identification of overlapped candidate genes

In this present study, mRNA expression profiling datasets by array (GSE87295, GSE49524 and GSE19649) were retrieved and obtained from the GEO database (Gene Expression Omnibus, https://www.ncbi.nlm.nih.gov/geo/). Of these, 5 HUVECs GDM samples and 5 HUVECs controls were enrolled in GSE87295 (platform: GPL10558); 3 Caucasian Gestational diabetes women and 3 Caucasian non diabetic women in GSE49524 (platform: GPL7020); 3 GDM and 2 control in GSE19649 (platform: GPL7350).

We identified differentially expressed genes (DEGs) for each of three datasets separately by using linear models for microarray (LIMMA) approach through GEO2R online tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/). The screening thresholds for promising DEGs were set at p < 0.05 and fold change (FC) > 1.5 or <0.667. Subsequently, DEGs and potential target genes of DEMs were overlapped to obtain candidate genes. Herein, those genes of intersection between two datasets would be considered as being included. The ggVennDiagram software package in R (https://cran.r-project.org/web/packages/ggVennDiagram/index.html) was used to plot the Venn diagrams.

Functional and pathway enrichment analysis

Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were conducted for candidate genes by DAVID v6.8 (Database for Annotation, Visualization and Integrated Discovery) (https://david.ncifcrf.gov/home.jsp). p < 0.05 for GO analysis, and p < 0.1 and count>2 for the KEGG were considered for further analysis.

Candidate small molecules discovery in CMap

We used Connectivity Map (CMap) to explore potential therapeutic agents related to GDM. CMap database (https://clue.io/query) is an open database that predict those potential small molecular compounds of altered expression of DEGs in cell lines, presenting a connectivity score from -100 to 100. Score closer to 100 indicates that gene list is more similar change to the molecule. Conversely, a negative score indicates that small molecular compounds express antagonism, which could be candidate molecules for the treatment of GDM.

Protein-protein interaction network construction and screening of hub genes

The investigation of protein-protein interaction (PPI) network is crucial in assessing the disease’s molecular process. The Search Tool for the Retrieval of Interacting Genes (STRING) online tool (http://stringdb.org/) was employed to construct a network of candidate genes. The node pairs with a combined score of less than 0.4 were selected for further exploration. The network was visualized using Cytoscape v3.9.1, and hub genes were screened according to degree using CytoHubba, a Cytoscape plugin. The Maximal Clique Centrality (MCC) method was used to choose the top 30 target genes.

Expression analysis of hub genes based on GSE103552

The GSE103552 database [platform GPL6244; 11 GDM and 8 normal foetoplacental arterial endothelial cells (AEC)] was used to examine the expression level of hub genes for external verification. The criteria for determining hub genes also definitely are consistent with the previous demonstration.

Statistical analysis

Continuous variables adhering to the normal distribution were represented as the mean ± standard deviation; otherwise, the interquartile range (P₂₅-P₇₅) was substituted. The difference of continuous variables was tested by the independent sample t-tests or the Mann-Whitney tests. Categorical variables were represented as n (proportion), and the difference were tested by the χ² tests. Body mass index (BMI) was categorized four parts according to the Working Group on Obesity in China recommended criteria (Underweight: BMI < 18.5; Normal: 18.5 ≤ BMI < 24; Overweight: 24 ≤ BMI < 28; and General obesity: BMI ≥ 28 kg/m²) (Zhou, 2002). A value of p < 0.05 was considered statistically significant. All analysis was performed with SPSS 20.0 (IBM Corp. Chicago, IL) and GraphPad Prime9.0 (GraphPad Software, Inc.).

Results

Figure 1 depicts a schematic representation of the study design. And, Table 1 summarizes the characteristics of 12 samples. There was no significant difference between GDM patients and healthy controls in terms of age, BMI, 2h-plasma glucose, systolic blood pressure (SBP) and diastolic blood pressure (DBP) (p > 0.05). GDM patients had a statistically significant higher fasting plasma glucose (FPG) and 1h-plasma glucose (p < 0.05) as compared to controls.

FIGURE 1

Identification and prediction of target genes of DEMs

In our sequencing dataset, a total of 20 upregulated and 30 downregulated miRNAs were explored and identified based on the selection standard of p < 0.05 and fold change (FC) > 1.5 or < 0.667. The volcano map of the database was presented in Figure 2. We utilized the web tools TargetScan, miRDB, and miRwalk to explore the overlapped target genes of miRNAs by Venn diagram analysis, as shown in Table2. There were 3,867 target genes for upregulated DEMs and 3,856 target genes for downregulated DEMs after removing those duplicates.

FIGURE 2

Identification of candidate genes

The mRNA expression profiling datasets by array (GSE87295, GSE49524 and GSE19649) were employed to screen DEGs via GRO2R and further analyzed the hub genes and pathways. Briefly, 324 DEGs were detected in HUVECs samples from 5 GDM and 5 control, including 180 upregulated and 144 downregulated genes in GSE87295. 334 DEGs were detected in 3 Caucasian GDM women and 3 healthy controls in GSE49524, including 206 upregulated genes and 128 downregulated genes. Meanwhile, 161 DEGs were evaluated in GSE104297, including 109 upregulated mRNAs and 52 downregulated mRNAs(3 GDM and 2 non-GDM).

After screening of internal 12 samples and 3 GEO databases, the comprehensive datasets shared 193 overlapped candidate genes (123 upregulated genes and 70 downregulated genes). The overlapped genes were displayed in Figures 3A,B. The DEMs-gene network was established to analyze the relationship between DEMs and genes intuitively (Figures 4A,B). Among these DEMs, there were 67 target genes for 14 upregulated miRNAs and 121 ones for 28 downregulated miRNAs. Owing to the unavailability of the intersection with DEGs for those remaining miRNAs’ target genes, miRNA-mRNA relationship have not been identified. Furthermore, with a view to 5 genes representing the intersection of two DEGs, miRNAs of these genes were unavailable [ANKRD16, STAT1 (upregulated genes) and IGFBP6, PLAT, PLAC9 (downregulated genes)].

FIGURE 3

FIGURE 4

Enrichment analysis of candidate genes

To explore biological features and enriched pathways of those candidate genes, GO and KEGG analysis were accomplished by DAVID online tools. GO enrichment results were shown in Figures 5A,B, respectively. The upregulated genes were primarily enriched in paracrine signaling, desmosome assembly, sequestering of TGF-β in extracellular matrix, etc., in the BP group; platelet alpha granule lumen, transcription factor complex, caveola, etc., in the CC group; platelet-derived growth factor binding, extracellular matrix structural constituent, type III transforming growth factor beta receptor binding, etc., in the MF group. While, the downregulated genes were mainly involved in glutamine transport, negative regulation of plasminogen activation, melanocyte proliferation, platelet-derived growth factor receptor-beta signaling pathway, etc., in the BP group; external side of apical plasma membrane, focal adhesion, lysosomal lumen, etc., in the CC group; L-glutamine transmembrane transporter activity, neutral amino acid transmembrane transporter activity, amino acid transmembrane transporter activity, etc., in the MF group.

FIGURE 5

Subsequently, KEGG were conducted to explore the enrichment analysis of these candidate genes. Upregulated genes were mostly associated with pancreatic cancer, TGF-β signaling pathway, AGE-RAGE signaling pathway in diabetic complications, growth hormone synthesis, secretion and action and phosphatidylinositol signaling system, whereas downregulated genes were mainly associated with proteoglycans in cancer, PI3K-Akt signaling pathway, fluid shear stress and atherosclerosis and regulation of actin cytoskeleton, as shown in Figures 6A,B.

FIGURE 6

Related small molecule compounds screening

Screening results of potential drugs of GDM therapy were downloaded from CMap, ranking based on connectivity scores. The top 10 small molecule compounds identified as potential options for GDM treatment were penicillic-acid, lacidipine, YC-1, RITA, ALW-II-49-7, SA-792709, isoliquiritigenin, VX-222, CNQX and WH-4023, respectively (Table 3).

TABLE 3

Chemical name	Chemical formula	Type	Score	Description
RITA		cp	−95.52	MDM inhibitor
penicillic-acid		cp	−94.11	other antibiotic
isoliquiritigenin		cp	−92.86	Guanylate cyclase activator
lacidipine		cp	−92.6	Calcium channel blocker
ALW-II-49–7		cp	−92.3	Ephrin inhibitor
SA-792709		cp	−86.22	Retinoid receptor agonist
VX-222		cp	−83.49	HCV inhibitor
WH-4023	—	cp	−83.19	SRC inhibitor
YC-1		cp	−81.88	Guanylyl cyclase activator
CNQX		cp	−80.91	Glutamate receptor antagonist

Ten compounds identified as potential GDM therapeutics.

Chemical formula were from Pubchem. cp denotes compound.

Construction of hub genes network

The PPI network was constructed by STRING and displayed by Cytoscape to identify hub genes. The top 30 hub genes for upregulated and downregulated genes were detected by the Maximal Clique Centrality (MCC) of CytoHubba, respectively (shown in Figures 7A,B). Then, a total of the top 20 hub genes were selected to validate their expression level using external dataset (Table 4).

FIGURE 7

Validation of hub genes expression

In this study, through GEO dataset GSE103552, 20 hub genes were validated, including COL3A1, TGFB2, FBN2, FBN1, THBS1, COL4A1, COL5A2, TGFBR2, TGFBI, PDGFB (upregulated genes) and PDGFRB, LOX, NOS3, IGF2, SLC1A5, PLAT, SLC7A5, PYCR1, SLC38A1, FAP (downregulated genes).

The expression of TGFB2, FBN2, FBN1, THBS1, COL4A1, COL5A2 and TGFBI presented a similar trend of upregulation with the bioinformatic analysis; SLC7A5 and SLC1A5 expression were similar to the previous findings. Whereas, the expression of FAP was incompatible with previous screening result, and the expression of other mRNAs were not significantly different (p > 0.05) (Figures 8, 9). Therefore, miR-29b-1-5p-TGFB2, miR-142-3p-TGFB2, miR-9-5p-FBN2, miR-212-5p-FBN2, miR-542-3p-FBN1, miR-9-5p-FBN1, miR-508-3p-FBN1, miR-493-5p-THBS1, miR-29b-3p-COL4A1, miR-432-5p-COL5A2, miR-9-5p-TGFBI, miR-486-3p-SLC7A5 and miR-6515-5p-SLC1A5 were revealed as thirteen possible regulating pathways in our study.

FIGURE 8

FIGURE 9

Discussion

As one of the most common complications of pregnancy, the biological mechanism of GDM remains to be clearly elucidated. Exploring the aetiology and progression of GDM, as well as supporting the development of disease-modifying treatments, is a clear and urgent need. miRNAs have potential function in essential biological activities and their dysregulation or dysfunction was revealed in metabolic researches regarding GDM, rendering them a potential role as biomarkers or therapeutic targets. By merging internal RNA-seq data and GEO datasets with integrated bioinformatic analysis, we were able to produce solid evidence to considerably intensify the likelihood of detecting candidate biological markers and tremendously improve the reliability of our findings. In our study, some core modules and critical signaling pathway were discovered and extracted from dysregulated ceRNA that shed a light on plausible etiology of GDM. Collectively, the present study highlighted the effect of dysregulated glycometabolism and hormone-related in GDM and revealed the complexity of miRNAs as important fine tune regulators in the biological processes and their potential as novel biomarkers and treatment targets in GDM.

Our study identified 193 candidate genes by overlapping DEGs and target genes of DEMs. The GO analysis exhibited that candidate genes were primarily involved in post-embryonic eye morphogenesis, paracrine signaling, sequestering of TGF-β in extracellular matrix, platelet-derived growth factor receptor-beta signaling pathway, melanocyte proliferation and glutamine transport. These result indicate that GDM is related to placental development and endocrine. Previously, a study reported that GDM could alter angiogenesis by modulating paracrine factors (Loegl et al., 2017). A study suggested that platelet derived growth factor receptor beta polypeptide (Pdgfrb) expression were decreased in maternal hyperglycemia compared to controls in animal model (Lehtoranta et al., 2016). Wang et al. found glutamine capacity significantly decreased in the fetuses of GDM (Wang et al., 2019). Interestingly, the upregulated genes were mostly found in TGF-β signaling pathway, AGE-RAGE signaling pathway in diabetic complications, Growth hormone synthesis, secretion and action and Pancreatic cancer via KEGG analysis, and these pathways were involved in inflammatory response. It suggests that inflammatory response are also involved in the pathology of GDM development.

Since there are no effective drugs for GDM, the online database was used to aid the prediction of some drugs. The analysis via CMap suggested that the top 10 small molecule compounds could have potential therapeutic effect on GDM. Isoliquiritigenin, as one of the most important chalcone compounds, presents the antidiabetic activity and plays a part in the suppression of inflammatory pathways (Gaur et al., 2014; Gupta et al., 2018; Zhao et al., 2019). Soluble guanylate cyclase (sGC) activator YC-1 mimicked Bradykinin (BK) enhance the uptake of insulin-stimulated glucose (Frigolet et al., 2017). However, more experimental studies are necessary to validate therapeutic effects of these potential drugs on GDM. And, their application to the clinical settings requires extensive basic research and clinical trials in the future.

The PPI network revealed the intercommunication of candidate miRNAs and target genes, among which 13 potential miRNA-target gene regulatory pathways were miR-29b-1-5p-TGFB2, miR-142-3p-TGFB2, miR-9-5p-FBN2, miR-212-5p-FBN2, miR-542-3p-FBN1, miR-9-5p-FBN1, miR-508-3p-FBN1, miR-493-5p-THBS1, miR-29b-3p-COL4A1, miR-432-5p-COL5A2, miR-9-5p-TGFBI, miR-486-3p-SLC7A5 and miR-6515-5p-SLC1A5, respectively. For those upregulated DEMs, plasma miR-486-3p was upregulated in T2DM patients compared to healthy controls (Meerson et al., 2019) and increased concentrations in prepubertal obesity (Prats-Puig et al., 2013), which indicate that miR-486-3p could participate in a metabolism-related mechanism of GDM. Among those downregulated DEMs, miR-9-5p was upregulated in endocrine pancreatic (EN) of pluripotent stem cells (hiPSCs) compared to undifferentiated hiPSCs, and the pancreas has a regulatory role in glycemia (Sebastiani et al., 2017). However, miR-9-5p had higher levels in GDM patients (Martínez-Ibarra et al., 2019). A study reported that miR-29b-3p expression level was decreased in mice obesity model induced by high-fat diet (Zhang et al., 2022), while miR-29b-3p showed higher levels in the early phase of placentation in GDM(6–15 weeks of gestation) (Gillet et al., 2019). These might be related to the differences of sample source and gestation periods. Additionally, miR-542-3p, miR-493-5p, miR-432-5p and miR-6515-5p were reported to participate in the progression of breast cancer (Jiang et al., 2022), diabetic osteoporosis (Zhai et al., 2020), pancreatic ductal adenocarcinoma (Shen et al., 2021), inflammation-related disorders (periodontitis) (Xu et al., 2020). Those critical miRNAs, in particular, play significant roles in regulating biological processes in GDM.

In our study, interestingly, the expression of 9 genes (TGFB2, FBN2, FBN1, THBS1, COL4A1, COL5A2, TGFBI, SLC7A5, SLC1A5), being consistent with the result in GSE103552, were further verified. TGFB2 encodes a secreted ligand of the transforming growth factor-beta (TGF-β) superfamily of proteins, which is involved in the pathogenesis of diabetes and complications. Zhou et al. reported that TGFB2 is an upregulated gene in the plasma of GDM patients and correlated with FBG levels in GDM patients (Zhou et al., 2021), which suggested that TGFB2 might play vital roles in the pathogenesis of GDM. FBN2 encodes a peptide hormone placensin, which stimulates cAMP-PKA signaling, glucose secretion and trophoblast invasion in human trophoblastic cells. During third trimester, serum placensin levels of GDM patients are increased to a bigger extent compared to healthy pregnant women (Yu et al., 2020). FBN1 can encode the protein hormone asprosin, which is involved in regulating glucose homeostasis (Goodarzi et al., 2021). Zhong et al. reported that asprosin was highly expressed in the plasma of GDM patients and their offspring (Zhong et al., 2020). The protein encoded by THBS1 is an adhesive glycoprotein and proinflammatory cytokine that promotes insulin resistance. A study suggested that THBS1 in obesity or type 2 diabetes mellitus was highly expressed than healthy controls (Tang et al., 2020). COL4A1 was involved in diabetic tubulointerstitial injury (Zeng et al., 2019) and COL5A2 participated in the progression of uterine fibroids (Giri et al., 2017) and colon adenocarcinoma (Wei et al., 2020), although their functions were not fully understood for GDM development. TGFBI, upregulated in the placenta and the plasma of GDM, was positively associated with FBG levels in GDM patients (Han et al., 2014; Zhou et al., 2021). SLC7A5 enables L-leucine/L-tryptophan transmembrane transporter activity and is a part of amino acid transport complex. One study reported the uptake of ¹⁴C-l-methionine by human trophoblasts derived from normal pregnancies is mainly mediated by L-type amino acid transporter 1 [LAT1 (L), SLC7A5]. However, the process exposed to the high glucose environment of GDM may alter the nature of transporters involved in the uptake process (Araújo et al., 2013). The study of HOLM et al. reported that SLC1A5 was involved in sphingolipid metabolism that contributes to genetic predisposition to type 1 diabetes, whereas there has been no direct evidence of the association between SLC1A5 and GDM risk (Holm et al., 2018).

Our study has obvious merits. We have systematically analyzed the miRNAs expression and their target genes between 6 GDM and 6 healthy controls integrating with public GEO dataset. The sample size of the integrated analysis is adequate to explore and verify those candidate miRNAs and key genes. Furthermore, we investigated the biological functions and potential therapeutic small molecule compounds of miRNA-target genes, as well as obtained pivotal genes through bioinformatics analysis. More importantly, when being compared with placenta, collecting blood samples is actually more accessible in clinical practice as early screening and diagnosis.

Also, there are some limitations that need be considered. First, due to the availability of GSE103552, our study did not analyze the association between clinical data (clinical parameters and prognosis) and genetic change. Additionally, as the statistical power based on its relatively larger sample size, GSE103552 were chosen as the external validation dataset. Second, we acknowledge the relatively low sample size of internal screening dataset as limitation of our study. We integrated internal samples and the 3 diverse public databases to comprehensively screen the critical miRNAs and key genes, and this could present a comprehensive screening results and predictions. However, external experimental verifications are further needed to elucidate the molecular mechanisms of GDM in the future.

In conclusion, this study demonstrated a series of aberrantly differentially expressed miRNAs and genes that are associated with epigenetic alternations of miRNAs in GDM by overlapping DEGs and targets of DEMs. A multitude of novel miRNA-targeting genes and regulatory pathways were identified that might serve as high priority targets for therapeutic interventions (miR-29b-1-5p-TGFB2, miR-142-3p-TGFB2, miR-9-5p-FBN2, miR-212-5p-FBN2, miR-542-3p-FBN1, miR-9-5p-FBN1, miR-508-3p-FBN1, miR-493-5p-THBS1, miR-29b-3p-COL4A1, miR-432-5p-COL5A2, miR-9-5p-TGFBI, miR-486-3p-SLC7A5 and miR-6515-5p-SLC1A5). Furthermore, ten chemicals were identified as putative therapeutic agents for GDM. Precise diagnosis and therapeutic targets of GDM would be further explored through putative genes in the future.

References

• 1

  Abu SamraN.JelinekH. F.AlsafarH.AsgharF.SeoudM.HusseinS. M.et al (2022). Genomics and Epigenomics of gestational diabetes mellitus: Understanding the molecular pathways of the disease pathogenesis. Int. J. Mol. Sci.23 (7), 3514. 10.3390/ijms23073514

  • CrossRef
  • Google Scholar

• 2

  AraújoJ. R.Correia-BrancoA.RamalhoC.GonçalvesP.PinhoM. J.KeatingE.et al (2013). L-Methionine placental uptake: Characterization and modulation in gestational diabetes mellitus. Reprod. Sci.20 (12), 1492–1507. 10.1177/1933719113488442

  • CrossRef
  • Google Scholar

• 3

  BrenneckeJ.StarkA.RussellR. B.CohenS. M. (2005). Principles of microRNA-target recognition. PLoS Biol.3 (3), e85. 10.1371/journal.pbio.0030085

  • CrossRef
  • Google Scholar

• 4

  DammP.Houshmand-OeregaardA.KelstrupL.LauenborgJ.MathiesenE. R.ClausenT. D. (2016). Gestational diabetes mellitus and long-term consequences for mother and offspring: A view from Denmark. Diabetologia59 (7), 1396–1399. 10.1007/s00125-016-3985-5

  • CrossRef
  • Google Scholar

• 5

  DuT.ZamoreP. D. (2005). microPrimer: the biogenesis and function of microRNA. Dev. Camb. Engl.132 (21), 4645–4652. 10.1242/dev.02070

  • CrossRef
  • Google Scholar

• 6

  FilardiT.CatanzaroG.GriecoG. E.SplendianiE.TrocchianesiS.SantangeloC.et al (2022). Identification and validation of miR-222-3p and miR-409-3p as plasma biomarkers in gestational diabetes mellitus sharing validated target genes involved in metabolic homeostasis. Int. J. Mol. Sci.23 (8), 4276. 10.3390/ijms23084276

  • CrossRef
  • Google Scholar

• 7

  FilardiT.CatanzaroG.MardenteS.ZicariA.SantangeloC.LenziA.et al (2020). Non-coding RNA: Role in gestational diabetes pathophysiology and complications. Int. J. Mol. Sci.21 (11), 4020. 10.3390/ijms21114020

  • CrossRef
  • Google Scholar

• 8

  FrigoletM. E.ThomasG.BeardK.LuH.LiuL.FantusI. G. (2017). The bradykinin-cGMP-PKG pathway augments insulin sensitivity via upregulation of MAPK phosphatase-5 and inhibition of JNK. Am. J. Physiol. Endocrinol. Metab.313 (3), E321–E334. 10.1152/ajpendo.00298.2016

  • CrossRef
  • Google Scholar

• 9

  GaurR.YadavK. S.VermaR. K.YadavN. P.BhakuniR. S. (2014). In vivo anti-diabetic activity of derivatives of isoliquiritigenin and liquiritigenin. Phytomedicine21 (4), 415–422. 10.1016/j.phymed.2013.10.015

  • CrossRef
  • Google Scholar

• 10

  GilletV.OuelletA.StepanovY.RodosthenousR. S.CroftE. K.BrennanK.et al (2019). miRNA profiles in extracellular vesicles from serum early in pregnancies complicated by gestational diabetes mellitus. J. Clin. Endocrinol. Metab.104 (11), 5157–5169. 10.1210/jc.2018-02693

  • CrossRef
  • Google Scholar

• 11

  GiriA.EdwardsT. L.HartmannK. E.TorstensonE. S.WellonsM.SchreinerP. J.et al (2017). African genetic ancestry interacts with body mass index to modify risk for uterine fibroids. PLoS Genet.13 (7), e1006871. 10.1371/journal.pgen.1006871

  • CrossRef
  • Google Scholar

• 12

  GoodarziG.SetayeshL.FadaeiR.KhamsehM. E.AliakbariF.HosseiniJ.et al (2021). Circulating levels of asprosin and its association with insulin resistance and renal function in patients with type 2 diabetes mellitus and diabetic nephropathy. Mol. Biol. Rep.48 (7), 5443–5450. 10.1007/s11033-021-06551-2

  • CrossRef
  • Google Scholar

• 13

  GuptaS. C.KunnumakkaraA. B.AggarwalS.AggarwalB. B. (2018). Inflammation, a double-edge sword for cancer and other age-related diseases. Front. Immunol.9, 2160. 10.3389/fimmu.2018.02160

  • CrossRef
  • Google Scholar

• 14

  HanB.LuoH.RaelsonJ.HuangJ.LiY.TremblayJ.et al (2014). TGFBI (βIG-H3) is a diabetes-risk gene based on mouse and human genetic studies. Hum. Mol. Genet.23 (17), 4597–4611. 10.1093/hmg/ddu173

  • CrossRef
  • Google Scholar

• 15

  HolmL. J.KrogvoldL.HasselbyJ. P.KaurS.ClaessensL. A.RussellM. A.et al (2018). Abnormal islet sphingolipid metabolism in type 1 diabetes. Diabetologia61 (7), 1650–1661. 10.1007/s00125-018-4614-2

  • CrossRef
  • Google Scholar

• 16

  IDA (2021). IDF Diabetes atlaa. Available at: https://diabetesatlasorg/atlas/tenth-edition/[Accessed April 18 2022].

  • Google Scholar

• 17

  JiangY.ZhangM.YuD.HouG.WuJ.LiF. (2022). CircRBM33 downregulation inhibits hypoxia-induced glycolysis and promotes apoptosis of breast cancer cells via a microRNA-542-3p/HIF-1α axis. Cell. Death Discov.8 (1), 126. 10.1038/s41420-022-00860-6

  • CrossRef
  • Google Scholar

• 18

  KwakS. H.KimS. H.ChoY. M.GoM. J.ChoY. S.ChoiS. H.et al (2012). A genome-wide association study of gestational diabetes mellitus in Korean women. Diabetes61 (2), 531–541. 10.2337/db11-1034

  • CrossRef
  • Google Scholar

• 19

  LangmeadB. (2010). Aligning short sequencing reads with Bowtie. Curr. Protoc. Bioinforma.11, Unit 11.7. Chapter 11:Unit 11.7. 10.1002/0471250953.bi1107s32

  • CrossRef
  • Google Scholar

• 20

  LehtorantaL.VuolteenahoO.LaineJ.PolariL.EkholmE.RäsänenJ. (2016). Placental structural abnormalities have detrimental hemodynamic consequences in a rat model of maternal hyperglycemia. Placenta44, 54–60. 10.1016/j.placenta.2016.06.002

  • CrossRef
  • Google Scholar

• 21

  LiaoY.SmythG. K.ShiW. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinforma. Oxf. Engl.30 (7), 923–930. 10.1093/bioinformatics/btt656

  • CrossRef
  • Google Scholar

• 22

  LoeglJ.NussbaumerE.CviticS.HuppertzB.DesoyeG.HidenU. (2017). GDM alters paracrine regulation of feto-placental angiogenesis via the trophoblast. Lab. Investig.97 (4), 409–418. 10.1038/labinvest.2016.149

  • CrossRef
  • Google Scholar

• 23

  Martínez-IbarraA.Martínez-RazoL. D.Vázquez-MartínezE. R.Martínez-CruzN.Flores-RamírezR.García-GómezE.et al (2019). Unhealthy levels of phthalates and bisphenol A in Mexican pregnant women with gestational diabetes and its association to altered expression of miRNAs involved with metabolic disease. Int. J. Mol. Sci.20 (13), 3343. 10.3390/ijms20133343

  • CrossRef
  • Google Scholar

• 24

  MeersonA.NajjarA.SaadE.SbeitW.BarhoumM.AssyN. (2019). Sex differences in plasma MicroRNA biomarkers of early and complicated diabetes mellitus in Israeli arab and jewish patients. Noncoding. RNA5 (2), 32. 10.3390/ncrna5020032

  • CrossRef
  • Google Scholar

• 25

  MetzgerB. E.GabbeS. G.PerssonB.BuchananT. A.CatalanoP. A.DammP.et al (2010). International association of diabetes and pregnancy study groups recommendations on the diagnosis and classification of hyperglycemia in pregnancy. Diabetes care33 (3), 676–682. 10.2337/dc09-1848

  • CrossRef
  • Google Scholar

• 26

  PlowsJ. F.StanleyJ. L.BakerP. N.ReynoldsC. M.VickersM. H. (2018). The pathophysiology of gestational diabetes mellitus. Int. J. Mol. Sci.19 (11), 3342. 10.3390/ijms19113342

  • CrossRef
  • Google Scholar

• 27

  Prats-PuigA.OrtegaF. J.MercaderJ. M.Moreno-NavarreteJ. M.MorenoM.BonetN.et al (2013). Changes in circulating microRNAs are associated with childhood obesity. J. Clin. Endocrinol. Metab.98 (10), E1655–E1660. 10.1210/jc.2013-1496

  • CrossRef
  • Google Scholar

• 28

  SebastianiG.ValentiniM.GriecoG. E.VentrigliaG.NigiL.MancarellaF.et al (2017). MicroRNA expression profiles of human iPSCs differentiation into insulin-producing cells. Acta Diabetol.54 (3), 265–281. 10.1007/s00592-016-0955-9

  • CrossRef
  • Google Scholar

• 29

  ShenP.YangT.ChenQ.YuanH.WuP.CaiB.et al (2021). CircNEIL3 regulatory loop promotes pancreatic ductal adenocarcinoma progression via miRNA sponging and A-to-I RNA-editing. Mol. Cancer20 (1), 51. 10.1186/s12943-021-01333-7

  • CrossRef
  • Google Scholar

• 30

  ShenY.JiaY.LiY.GuX.WanG.ZhangP.et al (2020). Genetic determinants of gestational diabetes mellitus: A case-control study in two independent populations. Acta Diabetol.57 (7), 843–852. 10.1007/s00592-020-01485-w

  • CrossRef
  • Google Scholar

• 31

  SliwinskaA.KasinskaM. A.DrzewoskiJ. (2017). MicroRNAs and metabolic disorders - where are we heading?Arch. Med. Sci.13 (4), 885–896. 10.5114/aoms.2017.65229

  • CrossRef
  • Google Scholar

• 32

  TagomaA.AlnekK.KirssA.UiboR.Haller-KikkataloK. (2018). MicroRNA profiling of second trimester maternal plasma shows upregulation of miR-195-5p in patients with gestational diabetes. Gene672, 137–142. 10.1016/j.gene.2018.06.004

  • CrossRef
  • Google Scholar

• 33

  TangX.MiaoY.LuoY.SriramK.QiZ.LinF. M.et al (2020). Suppression of endothelial AGO1 promotes adipose tissue browning and improves metabolic dysfunction. Circulation142 (4), 365–379. 10.1161/CIRCULATIONAHA.119.041231

  • CrossRef
  • Google Scholar

• 34

  WangJ.PanY.DaiF.WangF.QiuH.HuangX. (2020). Serum miR-195-5p is upregulated in gestational diabetes mellitus. J. Clin. Lab. Anal.34 (8), e23325. 10.1002/jcla.23325

  • CrossRef
  • Google Scholar

• 35

  WangY.BucherM.MyattL. (2019). Use of glucose, glutamine, and fatty acids for trophoblast respiration in lean women, women with obesity, and women with gestational diabetes. J. Clin. Endocrinol. Metab.104 (9), 4178–4187. 10.1210/jc.2019-00166

  • CrossRef
  • Google Scholar

• 36

  WeiS.ChenJ.HuangY.SunQ.WangH.LiangX.et al (2020). Identification of hub genes and construction of transcriptional regulatory network for the progression of colon adenocarcinoma hub genes and TF regulatory network of colon adenocarcinoma. J. Cell. Physiol.235 (3), 2037–2048. 10.1002/jcp.29067

  • CrossRef
  • Google Scholar

• 37

  XuX. Y.TianB. M.XiaY.XiaY. L.LiX.ZhouH.et al (2020). Exosomes derived from P2X7 receptor gene-modified cells rescue inflammation-compromised periodontal ligament stem cells from dysfunction. Stem Cells Transl. Med.9 (11), 1414–1430. 10.1002/sctm.19-0418

  • CrossRef
  • Google Scholar

• 38

  YuY.HeJ. H.HuL. L.JiangL. L.FangL.YaoG. D.et al (2020). Placensin is a glucogenic hormone secreted by human placenta. EMBO Rep.21 (6), e49530. 10.15252/embr.201949530

  • CrossRef
  • Google Scholar

• 39

  ZengM.LiuJ.YangW.ZhangS.LiuF.DongZ.et al (2019). Multiple-microarray analysis for identification of hub genes involved in tubulointerstial injury in diabetic nephropathy. J. Cell. Physiol.234, 16447–16462. 10.1002/jcp.28313

  • CrossRef
  • Google Scholar

• 40

  ZhaiZ.ChenW.HuQ.WangX.ZhaoQ.TuerxunyimingM. (2020). High glucose inhibits osteogenic differentiation of bone marrow mesenchymal stem cells via regulating miR-493-5p/ZEB2 signalling. J. Biochem.167 (6), 613–621. 10.1093/jb/mvaa011

  • CrossRef
  • Google Scholar

• 41

  ZhangC.RawalS.ChongY. S. (2016). Risk factors for gestational diabetes: Is prevention possible?Diabetologia59 (7), 1385–1390. 10.1007/s00125-016-3979-3

  • CrossRef
  • Google Scholar

• 42

  ZhangL.MaJ.PanX.ZhangM.HuangW.LiuY.et al (2022). LncRNA MIR99AHG enhances adipocyte differentiation by targeting miR-29b-3p to upregulate PPARγ. Mol. Cell. Endocrinol.550, 111648. 10.1016/j.mce.2022.111648

  • CrossRef
  • Google Scholar

• 43

  ZhaoT. T.XuY. Q.HuH. M.GongH. B.ZhuH. L. (2019). Isoliquiritigenin (ISL) and its formulations: Potential antitumor agents. Curr. Med. Chem.26 (37), 6786–6796. 10.2174/0929867325666181112091700

  • CrossRef
  • Google Scholar

• 44

  ZhongL.LongY.WangS.LianR.DengL.YeZ.et al (2020). Continuous elevation of plasma asprosin in pregnant women complicated with gestational diabetes mellitus: A nested case-control study. Placenta93, 17–22. 10.1016/j.placenta.2020.02.004

  • CrossRef
  • Google Scholar

• 45

  ZhouB. F. (2002). Predictive values of body mass index and waist circumference for risk factors of certain related diseases in Chinese adults--study on optimal cut-off points of body mass index and waist circumference in Chinese adults. Biomed. Environ. Sci.15 (1), 83–96.

  • Google Scholar

• 46

  ZhouH.ChenP.DaiF.WangJ. (2021). Up-regulation of TGFBI and TGFB2 in the plasma of gestational diabetes mellitus patients and its clinical significance. Ir. J. Med. Sci.191 (5), 2029–2033. 10.1007/s11845-021-02838-2

  • CrossRef
  • Google Scholar