Skip to main content
Download PDF

Transcriptomic analysis delineates preterm prelabor rupture of membranes from preterm labor in preterm fetal membranes

BMC Medical Genomics volume 17, Article number: 72 (2024) Cite this article

Abstract

Background

Globally, preterm birth remains the leading cause of death in children younger than 5 years old. Spontaneous preterm birth is comprised of two events that may or may not occur simultaneously: preterm labor and preterm prelabor rupture of membranes (PPROM). To further explore the concept that spontaneous preterm birth can result from the initializing of two separate but overlapping pathological events, we compared fetal membrane tissue from preterm labor deliveries to fetal tissue from preterm labor with PPROM deliveries. We hypothesized that the fetal membrane tissue from preterm labor with PPROM cases will have an RNA-seq profile divergent from the fetal membrane tissue from preterm labor controls.

Methods

Chorioamnion, separated into amnion and chorion, was collected from eight gestationally age-matched cases and controls within 15 min of birth, and analyzed using RNA sequencing. Pathway enrichment analyses and functional annotations of differentially expressed genes were performed using KEGG and Gene Ontogeny Pathway enrichment analyses.

Results

A total of 1466 genes were differentially expressed in the amnion, and 484 genes were differentially expressed in the chorion (log2 fold change > 1, FDR < 0.05) in cases (preterm labor with PPROM), versus controls (preterm labor only). In the amnion, the most significantly enriched (FDR < 0.01) KEGG pathway among down-regulated genes was the extracellular matrix receptor interaction pathway. Seven of the most significantly enriched pathways were comprised of multiple genes from the COL family, including COL1A, COL3A1, COL4A4, and COL4A6. In the chorion, the most significantly enriched KEGG pathways in up-regulated genes were chemokine, NOD receptor, Toll-like receptor, and cytokine-cytokine receptor signaling pathways. Similarly, KEGG pathway enrichment analysis for up-regulated genes in the amnion included three inflammatory pathways: cytokine-cytokine interaction, TNF signaling and the CXCL family. Six genes were significantly up regulated in chorionic tissue discriminated between cases (preterm labor with PPROM) and controls (preterm labor only) including GBP5, CXCL9, ALPL, S100A8, CASP5 and MMP25.

Conclusions

In our study, transcriptome analysis of preterm fetal membranes revealed distinct differentially expressed genes for PPROM, separate from preterm labor. This study is the first to report transcriptome data that reflects the individual pathophysiology of amnion and chorion tissue from PPROM deliveries.

Peer Review reports

Background

Preterm birth remains the leading cause of death in children younger than 5 years old and is defined by the World Health organization as births before 37 completed weeks of gestation or fewer than 259 days from the first date of a woman’s last menstrual period [1]. Most recent global estimates show approximately 10.6% of all live births are preterm, while in the U.S., the rate is 10.1% [2].

Multiple preterm birth classifications have been created to help guide research on risk determinants, pathology, and epidemiology of this syndrome. The etiology of preterm birth is multifactorial, with nutritional, environmental and socio-demographics all contributing to an increase in risk of occurrence [3]. It is believed that approximately 30% of preterm birth cases are caused by infection and inflammation [4]. More recently, it has been hypothesized that maternal vaginal, cervical, and intestinal microbiomes may contribute to the etiology of preterm birth [5,6,7].

Preterm birth can be classified as either spontaneous (spontaneous preterm birth), or provider initiated by cesarean or labor induction [4]. Spontaneous preterm birth is comprised of two events that may or may not occur simultaneously: preterm labor (PTL) and preterm prelabor rupture of membranes (PPROM). PTL occurs when regular contractions result in the opening of the cervix after week 20 but before week 37 of pregnancy. Although both term and preterm labor present with clinical similarities such as uterine contraction, and cervical dilation, it has been proposed that preterm labor is not merely labor that starts before 37 weeks of gestation, but a syndrome caused by multiple pathologies [8, 9].

PPROM is the spontaneous rupture of fetal membranes before labor before the 37th week of gestation. Fetal membranes are attached to the placenta, and are composed of two layers: the amnion, the layer directly in contact with amnionic fluid, and the chorion which is attached to the maternal decidua [10]. It has been hypothesized that PPROM is a disease of the fetal membranes initiated, in part, by the dysregulation of inflammatory pathways [11,12,13]. Connective tissue disorders are also associated with weakened fetal membranes and an increased incident of PPROM [14, 15].

Many studies have focused on the developmental events that differentiate the onset of spontaneous preterm birth from the onset of PPROM, yet we are still unable to elucidate the pathways that differentiate the two pathophysiologies. To further explore the concept that spontaneous preterm birth can result from the initializing of two separate but overlapping pathological events (PTL +/- PPROM), we compared fetal membrane tissue from PTL deliveries versus tissue from PTL w/ PPROM deliveries. In this study, we hypothesize that the fetal membrane tissue from PTL + PPROM cases will have an RNA-seq profile divergent from the fetal membrane tissue from PTL only controls.

