Next Article in Journal
Synthetic Melanin Acts as Efficient Peptide Carrier in Cancer Vaccine Strategy
Next Article in Special Issue
A Fatal Case of Metastatic Pulmonary Calcification during the Puerperium
Previous Article in Journal
Therapeutic Effects of Acer palmatum Thumb. Leaf Extract (KIOM-2015E) on Benzalkonium Chloride-Induced Dry Eye in a Mouse Model
Previous Article in Special Issue
IL-6 and IL-8: An Overview of Their Roles in Healthy and Pathological Pregnancies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 23
10.3390/ijms232314972
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

First Trimester Maternal Plasma Aberrant miRNA Expression Associated with Spontaneous Preterm Birth

by
Danai Mavreli
1,
Mariana Theodora
2,
Margaritis Avgeris
3,
Nikolas Papantoniou
2,
Panagiotis Antsaklis
2,
George Daskalakis
2 and
Aggeliki Kolialexi
4,*
1
Laboratory of Medical Genetics, School of Medicine, National and Kapodistrian University of Athens, 106 79 Athens, Greece
2
1st Department of Obstetrics and Gynecology, School of Medicine, National and Kapodistrian University of Athens, 106 79 Athens, Greece
3
Laboratory of Clinical Biochemistry–Molecular Diagnostics, Second Department of Pediatrics, School of Medicine, National and Kapodistrian University of Athens, “P. & A. Kyriakou” Children’s Hospital, 106 79 Athens, Greece
4
Department of Genetics, Institute of Child Health, 106 79 Athens, Greece
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(23), 14972; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms232314972
Submission received: 31 October 2022 / Revised: 17 November 2022 / Accepted: 23 November 2022 / Published: 29 November 2022
(This article belongs to the Special Issue Pathogenesis of Pregnancy-Related Complication)
Browse Figures
Versions Notes

Abstract

:
Spontaneous Preterm Delivery (sPTD) is one of the leading causes of perinatal mortality and morbidity worldwide. The present case–control study aims to detect miRNAs differentially expressed in the first trimester maternal plasma with the view to identify predictive biomarkers for sPTD, between 320/7 and 366/7 weeks, that will allow for timely interventions for this serious pregnancy complication. Small RNA sequencing (small RNA-seq) of five samples from women with a subsequent sPTD and their matched controls revealed significant down-regulation of miR-23b-5p and miR-125a-3p in sPTD cases compared to controls, whereas miR-4732-5p was significantly overexpressed. Results were confirmed by qRT-PCR in an independent cohort of 29 sPTD cases and 29 controls. Statistical analysis demonstrated that miR-125a is a promising early predictor for sPTL (AUC: 0.895; 95% CI: 0.814-0.972; p < 0.001), independent of the confounding factors tested, providing a useful basis for the development of a novel non-invasive predictive test to assist clinicians in estimating patient-specific risk.

1. Introduction

Spontaneous Preterm Delivery (sPTD), defined as delivery that occurs before the completion of 37 weeks of gestation, is a complex pregnancy-related complication with multiple etiologies affecting 5–18% of all pregnancies worldwide [1].
sPTD is a leading cause of perinatal mortality and morbidity accounting for approximately 16% of all deaths and 35% of deaths among newborns. [2]. Premature infants that survive are at higher risk of developmental and neurological dysfunction later in life [3]. It is therefore widely accepted that the foresight in women at risk for sPTD early in pregnancy may enable, through earlier intervention, prolonged pregnancy and improved neonatal outcome [4].
Despite extensive efforts, screening strategies to accurately predict sPTD are still unsatisfactory. Currently, the most effective predictive tool combines well-known risk factors including maternal age, maternal BMI, racial origin, maternal behavior, spontaneous or assisted conception and prior history of sPTD, resulting in the identification of ~18% of nulliparous and 38% of parous sPTD cases, with a 10% false positive rate [5]. The prediction rate has been shown to increase by 54.8% by combining a priori risk factors with the measurement of cervical length [6].
MicroRNAs (miRNAs) are a class of non-coding RNAs of 19–25nt long responsible for regulating gene expression at the posttranscriptional level, without altering the genetic code. Concerning mode of action, miRNAs target the 3′ untranslated region (UTR) of mRNAs, resulting in translational inhibition and finally in mRNA degradation [7]. miRNAs have been documented to modulate most biological processes, such as cell proliferation, survival, differentiation and apoptosis. Interestingly, miRNAs can be detected within the cells, but they are also released in a stable form, due to their small size, into body fluids, such as plasma and serum [8].
Often, altered expression levels of circulating miRNAs have been involved with the initiation and progression of various diseases, including pregnancy complications, suggesting their role as potential biomarkers.
Regarding sPTD, several studies have investigated circulating miRNAs as potential biomarkers for its prediction. However, most studies were restricted to miRNAs associated to cardiovascular disease, clusters from chromosome 14 (C14MC), chromosome 19 (C19MC) and miR-371-3 cluster or were performed using samples collected later during pregnancy, in the second or even the third trimester [9,10,11,12,13,14,15]. However, an ideal biomarker would be detectable during the first trimester of pregnancy, when routine screening for fetal aneuploidies is also offered, to reduce anxiety throughout pregnancy in low-risk women and allow for the close monitoring and timely interventions in the high-risk group [16].
The aim of the present study is to identify differentially expressed miRNAs in first trimester maternal plasma and to evaluate their clinical value as novel biomarkers for the early prediction of sPTD. Analysis was performed using small RNA next-generation sequencing (small RNA-seq) followed by quantitative real-time polymerase chain reaction (qRT-PCR) to confirm the results. The sPTD group consisted of women who experienced premature delivery between 320/7 and 366/7 weeks of gestation (moderate/late sPTD), which represent more than 8% of all premature deliveries [17]. Born at this gestational age, premature neonates, are at a relatively lower risk of mortality and morbidity than early preterm births. Still, the impact on healthcare worldwide is significant due to their increased risks compared to full-term births [18].

