Integrative analysis of circulating microRNAs and the placental transcriptome in recurrent pregnancy loss

Recurrent pregnancy loss (RPL) is a major type of pathological pregnancy that still lacks reliable early diagnosis and effective treatment. The placenta is critical to fetal development and pregnancy success because it participates in critical processes such as early embryo implantation, vascular remodeling, and immunological tolerance. RPL is associated with abnormalities in the biological behavior of placental villous trophoblasts, resulting in aberrant placental function. MicroRNAs (miRNAs) are increasingly being recognized as essential regulators of placental development, as well as potential biomarkers. In this study, plasma miRNAs and placental messenger RNAs (mRNAs) from RPL patients and normal pregnant (NP) controls were sequenced and analyzed. Compared to those in NP controls, 108 circulating miRNAs and 1199 placental mRNAs were differentially expressed in RPL samples. A total of 140 overlapping genes (overlapping between plasma miRNA target genes and actual placental disorder genes) were identified, and functional enrichment analysis showed that these genes were mainly related to cell proliferation, angiogenesis, and cell migration. The regulatory network among miRNAs, overlapping genes, and downstream biological processes was analyzed by protein–protein interactions and Cytoscape. Moreover, enriched mRNAs, which were predictive targets of the differentially expressed plasma miRNAs miR-766-5p, miR-1285-3p, and miR-520a-3p, were accordingly altered in the placenta. These results suggest that circulating miRNAs may be involved in the pathogenesis of RPL and are potential noninvasive biomarkers for RPL.

Introduction

Recurrent pregnancy loss (RPL) is one of the most complex and challenging scenarios in reproductive medicine and is defined as the failure of two or more pregnancies (Bolondi et al., 2016). The prevalence of RPL is difficult to estimate; however, most studies claim that 1–4% of women are affected (Rai and Regan, 2006; Rasmark Roepke et al., 2017; Magnus et al., 2019). RPL can be attributed to heredity, structure, infection, endocrine issues, immunity, or unknown causes (Rai and Regan, 2006). The placenta is the first organ to form during pregnancy and is critical for embryonic development and successful pregnancy because it participates in critical processes such as early embryo implantation, vascular remodeling, and immunological tolerance. Although the specific pathogenesis of RPL is unclear, it is commonly recognized that the placenta plays an important role in this process. The origin and development of the placenta are closely related to the proliferation, differentiation, and functional state of placental trophoblasts (Lala et al., 2021). Pathological pregnancies, such as preeclampsia (PE) and pregnancy loss, are thought to be mediated by placental trophoblast dysfunction (Huppertz, 2019; Lokeswara et al., 2021). Thus, it is crucial to understand the molecular mechanisms underlying trophoblast turnover and identify corresponding biological targets.

MiRNAs are endogenous, short noncoding single-stranded RNAs with an average length of 22 nucleotides. MiRNAs regulate gene expressions in a sequence-specific manner mainly through posttranscriptional inhibition or degradation of messenger RNAs. Approximately 600 miRNAs are expressed in the human placenta and play key roles in placental development and function by regulating the expressions of genes involved in trophoblast proliferation, differentiation, invasion, migration, apoptosis, and angiogenesis (Mouillet et al., 2015; Chen et al., 2021). A number of differentially expressed miRNAs have been confirmed in the placenta in pathological pregnancies, including RPL, PE, intrauterine growth restriction, and preterm delivery (Morales-Prieto et al., 2014; Tian et al., 2020; Tang et al., 2021). Besides their effects on intracellular silencing, placental-derived miRNAs are appealing because they can be released into maternal circulation by being packaged within extracellular vesicles (Luo et al., 2009; Morales-Prieto et al., 2020). It is possible to provide a hormone-like way for different maternal cells to communicate with each other and even between the mother and the fetus (Xu et al., 2021). Therefore, the differential expressions of miRNAs in maternal circulation can indicate potential placental pathology and can be used as biomarkers for the early detection of pregnancy complications (Zhao et al., 2013; Tian et al., 2020; Munjas et al., 2021).

Although several investigations have been performed on the expressions and functions of miRNAs in the placentas of patients with RPL, few researchers have conducted joint studies on the regulatory network of plasma miRNAs and placental mRNAs. In this study, we combined plasma miRNA sequencing with placental mRNA sequencing for the first time. The target genes of differentially expressed plasma miRNAs and placental mRNAs were cross-screened for further gene function analysis and protein–protein interaction (PPI) regulatory network construction. We focused on 10 hub genes and 7 central miRNAs. Functional analysis and literature reviews showed that they were related to pathological pregnancy by participating in cell migration, apoptosis, and angiogenesis. We aimed to elucidate the miRNA–mRNA regulatory network in RPL and provide a promising target for the diagnosis and treatment of RPL.

Materials and methods

Study population

Patients with RPL were recruited for this study at Shanghai Jiao Tong University School of Medicine International Peace Maternity and Child Health Hospital from December 2018 to December 2019. All participants were thoroughly examined to rule out any known causes of RPL. Patients with the following characteristics were excluded from the study: 1) genital malformation on pelvic examination and ultrasound, 2) abnormal immune function, 3) endocrine or metabolic disease symptoms, 4) infectious disease, 5) history of family genetic diseases, 6) pathogenic copy number variant (CNV) sequencing in chorionic villous tissue, and 6) other identified causes of miscarriage.

The normal pregnant (NP) group included normal early pregnant healthy women who underwent surgical abortion according to their wishes. We also excluded patients who could not sustain pregnancy due to their physical condition or patients whose fetuses had been identified with malformations or congenital illnesses.

The study was authorized by the Medical Ethics Committee of the International Peace Maternity and Child Health Hospital, China Welfare Institute, and all participants provided informed consent in accordance with ethical norms.

Biological samples

Placental chorionic villous tissue approximately 5 mm in size was collected during dilation and curettage. Under the dissecting microscope, chorionic villous samples were isolated from maternal decidua after being repeatedly washed with sterile saline. The clean tissue was immersed in RNAlater (ThermoFisher, AM7024) and immediately stored at −80°C.

Approximately 5 ml of venous blood was collected in an EDTA peripheral blood anticoagulant tube before the surgery. The blood was gently inverted and mixed 6–8 times after collection and centrifuged at 1600 × g for 10 min. The separated plasma was centrifuged at 12000 × g for 10 min. After two rounds of centrifugation, the plasma (approximately 2 ml) was stored at −80°C.

RNA isolation from villous samples