Methods

Tissue collection

Chorioamnion, was collected within 15 min of delivery from pregnancies complicated by preterm birth between the gestational ages of 27 to 34 weeks at Women and Infants Hospital, in Providence, RI from April 2018- February 2019. Pregnancies which included illicit drug use, IVF, or known congenital anomalies were excluded from the study. A 3 cm2 portion of the chorioamnionic membrane was taken from the place of rupture for the vaginal deliveries, and from the incision site for c-section deliveries. The chorioamnion membrane was then separated into amnion and chorion and immediately frozen in liquid nitrogen, or placed in RNAlater solution, and stored at -80 °C. Data obtained from clinical records included vaginal or cesarean section delivery, gestational age at delivery, and chronological initiation of PPROM or preterm labor, maternal age, race/ethnicity, and parity (See Additional File 3). s Informed consent for the study has been waived by our local ethics committee, The Care New England Institutional Review Board.The tissue being collected is considered residual pregnancy tissue; thus informed consent was not necessary.For this study, gestational age, and delivery type (vaginal or cesarean) were matched for 4 sets of cases and controls. Control samples were defined as those in which PTL presented first (PTL only), while cases were defined as those in which PPROM presented first, followed by PTL (PPROM + PTL). Both amnion and chorion samples were included in this study. The 16 tissue samples were shipped to BGI Americas (Hong Kong, China) for RNA extraction and analysis.

Total RNA was extracted from each sample using a phenol-chloroform extraction protocol to produce mRNA. Total concentration of RNA was evaluated on an Agilent Bioanalyzer 2100 (Santa Clara, CA). The RNA sequencing sample quality was determined by 28 S/18S and RIN (RNA Integrity Number) values. All 16 samples met the requirements of library construction and sequencing established by BGI.

RNA-Seq and data analysis

RNA-Seq library construction was conducted by BGI America. Double stranded PCR products were circularized by the splint oligo sequence, and the resulting single strand circle DNA (ssCir DNA) was formatted as the final library. The library was amplified with phi29 to make DNA nanoball (DNB) which has more than 300 copies of one molecule. The DNBs were loaded into a patterned nanoarray and paired-end 100 base reads were generated by Probe-Anchor Synthesis (cPAS).

Raw fastq files for all samples were checked for quality using FastQC v0.11.5 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and did not require quality or adapter trimming. Salmon v1.3.0 [16] was then used in mapping-based mode to align reads and quantify transcripts. The Gencode v36 GRCh38.p13 human genome and Gencode v36 transcriptome versions were used for alignment. The R package txImport v1.22.0 [17] was used to import Salmon quant files containing transcript counts with bias correction data into R. Finally, the DESeq2 v1.31.6 R package [18] was used with default settings to perform differential expression analysis using the bias-corrected transcript counts from Salmon. Differential expression analysis was performed with three replicates per experiment group and compared PPROM unaffected versus affected tissue expression profiles in amnion and chorion tissue samples separately.

Functional enrichment and annotation

We utilized the DAVID Database for Annotation, Visualization, and Integrated Discovery Functional Annotation Tool v6.8 [19, 20] to perform Gene Ontology (GO) and KEGG Pathway enrichment analyses on the genes which were differentially up-regulated or down-regulated by at least 2-fold. Genes with a p-value ≤ 0.05 and an absolute fold change ≥ 2 were significantly differentially expressed. GO terms and KEGG pathways were significantly enriched if they had a false discovery rate (FDR) corrected p-value ≤ 0.05.

Heatmaps of the DESeq2 functional enrichment results were prepared using the ggplot2 v3.3.6 R package [21].

RT-PCR validation of RNASeq gene expression results

To confirm differentially expressed (DE) gene data, six genes were chosen for RT-PCR analysis based on a log2-fold changes higher than 4.0 and a highly statistically significant padj: guanylate binding protein 5 (GBP5), alkaline phosphatase, biomineralization associated (ALPL), S100 calcium binding protein A8 (S100A8), c-x-c motif chemokine ligand 9 (CXCL9), matrix metallopeptidase 25 (MMP25), and caspase 5 (CASP5). 8 cases and 8 controls, separate from the samples used for RNA-Seq, were analyzed. Each sample was assessed in duplicate, and the %CV between the duplicates was < 3%. All primers were ‘off the shelf’ sequences purchased from ThermoFisher Scientific ( Asheville, NC). The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene was chosen as the reference genes to estimate relative quantification. The relative gene expression was calculated using the 2−∆∆Ct method.”