2. Results

2.1. Demographic and Clinical Characteristics

Table 1 describes the comparison between the demographic and the clinical characteristics of the cases and controls. Beside pre-pregnancy BMI, no significant differences were noted between the two groups.

2.2. Small RNA-seq Analysis

Differentially expressed miRNAs were screened using small RNA-seq in five pairs of plasma samples. A total of 1682 unique circulating miRNAs were identified across all first trimester maternal plasma samples analyzed. Among these, 387 miRNAs were up-regulated and 95 miRNAs were down-regulated in sPTD cases compared to controls. Three circulating miRNAs showed significant changes (p-Value < 0.05, FC >1.5) in the sPTD group as compared to the control group. Specifically, significantly decreased levels of miR-23b-5p (p < 0.047, FC = 0.60) and miR-125a-3p (p = 0.014, FC = 0.59) were detected during the first trimester of gestation in women who subsequently experienced sPTD between 320/7 and 366/7 weeks of gestation whereas miR-4732-5p was significantly overexpressed (p = 0.048, FC = 1.70) (Supplement S1).
A scatter plot is used to display differentially expressed circulating RNAs (Figure 1). Significantly dysregulated circulating miRNAs are presented in a volcano plot (Figure 2).

2.3. GO Analysis

The downstream targets of the significantly dysregulated miRNAs were identified commonly between Targetscan and miRDB databases and used for GO analysis. The significant biological processes associated with these miRNA target genes are presented in Table 2. KEGG pathway analysis revealed that the significantly underexpressed miRNAs target signaling pathways associated with the T cell receptor signaling pathway, Homologous recombination and Osteoclast differentiation (Table 3). No significant pathways were identified for the target genes of the up-regulated miRNAs.

2.4. qRT-PCR Assays

The expression levels miR-23b-5p and miR-125a-3p were further determined in an independent cohort of 58 first trimester maternal plasma samples (29 sPTD cases and 29 controls) by qRT-PCR in order to confirm the reliability of the small RNA-seq analysis. The selection of miRNAs for verification by qRT-PCR was based on p-Values (Supplement S1). When compared with the control group, miR-125a-3p was significantly down-regulated in sPTD cases (p < 0.001; Figure 3A); where for miR-23b-5p, no statistically significant changes were observed (p = 0.750; Figure 3B). The qRT-PCR data were further used to assess the ability of miRNAs to discriminate between women at risk for sPTD and controls using ROC analysis. As expected, miR-125a-3p was highlighted to differentiate sPTD cases from uncomplicated pregnancies (AUC: 0.895; 95% CI: 0.814–0.972; p < 0.001; Figure 3C). Binomial logistic regression analysis demonstrated that miR-125a-3p represents a significant predictor of pregnancy outcome independent of the confounding factors tested (Table 4).

3. Discussion

In the current study, we performed small RNA-seq in first trimester maternal plasma samples to investigate whether a differential miRNA expression profile is associated with subsequent sPTD and to provide a useful basis for the development of a novel non-invasive predictive test to assist clinicians in estimating patient-specific risks.
Compared to women who delivered at term, maternal plasma miRNA profiling of women who subsequently developed sPTD revealed significant alterations in the expression level of miR-23b-5p, mIR-125a-5p and mIR-4732-5p. It is noteworthy that none has previously been related to sPTD. Downstream targets, however, have been associated with biological processes with an important role in the pathogenesis of sPTD, e.g., immune response, inflammation and apoptosis.
More precisely, miR-4732-5p showed a significant over expression in the first trimester maternal plasma of women destined to deliver sPTD compared to those with uncomplicated pregnancies. miR4732-5p has a proven role in tumorigenesis [19,20]. However, our GO analysis revealed that its downstream targets include members of suppressors of the cytokine signaling (SOCS) family, namely SOCS1-SOCS7. The SOCS protein family is implicated in the negative regulation of cytokine signaling and therefore in the regulation of pregnancy and labor [21,22]. Blumenstein et al. reported overexpression of SOCS1, SOCS2 and SOCS3 proteins in human placenta in an attempt to either prevent the entrance of pathogens into the maternal circulation or suppress trophoblast cytokine production and inhibit their damaging effects on the fetus [23]. In the control group, i.e., women who delivered at term, a decrease in SOCS1 and SOCS3 was observed possibly reflecting the involvement of inflammatory cytokines in mechanisms leading to delivery. However, no alteration in SOCS protein levels was demonstrated in the sPTD group. Still, further investigation is required to uncover the relationship between circulating miR-4732-5p and SOCS proteins as the lack of change in SOCS levels in the sPTD group could be attributed to immaturity.
Moreover, 23b-5p, a pleiotropic miRNA demonstrating diverse effects on various biological processes, was significantly down-regulated in the first trimester maternal plasma in women who later delivered prematurely as compared to uncomplicated pregnancies. miR-23b is a member of the miR-23-27-24 family which consists of two paralogs with the miR-23a cluster (miR-23a-27a-24-2) found on chromosome 19 and the intragenic miR-23b cluster (miR-23b27b-24-1) located on chromosome 9 within the C9orf3 gene [24]. Previous research suggested that miR-23b plays vital roles in cancer development where it exerts either oncogenic or tumor suppressor activity [25,26]. It is also implicated in the regulation of angiogenesis and endothelial cells homeostasis and may serve as biomarkers for the early diagnosis of acute myocardial infarction [27,28]. Moreover, research findings suggest that miR-23b overexpression enhanced the expression of IL-10 which is vital for normal pregnancy, and low IL-10 levels are associated with pregnancy complications [29]. Hence, down-regulation of miR-23b has been reported in first trimester maternal plasma in women who later developed late onset pre-eclampsia as compared to uncomplicated pregnancies, indicating a possible involvement in the pathogenesis of the condition [24].
Finally, miR-125a-5p is involved in development and cell differentiation and therefore it is implicated in several malignancies including gastric, cervical, lung and ovarian cancers, retinoblastoma and neuroblastoma [30,31,32]. Recently, it was demonstrated that miR-125a suppresses cell proliferation and migration and inhibits angiogenesis by regulating its downstream target vascular endothelial growth factor A (VEGFA), indicating a potential role to the pathophysiology of pre-eclampsia [33].
Down-regulation of miR-125a was detected in the discovery cohort using small RNA-seq and was further confirmed by the verification cohort. Following statistical analysis, the AUC for miR-125a-5p was 89%, demonstrating promising diagnostic potential as a first trimester screening test for the prediction of subsequent sPTD at 33–36 weeks of gestation. Furthermore, logistic regression analysis revealed that miR-125a-5p is an early predictor of sPTD independently of well-known risk factors (previous sPTD, maternal smoking, maternal pre-pregnancy BMI, maternal age, fetal gender and mode of conception).
Hence, once validated in a large cohort, miR-125a-5p might be used to identify women at high risk for sPTD, allowing for close monitoring and/or clinical interventions such as cervical cerclage. They can also assist the development of novel therapeutics. To date, although the role of miRNA therapeutics has not yet been translated into clinical practice, miRNA mimics and miRNA suppressors have shown significant efficiency in various cancers [34].
To the best of our knowledge, this is the first study that aimed to identify novel biomarkers for the early prediction of women at risk for sPTD using a sensitive high-throughput method for miRNA plasma profiling followed by verification of the results obtained using RT-qPCR. We acknowledge, however, that the small sample size poses a significant limitation for our study.
Still, this pilot study aims to offer improved candidate biomarkers predictive of a subsequent sPTD to be validated in a larger series of pregnant women.

