Comprehensive Analysis of circRNAs, miRNAs, and mRNAs Expression Profiles and ceRNA Networks in Decidua of Unexplained Recurrent Spontaneous Abortion

The diagnosis and treatment of unexplained recurrent spontaneous abortion (URSA) are subject to debate, because the exact underlying mechanisms remain unclear. To address this issue, we elucidated the expression profiles of dysregulated circRNAs, miRNAs, and mRNAs and constructed circRNA-associated competitive endogenous RNA (ceRNA) networks by comparing the decidua of URSA with that of normal early pregnancy (NEP) using RNA-sequencing. In total, 550 mRNAs, 88 miRNAs, and 139 circRNAs were differentially expressed (DE) in decidua of URSA. Functional annotation revealed that DE mRNAs as well as potential target genes of DE miRNAs and DE circRNAs are mainly involved in immunologic function, such as antigen processing and presentation, allograft rejection, and T cell receptor signaling pathway. In addition, the top hub genes, including CCL4, DDX58, CXCL10, CXCL9, MX1, CD44, RPS2, SOCS3, RPS3A, and CXCL11, were identified. The mRNAs involved in ceRNA network were enriched in complement and coagulation cascades and protein processing in the endoplasmic reticulum. We found that circRNAs in the ceRNA network, which acted as decoys for hsa-miR-204-5p, were positively correlated with MFGE8 expression. Collectively, the results demonstrated that circRNAs, miRNAs, and mRNAs were aberrantly expressed in the decidua of patients with URSA and played a potential role in the development of URSA. Thus, the establishment of the ceRNA network may profoundly affect the diagnosis and therapy of URSA in the future.

Introduction

Recurrent spontaneous abortion (RSA), which refers to three or more spontaneous abortions with a single spouse, occurs in 1–5% of women of childbearing age and exerts a negative impact on human reproductive health (Medicine, 2008). Although the etiology of RSA has not yet been completely clarified owing to its complicated molecular and cellular regulation, known causes include chromosomal translocations, uterine abnormalities, thrombophilias, endocrine defects, autoimmune diseases, and infectious agents (Srinivas et al., 2006; Cohn et al., 2010; Daher et al., 2012; Rull et al., 2012; Asadpor et al., 2013; Pluchino et al., 2014; Dean et al., 2019). However, the cause and pathogeneses of approximately 40–50% of RSA cases remain unknown (Stephenson, 1996; Weimar et al., 2012); these are termed unexplained RSA (URSA). URSA poses a major obstacle for many couples who desire to conceive a child. Recent studies have shown a correlation between noncoding RNAs (ncRNAs) and the pathogenesis of URSA (Wang et al., 2012; Dong et al., 2014).

Circular RNAs (circRNAs) a class of ncRNAs, are characterized by a covalent bond linking the 3' and 5' ends generated by back splicing. With the advent of high-throughput RNA sequencing technology, circRNAs that were considered as nonfunctional (Capel et al., 1993), were found to be involved in many diseases, such as atherosclerosis, neurological diseases, diabetes, and cancer, all of which involve processes similar to those of embryo implantation (Chen L. et al., 2015; Qu et al., 2015). In recent years, some studies have investigated the role of circRNAs in URSA. Li et al. compared the decidua tissues of the early recurrent miscarriage group with those of the control group using circRNA microarray analysis and identified differentially expressed (DE) circRNAs (Li C. et al., 2020). Another study used an RSA mouse model to demonstrate a molecular mechanism involving migration and invasion of circ-ZUFSP regulating trophoblast cells (Li Z. et al., 2020). This provides a new indicator for early diagnosis and potential treatment of RSA. However, a comprehensive understanding of the processes involved in URSA is lacking, and the potential role of circRNAs, miRNAs, and mRNAs in URSA remains unresolved.

The purpose of the current study was to gain an overview of the expression profiles of circRNAs, miRNAs, and mRNAs in URSA via RNA sequencing. In addition, because the role of circRNAs as miRNA sponges that regulate placental behavior has been confirmed (Ou et al., 2019), we constructed a circRNA-miRNA-mRNA interaction network to investigate the potential functions of these DE circRNAs. Our findings provide new evidence that may help better understand the molecular mechanisms underlying the role of circRNAs, miRNAs, and mRNAs in the pathogenesis of URSA.

Materials and Methods

Sample Collection

All participants were recruited from the Outpatient Department of Guangdong Provincial Fertility Hospital (Guangzhou, China), during the period between June 2016 and December 2019. Patients included in the current study provided written consent before surgery, and the Ethical Committee of the Guangdong Provincial Fertility Hospital approved the study protocol. The decidua tissues were obtained from 9 women (age: 27 ± 6.38 years and gestational age: 8.36 ± 0.36 weeks) via dilation and curettage following spontaneous abortion. Normal decidua tissues were obtained from 9 women (age: 26.67 ± 3.64 years and gestational age: 7.88 ± 0.8 weeks) during early pregnancy (6–12 weeks) via dilatation and curettage after the termination of pregnancy. The clinical data are presented in Table 1. Inclusion criteria for URSA were three or more consecutive spontaneous abortions and no history of successful pregnancies. Exclusion criteria for URSA were infections, endocrine or metabolic disorders, anatomic abnormalities, autoimmune diseases, paternal or maternal chromosomal abnormalities, and abnormal karyotyping of the miscarriage product. Inclusion criteria for the normal controls were 1) elective termination of normal early pregnancy, 2) at least one live birth, and 3) no history of miscarriage, preeclampsia, ectopic pregnancy, preterm delivery, and systemic diseases. After dilation and curettage, decidua tissues were immediately separated from the products of conception, washed thoroughly with sterile normal saline, and stored in liquid nitrogen for future use.

