Bioinformatics Analysis of Differentially Expressed Genes, Methylated Genes, and miRNAs in Unexplained Recurrent Spontaneous Abortion

Approximately half of the unexplained recurrent spontaneous abortions remain unexplained (URSAs). We aimed to provide novel insights into the biological characteristics and related pathways of differentially expressed genes (DE-genes), DE-methylated genes, and DE-miRNAs in URSA, and construct a molecular miRNAs-mRNAs network. Four data sets (GSE22490, GSE121950, GSE73025, and GSE43256) were gained from GEO data sets. We identified the DE-genes, DE-methylated genes, and DE-miRNAs using the LIMMA package in R software. Function and enrichment analyses were conducted using DAVID. A protein-protein network was performed by STRING. We predicted the target genes of DE-miRNA using DIANA-microT-CDS. Then, we constructed miRNAs-mRNAs network. There were 137 genes that overlapped in two expression profile data sets (GSE121950 and GSE22490). We found 10 overlapping DE-methylated genes and DE-genes with opposite expression alteration trends. All those 10 genes were hypermethylated lowly expressed genes. Pathway analysis illustrated that DE-genes were enriched in osteoclast differentiation, leishmaniasis, NF-kappa B signaling pathway, Toll-like receptor signaling pathway, and tuberculosis. Based on protein-protein interaction analysis, TLR8, TLR2, CD86, TLR4, IL10, CD163, FCGR1A, CXCL8, FCGR3A, HCK, PLEK, and MNDA were identified as hub genes for DE-genes. We screened out 47 DE-miRNAs and 42 overlapping DE-genes between predicted target genes of DE-miRNAs and the 137 DE-genes. We then constructed miRNAs-mRNAs network. This study identified several genes and miRNAs involved in the development and progression of URSA, including FCGR1A, FCGR3A, CXCL8, HCK, PLEK, IL10, hsa-miR-498, and hsa-miR-4530. Although further in vivo and in vitro validations are required, our results may provide a theoretical basis for future studies.