4. Materials and Methods

4.1. Study Population

Maternal samples for this retrospective study for sPTD were collected from pregnant women during first trimester prenatal screening for fetal aneuploidies between March 2018 and December 2020 as previously described [35]. At the same time, maternal demographic characteristics and medical history were obtained and recorded in an electronic database. Gestational age was determined based on the last menstrual period and confirmed by ultrasound measurement of the fetal crown rump length (CRL). Pregnancy outcomes were made known after the completion of pregnancies from the maternity hospital files and were also recorded in the database.
The selection of samples for analysis was carried out using a case–control design. Cases were pregnant women with a subsequent sPTD, between 320/7 and 366/7 weeks of gestation. The control group consisted of participants with termly delivered neonates. Cases and controls were matched with respect to maternal age and duration of storage at −80°. Only European women with singleton pregnancies who delivered a phenotypically normal live born neonate were included in the study. Women with missing information or insufficient plasma for analysis, those with chronic diseases or pregnancy complications other than sPTD (e.g., pre-eclampsia, gestational diabetes, fetal growth restriction) as well as women with PPROM or signs of intra-uterine infection/inflammation at the time of admission were excluded from the study.
Over the study period, a total of 1809 plasma samples were collected and stored at −80 °C. Through the database search, 63 women were identified to be diagnosed with sPTD. Of those, 28 were excluded from the study because of fetal chromosomal abnormality/major fetal malformations (n = 4), presence of a multiple pregnancy (n = 6), miscarriage or fetal death before 24 weeks or termination (n = 5), lost to follow-up (n = 4) or due to inadequate plasma sample (n = 9). Plasma sample of a pregnant woman who subsequently delivered preterm before the 32nd week of gestation (early sPTD) was excluded to ensure increased homogeneity among cohort. Different aliquots from the same cohort have been used previously in a proteomic study [35].
Finally, 34 plasma samples from women who subsequently delivered prematurely, between 320/7 and 366/7 weeks, fulfilled the inclusion criteria and were retrieved for analyses along with their matched controls. Selected samples were analyzed in two phases: samples from five women with a subsequent sPTD between 320/7 and 366/7 weeks of gestation and with their matched controls were analyzed using small RNA-seq (discovery cohort) to identify miRNAs with significantly aberrant expression levels between the compared groups. Based on discovery cohort data, the remaining 58 samples (29 sPTL cases and 29 controls) were used in the second phase of the study to confirm the differential expression of selected miRNAs. None of the samples were previously thawed and refrozen.
Written informed consent to collect and use the biological samples and clinical information for research purposes was obtained from all participants prior to inclusion in the study. The study was approved by Alexandra’s Hospital ethics committee (P.N. 9/5-1-2018) and conducted according to the standards of the 1975 Declaration of Helsinki, as revised in 2008.

4.2. Clinical Definitions

sPTD is defined by the presence of regular uterine contractions (at least two uterine contractions every 10 min for 30 min, as confirmed by external tocometry) in combination with cervical changes occurring prior to 366/7 weeks that required hospitalization [36]. A term pregnancy is defined as a delivery from 37 completed weeks to less than 42 and was used to describe the optimal timing for a good outcome for the mother and baby [37].

4.3. Methods

4.3.1. Total RNA and miRN Isolation