The mean gestational age for both groups was 33.2 weeks Case and control data were compared using the Wilcoxon Rank sum test.

Total RNA was extracted from each sample using the RNeasy Plus mini kit (Qiagen, Germantown, MD), and the concentration and quality of the RNA was evaluated using a Nanodrop 2000 Spectrophotometer (ThermoFisher Scientific, Ashville, NC). cDNA was synthesized from the RNA using Superscript IV First-Strand Synthesis System (Thermo Scientific, Wilmington, DE)) using oligo primers, and performed in a SimpliAmp Thermal Cycler (Applied Biosystems, Asheville, NC) RT-PCR amplification was performed in Thermocycler 7500 using TaqMan Fast Universal Master Mix 2X Kit (ThermoFisher Scientific, Asheville, NC).

Results

Demographic and clinical characteristics of the cases and control groups are presented in Table 1.

Table 1 Demographic and clinical characteristics of cases and controls for RNA-Seq
Full size table

Evaluation of gene expression

Genome-wide gene expression in amnion and chorion tissue were evaluated using RNA-Seq technology for 4 cases (PPROM + PTL) and 4 control (PTL only) samples. Cases and controls were matched by gestational age and type of delivery (vaginal or cesarean).

A summary of the read alignments to the reference genome is presented in Additional File 1. Sixteen samples were sequenced using BGISEQ platform, averaging about 4.66Gb bases per sample. The average mapping ratio with reference genome was 95.98%, the average mapping ratio with gene was 77.97%; 18,844 genes were identified.

A principal component analysis (PCA) was performed on the gene expression log-transformed values for the amnion and for the chorion, resulting in one case and one control being identified as outliers. PCA was rerun with 6 samples, resulting in a large portion of the variance being assigned to PC1 and PC2 in both the amnion (61%, 16% respectively) and the chorion (65%, 19% respectively). The PCA plot for the amnion showed equal variability between cases and controls, while the plot of the chorion samples showed a tight grouping of the controls, and large variability in the cases.

A total of 1466 genes were differentially expressed (DE) in the amnion and 484 genes were differentially expressed in the chorion (log2fold change > 1, FDR < 0.05) (Fig. 1A, B). Among the 1466 DE genes in the amnion, 991 were upregulated and 475 were down regulated in cases versus controls (full list in Additional File 2a). Among the 484 DE genes in the chorion, 371 were upregulated and 113 were down regulated in cases versus controls (full list in Additional File 2b).

Fig. 1
figure 1

A. Heatmap of differentially expressed genes in PPROM affected amnion tissue. Heatmap of gene expression z-scores of 1466 differentially expressed genes clustered by row means between PPROM affected and unaffected amnion tissue. B. Heatmap of differentially expressed genes in PPROM affected chorion tissue. Heatmap of gene expression z-scores of 533 differentially expressed genes clustered by row means between PPROM affected and unaffected chorion tissue

Full size image

Enrichment analysis using gene ontology

Amnion

For KEGG pathways, the most significantly enriched (FDR < 0.01) pathway among the downregulated genes was the extracellular matrix (ECM) -receptor interaction pathway. Interestingly, 7 of the most significantly enriched pathways were comprised of multiple genes from the COL family, including COL1A1, COL3A1, COL4A4 and COL4A6 (Fig. 2). KEGG pathways that were significantly enriched (FDR < 0.01) for upregulated genes included the cytokine-cytokine receptor interaction pathway, TNF signaling pathway, and the CXCL family (Fig. 2).

Fig. 2
figure 2

Clustered heatmap of significantly altered KEGG pathways in PPROM affected amnion tissue. Heatmap of differentially expressed genes between PPROM affected and unaffected amnion tissue among significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways clustered by row and column means. The -log10(adjusted p-value) derived from differential expression analysis is illustrated by the Gene Score color bar. The false discovery rate (FDR)-corrected p-value derived from functional enrichment analysis using the DAVID Functional Annotation Tool is illustrated by the Enrichment Term Score color bar. The Log2FC color bar identifies sets of genes which are significantly up or downregulated. The significantly enriched KEGG pathways (FDR p-value < = 0.05) in which the differentially expressed genes are a member of are labeled on the lower axis

Full size image

A total of 89 GO terms were significantly enriched (FDR < = 0.01) among the DE genes for amnion. Within the Biological Process category, over 60 terms were significantly enriched for upregulated genes, while only the term extracellular matrix organization was significantly enriched in downregulated genes (Fig. 3A). In the Cellular Component category, the terms plasma membrane, extracellular region, and cytoplasm were among the 15 significantly enriched terms for upregulated genes, while plasma membrane, extracellular matrix extracellular exosome, and basement membrane were the only terms significantly enriched in downregulated genes (Fig. 3B). Lastly, in the Molecular Function category, signal transducer activity, receptor activity, protein binding and collagen binding were among the terms significantly enriched for the upregulated genes, while only the term extracellular matrix organization was significantly enriched for the downregulated genes (Fig. 3C).

