Transcriptome sequencing identified the ceRNA network associated with recurrent spontaneous abortion
BMC Medical Genomics volume 14, Article number: 278 (2021)
Recurrent spontaneous abortion (RSA) is one of the common complication of pregnancy, bringing heavy burden to the patients and their families. The study aimed to explore the lncRNA-miRNA-mRNA network associated with recurrent spontaneous abortion.
By transcriptome sequencing, we detected differences in lncRNA, miRNA and mRNA expression in villus tissue samples collected from 3 patients with RSA and 3 normal abortion patients. Differentially expressed lncRNAs, miRNAs and genes (DELs, DEMs and DEGs, respectively) were identified, and Geno Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to determine the functions of DELs and DEGs, which were analysed by Fisher’s test. We also observed the regulatory relationships between miRNA-mRNA and lncRNA-miRNA by Cytoscape 3.6.1.
The results showed that 1008 DELs (523 upregulated and 485 downregulated), 475 DEGs (201 upregulated and 274 downregulated) and 37 DEMs (15 upregulated and 22 downregulated) were identified. And we also constructed a novel lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 genes (NTNG2, GRIA1 and AQP1).
lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 mRNAs (NTNG2, GRIA1 and AQP1) was constructed. The results may provide a basic theory for elucidating the mechanism underlying RSA.
Recurrent spontaneous abortion (RSA), a complication of pregnancy, is defined as three or more consecutive spontaneous abortions with the same spouse, and gestational age < 20 weeks of spontaneous abortion . The incidence rate of that was about 5%. To date, chromosomal abnormalities of parents and embryos, anatomical factors, thrombosis, immunological factors, endocrine factors and environmental factors have been reported to be correlated with the occurrence of RSA [2,3,4,5]. RSA seriously endangers women's reproductive health and causes great physical and mental pain for patients and their families. In recent years, new risk factors have been gradually recognized, but the aetiology of many patients with RSA remains unknown. At present, the commonly used treatment methods include immunotherapy, endocrine therapy and anticoagulant therapy. However, the effects of these treatments are not satisfactory. Therefore, it is of great significance to study the aetiology of RSA.
Noncoding RNAs play an important role in the almost all pathological or pathological processes, such as embryonic development, cell proliferation, differentiation, apoptosis, infection and immune response, including RSA [6, 7]. Long noncoding RNAs (lncRNAs), highly conserved noncoding RNAs, have also been found to be involved in RSA related studies [8,9,10]. Gu et al. observed  that polymorphisms in lncRNA HULC may be related to the susceptibility to RSA in the Southern Chinese population. Xuan et al. also found that the lncRNA MALAT1 rs619586 G mutation reduced the risk of RSA . Che et al. found that lncRNA CCAT2 rs619586 G mutation may have a potential protective effect and reduce the risk of RSA in southern China . The results described above indicated that lncRNAs played a role in RSA. Furthermore, miRNAs have also been found to be indispensable for the pathogenesis of RSA [11, 12]. By assessing the influence of USP25 on trophoblasts, Ding et al. found that USP25 expression was negatively regulated by miR-27a-3p, and this effect contributed to the pathogenesis of RSA by suppressing the migration and invasion of trophoblasts . It was observed that the upregulation of miR-365 expression may promote the occurrence of RSA by reducing the expression of SGK1, suggesting that miR-365 may be used as a prognostic biomarker and therapeutic target for RSA reported by Zhao et al. . As a new model of gene expression regulation, the large regulatory network of ceRNAs is helpful for exploring the gene function and regulatory mechanisms at a deeper level and for more thoroughly and comprehensively understanding many biological phenomena. However, so far, lncRNA-miRNA interactions and lncRNA-miRNA-mRNA networks have not been reported in RSA.
In the study, we constructed a lncRNA-associated ceRNA network to explore the pathogenesis of RSA in 3 patients with RSA and 3 normal abortion patients, providing a theoretical basis for the elucidation and treatment of RSA in the future.
Material and methods
The villus tissue samples were collected from 3 patients with RSA and villus tissue samples from 3 normal abortion patients were served as controls. The fresh tissues were stored in liquid nitrogen tanks for subsequent use. The inclusion criteria for the RSA patients were as follows: (1) patients with RSA suffered three or more consecutive spontaneous abortions at a gestational age of < 20 weeks; (2) female patients with RSA who suffered primary abortion and had no previous history of live births; (3) RSA patients underwent routine examinations, including examination of maternal infection, chromosome aberration, endocrine dysfunction, anatomical factors and autoimmune diseases. Patients who did not meet these conditions were excluded. The controls had at least one childbirth and had no history of spontaneous abortion. Moreover, the controls had no pregnancy-related complications. All the subjects have signed an informed consent form. The study was approved by the Second Affiliated Hospital of Hainan Medical College (2018R005-F01).
Transcriptome sequencing data analysis
Using the TRIzol Reagent (Thermo Fisher Science, USA), we extracted total RNA from the villus tissue samples. Subsequently, we measured the RNA concentration and purity. We performed lncRNA, miRNA and mRNA sequencing with the Illumina transcriptome chip. FastQC software and the R package (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) were used to evaluate the quality of the original sequencing data. Using the Trimgalore method (http://www.bioinformatics.babraham.ac.uk/projects/trim_Galore/), we filtered raw reads to obtain clean reads for subsequent analysis. Besides, all the data were processed by quantile normalization.
Analysis of differentially expressed lncRNAs, miRNAs and genes
We used the Cuffdiff version 2.2.1 to identify differentially expressed lncRNAs, miRNAs and genes (DELs, DEMs and DEGs) in the villus tissue samples collected from 3 patients with RSA and 3 controls. p < 0:05 and |log2FC| > 1 were used as the screening criteria. We completed the heatmap analysis of DELs, DEMs and DEGs with the ComplexHeatmap in the R package.
In the present studie, Geno Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to determine the functions of DELs and DEGs, which were analysed by Fisher’s test using the clusterProfiler version 2.2.1. Biological process (BP), cell composition (CC) and molecular function (MF) were included in GO annotation analysis. KEGG enrichment analysis mainly focused on the related signaling pathways. p < 0.05 was regarded statistically significant.
Constructing the ceRNA network
We also observed the regulatory relationships between miRNA-mRNA and lncRNA-miRNA by Cytoscape version 3.6.1. The miRNA-mRNA-lncRNA network was constructed and visualized. According to the lncRNAs that directly interacted with mRNAs and regulated their activity as miRNA sponge, we analyze the data through the following three steps: (1) miRNAs targeted by DELs and the interaction between DELs and miRNAs were predicted using the LncTar software; (2) mRNAs targeted by DEMs and the interaction between DELs and miRNAs were predicted by the online tool (MiRDB, miRTarBase and Targetscan databases; (3) lncRNAs and miRNAs negatively regulated by mRNAs were integrated to construct a ceRNA network.
Identification of differentially expressed lncRNAs, miRNAs and mRNAs
According to the quality of the total RNA from 3 patients with RSA and 3 normal abortion personnel (Additional file 1: Figure S1), the transcriptome sequencing showed that 1008 DELs (523 upregulated and 485 downregulated), 475 DEGs (201 upregulated and 274 downregulated) and 37 DEMs (15 upregulated and 22 downregulated), which were shown in the heat map (Fig. 1) and the volcano map (Additional file 3: Figure S3, Additional file 4: Figure S4 and Additional file 5: Figure S5). Table 1 illustrated the DEGs and DELs (top 20), and Table 2 listed the top 15 DEMs. The thresholds of the screening data were p < 0.05 and |log2FC|> 1.
GO and pathway analysis of DELs
To further study the transcriptome differences between the two groups, we performed GO and KEGG pathway analyses of DELs. In the Table 3 and Fig. 2a, the results of the top 10 enriched GO pathways of DELs showed that the biological process (BP) changes were in the regulation of body fluid levels, embryonic skeletal system development, postsynapse organization, carbohydrate derivative transport, activation of JUN kinase activity, mammary gland epithelial cell proliferation, oxygen transport, gas transport, regulation of mammary gland epithelial cell proliferation and pericardium development. Additionally, the cell component (CC) changes of DELs were obviously enriched in transcription factor complex, axon part, postsynaptic specialization, histone methyltransferase complex, clathrin-coated pit, MLL1/2 complex, hemoglobin complex, MLL1 complex, haptoglobin-hemoglobin complex and exocyst. Moreover, molecular function (MF) changes were mainly enriched in DNA-binding transcription activator activity, RNA polymerase II-specific, enhancer sequence-specific DNA binding, enhancer binding, RNA polymerase II distal enhancer sequence-specific DNA binding, oxidoreductase activity, acting on NAD(P)H, molecular carrier activity, kinesin binding, laminin binding, haptoglobin binding and oxygen carrier activity. As shown in the Table 4 and Fig. 2b, the top 10 enriched KEGG pathways of DELs were in Alzheimer’s disease, Thermogenesis, Thyroid hormone signaling pathway, Hippo signaling pathway, Hepatocellular carcinoma, Adherens junction, Arrhythmogenic right ventricular cardiomyopathy (ARVC), Vibrio cholerae infection, Glycosphingolipid biosynthesis—lacto and neolacto series and Antifolate resistance.
GO and pathway analyses of DEGs
To further study the transcriptome differences between the two groups, we performed the GO and KEGG pathway analysis of DEGs. The results of the top 10 GO pathways of DEGs showed that changes in biological processes (BP) were mainly enriched in regulation of metal ion transport, leukocyte cell–cell adhesion, regulation of leukocyte proliferation, antigen processing and presentation, antibiotic catabolic process, gas transport, cellular extravasation, hydrogen peroxide catabolic process, oxygen transport and eosinophil migration. In addition, cell component (CC) changes were mainly concentrated in extracellular matrix, actin cytoskeleton, contractile fiber, contractile fiber part, myofibril, sarcomere, hemoglobin complex, haptoglobin-hemoglobin complex, MHC protein complex and MHC class II protein complex. Molecular function (MF) changes were mainly distributed in actin binding, actin filament binding, organic acid binding, molecular carrier activity, antioxidant activity, oxygen binding, peroxidase activity, oxidoreductase activity, acting on peroxide as acceptor, haptoglobin binding and oxygen carrier activity (Table 5 and Fig. 3a). The top 10 KEGG pathways of DEGs were mainly enriched in Cell adhesion molecules (CAMs), Chemokine signaling pathway, Staphylococcus aureus infection, Viral protein interaction with cytokine and cytokine receptor, Malaria, B cell receptor signaling pathway, Leishmaniasis, Asthma, African trypanosomiasis and Allograft rejection (Table 6 and Fig. 3b).
Construction of the lncRNA-miRNA-mRNA ceRNA network
First, we constructed a lncRNA-miRNA network and miRNA-mRNA network. The lncRNA-miRNA network included 4607 negative interactions (1945 downregulated lncRNAs-upregulated miRNAs and 2662 upregulated lncRNAs- downregulated miRNAs), and the miRNA-mRNA network included 15 negative interactions (6 downregulated miRNAs- upregulated mRNAs and 9 upregulated miRNAs-downregulated mRNAs). Then, we constructed the lncRNA-miRNA-mRNA ceRNA network to identify their relationships based on the lncRNA, miRNA, and mRNA expression profiles, and plotted them using Cytoscape version 3.6.1. First, based on the threshold values (r < 0 and p value < 0.05), we evaluated the relationship between downregulated lncRNAs and upregulated miRNAs shown in Fig. 4a, and the relationship between upregulated lncRNAs and downregulated miRNAs was displayed in Additional file 2: Figure S2. Additionally, the results of the miRNA-mRNA relationship showed a significant link between hsa-miR-210-5p and mRNAs (NTNG2, GRIA1 and AQP1), as shown in Fig. 4b. Besides, we constructed the ceRNA network between DELs, DEMs and DEGs by the Pearson correlation coefficient. Finally, the ceRNA network contained 31 lncRNAs (PSD2-AS1, lnc-ACAN-2, lnc-STON1-1, lnc-HPS4-8, lnc-SHC2-1, lnc-LMO7DN-6, lnc-TPTE-12, lnc-ARRDC3-5, lnc-CHPF-4, lnc-OR1J1-2, lnc-GPAT4-1, lnc-ARPC5L-1, LYPLAL1-DT, lnc-PIWIL4-1, lnc-CCR8-2, lnc-RHBDD3-3, lnc-PPP1R3G-9, RAMP2-AS1, LINC01771, lnc-SFRP4-3, lnc-C1QL3-2, lnc-C6orf223-1, lnc-IGFBP3-2, lnc-CUL2-3, lnc-SRGAP2C-5, PRKCQ-AS1, lnc-C11orf95-5, lnc-IGFBP1-1, lnc-CYP3A7-1, lnc-GPC6-7 and lnc-FUBP1-3), 1 miRNA (hsa-miR-210-5p) and 3 mRNAs (NTNG2, GRIA1 and AQP1) as displayed in Fig. 5 and Table 7 (top20), illustrating that these molecules may be involved in the development of RSA.
Recurrent spontaneous abortion is one of the common complications of pregnancy. In the past few decades, the disease has caused heavy psychological burden for couples who want to have children and their families. However, due to the high cost of treatment, many families have failed to realize their desire to have children. The present study showed that we found 1008 DELs, 475 DEGs and 37 DEMs in 3 patients with RSA and 3 normal abortion personnel by transcriptome sequencing of villous tissue samples. We also constructed a novel lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 mRNAs (NTNG2, GRIA1 and AQP1). The results may provide a theoretical basis for elucidating the mechanism of RSA.
NTNG2 (Netrin G2) the position of which on chromosome is 9q34.13 and encodes the protein NTNG2, a membrane anchor protein. It was found to promote the growth of axons and dendrites. Studies on the correlation of gene polymorphisms in schizophrenia revealed that NTNG1 and its paralogues for NTNG2 gene may be related to the pathophysiology of schizophrenia [13,14,15]. Another paper reported by Maroofian et al. illustrated that NTNG2 played a key role in neurotypical development . Therefore, we speculate that NTNG2 and NTNG1 may play a role in neurological disorders. In addition, based on the bioinformatics analysis of the pediatric onset of multiple sclerosis, genes such as NTNG2 were found to be nodes of the network, and the expression of some miRNAs were significantly correlated with brain volume . But to date, there is no report on the NTNG2-associated network in RSA.
GRIA1 (Glutamate Ionotropic Receptor AMPA Type Subunit 1) is located on the chromosome 5q33.2 and the encoded protein is the main excitatory neurotransmitter receptor in the mammalian brain. It is reported to be a ligand-gated ion channel, regulating the secretion of follicle-stimulating hormone and luteinizing hormone by controlling gonadotropin releasing hormone. Recently, Sugimoto et al. discovered that the gene was linked to the ovulation rate in cattle . Cushman et al. found the correlation between GRIA1 SNPs and cattle infertility . In addition, Sheikhha et al. also observed the relationship between GRIA1 variants and ovarian response to human menopausal gonadotropin in the group of Iranian women . The above studies show that GRIA1 plays an important role in diseases related to pregnancy in women. But so far, the role of GRIA1 in RSA has not been reported. In our study, we used Cytoscape software to construct a network to combine noncoding RNAs to explore its function in RSA.
AQP1 (Aquaporin 1), located on chromosome 7p14.3, contains 4 exons. Some reports have revealed that AQP1 plays an important role in acute lung injury caused by endotoxic shock, delaying the occurrence of renal cyst, and acute lung and brain injury [21, 22]. Su et al. used lipopolysaccharide (LPS)-induced murine model of acute lung injury to detect the function of AQP1, suggesting that AQP1 may be involved in the progression of acute lung injury . Also, noncoding RNAs interacting with AQP1, were involved in the development of acute lung injury. Long noncoding RNA CASC2 can reduce the apoptosis of lung epithelial cells and improve acute lung injury by regulating the miR-144-3p/AQP1 axis . Recent studies have shown that AQP1 participated in the occurrence of diseases through the ceRNA network [25, 26]. Tang et al. observed that lncRNA CASC2 acted as miR-144-3p, and directly interacted with AQP1 after LPS induced A549 cells . After lipopolysaccharide (LPS) induced sepsis, Fang et al. found that AQP1 has been reported to competitively bind to lncRNA H19 and regulated the expression of miRNA-874 . But so far, there has been no study on lncRNA-miRNA-AQP1 in RSA.
At present, some researchers have obtained some results about recurrent abortion through transcriptome sequencing [27, 28]. In this study, we firstly performed transcriptome sequencing analysis on the tissues of 3 patients with RSA and 3 patients with normal abortion, and found key molecules by constructing lncRNA-related ceRNA network, which is helpful to explore the pathogenic mechanism of RSA. However, there are some limitations: (1) the sample size was insufficient; (2) the lncRNAs-miRNAs linked with RSA were not verified; (3) the lncRNA-mediated ceRNA network in RSA was not verified. In the future, we will continue to collect a large number of samples for verification, and further analyze the ceRNA network in RSA by transcriptome analyses, and use molecular biology to verify this network, providing a theoretical basis for the elucidation and treatment of RSA.
In summary, a lncRNA-related ceRNA network containing 31 lncRNAs, 1 miRNA (hsa-miR-210-5p) and 3 mRNAs (NTNG2, GRIA1 and AQP1) was constructed. The results may provide the basic theory for elucidating the mechanism underlying RSA.
Availability of data and materials
The datasets generated and/or analysed during the current study are deposited in Genome Sequence Archive (GSA) for Human under analysis accession number "HRA001475" and "HRA001505".
Recurrent spontaneous abortion
Differentially expressed lncRNAs
Differentially expressed miRNAs
Differentially expressed genes
Long non coding RNAs
Kyoto encyclopedia of genes and genomes
Cell adhesion molecules
Glutamate Ionotropic Receptor AMPA Type Subunit 1
Cohn DM, Goddijn M, Middeldorp S, Korevaar JC, Dawood F, et al. Recurrent miscarriage and antiphospholipid antibodies: prognosis of subsequent pregnancy. J Thrombosis Haemostasis: JTH. 2010;8(10):2208–13.
Elkarhat Z, Kindil Z, Zarouf L, Razoki L, Aboulfaraj J, et al. Chromosomal abnormalities in couples with recurrent spontaneous miscarriage: a 21-year retrospective study, a report of a novel insertion, and a literature review. J Assist Reprod Genet. 2019;36(3):499–507.
Kaur R, Gupta K. Endocrine dysfunction and recurrent spontaneous abortion: an overview. Int J Appl Basic Med Res. 2016;6(2):79–83.
Hooker AB, Lemmers M, Thurkow AL, Heymans MW, Opmeer BC, et al. Systematic review and meta-analysis of intrauterine adhesions after miscarriage: prevalence, risk factors and long-term reproductive outcome. Hum Reprod Update. 2014;20(2):262–78.
Muyayalo KP, Li ZH, Mor G, Liao AH. Modulatory effect of intravenous immunoglobulin on Th17/Treg cell balance in women with unexplained recurrent spontaneous abortion. Am J Reprod Immunol (New York, NY: 1989). 2018;80(4):e13018.
Lou C, Goodier JL, Qiang R. A potential new mechanism for pregnancy loss: considering the role of LINE-1 retrotransposons in early spontaneous miscarriage. Reprod Biol Endocrinol: RB&E. 2020;18(1):6.
Jia N, Li J. Noncoding RNAs in unexplained recurrent spontaneous abortions and their diagnostic potential. Dis Mark. 2019;2019:7090767.
Fang Z, Yang Y, Xu Y, Mai H, Zheng W, et al. LncRNA HULC polymorphism is associated with recurrent spontaneous abortion susceptibility in the Southern Chinese population. Front Genet. 2019;10:918.
Che D, Yang Y, Xu Y, Fang Z, Pi L, et al. The lncRNA MALAT1 rs619586 G variant confers decreased susceptibility to recurrent miscarriage. Front Physiol. 2019;10:385.
Che D, Huang W, Fang Z, Li L, Wu H, et al. The lncRNA CCAT2 rs6983267 G allele is associated with decreased susceptibility to recurrent miscarriage. J Cell Physiol. 2019;234(11):20577–83.
Ding J, Cheng Y, Zhang Y, Liao S, Yin T, et al. The miR-27a-3p/USP25 axis participates in the pathogenesis of recurrent miscarriage by inhibiting trophoblast migration and invasion. J Cell Physiol. 2019;234(11):19951–63.
Zhao W, Shen WW, Cao XM, Ding WY, Yan LP, et al. Novel mechanism of miRNA-365-regulated trophoblast apoptosis in recurrent miscarriage. J Cell Mol Med. 2017;21(10):2412–25.
Aoki-Suzuki M, Yamada K, Meerabux J, Iwayama-Shigeno Y, Ohba H, et al. A family-based association study and gene expression analyses of netrin-G1 and -G2 genes in schizophrenia. Biol Psychiat. 2005;57(4):382–93.
Eastwood SL, Harrison PJ. Decreased mRNA expression of netrin-G1 and netrin-G2 in the temporal lobe in schizophrenia and bipolar disorder. Neuropsychopharmacology. 2008;33(4):933–45.
Zhu Y, Yang H, Bi Y, Zhang Y, Zhen C, et al. Positive association between NTNG1 and schizophrenia in Chinese Han population. J Genet. 2011;90(3):499–502.
Dias CM, Punetha J, Zheng C, Mazaheri N, Rad A, et al. Homozygous missense variants in NTNG2, encoding a presynaptic netrin-G2 adhesion protein, lead to a distinct neurodevelopmental disorder. Am J Hum Genet. 2019;105(5):1048–56.
Liguori M, Nuzziello N, Simone M, Amoroso N, Viterbo RG, et al. Association between miRNAs expression and cognitive performances of pediatric multiple sclerosis patients: a pilot study. Brain Behav. 2019;9(2):e01199.
Sugimoto M, Sasaki S, Watanabe T, Nishimura S, Ideta A, et al. Ionotropic glutamate receptor AMPA 1 is associated with ovulation rate. PLoS ONE. 2010;5(11):e13817.
Cushman RA, Miles JR, Rempel LA, McDaneld TG, Kuehn LA, et al. Identification of an ionotropic glutamate receptor AMPA1/GRIA1 polymorphism in crossbred beef cows differing in fertility. J Anim Sci. 2013;91(6):2640–6.
Golestanpour H, Javadi G, Sheikhha MH. Association of GRIA1 polymorphisms with ovarian response to human menopausal gonadotropin in Iranian women. Clin Exp Reprod Med. 2020;47(3):207–12.
Wang W, Li F, Sun Y, Lei L, Zhou H, et al. Aquaporin-1 retards renal cyst development in polycystic kidney disease by inhibition of Wnt signaling. FASEB J. 2015;29(4):1551–63.
Wang C, Yan M, Jiang H, Wang Q, Guan X, et al. Protective effects of puerarin on acute lung and cerebrum injury induced by hypobaric hypoxia via the regulation of aquaporin (AQP) via NF-κB signaling pathway. Int Immunopharmacol. 2016;40:300–9.
Su X, Song Y, Jiang J, Bai C. The role of aquaporin-1 (AQP1) expression in a murine model of lipopolysaccharide-induced acute lung injury. Respir Physiol Neurobiol. 2004;142(1):1–11.
Li H, Shi H, Gao M, Ma N, Sun R. Long non-coding RNA CASC2 improved acute lung injury by regulating miR-144-3p/AQP1 axis to reduce lung epithelial cell apoptosis. Cell Biosci. 2018;8:15.
Chen X, Hu J, Pan Y, Tang Z. Novel noncoding RNAs biomarkers in acute respiratory distress syndrome. Expert Rev Respir Med. 2020;14(3):299–306.
Fang Y, Hu J, Wang Z, Zong H, Zhang L, et al. LncRNA H19 functions as an aquaporin 1 competitive endogenous RNA to regulate microRNA-874 expression in LPS sepsis. Biomedicine Pharmacother = Biomedecine & pharmacotherapie. 2018;105:1183–91.
Mu Y, Yuan Y, Han W, Pian P, Li L, et al. Exploring the molecular mechanism and biomarker of recurrent spontaneous abortion based on RNA sequencing analysis. Clin Labor. 2020;66(10):2055–61.
Gu H, Li L, Du M, Xu H, Gao M, et al. Key gene and functional pathways identified in unexplained recurrent spontaneous abortion using targeted RNA sequencing and clinical analysis. Front Immunol. 2021;12:717832.
The authors thank all participants and volunteers of the Second Affiliated Hospital of Hainan Medical University in this study.
The project was supported by Major Science and Technology Program of Hainan Province (No. ZDKJ2017007). The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All participants have signed informed consent. And our study was approved by the Ethics committee of the Second Affiliated Hospital of Hainan Medical University and the 1964 Helsinki declaration.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Figure S1. The quality of the total RNA from 3 patients with RSA and 3 normal abortion patients
. Figure S2. The relationship between upregulated lncRNAs and downregulated miRNAs
. Figure S3. Volcano map of the differentially expressed lncRNAs
. Figure S4. Volcano map of the differentially expressed mRNAs
. Figure S5. Volcano map of the differentially expressed miRNAs. A represents the known miRNAs; B represents the novel miRNAs
About this article
Cite this article
Huang, Y., Hao, J., Liao, Y. et al. Transcriptome sequencing identified the ceRNA network associated with recurrent spontaneous abortion. BMC Med Genomics 14, 278 (2021). https://doi.org/10.1186/s12920-021-01125-4