Total RNA enriched in miRNAs was extracted from EDTA-preserved maternal plasma using the mirVana miRNA Isolation kit (Thermo Fisher, Waltham, MA, USA) according to the manufacturer’s instructions. The RNA integrity was determined with polyacrylamide gel electrophoresis. The concentration and purity of RNA were checked by NanoDrop ND-1000 spectrophotometer (Thermo Fisher, Wilmington, DE, USA) and 2100 Bioanalyzer Instruments (Agilent, Santa Clara, CA, USA).

4.3.2. Small RNA-seq

Profiling of circulating miRNA was performed by small RNA-seq in Illumina NextSeq 500 platform as previously described [24]. Briefly, small RNA libraries were generated using the TruSeq® miRNA Sample Prep kit v2 (Illumina, San Diego, CA, USA). RNA template of each sample was sequentially ligated to 3′ and 5′ adapters, reverse transcribed, PCR amplified and selected on agarose gels by size of ~130–150 bp (correspond to ~15–35nt small RNAs). Then, PCR amplified fragments were eluted from gel pieces, purified and quantified by Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). The completed libraries were diluted to a final concentration of 8 pM and processed for cluster generation on the Illumina cBOT using TruSeq Rapid SR cluster kit (Illumina, San Diego, CA, USA). Sequencing was performed on Illumina NextSeq 500 platform using TruSeq Rapid SRB kits (Illumina, San Diego, CA, USA).

4.3.3. Sequencing Data Analysis

The complete raw sequences from Illumina NextSeq 500 were generated as clean reads by real-time base calling and quality filtering. The clean reads were recorded in FASTQ format, containing the read information, sequences and quality encoding. The 3′ prime adapter sequences were trimmed and reads with lengths <15 nucleotides were removed. Those with length >15nt were aligned to all human miRNAs (miRBase v21) using NovoAlign software (NovoCraft, Selangor, Malaysia) allowing maximum one mismatch per sequence [38]. The miRNA read counts were used to estimate the expression level of each miRNA. All counts were normalized by reads per million [39]. Differential expression and statistical analysis of sequencing data were performed by the DESeq2 package in R.
Circulating miRNAs having p-Value ≤ 0.05 and FC ≥ 1.5 were considered significant. p-Values were adjusted for multiple hypotheses testing (<0.01 compared to control) using the method of Benjamini–Hochberg to establish a false discovery rate (FDR).

4.3.4. miRNAs Target Gene Prediction and Gene Ontology Analysis