Fig. 3
figure 3

A. Bubble plots of significantly altered biological process gene ontology terms in PPROM affected amnion tissue. Significantly up or downregulated biological process gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected amnion tissue. B. Bubble plots of significantly altered cellular components gene ontology terms in PPROM affected amnion tissue. Significantly up or downregulated cellular component gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected amnion tissue. C. Bubble plots of significantly altered molecular function gene ontology terms in PPROM affected amnion tissue. Significantly up or downregulated molecular function gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected amnion tissue

Full size image

Chorion

In the chorion, the most significantly enriched (FDR < 0.01) KEGG pathways in upregulated genes included chemokine signaling pathway, NOD-like receptor signaling pathway, cytokine-cytokine receptor interactions pathway, and Toll-like receptor signaling pathway (Fig. 4). There were no significantly enriched (FDR < 0.01) KEGG pathways for downregulated genes in the chorion (Fig. 4).

Fig. 4
figure 4

Clustered heatmap of significantly altered KEGG pathways in PPROM affected chorion tissue. Heatmap of differentially expressed genes between PPROM affected and unaffected chorion tissue among significantly enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways clustered by row and column means. The -log10(adjusted p-value) derived from differential expression analysis is illustrated by the Gene Score color bar. The false discovery rate (FDR)-corrected p-value derived from functional enrichment analysis using the DAVID Functional Annotation Tool is illustrated by the Enrichment Term Score color bar. The Log2FC color bar identifies sets of genes which are significantly up or downregulated. The significantly enriched KEGG pathways (FDR p-value < = 0.05) in which the differentially expressed genes are a member of are labeled on the lower axis.

Full size image

A total of 54 GO terms were significantly enriched (FDR < = 0.01) among the DE genes for chorion. Within the Biological Process category, 40 terms were significantly enriched for upregulated genes including cell adhesion, signal transduction and inflammatory response, while no terms were significantly enriched in downregulated genes (Fig. 5A). In the Cellular Component category, the terms plasma membrane, integral component of plasma membrane, and extracellular exosome were among the 9 significantly enriched terms for upregulated genes, while no terms were significantly enriched in downregulated genes (Fig. 5B). Lastly, in the Molecular Function category 5 terms were significantly enriched for upregulated genes, including rage receptor binding, carbohydrate binding, and oxygen transporter activity, while no terms were significantly enriched for downregulated genes (Fig. 5C).

Fig. 5
figure 5

A. Bubble plots of significantly altered biological process gene ontology terms in PPROM affected chorion tissue. Significantly upregulated biological process gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected chorion tissue. No downregulated terms were significantly enriched (FDR < = 0.01). B. Bubble plots of significantly altered cellular components gene ontology terms in PPROM affected chorion tissue. Significantly upregulated cellular component gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected chorion tissue. No downregulated terms were significantly enriched (FDR < = 0.01). C. Bubble plots of significantly altered molecular function gene ontology terms in PPROM affected chorion tissue. Significantly upregulated molecular function gene ontology (GO) terms which were significantly enriched (FDR < = 0.01) among differentially regulated genes in PPROM affected and unaffected chorion tissue. No downregulated terms were significantly enriched (FDR < = 0.01)

Full size image

Validation of gene expression by RT-PCR

To validate genes that were observed to be significantly upregulated in chorionic tissue of PPROM births during RNA-Seq analysis, the expression of 6 genes was assessed by PCR– S100A8, GBP5, MMP25, CXCL9, ALPL, and CASP5. All 6 genes showed upregulation between PPROM cases and controls, although the difference between the two groups did not reach statistical significance for the genes evaluated.

Discussion

Preterm birth is a major cause of morbidity and mortality in infants, yet minimal progress has been made in recent years to identify the pathophysiological pathways that lead to preterm birth. When studied as a single outcome, it is difficult to identify specific pathways or predictive biomarkers for preterm birth because it is a multifactorial syndrome. Recent studies of the transcriptome of gestational tissue have focused on elucidating the pathways involved in preterm birth by comparing tissue from term deliveries with tissue from preterm deliveries [22,23,24]. This approach leads to the comparison of DE genes from healthy control (term) tissue to pathological preterm tissue, which introduces a variability in the gestational age of the tissues being compared. Thus, by approaching preterm birth as the functional outcome of either PTL or PPROM, the gestational age variable is removed, and a more vigorous comparison of DE genes can occur. Our study is the first to report transcriptomic data for preterm birth amnion and chorionfrom PTL + PPROM cases compared to PTL only controls. We are also the first to separately analyze the fetal membrane layers to obtain a more complete profile of the fetal membrane transcriptome.