An mRNA isolation kit was used to isolate total RNA from placental tissues according to the manufacturer’s instructions (Ambion, Invitrogen, AM1561). The NanoDrop 2000 spectrophotometer was used to assess the purity and quantity of RNA (Thermo Scientific, United States). The Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, United States) was used to evaluate RNA integrity.

RNA isolation from maternal plasma samples

Total mRNA was extracted from 200 μl of plasma using the QIAGEN miRNeasy Serum/Plasma Advanced Kit (Qiagen, cat. # 217204) according to the manufacturer’s instructions. The NanoDrop 2000 was used for total RNA quantification (Thermo Fisher Scientific Inc., United States). The RNA integrity was determined using an Agilent 2100 Bioanalyzer (Agilent Technology, United States).

Library preparation and sequencing procedures

In the screening process, we investigated a pooled RNA sample extracted from villous and plasma samples of three RPL patients and three NP controls, matched for age and other clinical characteristics listed in Table 1.

TABLE 1

TABLE 1. Clinical characteristics of patients recruited.

The mRNA libraries of villous tissues were constructed using the TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, United States) according to the manufacturer’s instructions. The Illumina HiSeq X Ten platform was used to sequence the libraries, and 150-bp paired-end reads were produced. Each sample produced approximately 49.52 M raw readings. Trimmomatic was used to initially process the raw data (raw readings) in the fastq format, through the removal of low-quality reads to produce clean reads (Bolger et al., 2014). Then, for each sample, 48.57 M clean readings were kept for later analyses. HISAT2 was used to map the clean reads to the human genome (GRCh38) (Kim et al., 2015). The FPKM of each gene was computed using Cufflinks, and the HTSeq-count was used to determine the read counts of each gene (Trapnell et al., 2010; Roberts et al., 2011; Anders et al., 2015). All procedures were performed by OE Biotech Co., Ltd. (Shanghai, China).