Potential targets of the significantly dysregulated miRNAs were obtained using Targetscan7 (http://www.targetscan.org/vert_71/) and MirdbV5 (http://mirdb.org/miRDB/) algorithms accessed 21 January 2020. Only overlapping results between these databases were accepted as potential targets.
Targets gene were submitted to the WebGestalt web-tool http://bioinfo.vanderbilt.edu/webgestalt/ accessed on 21 January 2020 for Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling enrichment analyses [40,41]. GO terms and KEGG pathways having an adjusted p-Value ≤ 0.05 were considered significant.

4.3.5. Quantitative Real-Time Polymerase Chain Reaction Verification

The expression levels of representative miRNAs, differentially expressed in the first trimester maternal plasma of women who later experienced sPTD as compared to the control group were verified in an independent cohort consisting of 58 plasma samples (n = 29 sPTL and n = 29 controls) using qRT-PCR. cDNA synthesis and qRT-PCR were performed using a TaqMan miRNA Reverse Transcription kit (Applied Biosystems, Inc., Foster City, CA 94404, USA) and individual TaqMan MiRNA assays (Applied Biosystems, Inc., USA), following the manufacturer’s recommendations. All reactions were run in triplicate in a LC480 LightCycler system (Roche GmbH, Rotkreuz, Switzerland). The miRNA gene expression was determined using the 2−ΔΔCt method [42]. RNU44 (Applied Biosystems, Inc., Foster City, CA 94404, USA) was used for normalization purposes.

4.3.6. Statistical Analysis

Statistical analyses were conducted in IBM SPSS Statistics 20 software (IBM Corp., Armonk, NY, USA). Comparisons of maternal demographic and clinical characteristics between the two groups were compared using Pearson chi-square test for the evaluation of categorical variables or the Mann–Whitney U test for continuous variables. A two-tailed Fisher’s exact test was applied to test GO and pathway enrichment of the target genes of the significant differentially expressed miRNAs. Receiver Operating Characteristic curves (ROC) were applied to evaluate the diagnostic value of each miRNA using the qRT-qPCR data. Binomial logistic regression analyses were performed using the occurrence of sPTD as the dependent variable and miRNA expression levels, previous sPTD, maternal smoking, maternal pre-pregnancy BMI, maternal age, fetal gender and mode of conception as independent variables. p-Values < 0.05 were considered significant.

5. Conclusions

The present study demonstrated that small RNA profiling allows for the identification of novel biomarkers in the first trimester maternal plasma that can potentially be used in clinical practice for early, minimally invasive, prediction of sPTD. Analysis revealed for the first time that the expression levels of circulating, miR-125a-5p, miR-23b and miR-4732-5p are significantly different in plasma samples collected from women who subsequently experience sPTD, as compared to uncomplicated pregnancies, implying their critical roles in the pathology of sPTD. Further studies should be performed for a more complete understanding of the topic. More importantly, the results obtained in our cohort demonstrated that miR-125a-5p may be used as reliable independent biomarkers to predict poor pregnancy outcomes.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms232314972/s1. Supplement S1. MiRNAs with significant altered expression in the 1st trimester maternal plasma of women who later delivered preterm compared to term deliveries.

Author Contributions

D.M.: experimental work, data analysis, literature search and manuscript preparation, M.T.: sample and data collection, sample selection, literature search and manuscript preparation, M.A.: bioinformatics, statistical analyses and manuscript critical review, N.P.: sample selection, literature search and manuscript critical review, P.A.: sample selection, literature search and manuscript critical review, G.D.: sample and data collection, sample selection, and critical review A.K.: conceived and designed the experiments and provided critical review. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was approved by the Alexandra’s hospital ethics committees (P.N. 9/5-1-2018) and conducted according to the standards of the 1975 Declaration of Helsinki, as revised in 2008.

Informed Consent Statement

Written informed consent to collect and use the biological samples and clinical information for research purposes was obtained from all participants prior to inclusion in the study.

Data Availability Statement

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Please also refer to the Section 4 “Materials and Methods”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goldenberg, R.L.; Culhane, J.F.; Iams, J.D.; Romero, R. Epidemiology and causes of preterm birth. Epidemiology and causes of preterm birth. Lancet 2008, 371, 75–84. [Google Scholar] [CrossRef]
  2. Chawanpaiboon, S.; Vogel, J.P.; Moller, A.B.; Lumbiganon, P.; Petzold, M.; Hogan, D.; Landoulsi, S.; Jampathong, N.; Kongwattanakul, K.; Laopaiboon, M.; et al. Global, regional, and national estimates of levels of preterm birth in 2014: A systematic review and modelling analysis. Lancet Glob. Health 2019, 7, e37–e46. [Google Scholar] [CrossRef] [Green Version]
  3. Natarajan, G.; Shankaran, S. Short- and Long-Term Outcomes of Moderate and Late Preterm Infants. Am. J. Perinatol. 2016, 33, 305–317. [Google Scholar] [PubMed]
  4. Suff, N.; Story, L.; Shennan, A. The prediction of preterm delivery: What is new? Semin Fetal Neonatal Med. 2019, 24, 27–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Beta, J.; Akolekar, R.; Ventura, W.; Syngelaki, A.; Nicolaides, K.H. Prediction of spontaneous preterm delivery from maternal factors, obstetric history and placental perfusion and function at 11–13 weeks. Prenat. Diagn. 2011, 31, 75–83. [Google Scholar] [CrossRef] [PubMed]
  6. Greco, E.; Gupta, R.; Syngelaki, A.; Poon, L.C.; Nicolaides, K.H. First-trimester screening for spontaneous preterm delivery with maternal characteristics and cervical length. Fetal Diagn. Ther. 2012, 31, 154–161. [Google Scholar] [CrossRef]
  7. Mellios, N.; Sur, M. The Emerging Role of microRNAs in Schizophrenia and Autism Spectrum. Disord. Front. Psychiatry 2012, 25, 39. [Google Scholar] [CrossRef] [Green Version]
  8. Valadi, H.; Ekström, K.; Bossios, A.; Sjöstrand, M.; Lee, J.J.; Lötvall, J.O. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat. Cell Biol. 2007, 9, 654–659. [Google Scholar] [CrossRef] [Green Version]
  9. Hromadnikova, I.; Kotlabova, K.; Krofta, L. First Trimester Prediction of Preterm Delivery in the Absence of Other Pregnancy-Related Complications Using Cardiovascular-Disease Associated MicroRNA Biomarkers. Int. J. Mol. Sci. 2022, 23, 3951. [Google Scholar] [CrossRef]
  10. Winger, E.E.; Reed, J.L.; Ji, X. Early first trimester peripheral blood cell microRNA predicts risk of preterm delivery in pregnant women: Proof of concept. PLoS ONE 2017, 12, e0180124. [Google Scholar] [CrossRef]
  11. Winger, E.E.; Reed, J.L.; Ji, X.; Gomez-Lopez, N.; Pacora, P.; Romero, R. MicroRNAs isolated from peripheral blood in the first trimester predict spontaneous preterm birth. PLoS ONE 2020, 15, e0236805. [Google Scholar] [CrossRef] [PubMed]
  12. Ran, Y.; Yin, N.; Huang, D.; Zhao, Y.; Yang, J.; Zhang, H.; Qi, H. 
Identification and Characterization of Circular RNA as a Novel Regulator and Biomarker in Preterm Birth. Front. Bioeng. Biotechnol. 2020, 8, 566984. [Google Scholar] [CrossRef] [PubMed]
  13. Gray, C.; McCowan, L.M.; Patel, R.; Taylor, R.S.; Vickers, M.H. Maternal plasma miRNAs as biomarkers during mid-pregnancy to predict later spontaneous preterm birth: A pilot study. Sci. Rep. 2017, 7, 815. [Google Scholar] [CrossRef] [Green Version]
  14. Menon, R.; Debnath, C.; Lai, A.; Guanzon, D.; Bhatnagar, S.; Kshetrapal, P.K.; Sheller-Miller, S.; Salomon, C. Garbhini Study Team. Circulating Exosomal miRNA Profile During Term and Preterm Birth Pregnancies: A Longitudinal Study. Endocrinology 2019, 160, 249–275. [Google Scholar] [CrossRef] [Green Version]
  15. Enquobahrie, D.A.; Hensley, M.; Qiu, C.; Abetew, D.F.; Hevner, K.; Tadesse, M.G.; Williams, M.A. Candidate Gene and MicroRNA Expression in Fetal Membranes and Preterm Delivery Risk. Reprod Sci. 2016, 23, 731–737. [Google Scholar] [CrossRef] [Green Version]
  16. Mavreli, D.; Theodora, M.; Kolialexi, A. Known biomarkers for monitoring pregnancy complications. Expert Rev. Mol. Diagn. 2021, 21, 1115–1117. [Google Scholar] [CrossRef]
  17. Frey, H.A.; Klebanoff, M.A. The epidemiology, etiology, and costs of preterm birth. Semin Fetal Neonatal Med. 2016, 21, 68–73. [Google Scholar] [CrossRef]
  18. Di Renzo, G.C.; Cabero Roura, L.; Facchinetti, F.; Helmer, H.; Hubinont, C.; Jacobsson, B.; Jørgensen, J.S.; Lamont, R.F.; Mikhailov, A.; Papantoniou, N.; et al. Preterm Labor and Birth Management: Recommendations from the European Association of Perinatal Medicine. J. Matern Fetal Neonatal Med. 2017, 30, 2011–2030. [Google Scholar] [CrossRef]
  19. Qi, L.; Gao, C.; Feng, F.; Zhang, T.; Yao, Y.; Wang, X.; Liu, C.; Li, J.; Li, J.; Sun, C. MicroRNAs associated with lung squamous cell carcinoma: New prognostic biomarkers and therapeutic targets. J. Cell Biochem. 2019, 120, 18956–18966. [Google Scholar] [CrossRef] [PubMed]
  20. Chen, X.; Yu, L.; Hao, K.; Yin, X.; Tu, M.; Cai, L.; Zhang, L.; Pan, X.; Gao, Q.; Huang, Y. Fucosylated exosomal miRNAs as promising biomarkers for the diagnosis of early lung adenocarcinoma. Front. Oncol. 2022, 12, 935184. [Google Scholar] [CrossRef]
  21. Huang, S.; Liu, K.; Cheng, A.; Wang, M.; Cui, M.; Huang, J.; Zhu, D.; Chen, S.; Liu, M.; Zhao, X.; et al. SOCS Proteins Participate in the Regulation of Innate Immune Response Caused by Viruses. Front. Immunol. 2020, 11, 558341. [Google Scholar] [CrossRef] [PubMed]
  22. Blumenstein, M.; Bowen-Shauver, J.M.; Keelan, J.A.; Mitchell, M.D. Identification of suppressors of cytokine signaling (SOCS) proteins in human gestational tissues: Differential regulation is associated with the onset of labor. J. Clin. Endocrinol. Metab. 2002, 87, 1094–1097. [Google Scholar] [CrossRef] [PubMed]
  23. Blumenstein, M.; Keelan, J.A.; Bowen-Shauver, J.M.; Mitchell, M.D. Suppressors of cytokine signaling proteins in human preterm placental tissues. J. Mol. Endocrinol. 2005, 35, 165–175. [Google Scholar] [CrossRef] [PubMed]
  24. Mavreli, D.; Lykoudi, A.; Lambrou, G.; Papaioannou, G.; Vrachnis, N.; Kalantaridou, S.; Papantoniou, N.; Kolialexi, A. Deep Sequencing Identified Dysregulated Circulating MicroRNAs in Late Onset Preeclampsia. In Vivo 2020, 34, 2317–2324. [Google Scholar] [CrossRef] [PubMed]
  25. Ru, L.; Wang, X.M.; Niu, J.Q. The miR-23-27-24 cluster: An emerging target in NAFLD pathogenesis. Acta Pharmacol. Sin. 2022, 43, 1167–1179. [Google Scholar] [CrossRef] [PubMed]
  26. Guo, Y.X.; Wang, N.; Wu, W.C.; Li, C.Q.; Chen, R.H.; Zhang, Y.; Li, X. The Role of miR-23b in Cancer and Autoimmune Disease. J. Oncol. 2021, 2021, 6473038. [Google Scholar] [CrossRef] [PubMed]
  27. Bang, C.; Fiedler, J.; Thum, T. Cardiovascular importance of the microRNA-23/27/24 family. Microcirculation 2012, 19, 208–214. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, G.K.; Zhu, J.Q.; Zhang, J.T.; Li, Q.; Li, Y.; He, J.; Qin, Y.W.; Jing, Q. Circulating microRNA: A novel potential biomarker for early diagnosis of acute myocardial infarction in humans. Eur. Heart J. 2010, 31, 659–666. [Google Scholar] [CrossRef]
  29. Azizieh, F.Y.; Raghupathy, R. L-10 and pregnancy complications. Clin. Exp. Obstet. Gynecol. 2017, 44, 252–258. [Google Scholar] [CrossRef]
  30. Li, S.; Zheng, S. Down-Regulation of Circ_0032833 Sensitizes Colorectal Cancer to 5-Fluorouracil and Oxaliplatin Partly Depending on the Regulation of miR-125-5p and MSI1. Cancer Manag. Res. 2020, 12, 11257–11269. [Google Scholar] [CrossRef]
  31. Fu, K.; Zhang, L.; Liu, R.; Shi, Q.; Li, X.; Wang, M. MiR-125 inhibited cervical cancer progression by regulating VEGF and PI3K/AKT signaling pathway. World J. Surg. Oncol. 2020, 18, 115. [Google Scholar] [CrossRef]
  32. Szebeni, G.J.; Vizler, C.; Kitajka, K.; Puskas, L.G. Inflammation and Cancer: Extra- and Intracellular Determinants of Tumor-Associated Macrophages as Tumor Promoters. Mediat. Inflamm. 2017, 2017, 9294018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. 33 Wade, S.M.; Ohnesorge, N.; McLoughlin, H.; Biniecka, M.; Carter, S.P.; Trenkman, M.; Cunningham, C.C.; McGarry, T.; Canavan, M.; Kennedy, B.N.; et al. Dysregulated miR-125a promotes angiogenesis through enhanced glycolysis. EBioMedicine 2019, 47, 402–413. [Google Scholar] [CrossRef] [Green Version]
  34. Hanna, H.; Hossain, G.S.; Kocerha, J. The Potential for microRNA Therapeutics and Clinical Research. Front. Genet. 2019, 16, 478. [Google Scholar] [CrossRef] [Green Version]
  35. Mavreli, D.; Theodora, M.; Lambrou, G.; Avgeris, M.; Papantoniou, N.; Treager-Synodinos, J.; Daskalakis, G.; Kolialexi, A. First trimester maternal plasma proteomic changes predictive of spontaneous moderate/late preterm delivery. J. Matern Fetal Neonatal Med. 2017, submitted.
  36. Tsiartas, P.; Holst, R.M.; Wennerholm, U.B.; Hagberg, H.; Hougaard, D.M.; Skogstrand, K.; Pearce, B.D.; Thorsen, P.; Kacerovsky, M.; Jacobsson, B. Prediction of spontaneous preterm delivery in women with threatened preterm labour: A prospective cohort study of multiple proteins in maternal serum. BJOG 2012, 119, 866–873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. ACOG Committee Opinion No 579: Definition of term pregnancy. Obstet. Gynecol. 2013, 122, 1139–1140. [CrossRef]
  38. Kozomara, A.; Griffiths-Jones, S. miRBase: Annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014, 42, D68. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Bullard, J.H.; Purdom, E.; Hansen, K.D.; Dudoit, S. Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinform. 2010, 11, 94. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Zhang, B.; Kirov, S.; Snoddy, J. WebGestalt: An integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33, W741–W748. [Google Scholar] [CrossRef]
  41. Wang, J.; Duncan, D.; Shi, Z.; Zhang, B. WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): Update 2013. Nucleic Acids Res. 2013, 41, W77–W83. [Google Scholar] [CrossRef] [PubMed]
  42. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 14972 g001 550