The amnion layer of fetal membranes is composed of collagen fibrils, amnionic cells, and the non-cellular scaffolding ECM in which cells are embedded; thus, the amnionic tissue is believed to correlate with fetal membrane structure and integrity. Our KEGG pathways analysis of amnionic tissue revealed the ECM-receptor pathway to be the most significantly enriched among downregulated genes. Gene ontology term analysis also indicated ECM organization and ECM matrix pathways were enriched among down-regulated genes. Many of these pathways are comprised of multiple genes from the COL family. The downregulation of collagen genes may alter the composition ratio of collagen content in the amnion, which is known to affect tissue rigidity [25, 26]. This multi-pathway down-regulation reported in our study corresponds with previous reviews of the significant role ECM structure and metabolism plays in the pathophysiology of PPROM [27, 28].

The chorion layer of fetal membranes serves as a barrier to inflammatory, immune, and microbial invasion. In our study, the most significantly enriched KEGG pathways in up-regulated genes in the chorion revealed chemokine, NOD receptor, Toll-like receptor, and cytokine-cytokine receptor signaling pathways to be the most significantly enriched. Similarly, KEGG pathway enrichment analysis for upregulated genes in the amnion included three inflammatory pathways: cytokine-cytokine interaction, TNF signaling and the CXCL family. These findings are consistent with previous studies which have described dysregulation of cytokines and chemokines to be involved in PPROM [29,30,31], and confirm that both layers of the fetal membrane are targets of inflammatory pathway dysregulation during PPROM. Interestingly, neither KEGG pathways nor gene ontology enrichment analyses identified any pathways enriched for down-regulated genes in the chorion.

Fetal membrane weakening is a physiological process that occurs at term, but when dysregulated, can result in PPROM. Inflammation is an essential component of this process, with sterile and septic inflammation considered involved in PPROM. Several markers involved with inflammation have been reported to significantly increase in maternal blood, amnionic fluid, and fetal membrane in cases of PPROM [10, 11, 32, 33]. In our study, GBP5 was significantly upregulated in chorionic tissue from PTL + PPROM births compared to PTL only births. Guanylate-binding proteins (GBPs) are a family of GTPases that play an essential role in cell-autonomous immunity. Studies have also indicated that GPBs have essential functions in controlling inflammation and innate immunity through the activation of inflammasome complexes, with GBP5 promoting NLRP3 inflammasome assembly and immunity [34,35,36]. Inflammasomes have been shown to be important to host defense mechanisms in pathophysiological inflammatory processes in the chorioamniotic membranes that accompany labor [37].

CXCL9, a major chemokine that controls the migratory and recruitment patterns of T helper 1 cells, is thought to contribute to proinflammatory responses required for labor onset [38, 39]. Serum CXCL9 levels have been reported to be significantly lower in pregnancies complicated by PPROM when compared to normal pregnancies [40]. Fetal membrane CXCL9 expression has been reported to be significantly increased in patients with PTL compared to term labor [39]. Our study is the first to show that CXCL9 is significantly upregulated in the chorion layer of fetal membranes from patients with PTL + PPROM, versus only PTL.

ALPL, the gene that encodes the alkaline phosphatase enzyme, was also significantly upregulated in chorionic tissue from PPROM cases. Previously, studies have shown that serum ALPL rises in normal pregnancies resulting from an increased production of ALPL by the placenta, and it has often been suggested as a potential serum biomarker for preterm birth [41, 42]. Recently, a study by Paquette et al. determined that ALPL is strongly decreased in preterm birth placenta, but significantly increased in the chorioamniotic membranes in preterm birth cases [43]. Our study confirms the significant increase in ALPL in the chorioamniotic membranes, and further specifies the preterm birth pathophysiology to be PPROM.

S100A8, a low molecular weight proinflammatory protein, was significantly upregulated in the chorion of patients with PTL + PPROM compared to the chorion from patients with only PTL. Significantly elevated serum levels of S100A8 have previously been reported in patients with early pregnancy loss and other adverse pregnancy outcomes [44, 45]. Elevated S100A8 has also been reported in the uterine decidua of patients with recurrent early pregnancy loss [46].