Following the manufacturer’s guidelines, 1 ug of total RNA from each plasma sample was utilized for the small RNA library building using TruSeq Small RNA Sample Prep Kits (Cat. No. RS-200–0012, Illumina, United States). In a nutshell, the whole RNA was ligated to adapters at both ends. The adapter-ligated RNA was then reverse-transcribed to cDNA and PCR amplification was conducted. Small RNA libraries were created from PCR products ranging in size from 140 to 160 bp. The quality of the library was evaluated using DNA High Sensitivity Chips on the Agilent Bioanalyzer 2100 instrument. Base-calling was used to transform the raw reads from the basic reads into sequence data. The readings containing 5′ primer contamination and poly(A) were eliminated after low-quality reads were filtered. To acquire clean reads, raw data were filtered to remove reads with the 3′ adapter and insert tag as well as reads that were either 15 nt or 41 nt in length. For the primary analysis, the length distribution of the clean sequences in the reference genome was determined. Small nuclear RNAs, rRNAs, tRNAs, and other terms were used to annotate noncoding RNAs. Following the alignment of these RNAs, a Bowtie search was conducted against Rfam v.10.1 (http://www.sanger.ac.uk/software/Rfam) (Langmead et al., 2009). By matching the known miRNAs to the miRBase v22 database, the known miRNAs were located, and the known miRNA expression patterns in various samples were examined (Griffiths-Jones et al., 2008). Next, miRDeep2 was used to predict new miRNAs from unannotated readings (Friedländer et al., 2012). OE Biotech Co., Ltd., performed the small RNA sequencing and analysis (Shanghai, China).

Bioinformatic analysis

Differentially expressed mRNAs and miRNAs were used to identify the DESeq (2012) functions with estimateSizeFactors and nbinomTest. The threshold for significantly differential expression was established as a P or q < 0.05 and |log2-fold change| > 1.

Enrichment analysis

Target genes of miRNAs were predicted using the datasets collected from miRTarBase, miRWalk, and Diana Tools (Maragkakis et al., 2011; Chou et al., 2018; Sticht et al., 2018). The genes involved in overlapping between miRNA target genes and actual differentially expressed placental genes were referred to as overlapping genes. To identify the potential relevant functional annotation and pathway enrichment analysis of overlapping genes, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used. Analyses were conducted using the Database for Annotation, Visualization, and Integrated Discovery [DAVID3], with a significance level of p < 0.05 (Dennis et al., 2003).

Protein–protein interaction network

The PPI data were obtained from the STRING database (https://string-db.org/) (Szklarczyk et al., 2017). We built the interaction network with the STRING database and selected PPI pairs based on a combination score > 0.4 to assess the PPIs of overlapping genes. Then, the PPI network was analyzed by Cytoscape visualization (Shannon et al., 2003). The Cytoscape program CytoHubba was used to choose the hub genes.

Building a miRNA–mRNA regulatory network

A regulatory network of miRNA–mRNA interactions was constructed by Cytoscape to show the interaction between differentially expressed circulating miRNAs and overlapping genes.

Cumulative distribution plots

The fold change between target and nontarget genes of differentially expressed plasma miRNAs was compared using the cumulative distribution function. The Mann–Whitney U test was used to evaluate whether the log2fold change in target genes in the placentas of RPL patients was significantly changed compared with that of nontarget genes.

RT-qPCR validation analysis

To confirm the reliability of the sequencing results, crucial dysregulated villous mRNAs and plasma miRNAs were subsequently validated in all samples (20 RPL patients and 21 healthy controls). According to the manufacturer’s instructions, mRNA and miRNA were reverse-transcribed using PrimeScript RT reagent and Mir-X miRNA First-Strand Synthesis Kit (Takara, Japan). The TB Green qPCR Master Mix (Takara, Japan) was used to perform quantitative real-time PCR (qPCR) analyses according to the manufacturer’s instructions. The quantities of mRNA and miRNA transcripts were standardized to the housekeeping genes GAPDH and U6, respectively. Relative mRNA and miRNA expressions were plotted as 2−ΔΔct values. The primer sequences used for qPCR are shown in Table 2.

TABLE 2

TABLE 2. Sequences of the primers used in qRT-PCR experiments.

Statistical analysis

For statistical analysis, GraphPad Prism version 8 software (GraphPad Software, San Diego, United States) was used. Differences between groups of continuous variables were tested using a parametric two-tailed t test or Mann–Whitney U test when the normality assumption was not met. The chi-squared test was used for comparison of categorical variables. Statistical significance was accepted at p < 0.05.

Results

Clinical characteristics of the patients

In the current study, 20 RPL patients and 21 NP controls who fulfilled the inclusion criteria were included. The demographic information of the participants is shown in Table 1. There were no significant differences in maternal age, gestational age, maternal prepregnancy body mass index, smoking status, or cycle length between the two groups.

Differentially expressed circulating miRNAs and placental mRNAs

Considering a |log2-fold change| >1 and a q < 0.05, a total of 108 circulating miRNAs showed differences. A total of 77 miRNAs were upregulated and 31 were downregulated in the RPL group compared with the NP group (Figures 1A and B). Placental villous mRNA expression profiles with the same population were investigated. Considering a |log2-fold change| > 1 and p < 0.05, a total of 1199 differentially expressed transcripts were identified in placental tissue. Among them, 568 mRNAs were upregulated and 631 were downregulated (Figures 1C and D). The detailed miRNAs and mRNAs sequencing results are in Supplementary Table S1 and Supplementary Table S2.

FIGURE 1

FIGURE 1. Expression profiles of circulating miRNAs and placental mRNAs in the NP and RPL groups (n = 3 per group). (A) heatmap of circulating miRNAs. (B) volcano plot of circulating miRNAs. (C) heatmap of placental mRNAs. (D) volcano plot of placental mRNAs.

Identification of the candidate target genes of differentially expressed miRNAs

A total of 5607 target genes were identified, including 1890 targets of upregulated miRNAs and 2798 targets of downregulated miRNAs. There were overlaps between anticipated target genes and the measured mRNAs, adding to the evidence for the involved genes. We performed a joint analysis of the putative target genes of differentially expressed miRNAs and differentially expressed mRNAs. Eighty-one genes overlapped between the upregulated mRNAs and the target genes of the downregulated miRNAs were called upregulated overlapping genes. Fifty-nine genes overlapped between the downregulated mRNAs and the target genes of the upregulated miRNAs were called downregulated overlapping genes (Figure 2). A total of 140 overlapping genes were identified (Supplementary Table S3).

FIGURE 2

FIGURE 2. Identification of overlapping genes. (A) the intersection of upregulated mRNAs (red) and target genes of downregulated miRNAs (blue). (B) the intersection of downregulated mRNAs (red) and target genes of upregulated miRNAs (blue).

Functional annotation and pathway enrichment analysis

GO and KEGG pathway analyses were performed on the up/downregulated overlapping genes. GO consists of three categories: molecular cellular components, biological processes (BPs), and molecular functions (MFs). As shown in Figure 3A, BP-related results showed that downregulated overlapping genes were mainly involved in cell proliferation, angiogenesis, vasculogenesis, and cell migration, which are BPs that are closely related to abortion. In addition, the results of the KEGG pathway analysis showed that the downregulated overlapping genes were mainly enriched in the Wnt signaling pathway, focal adhesion, and axon guidance (Figure 3B). BP-related results showed that upregulated overlapping genes were particularly enriched in responses to estradiol, angiogenesis, and the cell surface receptor signaling pathway (Figure 3C). The KEGG pathway analysis results showed that the upregulated overlapping genes were mainly enriched in axon guidance (Figure 3D).

FIGURE 3

FIGURE 3. Functional analysis of overlapping genes: (A) gene ontology (GO) analysis of downregulated overlapping genes. (B) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of downregulated overlapping genes. (C) GO analysis of upregulated overlapping genes. (D) KEGG pathway analysis of upregulated overlapping genes.

Protein–protein interaction network

To explore the interaction network between target proteins and the core regulatory genes, the STRING method was used to predict the interactions of 140 overlapping genes at the protein level. In the PPI network of the overlapping genes, a total of 134 nodes and 99 edges were mapped (Figure 4A). We further identified the top 10 hub genes assessed by CytoHubba (Figure 4B). The top 10 hub genes were FYN, KDR, PECAM1, PDGFRB, S1PR1, KIT, COL1A1, EPHB2, EPHA7, and SULF1 (Table 3). The differential expressions of hub genes between the two groups are shown in Figure 5.

FIGURE 4

FIGURE 4. Protein–protein interaction (PPI) network. (A) PPI network of overlapping genes. The blue circle represents downregulated overlapping genes, and the orange circle represents upregulated overlapping genes. (B) PPI network of top 10 overlapping genes (hub genes). The darker the circle, the greater the connectivity score.

TABLE 3

TABLE 3. Top 10 overlapping genes (hub genes) with highest connectivities.

FIGURE 5

FIGURE 5. Expression changes of circulating central miRNAs and placental hub genes in RPL relative to the NP group.

miRNA–overlapping gene regulatory network establishment

Cytoscape was used to visualize the miRNA target gene regulatory network between overlapping genes and miRNAs (Figure 6A). Among the upregulated miRNAs, mir-200b-3p regulates 10 downregulated target genes and was the upregulated miRNA with the highest number of target genes. Both miR-145-5p and miR-493-5p regulate 9 downregulated overlapping genes. MiR-766-5p, miR-1285-3p, miR-4685-3p, and miR-520a-3p regulate 11, 11, 9, and 7 overlapping genes respectively and were the top four downregulated miRNAs with the largest number of target genes. We identified 7 central differentially expressed miRNAs (upregulated mir-200b-3p, miR-145-5p, and miR-493-5p; downregulated miR-766-5p, miR-1285-3p, miR-4685-3p, and miR-520a-3p) (Table 4). The differential expressions of central miRNAs in the plasma between the two groups are shown in Figure 5.

FIGURE 6

FIGURE 6. miRNA and mRNA network. (A) Construction of miRNA-overlapping genes regulatory network. The diamond represents differentially expressed circulating miRNAs, the blue circle represents downregulated overlapping genes, and the orange circle represents upregulated overlapping genes. (B) Predicted biological processes (BPs) of GO involving top 10 hub genes regulated by differentially expressed plasma miRNAs in RPL patients.

TABLE 4

TABLE 4. Seven central differentially expressed miRNAs and target overlapping genes.

Next, we established an interaction network among the differentially expressed miRNAs, hub genes, and the BPs of GO analysis (Figure 6B). Hub genes were extensively regulated by circulating miRNAs and are involved in BPs such as the ephrin receptor signaling pathway, angiogenesis, cell migration, and cell proliferation of placental villi.

Downregulation of circulating miRNAs leads to global upregulation of placental targets

We examined the cumulative distribution function of the mRNA targets of central differentially expressed miRNAs. As shown in Figure 7, the predicted targets of downregulated miRNAs (miR-766-5p, miR-1285-3p, and miR-520a-3p) were upregulated in the placenta (a righthand shift in the target mRNA compared with the nontarget mRNA indicates increased mRNA expression).

FIGURE 7

FIGURE 7. Cumulative frequency of mRNAs log2 fold change based on miRNA target status. mRNA targets of downregulated plasma miRNAs [miR-766-5p (A), miR-1285-3p (B), and miR-520a-3p (C)] are upregulated in the placental tissue.

Validation of differentially expressed miRNAs and mRNAs

We selected 4 hub genes and 4 central differentially expressed miRNAs and increased the sample size to 20 RPL patients and 21 healthy controls to validate our RNA-seq data using qRT-PCR (Figure 8). The experimental results showed that the levels of FYN and PDGFRB mRNA were elevated in the placenta, which was consistent with the RNA-seq data. In addition, the levels of KDR and PECAM1 mRNA showed the same trends as the RNA-seq results but were not significant. MiR-145-5p and miR-493-5p were significantly upregulated in the plasma of RPL patients. In addition, miR-1285-3p and miR-766-5p were downregulated in RPL patients.

FIGURE 8

FIGURE 8. Validation of placental mRNAs and circulating miRNAs in the NP and RPL groups. mRNAs (A) and miRNAs (B) were confirmed by qPCR. The expression trends of candidate mRNAs and miRNAs were consistent with the sequencing results. Data are presented as 2−ΔΔct values (mean ± SEM). *p < 0.05, ***p < 0.001, and ****p < 0.0001; ns, not significant.

Discussion

In this study, we characterized the circulating miRNA profile in the plasma and the mRNA profile in the placenta of patients with RPL. This study used a comprehensive approach to elucidate potential upstream regulators of key molecular pathways related to RPL by combining circulating miRNA and placental mRNA profiles. To our knowledge, this is the first study to integrate these two profiles from the same set of patients with RPL. In conclusion, our results show that the miRNA profile reflects the state of patients and is a potential noninvasive biomarker of RPL. Furthermore, alterations in plasma miRNA and placental mRNA expressions suggest the role of miRNAs in placental function in RPL.

In recent years, miRNAs have attracted extensive attention because of their high stability in circulating body fluids (including whole blood, plasma, serum, and saliva). Many ubiquitous and placenta-specific miRNAs are expressed and regulate trophoblast differentiation, proliferation, apoptosis, invasion/migration, and angiogenesis, suggesting that miRNAs play a significant role in placental development (Fu et al., 2013). These miRNAs can be produced by trophoblasts of the placenta and secreted into the maternal circulation. These miRNAs are potential diagnostic and prognostic biomarkers for pregnancy-related complications, as well as illness occurrence. We found 108 differentially expressed miRNAs in the plasma of RPL patients in our study, and there were 77 upregulated and 31 downregulated miRNAs. The tissues and cells from which circulating miRNAs are derived are often unknown. Therefore, we believe that the expressions of candidate diagnostic miRNAs in the circulation and placenta, as well as the evidence of important regulatory pathways involved in the pathophysiology of pregnancy loss, provide valuable indirect support for evaluating the diagnostic potential of these factors. Through combined placental RNA-seq analysis, we identified 7 central differentially expressed miRNAs that may play a broad regulatory role in the placenta. Then, potential biomarker miRNAs of RPL were further screened in larger samples (20 RPL patients and 21 NP controls) by real-time PCR. The results showed that in RPL patients, 2 miRNAs (miR-145-5p and miR-493-5p) were upregulated, and 2 miRNAs (miR-1285-3p and miR-766-5p) were downregulated. Our review of the literature revealed that these miRNAs, which had been replicated in more than one study, were expressed in the placenta (Lv et al., 2019; Inno et al., 2021; Smith et al., 2021). The four potential biomarker miRNAs identified in our study have never been linked to RPL before. The expression of miR-145-5p was increased in the placentas of patients with PE (Zhou et al., 2021). Mir-145-5p overexpression inhibited HTR-8/SVneo cell motility and invasion by negatively regulating Cyr61 expression through interactions with its 3′-untranslated region (UTR) (Wen et al., 2018). MiR-493-5p inhibits glioma progression by decreasing E2F3-mediated P53 and PI3K/AKT pathway dysfunction (Liu H. et al., 2022). Some studies have shown that mir-493-5p overexpression plays a negative regulatory role in esophageal cancer, osteosarcoma, cervical cancer, nonsmall cell lung cancer, and other malignancies by inhibiting cell proliferation/migration/invasion (Feng et al., 2021; Wang et al., 2021; Liu Q. et al., 2022; Huang et al., 2022). Downregulation of miR-1285-3p reversed the effects of weak expression of circRNA_400029 on the progression and apoptosis of cervical cancer cells (Ma et al., 2022). MiR-1285-3p overexpression promotes the development, metastasis, and invasion of a range of malignant tumors, including cervical cancer, lung adenocarcinoma, and colorectal cancer (Villanova et al., 2020; Zhang et al., 2021; Ma et al., 2022). Studies have shown that miR-766-5p can participate in the mitochondrial apoptotic pathway by targeting the 3′UTRs of BAX, BAK, and BOK, and the overexpression of miR-766-5p in SW480 cells has been shown to reduce cell death and improve viability (Dokanehiifard et al., 2020). Rapid proliferation, invasion, and migration of placental trophoblasts are required for effective embryo implantation and placental development (Brosens et al., 2022). These findings suggest that circulating miRNAs may be involved in RPL pathogenesis by targeting migration-, apoptosis-, and angiogenesis-related gene functions.

A total of 140 overlapping genes (overlapping between plasma miRNA target genes and actual placental disorder genes) were identified. This resulted in a total of 192 miRNA–mRNA interactions between 98 miRNAs and 140 target genes. In addition, using PPI analysis, we further identified the top 10 hub genes. In a consistent manner, the functional enrichment analysis of hub genes in placental villi showed that the regulated genes were mainly related to the regulation of cell proliferation, angiogenesis, cell migration, and focal adhesion, which are important processes involved in placental and fetal development. Studies have demonstrated that FYN expression, which is elevated at the fetomaternal interface of abortion-prone mice and RPL patients, plays a role in the regulation of fetomaternal immune tolerance by encouraging the expansion of Th17 cells and the expression of proinflammatory factors (Khan et al., 2014; Liu et al., 2016). A FYN inhibitor overcomes the aberrant inflammatory state and embryo resorption caused by LPS (Liu et al., 2016). KDR, also known as vascular endothelial growth factor receptor 2 (VEGFR2), was demonstrated to be elevated in the placenta of pathological pregnancy diseases (early pregnancy loss, PE) in many studies (Kumazaki et al., 2002; Plaisier et al., 2009; Fang et al., 2013; Sundrani et al., 2013). As a significant mediator of VEGF biological effects, KDR is important in uterine decidual and embryonic angiogenesis during early pregnancy (Douglas et al., 2009; Fang et al., 2013). In addition, a previous study has demonstrated that KDR participates in trophoblast migration, invasion, and proliferation by binding to decorin (Khan et al., 2011). PECAM-1 (platelet endothelial cell adhesion molecule-1), commonly known as CD31, is the most prevalent membrane glycoprotein that is expressed constitutively on the vascular endothelium (Chantraine et al., 2012; Caligiuri, 2020). In an interesting way, we found a trend toward higher levels of PECAM-1 expression in the placental villi of RPL patients. Consistent with our findings, another scientific study has shown an increased number of arterial, venous, and lymphatic vessels in placentas of RPL patients, and RPL is related to the decrease in venous and lymphatic invasion by extravillous trophoblasts rather than the total number of vessels (Windsperger et al., 2017). In an important way, our comprehensive approach showed that the downregulation of plasma miRNAs resulted in the upregulation of targeted mRNAs in placental tissue. Recent studies have identified correlations between circulating miRNAs and tissue gene expression, which lends support to our approach and findings. MiRNAs released from cells, such as in biological fluids like blood and breast milk, are thought to be present in vesicle-encapsulated, protein-bound forms, such as microvesicles, exosomes, and apoptotic bodies (Ouyang et al., 2014; Chang et al., 2017; Morales-Prieto et al., 2020). An emerging concept is that these miRNAs can silence neighboring or distant cells, allowing for hormonal-like intercellular communication associated with a variety of BP (Ouyang et al., 2014). Placental or other tissue-derived exosomal miRNAs are transferred to the maternal circulation and regulate neighboring or distant cells. Circulating exosomes released by placental syncytiotrophoblasts contain syncytin-1 and syncytin-2 on their surfaces, which can bind to MFSD2A and ASCT2 receptors on the surface of placental trophoblasts and contribute to exosome uptake (Vargas et al., 2014). In vivo and in vitro findings suggest that placental exosomal miR-15a-5p isolated from maternal plasma suppresses trophoblast proliferation, invasion, and apoptosis by downregulating CDK1 expression and impairing PI3K/AKT signaling, which is linked to the progression of PE (Wang et al., 2020a). Combined with our findings, key processes in the placenta may be affected by circulating dysregulated miRNAs in RPL.

In particular, our investigation discovered a number of previously unreported plasma biomarker miRNAs of RPL. This discovery may be due to the thorough control of the included population of pregnant women with normal fetal CNV, which is missing in other studies (Yang et al., 2018; Jairajpuri et al., 2021). From the genetic perspective, embryonic numerical and structural chromosomal abnormalities are known genetic causes of RPL, accounting for more than 50% of miscarriages (Wang et al., 2020b; Gu et al., 2022). Excluding the abortion caused by such genetic factors may more clearly identify the unknown pathogenesis of RPL. In addition, the presence of cellular contaminants such as platelets or erythrocytes and hemolysis both affect plasma miRNA detection results, and these factors may also be responsible for biases between conclusions (McDonald et al., 2011; Grasedieck et al., 2013). Therefore, we call for an effort to standardize the conclusions of different studies and accelerate their use in clinical practice by establishing a standardized process (monitoring and standardizing of sample collection, handling and storage, detection, and analysis).

Moreover, the limitations of this study must be noted. First, all RPL patients and healthy controls were Chinese. Our findings may not be applicable to patients of different ethnicities. Second, this study included a small number of samples from which miRNA characteristics were obtained. The most crucial evaluation requirements for circulating miRNA biomarkers are high sensitivity and specificity, but a single miRNA is less accurate because the levels in RPL patients and healthy controls may overlap, increasing the risk of false negative or positive diagnoses (Wang et al., 2018). Therefore, further research with more participants is needed to obtain a baseline of candidates and thus verify a panel of miRNAs as biomarkers. Third, placental villi tissue contains a variety of cells, mainly trophoblast cells but also immune cells such as monocytes, macrophages, T cells, and NK cells (Suryawanshi et al., 2018; Pique-Regi et al., 2019). This is our limitation because we chose mixed sample transcriptome sequencing, which did not allow us to deeply investigate gene expression analysis at the cellular level. In addition, understanding pathways linking miRNA binding, mRNA transcription, protein translation, and placental function may prove to be an ambitious task due to the complexity of the networks involved, which needs to be prioritized in future research to better understand the potential epigenetic mechanisms of RPL.

In conclusion, our study provides comprehensive data on circulating miRNAs and placental tissue mRNAs in RPL and shows that circulating miRNAs may become reliable biomarkers of RPL and may be involved in the pathogenesis of RPL. These results suggest a direction for the further study of RPL-related miRNAs and mRNAs and provide potential therapeutic targets for RPL.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI—PRJNA851566.

Ethics statement

The studies involving human participants were reviewed and approved by the Medical Ethics Committee of the International Peace Maternity and Child Health Hospital, China Welfare Institute. The patients/participants provided their written informed consent to participate in this study.

Author contributions

CX and SC designed the experiments. NX and XZ performed the experiments and wrote the manuscript. WS and MY assisted with the experiments. XC checked the English grammar and polished the English language in the manuscript. All the authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Key R&D Program of China (21YFC2701002 and 2018YFC1004900), the National Natural Science Foundation of China (Nos. 81971344, 81771638, and 82171677), the Shanghai Municipal Commission of Science and Technology Program (Nos. 21Y21901002 and 15411966700), the Shanghai Municipal Health Commission (202140110, GW-10.1-XK07), the Shanghai “Rising Stars of Medical Talent” Youth Development Program Clinical Laboratory Practitioners Program (201972), and the Shanghai Municipal Commission of Health and Family Planning (202140110 and 20215Y0216).

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/fphys.2022.893744/full#supplementary-material

References

Anders S., Pyl P. T., Huber W. (2015). HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics 31 (2), 166–169. doi:10.1093/bioinformatics/btu638

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolger A. M., Lohse M., Usadel B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30 (15), 2114–2120. doi:10.1093/bioinformatics/btu170

PubMed Abstract | CrossRef Full Text | Google Scholar

Bolondi G., Mocchegiani F., Montalti R., Nicolini D., Vivarelli M., De Pietri L., et al. (2016). Predictive factors of short term outcome after liver transplantation: a review. World J. Gastroenterol. 22 (26), 5936–5949. doi:10.3748/wjg.v22.i26.5936

PubMed Abstract | CrossRef Full Text | Google Scholar

Brosens J. J., Bennett P. R., Abrahams V. M., Ramhorst R., Coomarasamy A., Quenby S., et al. (2022). Maternal selection of human embryos in early gestation: insights from recurrent miscarriage. Semin. Cell Dev. Biol. doi:10.1016/j.semcdb.2022.01.007

CrossRef Full Text | Google Scholar

Caligiuri G. (2020). CD31 as a therapeutic target in atherosclerosis. Circ. Res. 126 (9), 1178–1189. doi:10.1161/circresaha.120.315935

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang G., Mouillet J. F., Mishima T., Chu T., Sadovsky E., Coyne C. B., et al. (2017). Expression and trafficking of placental microRNAs at the feto-maternal interface. FASEB J. 31 (7), 2760–2770. doi:10.1096/fj.201601146R

PubMed Abstract | CrossRef Full Text | Google Scholar

Chantraine F., Blacher S., Berndt S., Palacios-Jaraquemada J., Sarioglu N., Nisolle M., et al. (2012). Abnormal vascular architecture at the placental-maternal interface in placenta increta. Am. J. Obstet. Gynecol. 207 (3), 188. doi:10.1016/j.ajog.2012.06.083

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen X., Guo D. Y., Yin T. L., Yang J. (2021). Non-Coding RNAs regulate placental trophoblast function and participate in recurrent abortion. Front. Pharmacol. 12, 646521. doi:10.3389/fphar.2021.646521

PubMed Abstract | CrossRef Full Text | Google Scholar

Chou C. H., Shrestha S., Yang C. D., Chang N. W., Lin Y. L., Liao K. W., et al. (2018). miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 46 (D1), D296–D302. doi:10.1093/nar/gkx1067

PubMed Abstract | CrossRef Full Text | Google Scholar

Dennis G., Sherman B. T., Hosack D. A., Yang J., Gao W., Lane H. C., et al. (2003). DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 4 (5), R60. doi:10.1186/gb-2003-4-9-r60

CrossRef Full Text | Google Scholar

Dokanehiifard S., Soltani B. M., Ghiasi P., Baharvand H., Reza Ganjali M., Hosseinkhani S., et al. (2020). hsa-miR-766-5p as a new regulator of mitochondrial apoptosis pathway for discriminating of cell death from cardiac differentiation. Gene 736, 144448. doi:10.1016/j.gene.2020.144448

PubMed Abstract | CrossRef Full Text | Google Scholar

Douglas N. C., Tang H., Gomez R., Pytowski B., Hicklin D. J., Sauer C. M., et al. (2009). Vascular endothelial growth factor receptor 2 (VEGFR-2) functions to promote uterine decidual angiogenesis during early pregnancy in the mouse. Endocrinology 150 (8), 3845–3854. doi:10.1210/en.2008-1207

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang Y., Yu S., Ma Y., Sun P., Ma D., Ji C., et al. (2013). Association of Dll4/notch and HIF-1a -VEGF signaling in the angiogenesis of missed abortion. PLoS One 8 (8), e70667. doi:10.1371/journal.pone.0070667

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng N., Guo Z., Wu X., Tian Y., Li Y., Geng Y., et al. (2021). Circ_PIP5K1A regulates cisplatin resistance and malignant progression in non-small cell lung cancer cells and xenograft murine model via depending on miR-493-5p/ROCK1 axis. Respir. Res. 22 (1), 248. doi:10.1186/s12931-021-01840-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedländer M. R., Mackowiak S. D., Li N., Chen W., Rajewsky N. (2012). miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 40 (1), 37–52. doi:10.1093/nar/gkr688

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu G., Brkić J., Hayder H., Peng C. (2013). MicroRNAs in human placental development and pregnancy complications. Int. J. Mol. Sci. 14 (3), 5519–5544. doi:10.3390/ijms14035519

PubMed Abstract | CrossRef Full Text | Google Scholar

Grasedieck S., Sorrentino A., Langer C., Buske C., Döhner H., Mertens D., et al. (2013). Circulating microRNAs in hematological diseases: Principles, challenges, and perspectives. Blood 121 (25), 4977–4984. doi:10.1182/blood-2013-01-480079

PubMed Abstract | CrossRef Full Text | Google Scholar

Griffiths-Jones S., Saini H. K., van Dongen S., Enright A. J. (2008). miRBase: tools for microRNA genomics. Nucleic Acids Res. 36, D154–D158. doi:10.1093/nar/gkm952

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu C., Gao H., Li K., Dai X., Yang Z., Li R., et al. (2022). Copy number variation analysis of euploid pregnancy loss. Front. Genet. 13, 766492. doi:10.3389/fgene.2022.766492

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang Z., Wang C., Zhao X. (2022). circFIG 4 drives the carcinogenesis and metastasis of esophagus cancer via the miR-493-5p/E2F3 axis. Thorac. Cancer 13, 783–794. doi:10.1111/1759-7714.14321

PubMed Abstract | CrossRef Full Text | Google Scholar

Huppertz B. (2019). Traditional and new routes of trophoblast invasion and their implications for pregnancy diseases. Int. J. Mol. Sci. 21 (1), E289. doi:10.3390/ijms21010289

PubMed Abstract | CrossRef Full Text | Google Scholar

Inno R., Kikas T., Lillepea K., Laan M. (2021). Coordinated expressional landscape of the human placental miRNome and transcriptome. Front. Cell Dev. Biol. 9, 697947. doi:10.3389/fcell.2021.697947

PubMed Abstract | CrossRef Full Text | Google Scholar

Jairajpuri D. S., Malalla Z. H., Mahmood N., Khan F., Almawi W. Y. (2021). Differentially expressed circulating microRNAs associated with idiopathic recurrent pregnancy loss. Gene 768, 145334. doi:10.1016/j.gene.2020.145334

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan G. A., Girish G. V., Lala N., Di Guglielmo G. M., Lala P. K. (2011). Decorin is a novel VEGFR-2-binding antagonist for the human extravillous trophoblast. Mol. Endocrinol. 25 (8), 1431–1443. doi:10.1210/me.2010-0426

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan M. A., Manna S., Malhotra N., Sengupta J., Ghosh D. (2014). Expressional regulation of genes linked to immunity & programmed development in human early placental villi. Indian J. Med. Res. 139 (1), 125–140.

PubMed Abstract | Google Scholar

Kim D., Langmead B., Salzberg S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12 (4), 357–360. doi:10.1038/nmeth.3317

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumazaki K., Nakayama M., Suehara N., Wada Y. (2002). Expression of vascular endothelial growth factor, placental growth factor, and their receptors Flt-1 and KDR in human placenta under pathologic conditions. Hum. Pathol. 33 (11), 1069–1077. doi:10.1053/hupa.2002.129420

PubMed Abstract | CrossRef Full Text | Google Scholar

Lala P. K., Nandi P., Hadi A., Halari C. (2021). A crossroad between placental and tumor biology: what have we learnt? Placenta 116, 12–30. doi:10.1016/j.placenta.2021.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead B., Trapnell C., Pop M., Salzberg S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10 (3), R25. doi:10.1186/gb-2009-10-3-r25

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu Q., Tian F. J., Xie Q. Z., Zhang J., Liu L., Yang J., et al. (2016). Fyn plays a pivotal role in fetomaternal tolerance through regulation of Th17 cells. Am. J. Reprod. Immunol. 75 (5), 569–579. doi:10.1111/aji.12498

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu H., Li Z., Sun H. (2022a). MiR-493-5p inhibits the malignant development of gliomas via suppressing E2F3-mediated dysfunctions of P53 and PI3K/AKT pathways. Clin. Transl. Oncol. 24 (2), 363–370. doi:10.1007/s12094-021-02698-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu Q., Luo J., Wang H., Zhang L., Jin G. (2022b). SNHG1 functions as an oncogenic lncRNA and promotes osteosarcoma progression by up-regulating S100A6 via miR-493-5p. Acta Biochim. Biophys. Sin. 54 (1), 137–147. doi:10.3724/abbs.2021014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lokeswara A. W., Hiksas R., Irwinda R., Wibowo N. (2021). Preeclampsia: From cellular wellness to inappropriate cell death, and the roles of nutrition. Front. Cell Dev. Biol. 9, 726513. doi:10.3389/fcell.2021.726513

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo S. S., Ishibashi O., Ishikawa G., Ishikawa T., Katayama A., Mishima T., et al. (2009). Human villous trophoblasts express and secrete placenta-specific microRNAs into maternal circulation via exosomes. Biol. Reprod. 81 (4), 717–729. doi:10.1095/biolreprod.108.075481

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv Y., Lu X., Li C., Fan Y., Ji X., Long W., et al. (2019). miR-145-5p promotes trophoblast cell growth and invasion by targeting FLT1. Life Sci. 239, 117008. doi:10.1016/j.lfs.2019.117008

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma Y., Liu J., Yang Z., Chen P., Wang D. B. (2022). CircRNA_400029 promotes the aggressive behaviors of cervical cancer by regulation of miR-1285-3p/TLN1 axis. J. Cancer 13 (2), 541–553. doi:10.7150/jca.61437

PubMed Abstract | CrossRef Full Text | Google Scholar

Magnus M. C., Wilcox A. J., Morken N. H., Weinberg C. R., Håberg S. E. (2019). Role of maternal age and pregnancy history in risk of miscarriage: prospective register based study. BMJ 364, l869. doi:10.1136/bmj.l869

PubMed Abstract | CrossRef Full Text | Google Scholar

Maragkakis M., Vergoulis T., Alexiou P., Reczko M., Plomaritou K., Gousis M., et al. (2011). DIANA-microT web server upgrade supports fly and worm miRNA target prediction and bibliographic miRNA to disease association. Nucleic Acids Res. 39, W145–W148. doi:10.1093/nar/gkr294

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonald J. S., Milosevic D., Reddi H. V., Grebe S. K., Algeciras-Schimnich A. (2011). Analysis of circulating microRNA: preanalytical and analytical challenges. Clin. Chem. 57 (6), 833–840. doi:10.1373/clinchem.2010.157198

PubMed Abstract | CrossRef Full Text | Google Scholar

Morales-Prieto D. M., Ospina-Prieto S., Schmidt A., Chaiwangyen W., Markert U. R. (2014). Elsevier trophoblast research award lecture: origin, evolution and future of placenta miRNAs. Placenta 35, S39–S45. doi:10.1016/j.placenta.2013.11.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Morales-Prieto D. M., Favaro R. R., Markert U. R. (2020). Placental miRNAs in feto-maternal communication mediated by extracellular vesicles. Placenta 102, 27–33. doi:10.1016/j.placenta.2020.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouillet J. F., Ouyang Y., Coyne C. B., Sadovsky Y. (2015). MicroRNAs in placental health and disease. Am. J. Obstet. Gynecol. 213 (4), S163–S172. doi:10.1016/j.ajog.2015.05.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Munjas J., Sopić M., Stefanović A., Košir R., Ninić A., Joksić I., et al. (2021). Non-Coding RNAs in preeclampsia-molecular mechanisms and diagnostic potential. Int. J. Mol. Sci. 22 (19), 10652. doi:10.3390/ijms221910652

PubMed Abstract | CrossRef Full Text | Google Scholar

Ouyang Y., Mouillet J. F., Coyne C. B., Sadovsky Y. (2014). Review: placenta-specific microRNAs in exosomes - good things come in nano-packages. Placenta 35, S69–S73. doi:10.1016/j.placenta.2013.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Pique-Regi R., Romero R., Tarca A. L., Sendler E. D., Xu Y., Garcia-Flores V., et al. (2019). Single cell transcriptional signatures of the human placenta in term and preterm parturition. Elife 8, e52004. doi:10.7554/eLife.52004

PubMed Abstract | CrossRef Full Text | Google Scholar

Plaisier M., Dennert I., Rost E., Koolwijk P., van Hinsbergh V. W., Helmerhorst F. M., et al. (2009). Decidual vascularization and the expression of angiogenic growth factors and proteases in first trimester spontaneous abortions. Hum. Reprod. 24 (1), 185–197. doi:10.1093/humrep/den296

PubMed Abstract | CrossRef Full Text | Google Scholar

Rai R., Regan L. (2006). Recurrent miscarriage. Lancet 368 (9535), 601–611. doi:10.1016/s0140-6736(06)69204-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasmark Roepke E., Matthiesen L., Rylance R., Christiansen O. B. (2017). Is the incidence of recurrent pregnancy loss increasing? a retrospective register-based study in Sweden. Acta Obstet. Gynecol. Scand. 96 (11), 1365–1372. doi:10.1111/aogs.13210

PubMed Abstract | CrossRef Full Text | Google Scholar

Roberts A., Trapnell C., Donaghey J., Rinn J. L., Pachter L. (2011). Improving RNA-Seq expression estimates by correcting for fragment bias. Genome Biol. 12 (3), R22. doi:10.1186/gb-2011-12-3-r22

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon P., Markiel A., Ozier O., Baliga N. S., Wang J. T., Ramage D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith M. D., Pillman K., Jankovic-Karasoulos T., McAninch D., Wan Q., Bogias K. J., et al. (2021). Large-scale transcriptome-wide profiling of microRNAs in human placenta and maternal plasma at early to mid gestation. RNA Biol. 18 (1), 507–520. doi:10.1080/15476286.2021.1963105

PubMed Abstract | CrossRef Full Text | Google Scholar

Sticht C., De La Torre C., Parveen A., Gretz N. (2018). miRWalk: an online resource for prediction of microRNA binding sites. PLoS One 13 (10), e0206239. doi:10.1371/journal.pone.0206239

PubMed Abstract | CrossRef Full Text | Google Scholar

Sundrani D. P., Reddy U. S., Joshi A. A., Mehendale S. S., Chavan-Gautam P. M., Hardikar A. A., et al. (2013). Differential placental methylation and expression of VEGF, FLT-1 and KDR genes in human term and preterm preeclampsia. Clin. Epigenetics 5 (1), 6. doi:10.1186/1868-7083-5-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Suryawanshi H., Morozov P., Straus A., Sahasrabudhe N., Max K. E. A., Garzia A., et al. (2018). A single-cell survey of the human first-trimester placenta and decidua. Sci. Adv. 4 (10), eaau4788. doi:10.1126/sciadv.aau4788

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk D., Morris J. H., Cook H., Kuhn M., Wyder S., Simonovic M., et al. (2017). The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 45 (D1), D362–D368. doi:10.1093/nar/gkw937

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang M., Bai L., Wan Z., Wan S., Xiang Y., Qian Y., et al. (2021). circRNA-DURSA regulates trophoblast apoptosis via miR-760-HIST1H2BE axis in unexplained recurrent spontaneous abortion. Mol. Ther. Nucleic Acids 26, 1433–1445. doi:10.1016/j.omtn.2021.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian S., Yu J., Zhang Y., Bian Y., Ma J., Yan J., et al. (2020). Overexpression of PTEN regulated by miR-19b and miR-494 in the villous of recurrent spontaneous abortion patients. J. Reprod. Immunol. 140, 103133. doi:10.1016/j.jri.2020.103133

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell C., Williams B. A., Pertea G., Mortazavi A., Kwan G., van Baren M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28 (5), 511–515. doi:10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

Vargas A., Zhou S., Éthier-Chiasson M., Flipo D., Lafond J., Gilbert C., et al. (2014). Syncytin proteins incorporated in placenta exosomes are important for cell uptake and show variation in abundance in serum exosomes from patients with preeclampsia. FASEB J. 28 (8), 3703–3719. doi:10.1096/fj.13-239053

PubMed Abstract | CrossRef Full Text | Google Scholar

Villanova L., Barbini C., Piccolo C., Boe A., De Maria R., Fiori M. E., et al. (2020). miR-1285-3p controls colorectal cancer proliferation and escape from apoptosis through DAPK2. Int. J. Mol. Sci. 21 (7), E2423. doi:10.3390/ijms21072423

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang H., Peng R., Wang J., Qin Z., Xue L. (2018). Circulating microRNAs as potential cancer biomarkers: the advantage and disadvantage. Clin. Epigenetics 10, 59. doi:10.1186/s13148-018-0492-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang Y., Du X., Wang J. (2020a). Transfer of miR-15a-5p by placental exosomes promotes pre-eclampsia progression by regulating PI3K/AKT signaling pathway via CDK1. Mol. Immunol. 128, 277–286. doi:10.1016/j.molimm.2020.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang Y., Li Y., Chen Y., Zhou R., Sang Z., Meng L., et al. (2020b). Systematic analysis of copy-number variations associated with early pregnancy loss. Ultrasound Obstet. Gynecol. 55 (1), 96–104. doi:10.1002/uog.20412

CrossRef Full Text | Google Scholar

Wang J., Liu Y., Cai H., Jiang H., Li W., Shi Y., et al. (2021). Long coding RNA CCAT2 enhances the proliferation and epithelial-mesenchymal transition of cervical carcinoma cells via the microRNA-493-5p/CREB1 axis. Bioengineered 12 (1), 6264–6274. doi:10.1080/21655979.2021.1969834

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen Z., Chen Y., Long Y., Yu J., Li M. (2018). Tumor necrosis factor-alpha suppresses the invasion of HTR-8/SVneo trophoblast cells through microRNA-145-5p-mediated downregulation of Cyr61. Life Sci. 209, 132–139. doi:10.1016/j.lfs.2018.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Windsperger K., Dekan S., Pils S., Golletz C., Kunihs V., Fiala C., et al. (2017). Extravillous trophoblast invasion of venous as well as lymphatic vessels is altered in idiopathic, recurrent, spontaneous abortions. Hum. Reprod. 32 (6), 1208–1217. doi:10.1093/humrep/dex058

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu P., Ma Y., Wu H., Wang Y. L. (2021). Placenta-derived MicroRNAs in the pathophysiology of human pregnancy. Front. Cell Dev. Biol. 9, 646326. doi:10.3389/fcell.2021.646326

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang Q., Gu W. W., Gu Y., Yan N. N., Mao Y. Y., Zhen X. X., et al. (2018). Association of the peripheral blood levels of circulating microRNAs with both recurrent miscarriage and the outcomes of embryo transfer in an in vitro fertilization process. J. Transl. Med. 16 (1), 186. doi:10.1186/s12967-018-1556-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang R. N., Wu D. M., Wu L. P., Gao G. W. (2021). LncRNA LINC00337 sponges mir-1285-3p to promote proliferation and metastasis of lung adenocarcinoma cells by upregulating YTHDF1. Cancer Cell Int. 21 (1), 550. doi:10.1186/s12935-021-02253-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao Z., Moley K. H., Gronowski A. M. (2013). Diagnostic potential for miRNAs as biomarkers for pregnancy-specific diseases. Clin. Biochem. 46 (10-11), 953–960. doi:10.1016/j.clinbiochem.2013.01.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou Q., Li H., Zhang Y., Peng W., Hou H., Gu M., et al. (2021). MicroRNA-513c-5p is involved in the pathogenesis of preeclampsia by regulating of low-density lipoprotein receptor-associated protein 6. BMC Pregnancy Childbirth 21 (1), 837. doi:10.1186/s12884-021-04069-w

PubMed Abstract | CrossRef Full Text | Google Scholar