Figure 1. Scatter plot of circulating miRNA levels in the first trimester maternal plasma from women who subsequently experienced sPTD when compared to term ones. Red points show up-regulated miRNAs and green points show down-regulated miRNAs. Data from the sPTD group are plotted on the vertical axis (y-axis). Data from the term deliveries (controls) are plotted on the horizontal axis (x-axis).
Figure 1. Scatter plot of circulating miRNA levels in the first trimester maternal plasma from women who subsequently experienced sPTD when compared to term ones. Red points show up-regulated miRNAs and green points show down-regulated miRNAs. Data from the sPTD group are plotted on the vertical axis (y-axis). Data from the term deliveries (controls) are plotted on the horizontal axis (x-axis).
Ijms 23 14972 g001
Ijms 23 14972 g002 550
Figure 2. Volcano plot illustrates significantly differentially abundant circulating miRNAs. The -log10 (Benjamini–Hochberg corrected p-Value) is plotted against the log2 (fold change: sPTD cases/controls). The non-axial vertical lines denote ±1.5-fold change while the non-axial horizontal line denotes p-Value = 0.05, which is our significance threshold (prior to logarithmic transformation).
Figure 2. Volcano plot illustrates significantly differentially abundant circulating miRNAs. The -log10 (Benjamini–Hochberg corrected p-Value) is plotted against the log2 (fold change: sPTD cases/controls). The non-axial vertical lines denote ±1.5-fold change while the non-axial horizontal line denotes p-Value = 0.05, which is our significance threshold (prior to logarithmic transformation).
Ijms 23 14972 g002
Ijms 23 14972 g003 550
Figure 3. Statistical analysis of the data obtained using qRT-PCR. (A,B): Box-plots of the expression levels of miR-23b and miR-125a in the first trimester maternal plasma in women with sPTD versus uncomplicated pregnancies based on data derived from qRT-qPCR. p-Values calculated by Mann–Whitney U test. (C): ROC curve analysis for the prediction of sPTD based on miRNAs levels. The diagram is a plot of the % sensitivity (true-positive rate) vs. 100% specificity (false-positive rate).
Figure 3. Statistical analysis of the data obtained using qRT-PCR. (A,B): Box-plots of the expression levels of miR-23b and miR-125a in the first trimester maternal plasma in women with sPTD versus uncomplicated pregnancies based on data derived from qRT-qPCR. p-Values calculated by Mann–Whitney U test. (C): ROC curve analysis for the prediction of sPTD based on miRNAs levels. The diagram is a plot of the % sensitivity (true-positive rate) vs. 100% specificity (false-positive rate).
Ijms 23 14972 g003
Table
Table 1. Maternal and neonatal characteristics of moderate/late sPTD cases and controls included in the study.
Table 1. Maternal and neonatal characteristics of moderate/late sPTD cases and controls included in the study.
CharacteristicControls
(n = 34)		sPTD Cases
(n = 34)		p-Value
Maternal age (y)
Median (min–max)		31.00
(24.80–38.80)		30.65
(24.80–38.20)		0.898
Pre-pregnancy BMI (kg/m2) Median (min–max)25.95
(22.00–31.20)		24.65
(20.90–31.90)		0.028
Cigarette smoker
  • Yes
  • No