MMP25 is a glycosyl-phosphatidylinositol protease, which degrades ECM proteins, including type 4 collagen, gelatin, fibronectin, and fibrin [47]. Fortunato and Menon have shown that human fetal membranes from PPROM pregnancies have fully functioning MMP systems, with several MMPs expressed in the amniochorion during PPROM, but not MMP25 [48, 49]. Our study showed MMP25 to have the most significant upregulation in chorionic tissue of PPROM pregnancies than any of the other MMPs in the transcriptomic profile. The difference in reported results may be based in the fact that the fetal membrane samples from their study were from PPROM patients not in labor, whereas our fetal membrane samples were from PPROM patients with PTL. Additionally, their studies extracted RNA from intact chorioamniotic tissue, whereas we measured RNA from the amnion and chorion separately.

A caveat to the findings described in this study is the small sample size used in the study. Since our study only collected tissue from preterm births, this limited the number of eligible births to 10% of the total births. In a normally distributed population, PPROM is the cause of 1/3 of preterm births, thus decreasing the occurrence even further. To minimize the potential for confounding variables in this small sample size, we carefully matched the samples by gestational age and mode of delivery. Additionally we were unable to extract sufficient amounts of RNA from the 4 cases and 4 control samples to perform both RNA-Seq and RT-PCR on each sample. We chose to run RNA-Seq on the gestationally age-matched and delivery mode-matched samples and run the RT-PCR controls with the same mean gestational age (33w2d), but not delivery-mode matched. Lastly, it would be informative to perform immunohistochemical localization on the fetal membranes to provide further insight into the contribution of cellular heterogeneity to the transcriptomic differences presented in our study.

Conclusion

In our study, transcriptome analysis of preterm fetal membranes revealed previously reported as well as distinct DE pathways and genes for PPROM, separate from PTL. We are the first to report transcriptome data that reflects the individual pathophysiology of amnion and chorion tissue in PPROM. Although the fetal membrane sample was small, we are certain that we have presented riveting data to support our hypothesis that PPROM and PTL are two separate pathologies of preterm birth. These findings will need to be verified in a larger sample of fetal membranes.

Data availability

The data sets generated and/or analyzed during the current study are available in the GEO repository and are accessible through GEO series accession number GSE243831.

https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE243831.

Abbreviations

ALPL:

alkaline phosphatase, biomineralization associated

CASP5:

caspase 5

CXCL9:

c-x-c motif chemokine ligand 9

DE:

differentially expressed

DNB:

DNA nanoball

ECM:

extracellular matrix

FDR:

false discovery rate

GBP5:

guanylate binding protein 5

GO:

gene ontology

S100A8:

S100 calcium binding protein A8

MMP25:

matrix metallopeptidase 25

PCA:

principal component analysis

PPROM:

preterm prelabor rupture of membranes

PTL:

preterm labor

RIN:

RNA integrity number

RNA-Seq:

RNA sequencing

RT-PCR:

ssCIR DNA:single strand circular DNA

