Abstract
Preterm prelabour rupture of membranes beyond the 34th week of gestation (late PPROM) is frequently associated with the risk of the microbial invasion of the amniotic fluid (MIAC) and histological chorioamnionitis (HCA). Hence, we employed a Tandem Mass Tag-based approach to uncover amniotic fluid proteome response to the presence of MIAC and HCA in late PPROM. Protein dysregulation was associated with only five cases in the group of 15 women with confirmed MIAC and HCA. Altogether, 138 amniotic fluid proteins were changed in these five cases exclusively. These proteins were particularly associated with excessive neutrophil responses to infection, such as neutrophil degranulation and extracellular trap formation. We believe that the quantification of these proteins in amniotic fluid may assist in revealing women with the highest risk of excessive inflammatory response in late PPROM.
Similar content being viewed by others
Introduction
Preterm prelabour rupture of membranes (PPROM) is defined as leakage of amniotic fluid prior to 37 weeks of gestation and the leakage of amniotic fluid happening at least two hours prior to the onset of regular uterine contractions1. Late PPROM, occurring between 34 and 37 weeks of gestation, represents more than 50% of cases2. Despite being considered as “near-term,” late PPROM is associated with significant short-term neonatal morbidity, mainly respiratory distress syndrome and a higher risk of early-onset neonatal sepsis3. Concerns about the long-term consequences have also been raised because children born late preterm are at a higher risk of adverse developmental outcomes up to 7 years of age4,5.
Women with late PPROM can be managed actively with an early planned delivery or expectantly through an attempt to prolong the pregnancy. The aim of the active approach is to avoid prolonged fetal exposure to a potential intra-amniotic infection. On the other hand, expectant management may reduce the incidence of neonatal morbidity (especially respiratory distress syndrome) as the neonates are more matured6. A recent meta-analysis did not reveal evidence of a higher incidence of early-onset neonatal sepsis in the group with expectant management7. In addition, neonates in the active-management group were more prone to develop the respiratory distress syndrome7. The authors therefore suggest that expectant management is an acceptable alternative to immediate delivery even in the near-term pregnancy.
Nonetheless, expectant management for all women is still controversial and most national guidelines recommend planned delivery for PPROM shortly after 34 weeks of gestation8,9. The main rationale for the guidelines is the possible presence of subclinical complications such as microbial invasion of the amniotic cavity (MIAC) and intra-amniotic inflammation such as acute histological chorioamnionitis (HCA). These complications are related to the highest incidence of fetal inflammatory response syndrome and short-term neonatal morbidity10,11,12. Moreover, MIAC and HCA are frequently asymptomatic and therefore difficult to being detected. The current dilemma with the delivery timing suggests that a universal approach may not offer an apparent benefit in all cases. Early identification of cases with intra-amniotic and intrauterine complications and personalized clinical management might improve the outcomes of late PPROM pregnancies.
Our group has previously shown that the presence of both MIAC and HCA in women with late PPROM was associated with an increase of amniotic fluid interleukin-6 (IL-6) in comparison with women with either HCA alone, MIAC alone or with absence of both conditions13. However, little is known about the intra-amniotic environment in terms of other proteins, when MIAC and HCA are present in late PPROM. Moreover, the concentration of amniotic fluid IL-6 in the presence of MIAC and HCA may vary from hundreds to hundreds of thousands pg/mL14. Thus, the heterogeneity of individual cases may hamper proper assessment of how severe the intra-amniotic inflammation is.
Despite these facts, a complex proteomic study focused solely on late PPROM is lacking. Liquid chromatography separation coupled to tandem mass spectrometry (LC–MS/MS) enables identification of thousands of proteins in different kinds of biological samples15. One of the main and widely used approaches of protein quantification is multiplexed isobaric tagging allowing simultaneous protein quantification in multiple samples. In this approach, peptides processed from previous protein digestion are labeled with a chemical tag (reporter) in each sample. Differently-labeled samples are then combined before LC–MS/MS analysis. The fragmentation of peptide precursors separates the reporter. Subsequently, the difference in peptide level among samples is assessed according to the relative abundance of reporters in the fragmentation spectra16. Finally, protein quantification is based on the relative quantification of tagged peptides, unique to identified proteins (Fig. 1A).
The aim of this study was to take advantage of the modern methodology to assess how the amniotic fluid proteome responds to the presence of MIAC and/or HCA in women with late PPROM. We believed that these findings might improve our understanding of the molecular processes taking place due to MIAC and HCA in late PPROM pregnancies, and hence facilitate personalized clinical management.
Methods
Reagents and Material
Unless otherwise indicated, chemicals and reagents were purchased from Sigma-Aldrich (St. Louis, MO, USA) in the highest available grade.
Ethical approval
This study was approved by the Institutional Review Board—Ethics Committee, University Hospital Hradec Kralove (March 19, 2008, No. 200804 SO1P; and July 2014, No. 201407 S14P). All research was performed in accordance with relevant guidelines and regulations. Informed consent was obtained from all participants.
Design of late PPROM cohort
The late PPROM cohort was completed considering the following criteria as described previously17. A total of 60 samples of amniotic fluid were collected by transabdominal amniocentesis. Women older than 18 years with a singleton pregnancy, with gestational age between 34 and 37th week, and admitted with a diagnosis of PPROM were included. The diagnosis of PPROM was conditioned by pooling of amniotic fluid in the vagina and detection of insulin-like growth factor-binding protein 1 (ACTIM PROM test; MedixBiochemica, Kauniainen, Finland) in the vaginal fluid. Pregnancy complications such as intrauterine growth restriction, structural malformation or chromosomal abnormality of the fetus, maternal hypertension, diabetes mellitus, pre-eclampsia, thyroid disease, vaginal bleeding, and fetal hypoxia were exclusion criteria. All clinical characteristics are listed in Table 1. Potential differences among the subgroups of late PPROM were tested by nonparametric Kruskal–Wallis test for continuous variables and are presented in Table 1 as median with interquartile range. Categorical variables were compared by the χ2 test and are presented as numbers (%).
Diagnosis of MIAC and HCA
MIAC was established based on a qPCR detection of Ureaplasma spp., Mycoplasma hominis, Chlamydia trachomatis, and positivity for bacterial DNA identified by PCR targeting 16S rRNA and analyzed by sequencing, and aerobic/anaerobic culture of amniotic fluid. Detection of DNA of Ureaplasma spp., M. hominis, C. trachomatis and other bacteria and aerobic and anaerobic cultures of amniotic fluid was carried out as in previously published work18. HCA was diagnosed on the basis of the degree of neutrophil infiltration that was inspected separately in the amnion, chorion-decidua, chorionic plate, and umbilical cord. The degree of infiltration was assessed according to the criteria described by Salafia19. A diagnosis of HCA was conditioned by the grades of infiltration 3–4 for the chorion-decidua, 3–4 for the chorionic plate, 1–4 for the umbilical cord, and/or 1–4 for the amnion. Samples were divided according to the MIAC and HCA status and were distributed into four subgroups: subgroup with the presence of both MIAC and HCA, subgroup with MIAC alone, subgroup with HCA alone, and subgroup with absence of both MIAC and HCA (Table 1).
Amniotic fluid sample preparation
Amniotic fluid samples were prepared according to the workflow successfully applied in our previous reports17,20. Briefly, ultrasound-guided transabdominal amniocentesis was carried out upon admission before the administration of corticosteroids, antibiotics (penicillin G, clindamycin, azithromycin, or clarithromycin), or tocolytics. Up to five mL of amniotic fluid was supplemented with a protease inhibitor cocktail (Complete Mini, EDTA free, Roche Diagnostics, Basel, Switzerland). The treated samples were then centrifuged at 300 × g for 15 min and passed through a 0.22 μm syringe filter (TPP, Trasadingen, Switzerland) to remove cells and debris. The total protein concentration was determined by Bicinchoninic acid assay and ballast proteins were removed as described in the Supplementary Document.
Protein digestion
Enzymatic digestion of proteins was processed in 200 mM triethylammonium bicarbonate buffer, pH 8.5 (TEAB, Thermo Scientific, Rockford, IL, USA) with 1% of sodium deoxycholate. Disulfide bonds were reduced with 5 mM tris(2-carboxyethyl) phosphine hydrochloride at 60 °C for 1 hour (h), and the free thiol groups were blocked using 10 mM methyl methanethiosulfonate at room temperature for 30 min. Amniotic fluid proteins were digested with lysyl endopeptidase (Waco Pure Chemical Industries, Osaka, Japan) at 37 °C for 4 h, followed by digestion with trypsin (Promega, Madison, WI, USA) at 37 °C overnight. Both enzymes were added at a 1:50 enzyme to protein ratio (w/w).
Tandem Mass Tag 10plex isobaric labeling and multiplex strategy
Tandem Mass Tag (TMT) 10plex isobaric label reagents (Thermo Scientific, Rockford, IL, USA) were dissolved according to the manufacturer’s instructions. All the amniotic fluid samples were evaporated to dryness and re-suspended in 8.57 µL of 100 mM TEAB, and 1.57 µL was taken from each sample and combined to prepare a global internal standard (GIS). The total volume of the GIS was split into 10 aliquots. Individually digested samples and the GIS aliquots, both at concentration 1 µg/µL, were incubated with 10 μL of the TMT reagents at room temperature for 1 h. The reaction was terminated by adding 5% hydroxylamine. Nine amniotic fluid samples and one GIS aliquot were mixed in 1:1 ratio to create seven multiplexes (Fig. 1A). To complete the seventh multiplex, the last three GIS aliquots were labeled by 126, 128C, 129C channels. However, these labeled GIS aliquots were not considered in the data evaluation. Prepared multiplexes were cleaned by peptide extraction on solid phase as described in the Supplementary Document.
Liquid chromatography coupled to tandem mass spectrometry (LC–MS/MS)
The multiplexes were redissolved in 0.1% trifluoroacetic acid (TFA), 2% acetonitrile (AcN) to a concentration of 1 µg/µL and injected twice on an UltiMate 3000 RSLCnano System (Thermo Scientific, Bremen, Germany) in three technical replicates. The analytical system consisted of a PepMap100 C18, 3 µm, 100 Å, 75 µm × 20 mm trap column and a PepMap RSLC C18, 2 µm, 100 Å, 75 µm × 500 mm analytical column (both from Thermo Scientific, Bremen, Germany) and was previously used in our proteomic analyses21. The samples were loaded onto the trap column at a flow rate 8 µl/min of 0.1% TFA, 2% AcN for 3 min. Tryptic peptides were separated via a linear gradient running from 2 to 45% of 0.1% formic acid, 80% AcN at a flow rate 200 nL/min for 240 min. Eluted peptides were sprayed into a Q Exactive Plus mass spectrometer using a Nanospray Flex ion source (Thermo Scientific, Bremen, Germany) at 1.8 kV spray voltage. Positive ion full scan MS spectra were acquired in the range of m/z 350–1600, with 3 × 106 AGC target in the Orbitrap at 70,000 resolution with a maximum ion time of 50 ms. A lock mass of m/z 445.12003 ([C2H6SiO]6) was used for internal calibration of the mass spectra. The fragmentation (MS/MS) spectra were acquired for the 15 most intense precursors. Other parameter settings were assumed from the literature22. Isolation window of 1.2 m/z and normalized collision energy of 33 was used. Each fragmentation spectrum was acquired at a resolution of 35,000, with a 2 × 105 AGC target and a maximum injection time of 100 ms. The first mass was fixed to 120 m/z.
Data treatment
Recorded spectra were processed and searched in MaxQuant 1.5.5.1 against a reviewed SwissProt human protein database including contaminants downloaded in 2016. Trypsin digestion specificity with up to two missed cleavages allowed was used. Cysteine thiomethylation was set as a fixed modification and oxidation of methionine was selected as a variable modification. The mass tolerance in MS and MS/MS mode was left at the default value for the initial search and 6 ppm was set as mass tolerance in MS mode for the main search. The false discovery rate for protein identification was left at the default value. Output files from MaxQuant were processed in Perseus23 and R24. Potential contaminants and proteins identified by site and by reverse sequence were removed. The reporter intensity of proteins in each sample was divided by reporter intensity of proteins in the GIS, both adjusted to correction factors. Resulting ratios were log2 transformed and normalized to the median. One-way analysis of variance (ANOVA) with Benjamin–Hochberg false discovery rate (FDR) correction was applied to compare normalized protein ratios in the subgroups of late PPROM. Perseus Tukey’s significant honest difference post-hoc test on the significance level of 0.05 was used on the ANOVA significant hits. Gene ontology (GO) analysis was carried out using StringApp plugin25 in the Cytoscape platform26 (enrichment p value < 0.01). Violin plots, principal component analysis (PCA) plot and bubble plots were designed in the R software. Bar graphs of GO terms were prepared in GraphPad Prism 6.
Amniotic fluid IL-6 and histone H1.5 level determination by ELISA
ELISA kit for human IL-6 (R&D systems, MN, USA) and human histone H1.5 (MyBiosource, CA, USA) were used to determine IL-6 and H1.5 levels in amniotic fluid. Samples with IL-6 concentration above the top calibration point were diluted to fit the calibration range. Resulting absorbance was read at 450 nm with correction to 540 nm. Data was evaluated in GraphPad Prism 6.
Results
Clinical characteristics of the patient cohort
A late PPROM cohort of 60 amniotic fluid samples with the presence or absence of MIAC and/or HCA was designed (Fig. 1A) and important clinical characteristics were compared among the subgroups (Table 1). Gestational age at delivery was considered as important clinical confounder of protein level due to a substantial increase in protein turnover during pregnancy27. Hence, the late PPROM cohort was carefully designed to include women in similar gestational age at admission and delivery. The median time of latency between late PPROM and amniocentesis was 4 h. We found a significant difference only between the subgroup with the presence of both MIAC and HCA (median = 7 h; IQR 4–11 h) compared to the subgroup with the absence of both MIAC and HCA (median = 3 h; IQR 2–6 h; p = 0.04). The median time of latency between amniocentesis and delivery was 10 h and the median time of latency between late PPROM and delivery was 16 h. There was not any significant difference in the duration between amniocentesis/late PPROM and delivery (p = 0.4 and p = 0.15, respectively). Based on these results, we infer that the changes of protein levels were not due to differences in the gestational age.
Tandem Mass Tag 10plex LC–MS/MS analysis and data treatment
The wide dynamic range of protein concentration represents a major obstacle in mass spectrometry-based protein analysis in human biological fluids such as plasma or amniotic fluid. Human serum albumin (HSA) is the most abundant protein in amniotic fluid and may smother information on important but far less abundant proteins28. The immunoaffinity depletion of abundant ballast proteins allowed us to detect a total of 1132 amniotic fluid proteins. The protein quantity was expressed as the ratio of TMT reporter ion intensity of each identified protein to the TMT reporter ion intensity of the global internal standard (GIS) in the respective multiplex (Fig. 1B). This allowed us to determine protein level changes in all 60 samples by the implementation of ten isobaric tags. Subsequent filtration of protein contaminants and proteins with a high number of missing values reduced the dataset to 816 quantified amniotic fluid proteins. Out of those, a total of 716 proteins were quantified in at least two from three technical replicates in all samples and 100 proteins were quantified in at least two technical replicates in eight samples in each late PPROM group (Fig. 1B). Thus, 72% of the identified proteins were kept for further consideration.
Comparative analysis of amniotic fluid proteins in the presence or absence of MIAC and/or HCA
The initial comparison was based on determination of the protein difference in the subgroups divided according to MIAC and/or HCA (the presence of both MIAC and HCA, MIAC alone, HCA alone, and the absence of both MIAC and HCA). Firstly, the multiple comparisons of protein levels were adjusted by Benjamin-Hochberg correction with the FDR threshold 0.05. However, no significant protein changes related to MIAC and/or HCA were found. The histone H1.5 (H1.5) was found as the first protein above the FDR threshold (FDR q value = 0.0625). Thus, violin plots based on the level of H1.5 in amniotic fluid were constructed to explain the uniformity of the late PPROM subgroups despite distinct MIAC and HCA status (Fig. 2A). Although the mean of the subgroup with the presence of both MIAC and HCA was higher than the mean of other subgroups, we observed a distinctly higher variance of the subgroup with the presence of both MIAC and HCA. Subsequently, a PCA plot of all individual samples was constructed and all 816 proteins were used as variables (Fig. 2B). The PCA plot also showed a distinct distribution of the subgroup with the presence of both MIAC and HCA compared to other late PPROM subgroups. The PCA plot illustrating the distribution of all individual samples revealed a cluster of five samples from the subgroup with both MIAC and HCA. The distribution of these five samples was dissimilar compared to the distribution of the remaining samples in the PCA plot (Fig. 2B). Violin and PCA plots together illustrated substantial heterogeneity of protein response to MIAC and HCA and a unique subgroup of five women was identified. Thus, we set those five samples as an individual group and conducted a subsequent analysis of variance of the following late PPROM subgroups: five outlying samples with the presence of both MIAC and HCA, the remaining ten samples with the presence of both MIAC and HCA, MIAC alone, HCA alone, or the absence of both MIAC and HCA. We applied the Benjamin-Hochberg correction with the FDR threshold 0.01 in this case.
The subsequent comparison by ANOVA and computed Tukey’s honestly significant difference uncovered that 138 proteins were significantly up- or downregulated in the five outlying samples compared to all other subgroups. In total, 112 proteins were upregulated and 26 proteins were downregulated (Supplementary table). To further characterize these proteins and explain the distinction of the five samples from the remaining samples with the presence of both MIAC and HCA, we conducted gene ontology analysis using the StringApp plugin25. All over-represented biological processes, cellular localization, and pathway terms are summed on Supplementary table and the most significant GO terms are depicted in Fig. 3A. Most of the upregulated proteins were localized in neutrophil granules (Fig. 3B). In terms of neutrophil activity, the protein upregulation was associated with platelet activation, glycolysis, the pentose phosphate pathway, and redox processes (Fig. 3C). By contrast, downregulated proteins were annotated as regulatory proteins (Fig. 3A).
Correlation of the amniotic fluid proteome with the level of interleukin-6, HCA and MIAC
Amniotic fluid IL-6 is currently the most relevant and discussed protein for evidence of intra-amniotic inflammation in PPROM29. Accordingly, we determined also the concentration of IL-6 in our samples and examined the correlation of IL-6 with the protein upregulation discovered through TMT 10plex LC–MS/MS analysis. We sought to assess whether proteome differences in the presence of MIAC and/or HCA correspond with the level of IL-6 in amniotic fluid. To this end, we constructed bubble plots as a correlation of H1.5 level with that of two other proteins, neutrophil defensin 3 (DEFA3) and protein S100A12. The size of bubbles are set according to the concentration of IL-6 and the color of bubbles represents the presence or absence of HCA and funisitis, the most severe form of HCA characterized by the presence of neutrophils in the umbilical cord (Fig. 5A,B). Histone H1.5, DEFA3, and S100A12 were found as highly upregulated proteins in the five outlying samples with the presence of both MIAC and HCA (Supplementary table). Moreover, DEFA3 and S100A12 were classified as proteins enriched in the antimicrobial immune response, neutrophil degranulation, and cell secretion (Fig. 3B). DNA-associated histones are frequently discussed as forming neutrophil extracellular traps upon infectious stimuli30 (Fig. 4).
Amniotic fluid levels of H1.5, DEFA3 and S100A12 showed substantial enhancement in the case of the five samples with the presence of both MIAC and HCA and also positive correlation with the highest levels of IL-6 (Fig. 5). Together, funisitis was confirmed in these five samples (Fig. 5A,B). These five samples were also found to have a high microbial load of Ureaplasma spp. and M. hominis DNA (threshold cycle values < 22, responding quantity more than 105 copies/mL) or the presence of Haemophilus influenzae and Sneathia sanquinegens was confirmed (Fig. 5C,D and Table 2).
Validation of amniotic fluid histone H1.5
As aforementioned, amniotic fluid histone H1.5 level was found to be significantly higher in the five samples from the subgroup with the presence of MIAC and HCA (Supplementary table). Histone H1.5 is an important intracellular protein related to excessive inflammatory response31. Thus, we determined the amniotic fluid level of histone H1.5 in the subgroups with the presence or absence of MIAC and/or HCA using ELISA. All but six samples had the level of H1.5 on concentration below the limit of detection of the ELISA kit (Fig. 6). The H1.5 level of the five outlying samples from the subgroup with both MIAC and HCA was in the range from 1.44 to 4.56 ng/mL. The sixth sample with H1.5 level of 2.11 ng/mL belonged to the group with MIAC alone but the IL-6 level of this sample was more than 15,000 pg/mL despite the absence of chorioamnionitis or funisitis (Fig. 6).
Discussion
In this study, we utilized a quantitative proteomic approach based on isobaric labeling of peptides to reveal proteome differences related to the MIAC and HCA in women with late PPROM. The principal findings of this study are: (1) The levels of amniotic fluid proteins in the subgroup with both MIAC and HCA were surprisingly uniform compared to the other late PPROM subgroups where both conditions weren’t confirmed. We found only five cases of late PPROM women with a substantial protein dysregulation related to both MIAC and HCA. We speculate that particularly these women can profit from an active management beyond the 34th week of gestation. (2) Proteome change was associated with a high concentration of IL-6 in amniotic fluid. (3) The positive correlation of proteomic findings with IL-6 concentration identified five women with strong intra-amniotic inflammation in late PPROM.
In our previous reports, LC–MS/MS enabled us to reveal changes in amniotic fluid proteome associated with the presence of MIAC and HCA in PPROM from 24 to 36th week of gestation17,20,32. However, to the best of our knowledge, a global proteomic study focused on PPROM beyond the 34th week has not been carried out yet. In TMT 10plex-based quantification, peptides of each sample are tagged with one different isotopic variant from a set of ten mass tags16,33. However, since the number of patients in our late PPROM cohort exceeded the number of TMT channels, we designed seven multiplexes and analyzed them in different LC–MS/MS runs (Fig. 1A). Normalization of mass tag information to the GIS removed the eventual intensity fluctuation during LC–MS/MS analysis of a particular multiplex, e.g. due to instability of electrospray (Fig. 1B).
Protein abundance may vary by orders of magnitude in complex samples. To make this variance similar across the global amniotic fluid proteome, the data were log2 transformed and normalized to the median34. Also, a high number of missing values may debase mass tag-based quantification and provide incorrect information about protein levels35. Perseus software replaces the missing values by random numbers that are estimated from normal distributions23, but that would mean working with a dataset containing artificial values. Hence, we rather decided to exclude proteins with a high number of missing values in technical replicates of analyzed samples (Fig. 1B).
This study was inspired by published results which indicated that a change of amniotic fluid proteome is associated with the presence of both MIAC and HCA in PPROM20,32. The authors applied a CysTRAQ approach, consisting of iTRAQ labeling to quantify selectively the fraction of peptides with cysteine residue and the fraction of non-cysteinyl peptides and found more than 80 amniotic fluid proteins significantly dysregulated in the presence of MIAC and HCA in PPROM pregnancies. Histone H4 together with other antimicrobial proteins (cathelicidin, neutrophil defensin 1, myeloperoxidase, etc.) and proteins integrated in neutrophil extracellular traps (myosin-9, plastin-2, matrix metalloproteinase-9, etc.) were shown to be considerably upregulated in the presence of both MIAC and HCA. Besides the CysTRAQ article, recent papers have discussed the elevation of IL-6 in the presence of both MIAC and HCA in cases of both PPROM and late PPROM13,36. The authors also concluded that the presence of both MIAC and HCA and a higher level of IL-6 were related to early-onset sepsis and neonatal morbidity.
In the CysTRAQ study, a pooled sample strategy was used. The pooling of samples is a time-efficient strategy compensating for limited amounts of the samples or high biological variation37. However, protein inflammatory response to MIAC and HCA is not a homogeneous phenomenon and a pooling strategy may disguise the variance of inflammatory protein concentration in individual samples. Keeping in mind the heterogeneity of MIAC and HCA, we required information about levels of proteins in each individual sample of the patient cohort. Moreover, we included also the subgroups with MIAC alone and HCA alone. Histological chorioamnionitis itself may evoke the pro-inflammatory protein response, but only the presence of both MIAC and HCA is apparently related to a higher number of cases of early-onset sepsis and fetal inflammatory response syndrome13,38. Hence, it is desirable to distinguish women with both MIAC and HCA from those with HCA alone.
Analysis of data achieved through TMT 10plex LC–MS/MS and correlation of these data with the IL-6 concentration in amniotic fluid, diagnosed funisitis, and MIAC revealed the substantial disparity in the subgroup with the presence of both MIAC and HCA. Dysregulation of amniotic fluid proteins occurred only in five cases in the subgroup with the presence of both MIAC and HCA whereas no significant distinction of protein level was observed in the remaining 55 entities regardless of MIAC and/or HCA (Fig. 2). These five samples were also affected by a high DNA load of Ureaplasma spp. and/or M. hominis and/or the presence of Haemophilus influenzae and Sneathia sanquinegens (Fig. 5C,D, Table 2). Nevertheless, some of the samples from women with MIAC alone with a low level of protein determined by LC–MS/MS as well as a low level of IL-6 were also affected by a high load of bacterial DNA (Fig. 5C,D). This confirmed previous findings that MIAC can occur as isolated bacterial colonization without any inflammation-associated change of proteome2. These five “outliers” were also diagnosed as pregnancies with funisitis (Fig. 5A,B). Funisitis is defined as the presence of neutrophils in the wall of umbilical vessels and/or Wharton’s jelly39. It seems that these five cases of the subgroup with both MIAC and HCA are the most affected cases with potential severe fetal inflammatory involvement.
A total of 138 proteins were significantly up- or downregulated in the five samples of the group with the presence of both MIAC and HCA compared to the other samples from women with late PPROM (Supplementary table). Substantial portion of upregulated proteins was annotated as secreted proteins present in neutrophil granules, proteins participating in glycolysis and the pentose phosphate pathway, in the regulation of cytoskeletal actin, and redox process (Fig. 3). The downregulated proteins were classified as regulators of acute inflammatory response. Neutrophil degranulation and secretion is host defense mechanism against potential intra-amniotic infection40. Upon cell activation, antimicrobial proteins such as myeloperoxidase (MPO), DEFA3, cathelicidin are released from neutrophil granules20,41. The lumen of neutrophil granules was the most over-represented cellular component in our GO analysis (Fig. 3A).
Besides the granular proteins mentioned above, we found other 34 proteins present in neutrophil granules and most of them were highly upregulated in the five discussed samples (Fig. 3B). Upregulated calcium-binding S100 proteins are small secreted proteins crucial for neutrophil migration and they may act as alarmins in the pathology of funisitis and early neonatal sepsis42. In addition, protein S100-A12 has a function as a ligand for the receptor for advanced glycation end‐products (RAGE)43. The ligand-RAGE interplay is a key mechanism of reactive oxygen species (ROS) generation and oxidative burst44. Other target molecule of S100 proteins is glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a key enzyme of glycolysis45. Glycolysis and the pentose phosphate pathway (PPP) are the key metabolic pathways in neutrophils and we found both of them as significant pathways in GO enrichment (Fig. 3A). We observed an expressive upregulation of enzymes participating in glycolysis and pyruvate generation in the five samples of the subgroup with the presence of both MIAC and HCA (Fig. 3C).
Besides energy generation, glycolysis and PPP directly affect also the formation of ROS by providing NADPH. A total of 19 proteins related to the redox process were upregulated in the five samples. We found upregulation of enzymes important for formation and elimination of ROS, including MPO, mitochondrial superoxide dismutase, catalase, glutaredoxin-1, peroxiredoxin-1, etc. (Supplementary table, Fig. 3C). Glycolysis, PPP and MPO- induced ROS may directly affect releasing of neutrophil extracellular traps (NETs)46,47,48. The presence of NETs in intra-amniotic infection and inflammation has been already established49. Unambiguous confirmation of NETs in amniotic fluid for instance using immunofluorescence microscopy50 was beyond the scope of this study. Nevertheless, we found alterations in the redox homeostasis affecting the neutrophil activity towards to NETs formation as well as the substantial upregulation of proteins being a part of NETs defense mechanisms such as MPO, DEFA3, MMP9 and histone isoforms representing strong indicia of the NETs presence (Fig. 4)51.
Histones are frequently discussed as key part of NETs and potential biomarkers of MIAC and HCA in PPROM20,32,49. In the five outlying samples from the subgroup with both MIAC and HCA, histones H1.5, H2B, H4 and H3.2 were highly upregulated. The level of H1.5 in amniotic fluid determined by ELISA also showed a significant increase in the five outlying samples with the presence of both MIAC and HCA compared to the other 55 samples (Fig. 6). Beyond the direct antimicrobial function, extracellular histones in released NETs bind platelets and induce platelet aggregation52. NETs have been even signed as a link between inflammation and thrombosis53. We found twenty up- or downregulated proteins related to platelet activation, signaling and aggregation (Fig. 3A, Supplementary table). From those, cofilin-1 and actinins were the most upregulated. Cofilin-1 is an actin-binding protein responsible for actin re-modelling preceding platelet aggregation54,55. On the other hand, another platelet protein thrombospondin-1 was downregulated in the five outlying samples with the presence of both MIAC and HCA compared to the other 55 samples. Thrombospondin-1 inhibits the activity of neutrophil elastase and cathepsin G and mitigates the inflammatory response56. We suggest that downregulation of thrombospondin-1 may promote the activity of neutrophil proteases as a consequence of excessive intra-amniotic inflammation in the specific subgroup with the presence of both MIAC and HCA in late PPROM.
Conclusion
In this study, we analyzed a unique collection of late PPROM amniotic fluid samples without a significant difference in gestational age at PPROM, amniocentesis, and delivery. We revealed the dysregulation of 138 proteins along with the high level of amniotic fluid IL-6, the presence of funisitis and MIAC. This study was initially aimed at discovery of protein changes associated with MIAC and HCA. We had presumed that the presence of neutrophils in the placenta, umbilical cord, and fetal membranes would be reflected by the change in the amniotic fluid proteome. However, we only detected five samples in the subgroup with the presence of both MIAC and HCA with substantial change in proteome composition. Consequently, the original stratification of the late PPROM cohort based on the status of HCA turned out to be a limitation of this study because the number of MIAC/HCA women with a strong proteome response was very low. However, these women can particularly profit from active management compared to other women with late PPROM. With respect to the achieved results, we recommend diagnosing the intra-amniotic inflammation rather based on the level of amniotic fluid IL-6 than as the diagnosed HCA. Nevertheless, the threshold values of intra-amniotic inflammation in PPROM vary from 2.6 to 11.3 ng/mL in previous publications57, and the conclusive cut-off value has not been established. Future research should thus unambiguously define the concentration of IL-6 representing the threshold for the confirmation of intra-amniotic inflammation.
Data availability
The mass spectrometry data together with the MaxQuant output files have been deposited in the ProteomeXchange Consortium via the PRIDE partner58 repository with the dataset identifier PXD020329.
References
Mercer, B. M. Preterm premature rupture of the membranes: current approaches to evaluation and management. Obstet. Gynecol. Clin. North Am. 32, 411–428 (2005).
Musilova, I. et al. Intraamniotic inflammation in women with preterm prelabor rupture of membranes. PLoS ONE 10, e0133929 (2015).
McIntire, D. D. & Leveno, K. J. Neonatal mortality and morbidity rates in late preterm births compared with births at term. Obstet. Gynecol 111, 35–41 (2008).
McGowan, J. E., Alderdice, F. A., Holmes, V. A. & Johnston, L. Early childhood development of late-preterm infants: a systematic review. Pediatrics 127, 1111–1124 (2011).
Morse, S. B., Zheng, H., Tang, Y. & Roth, J. Early school-age outcomes of late preterm infants. Pediatrics 123, e622-629 (2009).
Bond, D. M. et al. Planned early birth versus expectant management for women with preterm prelabour rupture of membranes prior to 37 weeks’ gestation for improving pregnancy outcome. Cochrane Database Syst. Rev.3, CD004735 (2017).
Quist-Nelson, J. et al. Immediate delivery compared with expectant management in late preterm prelabor rupture of membranes: an individual participant data meta-analysis. Obstet. Gynecol. 131, 269–279 (2018).
Committee on Practice Bulletins-Obstetrics. ACOG practice bulletin no. 188: prelabor rupture of membranes. Obstet. Gynecol. 131, e1–e14 (2018).
Tsakiridis, I., Mamopoulos, A., Chalkia-Prapa, E.-M., Athanasiadis, A. & Dagklis, T. Preterm premature rupture of membranes: a review of 3 national guidelines. Obstet. Gynecol. Surv. 73, 368–375 (2018).
DiGiulio, D. B. et al. Prevalence and diversity of microbes in the amniotic fluid, the fetal inflammatory response, and pregnancy outcome in women with preterm pre-labor rupture of membranes. Am. J. Reprod. Immunol. 64, 38–57 (2010).
Lee, S. M. et al. Acute histologic chorioamnionitis is a risk factor for adverse neonatal outcome in late preterm birth after preterm premature rupture of membranes. PLoS ONE 8, e79941 (2013).
Ofman, G., Vasco, N. & Cantey, J. B. Risk of early-onset sepsis following preterm, prolonged rupture of membranes with or without chorioamnionitis. Am. J. Perinatol. 33, 339–342 (2016).
Kacerovsky, M. et al. Prelabor rupture of membranes between 34 and 37 weeks: the intraamniotic inflammatory response and neonatal outcomes. Am. J. Obstet. Gynecol. 210, 325.e1–325.e10 (2014).
Musilova, I. et al. Late preterm prelabor rupture of fetal membranes: fetal inflammatory response and neonatal outcome. Pediatr. Res. 83, 630–637 (2018).
Michalski, A., Cox, J. & Mann, M. More than 100,000 detectable peptide species elute in single shotgun proteomics runs but the majority is inaccessible to data-dependent LC-MS/MS. J. Proteome Res. 10, 1785–1793 (2011).
Rauniyar, N. & Yates, J. R. Isobaric labeling-based relative quantification in shotgun proteomics. J. Proteome Res. 13, 5293–5309 (2014).
Vajrychova, M., Kacerovsky, M., Tambor, V., Hornychova, H. & Lenco, J. Microbial invasion and histological chorioamnionitis upregulate neutrophil-gelatinase associated lipocalin in preterm prelabor rupture of membranes. J. Matern. Fetal. Neonatal. Med. 29, 12–21 (2016).
Kacerovsky, M. et al. Lactobacilli-dominated cervical microbiota in women with preterm prelabor rupture of membranes. Pediatr. Res. 87, 952–960 (2020).
Salafia, C. M., Weigl, C. & Silberman, L. The prevalence and distribution of acute placental inflammation in uncomplicated term pregnancies. Obstet Gynecol 73, 383–389 (1989).
Tambor, V. et al. Amniotic fluid cathelicidin in PPROM pregnancies: from proteomic discovery to assessing its potential in inflammatory complications diagnosis. PLoS ONE 7, e41164 (2012).
Vajrychova, M. et al. Quantification of cellular protein and redox imbalance using SILAC-iodoTMT methodology. Redox Biol 24, 101227 (2019).
Werner, T. et al. Ion coalescence of neutron encoded TMT 10-plex reporter ions. Anal. Chem. 86, 3594–3601 (2014).
Tyanova, S. et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 13, 731–740 (2016).
R Core Team. R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, 2017).
Doncheva, N. T., Morris, J. H., Gorodkin, J. & Jensen, L. J. Cytoscape StringApp: network analysis and visualization of proteomics data. J. Proteome Res. 18, 623–632 (2019).
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003).
Thompson, G. N. & Halliday, D. Protein turnover in pregnancy. Eur J Clin Nutr 46, 411–417 (1992).
Perluigi, M. et al. Proteomic analysis for the study of amniotic fluid protein composition. J Prenat Med 3, 39–41 (2009).
Chaemsaithong, P. et al. A point of care test for interleukin-6 in amniotic fluid in preterm prelabor rupture of membranes: a step toward the early treatment of acute intra-amniotic inflammation/infection. J. Matern. Fetal. Neonatal. Med. 29, 360–367 (2016).
Knopf, J., Leppkes, M., Schett, G., Herrmann, M. & Muñoz, L. E. Aggregated NETs sequester and detoxify extracellular histones. Front .Immunol. 10, 2176 (2019).
Silk, E., Zhao, H., Weng, H. & Ma, D. The role of extracellular histone in organ injury. Cell Death Dis. 8, e2812 (2017).
Tambor, V. et al. CysTRAQ - A combination of iTRAQ and enrichment of cysteinyl peptides for uncovering and quantifying hidden proteomes. J. Proteomics https://doi.org/10.1016/j.jprot.2011.09.027 (2011).
Bai, B. et al. Deep profiling of proteome and phosphoproteome by isobaric labeling, extensive liquid chromatography, and mass spectrometry. Methods Enzymol. 585, 377–395 (2017).
Kammers, K., Cole, R. N., Tiengwe, C. & Ruczinski, I. Detecting significant changes in protein abundance. EuPA Open Proteomics 7, 11–19 (2015).
Palstrøm, N. B., Matthiesen, R. & Beck, H. C. Data imputation in merged isobaric labeling-based relative quantification datasets. Methods Mol. Biol. 2051, 297–308 (2020).
Kacerovsky, M. et al. Bedside assessment of amniotic fluid interleukin-6 in preterm prelabor rupture of membranes. Am. J. Obstet. Gynecol. 211(385), e1-9 (2014).
Kusonmano, K. et al. Effects of pooling samples on the performance of classification algorithms: a comparative study. ScientificWorldJournal 2012, 278352 (2012).
Kacerovsky, M. et al. The fetal inflammatory response in subgroups of women with preterm prelabor rupture of the membranes. J. Matern. Fetal. Neonatal. Med. 26, 795–801 (2013).
Pacora, P. et al. Funisitis and chorionic vasculitis: the histological counterpart of the fetal inflammatory response syndrome. J. Matern. Fetal. Neonatal. Med. 11, 18–25 (2002).
Gomez-Lopez, N. et al. Are amniotic fluid neutrophils in women with intraamniotic infection and/or inflammation of fetal or maternal origin?. Am. J. Obstet. Gynecol. 217(693), e1-693.e16 (2017).
Faurschou, M., Sørensen, O. E., Johnsen, A. H., Askaa, J. & Borregaard, N. Defensin-rich granules of human neutrophils: characterization of secretory properties. Biochim. Biophys. Acta 1591, 29–35 (2002).
Buhimschi, I. A. & Buhimschi, C. S. The role of proteomics in the diagnosis of chorioamnionitis and early-onset neonatal sepsis. Clin. Perinatol. 37, 355–374 (2010).
Leclerc, E. & Vetter, S. W. The role of S100 proteins and their receptor RAGE in pancreatic cancer. Biochim. Biophys. Acta 1852, 2706–2711 (2015).
Daffu, G. et al. Radical roles for RAGE in the pathogenesis of oxidative stress in cardiovascular diseases and beyond. Int. J. Mol. Sci. 14, 19891–19910 (2013).
Hatakeyama, T., Okada, M., Shimamoto, S., Kubota, Y. & Kobayashi, R. Identification of intracellular target proteins of the calcium-signaling protein S100A12. Eur. J. Biochem. 271, 3765–3775 (2004).
Azevedo, E. P. et al. A metabolic shift toward pentose phosphate pathway is necessary for amyloid fibril- and phorbol 12-myristate 13-acetate-induced neutrophil extracellular trap (NET) formation. J. Biol. Chem. 290, 22174–22183 (2015).
Rodríguez-Espinosa, O., Rojas-Espinosa, O., Moreno-Altamirano, M. M. B., López-Villegas, E. O. & Sánchez-García, F. J. Metabolic requirements for neutrophil extracellular traps formation. Immunology 145, 213–224 (2015).
Kirchner, T. et al. The impact of various reactive oxygen species on the formation of neutrophil extracellular traps. Mediators Inflamm. 2012, 849136 (2012).
Gomez-Lopez, N. et al. Neutrophil extracellular traps in the amniotic cavity of women with intra-amniotic infection: new mechanism of host defense. Reprod. Sci. 24, 1139–1153 (2017).
de Buhr, N. & von Köckritz-Blickwede, M. How neutrophil extracellular traps become visible. J. Immunol. Res. 2016, 4604713 (2016).
Dąbrowska, D., Jabłońska, E., Garley, M., Ratajczak-Wrona, W. & Iwaniuk, A. New aspects of the biology of neutrophil extracellular traps. Scand. J. Immunol. 84, 317–322 (2016).
Fuchs, T. A., Bhandari, A. A. & Wagner, D. D. Histones induce rapid and profound thrombocytopenia in mice. Blood 118, 3708–3714 (2011).
Fuchs, T. A. et al. Extracellular DNA traps promote thrombosis. Proc. Natl. Acad. Sci. U.S.A. 107, 15880–15885 (2010).
Garg, P. et al. Actin-depolymerizing factor cofilin-1 is necessary in maintaining mature podocyte architecture. J. Biol. Chem. 285, 22676–22688 (2010).
Shin, E.-K., Park, H., Noh, J.-Y., Lim, K.-M. & Chung, J.-H. Platelet shape changes and cytoskeleton dynamics as novel therapeutic targets for anti-thrombotic drugs. Biomol. Ther. (Seoul) 25, 223–230 (2017).
Zhao, Y. et al. Thrombospondin-1 restrains neutrophil granule serine protease function and regulates the innate immune response during Klebsiella pneumoniae infection. Mucosal Immunol. 8, 896–905 (2015).
Combs, C. A. et al. Amniotic fluid infection, inflammation, and colonization in preterm labor with intact membranes. Am. J. Obstet. Gynecol. 210(125), e1-125.e15 (2014).
Perez-Riverol, Y. et al. The PRIDE database and related tools and resources in 2019: improving support for quantification data. Nucleic Acids Res. 47, D442–D450 (2019).
Acknowledgements
This work was supported by PERSONMED (Reg. Nr. CZ.02.1.01/0.0/0.0/17_048/0007441) and MH CZ—DRO (UHHK, 00179906). We wish to thank Dr. Barbora Šalovská from Institute of Molecular Genetics of ASCR (Prague, Czech Republic) for her valued help with the design of R scripts. The authors are grateful to Ian McColl MD, PhD for assistance with the manuscript.
Author information
Authors and Affiliations
Contributions
M.V. carried out ELISA experiments, analyzed protein data, and prepared figures. M.V, J.S., M.B., and B.J. wrote the manuscript. M.K. and J.S. designed the patient cohort and collected amniotic fluid samples. K.P. and V.T. designed the TMT 10plex experiment, prepared samples, and carried out the LC–MS/MS analyses. R.K. analyzed the MIAC data. J.L. participated in data evaluation. P.Z. carried out immunoaffinity depletion. R. B. and L.P. carried out qPCR detection of bacteria in amniotic fluid. H.H. investigated HCA in fetal membranes, umbilical cord, and placenta. C.A. assisted with ELISA. M.K. and M.V. got funding. All authors contributed to the revision of the manuscript.
Corresponding authors
Ethics declarations
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.
Supplementary information
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/.
About this article
Cite this article
Vajrychová, M., Stráník, J., Pimková, K. et al. Comprehensive proteomic investigation of infectious and inflammatory changes in late preterm prelabour rupture of membranes. Sci Rep 10, 17696 (2020). https://doi.org/10.1038/s41598-020-74756-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-74756-9
This article is cited by
-
Transcriptomic analysis delineates preterm prelabor rupture of membranes from preterm labor in preterm fetal membranes
BMC Medical Genomics (2024)
-
Proteomic analysis of plasma to identify novel biomarkers for intra-amniotic infection and/or inflammation in preterm premature rupture of membranes
Scientific Reports (2023)
-
Proteomic profiling of human amnion for preterm birth biomarker discovery
Scientific Reports (2021)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.