12
22
8
26
0.40

Mode of conception
  • Spontaneous
  • In vitro fertilization

28 (82.4%)
6 (17.6%)
26 (76.4%)
8 (23.6%)
0.750

Parity
  • Nulliparous
  • Parous

17 (50%)
17 (50%)
14 (41.2%)
20 (58.8%)
0.616

Previous PTD
  • No
  • Yes

32 (94.1%)
2 (5.88%)
27 (79.4%)
7 (20.58%)
0.104

Mode of Delivery
  • Vaginal
  • Caesarean section

29 (85.2%)
5 (14.7%)
22 (64.7%)
12 (35.2%)
0.075

Neonatal Gender
  • Male
  • Female

18 (52.95%)
16 (47.05%)
16 (47.05%)
18 (52.95%)
0.95
Table
Table 2. GO annotations of miRNA’s target genes.
Table 2. GO annotations of miRNA’s target genes.
Down-Regulated miRNAs
GO.IDTermp-ValueGenes
GO:0048585Negative regulation of response to stimulus0.000RGS4//IL10//NDFIP1//GPRC5A//STAP1//MTMR4//TMEM127//BRCA1//GREM2//HGS//STK38//FAP//TPT1
GO:0014854Response to inactivity0.001MTMR4//IL10
GO:1900119Positive regulation of execution phase of apoptosis0.001FAP//BOK
GO:0009968Negative regulation of signal transduction0.001RGS4//GPRC5A//IL10//MTMR4//TMEM127//BRCA1//GREM2//HGS//STK38//STAP1//TPT1
GO:0043032Positive regulation of macrophage activation0.001STAP1//IL10
GO:0010324Membrane invagination0.001ABCA1//STAP1//HGS
GO:1900120Regulation of receptor binding0.002GREM2//IL10
GO:0010648Negative regulation of cell communication0.002RGS4//GPRC5A//IL10//MTMR4//TMEM127//BRCA1//GREM2//HGS//STK38//STAP1//TPT1
GO:0023057Negative regulation of signaling0.002RGS4//GPRC5A//IL10//MTMR4//TMEM127//BRCA1//GREM2//HGS//STK38//STAP1//TPT1
GO:0070230Positive regulation of lymphocyte apoptotic process0.002IL10//PDCD1
Up-Regulated miRNAs
GO.IDTermp-ValueGENES
GO:0006397mRNA processing0.008RAVER2//LUC7L3//ALKBH5
GO:0045625Regulation of T-helper 1 cell differentiation0.009SOCS5
GO:0071071Regulation of phospholipid biosynthetic process0.009HTR2A
GO:0002829Negative regulation of type 2 immune response0.010SOCS5
GO:0006054N-acetylneuraminate metabolic process0.010CMAS
GO:0007175Negative regulation of epidermal growth factor-activated receptor activity0.010SOCS5
GO:0035970Peptidyl-threonine dephosphorylation0.010PPM1E
GO:0042118Endothelial cell activation0.010SOCS5
GO:0045623Negative regulation of T-helper cell differentiation0.010SOCS5
GO:0043371Negative regulation of CD4-positive, alpha-beta T cell differentiation0.011SOCS5
Table
Table 3. Pathway enrichment analysis of the target genes of down-regulated miRNAs.
Table 3. Pathway enrichment analysis of the target genes of down-regulated miRNAs.
Pathway IDDefinitionp-ValueFDREnricment
Score		Gene
Ratio		Genes
hsa04660T cell receptor signaling0.00210.6582.6716810.16IL10//PDCD1//TEC
hsa03440Homologous recombination0.00460.7132.3355660.11BRCA1//FAM175A
hsa04380Osteoclast differentiation0.042511.371060.11SIRPA//TEC
Table
Table 4. Logistic regression analysis for the prediction of sPTD patients according to maternal. plasma miRNAs expression levels.
Table 4. Logistic regression analysis for the prediction of sPTD patients according to maternal. plasma miRNAs expression levels.
Univariate Analysis
CovariantOR a95% CI bp-Value c
miR-23b-5p1.0570.943–1.1840.343
miR-125a-3p0.6540.490–0.8720.004
Previous PTD7.1540.809–63.300.077
Smoking0.5350.176–1.6240.269
Maternal pre-pregnancy
BMI		0.8080.650–1.0030.053
Mode of conception
Natural
IVF
1.00
1.512