References

  1. Chawanpaiboon S, Moller VJ, Lumbiganon AB, Petzold P, Hogan M. Global, regional, and mational estimates of levels of preterm birth in 2014: a systemic review and modelling analysis. Lancet Glob Health. 2019. Jan;7(1):e37–e46.

  2. Howson CP, McDougall KM, Lawn L. Born too soon: preterm birth matters. Reproductive Health. 2013;10:S1.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Vogel JP, Moller CS, Watananirun A-B, Bonet K, Luumbiganon M. P., The global epidemiology of preterm birth Best practice & research clinical obstetrics and gynaecology, 2018. 52: p. 3–12.

  4. Goldenberg RI, Iams CJ, Romero JD. Preterm Birth 1: epidemiology and causes of preterm birth. Lancet Glob Health. 2008;371:75–84.

    Google Scholar 

  5. Parnell LA, B.C., and, Mysorekar IU. Maternal microbiomes in preterm birth: recent progress and analytical pipelines. Semin Perinat. 2017;41:392–400.

    Article  Google Scholar 

  6. Romero R, Gajer HS, Tarca P, Fadrosh AL, Nikita DW. L The composition and stability of the vaginal microbiota of normal women is different from that of non-pregnant women. Micr Obiome, 2014. 2(4).

  7. Bennett PR, Macintyre BR. Vaginal microbiome in Preterm rupture of membranes. Obstet Gynecol Clin North Am. 2020;47(4):503–21.

    Article  PubMed  Google Scholar 

  8. Romero R, D.S., and, Fisher Preterm SJ. Labor: one syndrome, many causes. Science. 2014;345(6198):760–5.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  9. Gomez-Lopez N, Miller GJ, Farias-Jofre D. M The immunobiology of preterm labor and birth: intra-amniotic inflammation or breakdown of maternal-fetal homeostasis. Reproduction 2022. 164(2).

  10. Choltus H, De Sousa Do Outeiro LM, Coste C, Belville K, Blanchon C. L, and V sapin., Pathophysiological Implication of Pattern Recognition Receptors in Featl Membranes Rupture: RAGE and NLRP Inflammasome Biomedicines 2021. 9: p. 1123.

  11. Romero RMJ, Chaiworaponga CP. Sterile and microbial-associated intra-amnionic inflammation in preterm premature rupture of membrane. J Matern Fetal Neonatal Med. 2015;28(12):1394–409.

    Article  CAS  PubMed  Google Scholar 

  12. Menon R. Preterm prelabor rupture of the membranes: a disease of the fetal membrane. Semin Perinatol. 2017;41:409–17.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Menon R, R.L., and, Lappas M. Fetal membrane architecture, aging and inflammation in pregnancy and parturition. Placenta. 2019;79:40–5.

    Article  CAS  PubMed  Google Scholar 

  14. Barabas AP. Ehlers-Danlos syndrome:associated with prematurity and premature rupture of foetal membranes; possible increase in incidence. Br Med J. 1966;2(5515):682–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Underhill LA, Collis BC, Tucker S R, and, Lechner BE. Association of maternal versus fetal Ehlers-Danlos syndrome status with poor pregnancy outcomes. Reprod Sci, 2022. 29.

  16. Patro R, Love DG, Irizarry M, Kingsford RA. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Soneson C, Robinson LM. MD., Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences F1000Res 2015. 4: p. 1521.

  18. Love MI, Anders HW. Moderated estimation of Fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Huang DW, Lempicki SB. Systemic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  PubMed  Google Scholar 

  20. Huang DW, Lempicki SB. RA., Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res, 2009. 37(1).

  21. H., W., ggplot2: elegant graphics for data analysis. 2016, Springer-: New York.

  22. Eidem HR, McGary AW, Abbot KL, Rokas P. A., Gestational tissue transciptomics in term and preterm human pregnancies: a sytemic review and meta-analysis. BMC Med Genom, 2015. 8(27).

  23. Pereyra S, Sosa C, Bertoni B, Sapiro R. Transcriptomic analysis of fetal membranes reveals pathways involved in preterm birth. BMC Medical Genomics; 2019. p. 12.

  24. Yong HEJ. Current approaches and developments in transcript profiling of the human placenta. Hum Reprod Update. 2020;26(6):799–840.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Oyen ML, C.S., and, Landers DV. Premature rupture of the fetal membranes: is the amnion the major determinant? Am J Obstet Gynecol. 2005;195:510–5.

    Article  Google Scholar 

  26. Shi J-W, Yang LZ-Z. Collagen at the maternal-fetal interface in human pregnancy. Int J Biol Sci. 2020;16(12):2220–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Strauss J. Extracellular matrix dynamics and fetal membrane rupture. Reprod Sci, 2013. 20(2).

  28. Lannon SMR, Eschenbach VJ, Gravett DA MG, and, Adams Waldorf KM. Synergy and interactions among biological pathways leading to preterm premature rupture of membranes. Reprod Sci. 2014;21(10):1215–27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kumar D, Mercer MR, Mansour BM, Redline JM RW and, Moore JJ. The physiology of fetal membrane weakening and rupture; insights gained from the determination of physical properties revisited. Placenta. 2016;42:59–73.

    Article  CAS  PubMed  Google Scholar 

  30. Li W, Li ZX, Chen S, Cui X. H Upregulation of TNF and IL-6 induces preterm premature rupture of membranes by activation of ADAMTS-9 in embryonic membrane cells. Life Sci, 2020. 260.

  31. Zhu J, Ma HM, Peng C, Su F Y, and, Huang L. Expression and clinical significance of NOD-like receptor protein 3 (NLRP3) and caspase-1 in fetal membrane and placental tissues of patients with premature rupture of membrane. Med Sci Monit. 2018;24:1560–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Underhill LA, Tucker AN, Zhand R, Masserlian Z G, and, Lechner BE. Serum decorin and biglycan as potential biomarkers to predict PPROM in early gestation. Reprod Sci. 2019;27(8):1620–6.

    Article  Google Scholar 

  33. Vajrychova M, Pimkova SJ, Barman K, Kukla M. R Comprehensive proteomic investigation of infectious and inflammatory changes in late preterm prelabour rupture of membranes. Sci Rep 2020. 10.

  34. Meunier E, Broz P. Interferon-inducible GTPases in cell autonomous and innate immunity. Cell Microbiol. 2016;18(2):168–80.

    Article  CAS  PubMed  Google Scholar 

  35. Man SIaK. Regulation of inflammasome activation. Immunol Rev. 2015;265(1):6–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Shenoy AR, Kumar WD, Kassa P. H GBP5 promotes NLRP3 inflammasome assembly and immunity in mammals. Sci Express 2012. 29 March.

  37. Modi BP, Pearson TM, Parikh LN. Mutations in fetal genes involved in innate immunity and host defense against microbes increase risk of preterm premature rupture of membranes (PPROM). Mol Genet Genom Med. 2017;5(6):720–9.

    Article  CAS  Google Scholar 

  38. Sykes L, Yap MD, Teoh XJ TG, and, Bennett PR. The Th1:Th2 dichotomy of pregnancy and preterm labor. Mediato rs of Inflammation; 2012.

  39. Mei Y, Liu RY, Zhou Z. IL-27 mediates Th1 cells infiltration in fetal membranes in preterm labor. Reproductive Sci. 2022;29:1764–75.

    Article  CAS  Google Scholar 

  40. Ronzoni S, D’Souza SV, Murphy R. Cytokine changes in maternal peripheral blood correlate with time-to-delivery in pregnancies complicated by premature prelabor rupture of membranes. Reprod Sci. 2018;26(9):1266–76.

    Article  PubMed  Google Scholar 

  41. Moawad AH, Mercer GR, Meis B. The Preterm Prediction Study: the value of serum alkaline phosphatase, alpha-fetoprotein, plasma corticotropin-releasing hormone, and other serum markers for the prediction of spontaneous preterm birth. Am J Obstet Gynecol. 2002;186(5):990–6.

    Article  CAS  PubMed  Google Scholar 

  42. Meyers RE, Addy TS, Garrison CI, Best CZ. Maternal serum placental alkaline phosphatase level and risk for preterm delivery. Am J Obstet Gynecol. 1995;173(1):181–6.

    Article  Google Scholar 

  43. Paquette AG, Bammier MJ, Day T. Placental transcriptomic signatures of spontaneous preterm birth. Am J Obstet Gynecol. 2023;228(73):e1–18.

    Google Scholar 

  44. Nair RR, K.A., and, Singh K. Association of increased S100A8 serum protein with early pregnancy loss. Am J Reprod Immunol. 2014;73:91–4.

    Article  PubMed  Google Scholar 

  45. DelVecchio G, Li LQ, Thamotharan W. Cell-free DNA methylation and transcriptomic signature prediction of pregnancies with adverse outcomes. Epigenetics. 2021;16(6):642–61.

    Article  Google Scholar 

  46. Nair RR, K.A., and, Singh K. Role of inflammatory proteins S100A8 and S100A9 in pathophysiology of recurrent early pregnancy loss. Placenta. 2013;34:824–7.

    Article  CAS  PubMed  Google Scholar 

  47. English WR, Stracke VG, Knauper JO V, and, Murphy G. Catalytic activities of membrane-type 6 matrix metalloproeinase (MMP25). FEBS Lett. 2001;491:137–42.

    Article  CAS  PubMed  Google Scholar 

  48. Fortunato S. a.M.R., Screening of Novel Matrix metalloproteinases (MMPs) in human fetal membranes. J Assist Reprod Genet 2002. 19(10).

  49. Menon R. The role of matrix degrading enzymes and apoptosis in ruptures of membranes. J Soc Gynecol Investig. 2004;11(7):427–37.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

The author(s) disclose receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Institute of General Medical Sciences COBRE for Reproductive Health [P20GM121298].

Author information

Authors and Affiliations

  1. Warren Alpert Medical School at Brown University, Providence, RI, USA

    Lori A. Underhill, J. M. Mennella, A. Uzun & B. E. Lechner

  2. Department of Pediatrics, Women and Infants Hospital, Providence, RI, USA

    Lori A. Underhill, J. M. Mennella, G. A. Tollefson, A. Uzun & B. E. Lechner

  3. Women and Infants Hospital, 101 Dudley St, 02905, Providence, RI, USA

    Lori A. Underhill

Authors
  1. Lori A. Underhill
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. J. M. Mennella
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. G. A. Tollefson
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. A. Uzun
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. B. E. Lechner
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

J.M., B.L., and L.U. conceived and designed the study. J.M. assisted in sample collection. G.T. and A.U. analyzed the data, and created all figures. L.U., G.T., A.U., and B.L. discussed the results. L.U., G.T. and A.U. wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lori A. Underhill.

Ethics declarations

Ethics approval and consent to participate

This prospective research study was conducted with approval by the Care New England Institutional Review Board. Informed consent for this study has been waived by the Care New England Institutional Review Board.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Underhill, L.A., Mennella, J.M., Tollefson, G.A. et al. Transcriptomic analysis delineates preterm prelabor rupture of membranes from preterm labor in preterm fetal membranes. BMC Med Genomics 17, 72 (2024). https://doi.org/10.1186/s12920-024-01841-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12920-024-01841-7

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us