TABLE 1

TABLE 1. Clinical characteristics of the participants.

RNA Library Preparation and Sequencing

A total amount of 5 μg RNA per sample was used as input material for the RNA sample preparations. Firstly, ribosomal RNA was removed by Epicentre Ribo-zero™ rRNA Removal Kit (Epicentre, United States), and rRNA free residue was cleaned up by ethanol precipitation. The sequencing libraries were generated by NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina® (NEB, United States) following manufacturer’s recommendations. In addition, A total amount of 3 μg total RNA per sample was used as input material for the small RNA library. Sequencing libraries were generated using NEBNext ® Multiplex Small RNA Library Prep Set for Illumina ® (NEB, United States) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. The clustering of the index coded samples was performed on a cBot Cluster Generation System using TruSeq SR Cluster Kit v3 cBot HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq 2500/2000 platform and 50bp single-end reads were generated. The raw sequencing data presented in the study has been deposited in the Sequence Read Archive (SRA) database (SRA: PRJNA819210, PRJNA819201).

Protein-Protein Interaction Network Construction, Hub Gene Selection, and Module Analysis

The protein level interactions of DE mRNAs were analyzed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database (https://string-db.org/cgi/input.pl). Interactions with a combined score ≥0.4 were considered statistically significant. Cytoscape software (version 3.7) was applied to visualize the PPI network. The proteins in the central node (hub proteins) that were highly connected to a series of genes in the network were considered most closely associated with pathogenesis. The top 10 hub genes, ranked by degree of connectivity, were identified using the cytoHubba APP plug-in of Cytoscape. Significant modules in the PPI networks were determined using the Molecular Complex Detection (MCODE) APP (v1.5.1) plug-in of Cytoscape.

Real-Time Quantitative PCR (RT-qPCR) Validation

Total RNAs were extracted using a Qiagen RNeasy mini kit (Qiagen, Hilden, Germany) and reverse transcribed using SuperScript III First-Strand Synthesis system (Life Technologies, Carlsbad, CA, USA). Real-time quantitative PCR was carried out using TB Green Premix EX Taq II (Takara, Beijing, China) and analyzed using the AB Step One Plus System (Applied Biosystems, AB, United States). Convergent primers were designed for linear mRNA production while divergent primers targeting the predicted back-spliced region were designed for the circular form. The specific quantitative primers designed using oligo7 software are listed in Supplementary Table S1. Primers of β-actin (for circRNA and mRNA) and U6 (for miRNA) were designed as an endogenous control. Each experiment was performed in triplicate.

GO and KEGG Analyses and GSEA

GO, which describes genes from any organism (http://www.geneontology.org), covering the domains of Biological Process (BP), Cellular Component (CC), and Molecular Function (MF), was used. KEGG pathway analysis (http://www.genome.jp/kegg) was used to analyze the functions of mapped genes. Fisher’s exact test, the chi-square test, and the false discovery rate (FDR) were used for level of significance detection, where the p-value denoted the significance of the GO term and pathways correlated with the condition. A small FDR indicated a small error in judging the p-value. Gene Set Enrichment Analysis (GSEA) was performed using GSEA-P software to evaluate the RNA-Seq data at the gene level, with gene sets from the Molecular Signatures Database (MSigDB).

Construction of a circRNA-miRNA-mRNA Network

The expression levels of circRNAs, miRNAs, and mRNAs that showed significant differences between URSA and NEP were analyzed. We used miRanda (http://www.microrna.org/microrna/) to predict miRNA binding seed sequence sites, where overlapping of the same miRNA binding sites on both circRNAs and mRNAs indicated circRNA-miRNA-mRNA interactions. Finally, we constructed the network using Cytoscape for Windows (National Institute of General Medical Sciences) software (Supplementary Figure S1).

Western Blot Analysis

Briefly, whole protein samples were extracted, and their concentrations were measured using a Pierce BCA Protein Assay Kit (Thermo Fisher Scientific), following which 20 μg of protein was loaded and separated on a 10% sodium dodecyl sulfate (SDS) polyacrylamide gel. The membranes were incubated overnight at 4°C with specific primary antibodies MFGE8 (Abcam, dilution in 1:500) and GAPDH (Proteintech, dilution in 1:1000), followed by incubation with goat anti-rabbit IgG conjugated to horseradish peroxidase (Santa Cruz, dilution in 1:5000) for 1 h at room temperature. The membranes were then visualized via a film following incubation with ECL solution (Cell Signaling Technology).

Statistical Analysis

Statistical analyses were performed using SPSS 20.0 software. All data were expressed as mean ± SD. Statistical significance was set at p < 0.05. The Student’s t test was used to compare RT-qPCR results.

Results

Identification of DE circRNAs, miRNAs, and mRNAs

The expression of circRNAs, miRNAs, and mRNAs in the URSA and NEP control groups were analyzed using RNA sequencing. After processing the raw data, we identified 5126 circRNAs, 1540 miRNAs, and 16065 mRNAs that were aberrantly expressed in URSA. Subsequently, we used p < 0.05 as the selection criteria and determined that 139 circRNAs were DE in the decidua of RSA patients, of which 65 DE circRNAs were upregulated and 74 DE circRNAs were downregulated (Figure 1A; Table 2). Further, differences in miRNA profiles of URSA and control groups in response to changes in the expression profiles of circRNAs were investigated. Of these, 88 were DE miRNAs, including 51 upregulated DE miRNAs and 37 downregulated DE miRNAs (Figure 1B; Table 3). Besides these, 355 and 195 DE mRNAs were upregulated and downregulated, respectively (Figure 1C; Table 4). Results of the cluster analyses of the expression levels of DE circRNAs, microRNAs, and mRNAs are displayed in Figures 1D–F. Six DE RNAs (hsa_circ_0003653 and hsa_circ_0003234, hsa-miR-200c-3p, hsa-miR-141-3p, ITGVA, TP53BP2) were further verified by RT-qPCR, showing a significant difference between URSA and NEP, which was consistent with the RNA-sequencing results (Figure 1G).

FIGURE 1

FIGURE 1. The expression profiles of global circRNAs, miRNAs and mRNAs in decidua of URSA and NEP. (A–C) Volcano plots of differentially expressed circRNAs (A), miRNAs (B), and mRNAs (C). Red dots indicate upregulated and blue dots indicate downregulated RNAs. Heat map of differentially expressed circRNAs (D) miRNAs (E) and mRNAs (F) generated with RNA-seq in decidua of URSA patients and NEP. (G) The relative expression levels of six RNAs in URSA and NEP verified using RT-qPCR. URSA, unexplained recurrent spontaneous abortion; NEP, normal early pregnancy, *p < 0.05, **p < 0.01.

TABLE 2

TABLE 2. Top 10 differentially expressed circRNAs in decidua of URSA patients.

TABLE 3

TABLE 3. Top 10 differentially expressed microRNAs in decidua of URSA patients.

TABLE 4

TABLE 4. Top 10 differentially expressed mRNAs in decidua of URSA patients.

GO and KEGG Pathway Analysis of DE mRNA Combined With GSEA

GO, and KEGG analyses were conducted to predict the potential functions of the DE mRNAs. DE mRNAs in the BP group were mainly enriched in immune response, immune system process, and regulation of response to stimuli (Figure 2A). Furthermore, evaluation of RNA-seq data at the genome-wide level via GSEA confirmed and supplemented the GO results. Using BP GO annotations, GSEA showed that the activation of immune response, positive regulation of interleukin 2 production and positive regulation of CD4 positive alpha-beta T cell differentiation were upregulated in the decidua of URSA (Figure 2C). To determine the signaling pathways associated with DE-mRNAs, KEGG pathway analysis was performed. The DE mRNAs were primarily enriched in antigen processing and presentation, Toll-like receptor signaling pathway, and Natural killer cell mediated cytotoxicity (Figure 2B). GSEA results annotated via the KEGG pathway further showed that gene sets regulating antigen processing and presentation, the MAPK signaling pathway, and allograft rejection were upregulated in the decidua of URSA (Figure 2D).

FIGURE 2

FIGURE 2. GO and KEGG pathway enrichment analyses and GSEA. GO analysis (A) and KEGG (B) of dysregulated mRNAs. (C) GSEA showed DE mRNAs involved in the BP GO annotations. (D) GSEA showed DE mRNAs in KEGG.

PPI Analysis of DE mRNAs

To further analyze the functions of DE mRNAs at the protein level and also to explore the core mRNAs involved in the cellular processes of URSA, we used the STRING database to screen for functional genes and constructed an intuitive network for the annotation of PPIs (Figure 3A). The core mRNAs (known as hub genes) that corresponded to the DEGs included CCL4, DDX58, CXCL10, CXCL9, MX1, CD44, RPS2, SOCS3, RPS3A, and CXCL11 (Figure 3B). Module analysis revealed two closely connected regions with the highest scores in the PPI network (Figures 3C,D). CCL4, DDX58, and CXCL10 in the PPI network were upregulated in decidua of patients with URSA (Figure 3E).

FIGURE 3

FIGURE 3. PPI network analysis of DE mRNAs and validation of mRNA-Seq data. (A) The PPI network of DE mRNAs. (B) The subnet involving the top 10 DEGs (hub genes) with high connectivity degrees obtained from the PPI network. (C,D) The top 2 significant modules were obtained from the PPI network. (E) The relative expression levels of three mRNAs in the PPI network were validated by RT-qPCR. *p < 0.05, **p < 0.01.

Distinct circRNA and miRNA Expression Profiles in Decidua Tissues of Patients With URSA

GO analysis showed that the target genes of DE circRNAs in the biological processes mainly regulated response to DNA damage stimuli, T cell activation via T cell receptor contact with antigen bound to MHC molecules on antigen presenting cells, and the Toll-like receptor 7 signaling pathway (Figure 4A). The DE miRNAs were significantly enriched in cellular developmental processes, cellular component assembly, and developmental processes (Figure 4C).

FIGURE 4

FIGURE 4. Gene ontology and pathway analysis of the potential target genes (PTGs) of DE miRNAs and circRNAs. (A,C) GO analysis of miRNA and circRNA, respectively. (B,D) KEGG of miRNA and circRNA, respectively.

KEGG pathway enrichment analysis revealed that the DE circRNAs were mainly enriched in the Rap1 signaling pathway, the TGF-beta signaling pathway, and the tight junction and T cell receptor signaling pathway (Figure 4B). The pathways enriched in miRNA included the MAPK signaling pathway, the Ras signaling pathway, and gap junction and focal adhesion (Figure 4D).

Construction of circRNA-miRNA-mRNA ceRNA Network

To identify the ceRNA network in URSA, we extracted miRNAs, mRNAs, and circRNAs that were abnormally expressed in the triple network and constructed a network based on the ceRNA hypothesis using Cytoscape (Figures 5A,B). The associated mRNA correlated pathways were complement and coagulation cascades, protein processing in the endoplasmic reticulum, and transcriptional misregulation in cancer (Figure 5C). We selectively analyzed pairs in which the circRNAs, miRNAs, and target genes were possibly associated with URSA. A subnetwork consisting of circRNAs (hsa_circ_0004558, novel_circ_0006147, hsa_circ_0008546, etc.) included those acting as ceRNAs of hsa-miR-204-5p targeting MFGE8 (Figure 5B). Subsequently, hsa-miR-204-5p and MFGE8 were selected for RT-qPCR validation, which revealed that the expression trends and amplitudes were consistent with those of RNA-seq (Figure 5D). To build on these observations, we further investigated MFGE8 protein expression in the decidua of URSA, and the result was in concordance with those of RT-qPCR (Figure 5E). Additional results are listed in Supplementary Table S2.

FIGURE 5

FIGURE 5. Establishment of a ceRNA network and validation of the expression of has-miR-204-5p and MFGE8. (A) circRNA (down in URSA)-miRNA (up in URSA)-mRNA (down in URSA).(B) circRNA (up in URSA)-miRNA (down in URSA)-mRNA (up in URSA), RNAs marked in red are most likely to participate in the pathogenesis of URSA. (C) The KEGG pathway of mRNA involved in ceRNA. (D) The expression of has-miR-204-5p/MFGE8 in decidua of URSA and NEP by RT-qPCR. (E) The expression of MFGE8 in decidua of URSA and NEP by western blot analysis. The rectangles indicate miRNAs, ellipses represent mRNAs, and diamonds represent circRNAs, *p < 0.05, **p < 0.01.

Discussion

The current study revealed significant differences between the RNA expression profiles in the decidua tissues of URSA and NEP groups. The pathways of DE mRNAs as well as the potential target genes of DE miRNAs, and DE circRNAs were mainly involved in antigen processing and presentation, allograft rejection, the T cell receptor signaling pathway, and junction related pathways. A ceRNA network was constructed, which revealed that the pathways associated with the target mRNAs were complement and coagulation cascades and protein processing in the endoplasmic reticulum. Additionally, strict limitative requirements were applied to select possible circRNA-associated-ceRNA networks that were likely involved in the occurrence and development of URSA. Moreover, the potential mechanisms via which circRNAs binding to hsa-miR-204-5p enhanced MFGE8 expression, thereby contributing to URSA were verified.

Our results showed that mRNA regulated URSA via multiple signaling pathways, including antigen processing and presentation, the Toll-like receptor signaling pathway, and Natural killer cell mediated cytotoxicity. Classically, activated macrophages in decidua are mainly involved in antigen processing and presentation (Li et al., 2021) which has been shown to be associated with URSA (Huang et al., 2021). Toll like receptors are a family of pattern recognition receptors, which play an essential role in maintaining innate host immunity and Th1/Th2 balance (Koga and Mor, 2010; Chen et al., 2016), thereby permitting and accepting implanted and growing semi allogenic embryos. Uterine NK (uNK) cells, which form an important part of the innate endometrial immune system, make up approximately 70% of all uterine lymphocytes during pregnancy and play a critical role in establishing and maintaining pregnancy. Numerous studies have investigated the role of uNK cells in adverse pregnancy outcomes, such as recurrent miscarriage (RM) or recurrent implantation failure (RIF), where many fertility clinics have provided immunotherapy to patient groups based on abnormal levels or activity of NK cells (Tenger et al., 2005; Perricone et al., 2006; Moffett and Shreeve, 2015; Mekinian et al., 2016; Plaçais et al., 2020). Using the PPI network, we identified ten hub genes, CCL4, DDX58, CXCL10, CXCL9, MX1, CD44, RPS2, SOCS3, RPS3A, and CXCL11. According to reports, these genes, which are mainly involved in inflammatory functions (Camacho-Arroyo et al., 2021; Jiang et al., 2021), and innate immunity (Nizyaeva et al., 2019), are essential for maintaining normal pregnancy.

Previous studies have indicated that ncRNAs are involved in several stages of implantation and contribute to the dysfunction of trophoblasts in pregnancy-related disorders (He et al., 2013; Chen H. et al., 2015). A study conducted by Dong et al., revealed that miRNAs that were abnormally and significantly expressed in villi and decidua, were involved in several biological processes, including the regulation of transcription and protein amino acid phosphorylation, while KEGG pathway analysis indicated that the predicted target genes were related to adherence junctions, apoptosis, and the T cell receptor signaling pathway (Dong et al., 2014). Another study revealed that the target genes of these miRNAs were closely associated with the mitogen-activated protein kinase (MAPK) pathway, as well as with the B cell receptor, and T cell receptor signaling pathways in the URSA development process (Tang et al., 2016). Similar to these results, the results of the current study indicated that 139 circRNAs were significantly differentially expressed in human decidua, compared to NEP, and that these were primarily involved in Rap1 signaling pathways, the TGF-beta signaling pathway, and the tight junction and T cell receptor signaling pathways. Rap1 signaling pathway is implicated in several basic cellular functions, including cell-cell interactions, cell–matrix adhesion, proliferation, and regulation of cellular polarity during fetal development. Magdalena et al., highlighted the key role of Rap1 signaling in maintaining epithelial and endothelial cell junction integrity that helps maintain fetus viability (Chrzanowska-Wodnicka et al., 2015). TGF-beta, which is important in ensuring maternal support for embryo development, modulates the extent of decidua invasion (Jones et al., 2006), which finding was substantiated by our circRNA pathway analysis. Both circRNAs and miRNAs were enriched in junction related pathways, including tight junction, gap junction and focal junction, which are essential for normal pregnancy (Fukuda and Sugihara, 2008; 2012), and in T cell adhesion to extracellular matrix proteins, which have been described as determinants of successful or failing pregnancy (Chaiworapongsa et al., 2002). Tang et al., suggested that downregulation of miR-181d (Tang et al., 2016) which contains relevant binding sites for upregulated circRNAs, may contribute to the inhibition of adhesion and cell proliferation in RSA NK cells (Qian et al., 2017). But, hsa-miR-141-3p which was top 10 of differentially expressed microRNAs was low expression in the present study, while high expression in NK cells isolated from the decidua of URSA (Li and Li, 2016). A previous study in mouse showed that embryo implantation sites were significantly decreased (Liu et al., 2013), following uteri horn injection of inhibitors or mimic of miR-141. The inconsistency between the results of the present study and previous studies may be due to the existence of a balance between the differentiation and apoptosis functions performed by hsa-miR-141-3p. Considered together, our findings further demonstrated the vital role of ncRNAs in URSA development.

The functions of circRNAs include sponging miRNAs, competing with endogenous RNAs, regulating gene transcription, and interacting with RNA-binding proteins (RBP) (Hansen et al., 2013). In particular, our results confirmed that the upregulated circRNA (has_circ_0004558, novel_circ_0006147, etc.) acted as ceRNAs of hsa-miR-204-5p, targeting MFGE8 in patients with URSA. MFGE8 is a glycoprotein which has been widely studied by many engaged in physiological and pathological research, especially regarding its role in the immune system (Matsuda et al., 2011; Oba et al., 2011). Previous studies have suggested that it participates in the implantation process, because it was found upregulated during the window of implantation (Mirkin et al., 2005; Franchi et al., 2011; Schmitz et al., 2014). On the other hand, MFGE8 is also known to be involved in the inflammatory process (Komura et al., 2009), and is regulated by TNF-α in the human endometrium (Yu et al., 2014). It has been reported that overexpression of MFGE8 may impair physiological processes (Yamamoto et al., 2014; Zhao et al., 2015). In keeping with this notion, Schmitz et al.(Schmitz et al., 2017) reported that MFGE8 expression in the glandular epithelium of patients with endometriosis and infertility was significantly increased, when compared to that of healthy fertile patients. The present study also demonstrated that the expression levels of MFGE8 mRNA and protein were increased in women with URSA. In summary, these results indicated that circRNAs played a crucial role in the occurrence of URSA, and that the sponging function of miRNAs may be the main pathway involved. Our study was affected by several limitations which should be acknowledged. Firstly, the small cohort studied was a major limitation that prevented us from drawing clear conclusions from our results. Secondly, the present study was unable to substantiate direct interaction between ceRNAs.

Collectively, we found that the expression levels of circRNAs, miRNAs, and mRNAs in the URSA differed significantly from those in the NEP group. In addition, our results enabled us to propose probable pathways via which circRNAs function in URSA. Furthermore, circRNA-associated ceRNAs were constructed and MFGE8 was identified as a critical component of the circRNA-miRNA network. These novel networks may lead to the discovery of potential biomarkers or therapeutic targets in URSA.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: https://dataview.ncbi.nlm.nih.gov/object/PRJNA819201 and https://dataview.ncbi.nlm.nih.gov/object/PRJNA819210.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethical Committee of the Guangdong Provincial Fertility Hospital. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual (s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

WQ and LH designed the study and coordinated and participated in every part of the study. XL and JW participated in the design of the study and helped in the study and drafting of the manuscript. GS participated in the design and statistical analyses. XZ and HN participated in the design of the study and performed the PCR analyses. All of the authors read and approved the final manuscript.

Funding

This work was supported by the Natural Science Foundation of Guangdong Province (No. 2017A030313789, 2018A030313528, 2019A1515011984), the Science and Technology Planning Foundation of Guangzhou City (201904010017, 202102080102), Guangdong Province Medical Research Funding (No. A2021269, A2022326, A2022314, A2020614), the Family Planning Research Institute Innovation Team of Guangdong Province grants (C-03), and the Family Planning Research Institute of Guangdong Province grants (S2018010, S2018007).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.858641/full#supplementary-material

References

Asadpor, U., Totonchi, M., Sabbaghian, M., Hoseinifar, H., Akhound, M. R., Zari Moradi, S., et al. (2013). Ubiquitin-specific Protease (USP26) Gene Alterations Associated with Male Infertility and Recurrent Pregnancy Loss (RPL) in Iranian Infertile Patients. J. Assist. Reprod. Genet. 30 (7), 923–931. doi:10.1007/s10815-013-0027-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Camacho-Arroyo, I., Flores-Ramos, M., Mancilla-Herrera, I., Cruz, F. M. C., Hernández-Ruiz, J., Diaz, G. P., et al. (2021). Chemokine Profile in Women with Moderate to Severe Anxiety and Depression during Pregnancy. BMC Pregnancy Childbirth 21 (1), 807. doi:10.1186/s12884-021-04225-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Capel, B., Swain, A., Nicolis, S., Hacker, A., Walter, M., Koopman, P., et al. (1993). Circular Transcripts of the Testis-Determining Gene Sry in Adult Mouse Testis. Cell 73 (5), 1019–1030. doi:10.1016/0092-8674(93)90279-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaiworapongsa, T., Romero, R., Yoshimatsu, J., Espinoza, J., Kim, Y. M., Park, K., et al. (2002). Soluble Adhesion Molecule Profile in Normal Pregnancy and Pre-eclampsia. J. Maternal-Fetal Neonatal Med. 12 (1), 19–27. doi:10.1080/jmf.12.1.19.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H., Meng, T., Liu, X., Sun, M., Tong, C., Liu, J., et al. (2015a). Long Non-coding RNA MALAT-1 Is Downregulated in Preeclampsia and Regulates Proliferation, Apoptosis, Migration and Invasion of JEG-3 Trophoblast Cells. Int. J. Clin. Exp. Pathol. 8 (10), 12718–12727.

PubMed Abstract | Google Scholar

Chen, J.-Q., Szodoray, P., and Zeher, M. (2016). Toll-like Receptor Pathways in Autoimmune Diseases. Clin. Rev. Allerg. Immunol. 50 (1), 1–17. doi:10.1007/s12016-015-8473-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Huang, C., Wang, X., and Shan, G. (2015b). Circular RNAs in Eukaryotic Cells. Cg 16 (5), 312–318. doi:10.2174/1389202916666150707161554

PubMed Abstract | CrossRef Full Text | Google Scholar

Chrzanowska-Wodnicka, M., White, G. C., Quilliam, L. A., and Whitehead, K. J. (2015). Small GTPase Rap1 Is Essential for Mouse Development and Formation of Functional Vasculature. PLoS One 10 (12), e0145689. doi:10.1371/journal.pone.0145689

PubMed Abstract | CrossRef Full Text | Google Scholar

Cohn, D. M., Goddijn, M., Middeldorp, S., Korevaar, J. C., Dawood, F., and Farquharson, R. G. (2010). Recurrent Miscarriage and Antiphospholipid Antibodies: Prognosis of Subsequent Pregnancy. J. Thromb. Haemost. 8 (10), 2208–2213. doi:10.1111/j.1538-7836.2010.04015.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Daher, S., Mattar, R., Gueuvoghlanian-Silva, B. Y., and Torloni, M. R. (2012). Genetic Polymorphisms and Recurrent Spontaneous Abortions: an Overview of Current Knowledge. Am. J. Reprod. Immunol. 67 (4), 341–347. doi:10.1111/j.1600-0897.2012.01123.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dean, D. D., Agarwal, S., and Muthuswamy, S. (2019). Defining the Role of FMR1 Gene in Unexplained Recurrent Spontaneous Abortion. J. Assist. Reprod. Genet. 36 (11), 2245–2250. doi:10.1007/s10815-019-01591-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, F., Zhang, Y., Xia, F., Yang, Y., Xiong, S., Jin, L., et al. (2014). Genome-wide miRNA Profiling of Villus and Decidua of Recurrent Spontaneous Abortion Patients. Reproduction 148 (1), 33–41. doi:10.1530/rep-14-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

Franchi, A., Bocca, S., Anderson, S., Riggs, R., and Oehninger, S. (2011). Expression of Milk Fat Globule EGF-Factor 8 (MFG-E8) mRNA and Protein in the Human Endometrium and its Regulation by Prolactin. Mol. Hum. Reprod. 17 (6), 360–371. doi:10.1093/molehr/gaq102

PubMed Abstract | CrossRef Full Text | Google Scholar

Fukuda, M. N., and Sugihara, K. (2012). Cell Adhesion Molecules in Human Embryo Implantation. Sheng Li Xue Bao 64 (3), 247–258.

PubMed Abstract | Google Scholar

Fukuda, M. N., and Sugihara, K. (2008). An Integrated View of L-Selectin and Trophinin Function in Human Embryo Implantation. J. Obstet. Gynaecol. Res. 34 (2), 129–136. doi:10.1111/j.1447-0756.2008.00776.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., et al. (2013). Natural RNA Circles Function as Efficient microRNA Sponges. Nature 495 (7441), 384–388. doi:10.1038/nature11993

PubMed Abstract | CrossRef Full Text | Google Scholar

He, X., He, Y., Xi, B., Zheng, J., Zeng, X., Cai, Q., et al. (2013). LncRNAs Expression in Preeclampsia Placenta Reveals the Potential Role of LncRNAs Contributing to Preeclampsia Pathogenesis. PLoS One 8 (11), e81437. doi:10.1371/journal.pone.0081437

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y., Hao, J., Liao, Y., Zhou, L., Wang, K., Zou, H., et al. (2021). Transcriptome Sequencing Identified the ceRNA Network Associated with Recurrent Spontaneous Abortion. BMC Med. Genomics 14 (1), 278. doi:10.1186/s12920-021-01125-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Huang, F., Chai, X., Yuan, W., Ding, H., and Wu, X. (2021). The Role of IP-10 and its Receptor CXCR3 in Early Pregnancy. Mol. Immunol. 140, 59–69. doi:10.1016/j.molimm.2021.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, R. L., Stoikos, C., Findlay, J. K., and Salamonsen, L. A. (2006). TGF-β Superfamily Expression and Actions in the Endometrium and Placenta. Reproduction 132 (2), 217–232. doi:10.1530/rep.1.01076

PubMed Abstract | CrossRef Full Text | Google Scholar

Koga, K., and Mor, G. (2010). Toll-Like Receptors at the Maternal-Fetal Interface in Normal Pregnancy and Pregnancy Disorders. Am. J. reproductive Immunol. 63 (6), 587–600. doi:10.1111/j.1600-0897.2010.00848.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Komura, H., Miksa, M., Wu, R., Goyert, S. M., and Wang, P. (2009). Milk Fat Globule Epidermal Growth Factor-Factor VIII Is Down-Regulated in Sepsis via the Lipopolysaccharide-CD14 Pathway. J. Immunol. 182 (1), 581–587. doi:10.4049/jimmunol.182.1.581

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Chen, X., Liu, X., Liu, X., He, J., Ding, Y., et al. (2020a). CircRNA Expression Profiles in Decidual Tissue of Patients with Early Recurrent Miscarriage. Genes & Dis. 7 (3), 414–423. doi:10.1016/j.gendis.2019.06.003

CrossRef Full Text | Google Scholar

Li, D., and Li, J. (2016). Association of miR-34a-3p/5p, miR-141-3p/5p, and miR-24 in Decidual Natural Killer Cells with Unexplained Recurrent Spontaneous Abortion. Med. Sci. Monit. 22, 922–929. doi:10.12659/msm.895459

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Zheng, L., Zhao, D., Xu, Y., and Wang, Y. (2021). The Role of Immune Cells in Recurrent Spontaneous Abortion. Reprod. Sci. 28 (12), 3303–3315. doi:10.1007/s43032-021-00599-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Zhou, G., Tao, F., Cao, Y., Han, W., and Li, Q. (2020b). Circ-ZUFSP Regulates Trophoblasts Migration and Invasion through Sponging miR-203 to Regulate STOX1 Expression. Biochem. Biophys. Res. Commun. 531 (4), 472–479. doi:10.1016/j.bbrc.2020.06.117

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Gao, R., Chen, X., Zhang, H., Zheng, A., Yang, D., et al. (2013). Possible Roles of Mmu-miR-141 in the Endometrium of Mice in Early Pregnancy Following Embryo Implantation. PLoS One 8 (6), e67382. doi:10.1371/journal.pone.0067382

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuda, A., Jacob, A., Wu, R., Zhou, M., Nicastro, J. M., Coppa, G. F., et al. (2011). Milk Fat Globule-EGF Factor VIII in Sepsis and Ischemia-Reperfusion Injury. Mol. Med. 17 (1-2), 126–133. doi:10.2119/molmed.2010.00135

PubMed Abstract | CrossRef Full Text | Google Scholar

Medicine, P. C. o. t. A. S. f. R. (2008). Definitions of Infertility and Recurrent Pregnancy Loss. Fertil. Steril. 89 (6), 1603. doi:10.1016/j.fertnstert.2008.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Mekinian, A., Cohen, J., Alijotas-Reig, J., Carbillon, L., Nicaise-Roland, P., Kayem, G., et al. (2016). Unexplained Recurrent Miscarriage and Recurrent Implantation Failure: Is There a Place for Immunomodulation? Am. J. Reprod. Immunol. 76 (1), 8–28. doi:10.1111/aji.12493

PubMed Abstract | CrossRef Full Text | Google Scholar

Mirkin, S., Arslan, M., Churikov, D., Corica, A., Diaz, J. I., Williams, S., et al. (2005). In Search of Candidate Genes Critically Expressed in the Human Endometrium during the Window of Implantation. Hum. Reprod. 20 (8), 2104–2117. doi:10.1093/humrep/dei051

PubMed Abstract | CrossRef Full Text | Google Scholar

Moffett, A., and Shreeve, N. (2015). First Do No Harm: Uterine Natural Killer (NK) Cells in Assisted Reproduction. Hum. Reprod. 30 (7), 1519–1525. doi:10.1093/humrep/dev098

PubMed Abstract | CrossRef Full Text | Google Scholar

Nizyaeva, N. V., Lomova, N. A., Amiraslanov, E. Y., Kan, N. E., Nagovitsyna, M. N., and Shchegolev, A. I. (2019). Peculiarities of RIG-1 Expression in Placental Villi in Preeclampsia. Bull. Exp. Biol. Med. 167 (6), 791–794. doi:10.1007/s10517-019-04624-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Oba, J., Moroi, Y., Nakahara, T., Abe, T., Hagihara, A., and Furue, M. (2011). Expression of Milk Fat Globule Epidermal Growth Factor-VIII May Be an Indicator of Poor Prognosis in Malignant Melanoma. Br. J. Dermatol 165 (3), 506–512. doi:10.1111/j.1365-2133.2011.10409.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ou, Y., Liu, M., Zhu, L., Deng, K., Chen, M., Chen, H., et al. (2019). The Expression Profile of circRNA and its Potential Regulatory Targets in the Placentas of Severe Pre-eclampsia. Taiwan. J. Obstetrics Gynecol. 58 (6), 769–777. doi:10.1016/j.tjog.2019.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Perricone, R., Muzio, G. D., Perricone, C., Giacomelli, R., Nardo, D. D., Fontana, L., et al. (2006). High Levels of Peripheral Blood NK Cells in Women Suffering from Recurrent Spontaneous Abortion Are Reverted from High-Dose Intravenous Immunoglobulins. Am. J. Reprod. Immunol. 55 (3), 232–239. doi:10.1111/j.1600-0897.2005.00356.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Plaçais, L., Kolanska, K., Kraiem, Y. B., Cohen, J., Suner, L., Bornes, M., et al. (2020). Intralipid Therapy for Unexplained Recurrent Miscarriage and Implantation Failure: Case-Series and Literature Review. Eur. J. Obstetrics Gynecol. Reproductive Biol. 252, 100–104. doi:10.1016/j.ejogrb.2020.06.017

CrossRef Full Text | Google Scholar

Pluchino, N., Drakopoulos, P., Wenger, J. M., Petignat, P., Streuli, I., and Genazzani, A. R. (2014). Hormonal Causes of Recurrent Pregnancy Loss (RPL). Hj 13 (3), 314–322. doi:10.14310/horm.2002.1505

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, Y., Wang, X., Ruan, H., Rui, C., Mao, P., Cheng, Q., et al. (2017). Circular RNAs Expressed in Chorionic Villi Are Probably Involved in the Occurrence of Recurrent Spontaneous Abortion. Biomed. Pharmacother. 88, 1154–1162. doi:10.1016/j.biopha.2017.01.172

CrossRef Full Text | Google Scholar

Qu, S., Yang, X., Li, X., Wang, J., Gao, Y., Shang, R., et al. (2015). Circular RNA: a New Star of Noncoding RNAs. Cancer Lett. 365 (2), 141–148. doi:10.1016/j.canlet.2015.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Rull, K., Nagirnaja, L., and Laan, M. (2012). Genetics of Recurrent Miscarriage: Challenges, Current Knowledge, Future Directions. Front. Gene. 3, 34. doi:10.3389/fgene.2012.00034

CrossRef Full Text | Google Scholar

Schmitz, C. R., Oehninger, S., Genro, V. K., Chandra, N., Lattanzio, F., Yu, L., et al. (2017). Alterations in Expression of Endometrial Milk Fat Globule-EGF Factor 8 (MFG-E8) and Leukemia Inhibitory Factor (LIF) in Patients with Infertility and Endometriosis. JBRA Assist. Reprod. 21 (4), 313–320. doi:10.5935/1518-0557.20170056

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmitz, C., Yu, L., Bocca, S., Anderson, S., Cunha-Filho, J. S., Rhavi, B. S., et al. (2014). Role for the Endometrial Epithelial Protein MFG-E8 and its Receptor Integrin αvβ3 in Human Implantation: Results of an In Vitro Trophoblast Attachment Study Using Established Human Cell Lines. Fertil. Steril. 101 (3), 874–882. doi:10.1016/j.fertnstert.2013.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Srinivas, S. K., Ma, Y., Sammel, M. D., Chou, D., McGrath, C., Parry, S., et al. (2006). Placental Inflammation and Viral Infection Are Implicated in Second Trimester Pregnancy Loss. Am. J. Obstetrics Gynecol. 195 (3), 797–802. doi:10.1016/j.ajog.2006.05.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Stephenson, M. D. (1996). Frequency of Factors Associated with Habitual Abortion in 197 Couples. Fertil. Steril. 66 (1), 24–29.

PubMed Abstract | Google Scholar

Tang, L., Gao, C., Gao, L., Cui, Y., and Liu, J. (2016). Expression Profile of Micro-RNAs and Functional Annotation Analysis of Their Targets in Human Chorionic Villi from Early Recurrent Miscarriage. Gene 576 (1 Pt 2), 366–371. doi:10.1016/j.gene.2015.10.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Tenger, C., Sundborger, A., Jawien, J., and Zhou, X. (2005). IL-18 Accelerates Atherosclerosis Accompanied by Elevation of IFN-γ and CXCL16 Expression Independently of T Cells. Atvb 25 (4), 791–796. doi:10.1161/01.Atv.0000153516.02782.65

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Li, B., Wang, J., Lei, J., Liu, C., Ma, Y., et al. (2012). Evidence that miR-133a Causes Recurrent Spontaneous Abortion by Reducing HLA-G Expression. Reprod. Biomed. Online 25 (4), 415–424. doi:10.1016/j.rbmo.2012.06.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Weimar, C. H. E., Kavelaars, A., Brosens, J. J., Gellersen, B., de Vreeden-Elbertse, J. M. T., Heijnen, C. J., et al. (2012). Endometrial Stromal Cells of Women with Recurrent Miscarriage Fail to Discriminate between High- and Low-Quality Human Embryos. PLoS One 7 (7), e41424. doi:10.1371/journal.pone.0041424

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamoto, N., Yamaguchi, H., Ohmura, K., Yokoyama, T., Yoshifuji, H., Yukawa, N., et al. (2014). Serum Milk Fat Globule Epidermal Growth Factor 8 Elevation May Subdivide Systemic Lupus Erythematosus into Two Pathophysiologically Distinct Subsets. Lupus 23 (4), 386–394. doi:10.1177/0961203314523870

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, L., Anderson, S., Oehninger, S., and Bocca, S. (2014). Tumor Necrosis Factor α Up-Regulates Endometrial Milk Fat Globule-Epidermal Growth Factor 8 Protein Production via Nuclear Factor κB Activation, Resulting in Cell Migration of Epithelial Cells. Fertil. Steril. 101 (2), 552–559. doi:10.1016/j.fertnstert.2013.10.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Gong, P., Chen, Y., Nwachukwu, J. C., Srinivasan, S., Ko, C., et al. (2015). Dual Suppression of Estrogenic and Inflammatory Activities for Targeting of Endometriosis. Sci. Transl. Med. 7 (271), 271ra279. doi:10.1126/scitranslmed.3010626

CrossRef Full Text | Google Scholar