0.425–5.384

0.523
Fetal gender
Male
Female
1.00
1.319

0.470–3.703

0.600
Maternal age1.0160.885–1.1670.819
Multivariate analysis d
CovariantOR a95% CI bp-Value c
miR125a-3p0.4910.295–0.8170.006
a Odds Ratio. b 95% Confidence Interval of the estimated OR. c calculated by test for trend. d Multivariate models adjusted for tested miRNAs, previous PTD, smoking, BMI, mode of conception, fetal gender and age.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mavreli, D.; Theodora, M.; Avgeris, M.; Papantoniou, N.; Antsaklis, P.; Daskalakis, G.; Kolialexi, A. First Trimester Maternal Plasma Aberrant miRNA Expression Associated with Spontaneous Preterm Birth. Int. J. Mol. Sci. 2022, 23, 14972. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms232314972

AMA Style

Mavreli D, Theodora M, Avgeris M, Papantoniou N, Antsaklis P, Daskalakis G, Kolialexi A. First Trimester Maternal Plasma Aberrant miRNA Expression Associated with Spontaneous Preterm Birth. International Journal of Molecular Sciences. 2022; 23(23):14972. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms232314972

Chicago/Turabian Style

Mavreli, Danai, Mariana Theodora, Margaritis Avgeris, Nikolas Papantoniou, Panagiotis Antsaklis, George Daskalakis, and Aggeliki Kolialexi. 2022. "First Trimester Maternal Plasma Aberrant miRNA Expression Associated with Spontaneous Preterm Birth" International Journal of Molecular Sciences 23, no. 23: 14972. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms232314972

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop