Next Article in Journal
Elicitor-Driven Defense Mechanisms: Shielding Cotton Plants against the Onslaught of Cotton Leaf Curl Multan Virus (CLCuMuV) Disease
Next Article in Special Issue
Once upon a Time Oral Microbiota: A Cinderella or a Protagonist in Autism Spectrum Disorder?
Previous Article in Journal
Sphingolipids in Childhood Asthma and Obesity (SOAP Study): A Protocol of a Cross-Sectional Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Profiling of Second-Trimester Amniotic Fluid for Predicting Preterm Delivery: Insights from NMR Analysis

by
Charalampos Kolvatzis
1,*,
Paris Christodoulou
2,
Ioannis Kalogiannidis
1,
Konstantinos Tsiantas
2,
Ioannis Tsakiridis
1,
Charikleia Kyrkou
3,
Antigoni Cheilari
4,
Nikolaos S. Thomaidis
5,
Panagiotis Zoumpoulakis
2,
Apostolos Athanasiadis
1 and
Alexandra-Maria Michaelidou
3
1
Third Department of Obstetrics and Gynecology, School of Medicine, Faculty of Health Sciences, Aristotle University of Thessaloniki, 54642 Thessaloniki, Greece
2
Department of Food Science and Technology, University of West Attica, Ag. Spyridonos, 12243 Egaleo, Greece
3
Department of Food Science and Technology, Faculty of Agriculture, Forestry and Natural Environment, School of Agriculture, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
4
Department of Pharmacognosy and Natural Products Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, 17551 Athens, Greece
5
Laboratory of Analytical Chemistry, Department of Chemistry, National and Kapodistrian University of Athens, Panepistimiopolis Zografou, 15771 Athens, Greece
*
Author to whom correspondence should be addressed.
Metabolites 2023, 13(11), 1147; https://doi.org/10.3390/metabo13111147
Submission received: 23 October 2023 / Revised: 9 November 2023 / Accepted: 9 November 2023 / Published: 12 November 2023
(This article belongs to the Special Issue Preclinical and Clinical Application of Metabolomics in Medicine)

Abstract

:
Preterm delivery (PTD) is a notable pregnancy complication, affecting one out of every ten births. This study set out to investigate whether analyzing the metabolic composition of amniotic fluid (AF) collected from pregnant women during the second trimester of pregnancy could offer valuable insights into prematurity. The research employed 1H–NMR metabolomics to examine AF samples obtained from 17 women who gave birth prematurely (between 29+0 and 36+5 weeks of gestation) and 43 women who delivered at full term. The application of multivariate analysis revealed metabolites (dimethylglycine, glucose, myo-inositol, and succinate) that can serve as possible biomarkers for the prognosis and early diagnosis of preterm delivery. Additionally, pathway analysis unveiled the most critical metabolic pathways relevant to our research hypothesis. In summary, these findings suggest that the metabolic composition of AF in the second trimester can be a potential indicator for identifying biomarkers associated with the risk of PTD.

1. Introduction

Spontaneous preterm delivery (PTD), affecting one out of every ten births, is recognized as a syndrome influenced by multiple contributing factors [1]. Among the spectrum of suspected causes of PTD, infection and/or inflammation characterized as the body’s response to signals of microbial or non-microbial danger stand out as the only pathological processes for which a confirmed causal connection with PTD has been established, along with a clearly defined molecular pathophysiology [2].
Infants born prematurely, particularly those born before 34 weeks of gestation, have an elevated risk of mortality and health problems. Furthermore, infants born during the late preterm period, i.e., 34–37 weeks of gestation, face increased health complications and a higher probability of developing health conditions like obesity, metabolic syndrome, hypertension, and type 2 diabetes later in life [3].
Prediction and early diagnosis of PTD are often challenging because of their complexity [4]. Hence, it is not surprising that metabolomics, utilizing advanced techniques such as Nuclear Magnetic Resonance Spectroscopy (NMR), Gas Chromatography–Mass Spectrometry (GC-MS), and Liquid Chromatography–Mass Spectrometry (LC-MS/MS), have gained prominence in identifying biomarkers and pathways crucial to the development and progression of this syndrome [4,5,6,7,8,9,10].
Among the various options, amniotic fluid (AF) emerges as a particularly promising biofluid; it serves as a dynamic repository, reflecting the metabolic profile of the developing fetus. Research studies [11,12,13,14,15] suggest that conducting a metabolomic analysis on AF obtained during the prenatal period holds the potential to identify metabolic deviations before PTD occurs. However, a recent study using untargeted LC-MS mid-trimester amniotic fluid metabolic profiling on two groups of 37 pregnant women (full term/preterm) that underwent amniocentesis showed no evidence of metabolite differentiation for spontaneous PTD [16].
Studies in this area have employed diverse metabolomics techniques, and the consensus on critical metabolites serving as biomarkers remains elusive [10]. The apparent inconsistency highlights the intricate nature of the PTD syndrome and implies the presence of unexplored aspects. Within this context, our study aims to bridge this gap by identifying potential predictive biomarkers for spontaneous PTD within second-trimester AF using NMR analysis.

2. Materials and Methods

2.1. Study Design and Population

This prospective cohort study investigated AF samples between 2013 and 2014 at the Third Department of Obstetrics and Gynecology, School of Medicine, Faculty of Health Sciences, Aristotle University of Thessaloniki, Greece. The study was approved by the Aristotle University of Thessaloniki Research Ethics Committee (Prot. No. 1.662/21 November 2018) and conducted in compliance with the declaration of Helsinki. No incentives were provided, and all the participants signed a relevant informed written consent.

Data Collection and Eligibility Criteria

All the pregnant women in singleton pregnancies who underwent amniocentesis (between 16 and 22 weeks of gestation) for prenatal screening were eligible to participate in the study; a small amount of AF was donated and kept at −80 °C. Detailed obstetric and medical histories were recorded for each woman in the study.
Exclusion criteria were (1) multiple pregnancies, (2) short cervical length at second-trimester ultrasound (<25 mm), and (3) preterm prelabor rupture of membranes, placental ischemic disease, and iatrogenic PTD. In addition, pregnant women who had a miscarriage before 24 weeks or PTD within the first 3 weeks of amniocentesis were excluded from the analysis.
With regards to the indications of amniocentesis, increased risk of fetal genetic abnormality after combined first-trimester ultrasound and biochemical screening, detection of ultrasound markers at first- and second-trimester screenings associated with an increased likelihood of chromosomal abnormalities (e.g., nasal bone hypoplasia, short femur, and hyperechoic bowel), maternal desire, and high risk of vertical transmission of congenital infection were the most common ones, as recommended by the majority of the guidelines [17].
Pregnancy outcomes were collected to determine the gestational age at delivery (preterm <37 weeks or term >37 weeks), mode of delivery (vaginal/cesarean), birthweight, and possible complications during pregnancy. The metabolomic profile of the AF was checked and compared according to the week of delivery and birthweight.

2.2. NMR Metabolomics Analysis

2.2.1. Sample Preparation

Amniotic fluid samples were thawed at room temperature (25 °C) and extracted according to a common method established in the literature [13]. Specifically, 10 mL of AF was centrifuged (14,000 rpm, 4 °C, 10 min), and 1 mL of the supernatants was lyophilized overnight until dry. The dry residues were then reconstituted into 540 µL of phosphate buffer (0.2 M, Na2HPO4 2H2O, and NaH2PO4, pH = 7.0) in D2O and 60 µL (5 mM) of d6- trimethylsilyl propionic acid sodium salt (TSP) as the internal standard.

2.2.2. 1H–NMR Analysis

For 1H–NMR measurement, 600 μL of the final samples was transferred to 5 mm NMR tubes (LabScape, Bruker, Germany). The NMR spectra were acquired at 300 K, after a 5 min resting period for temperature stabilization, on a Bruker Ascend 500 MHz NMR spectrometer equipped with a 5 mm double resonance broadband inverse (BBI) detection probe. Experiments were performed in automation mode, using a SampleCase-24 sample changer operated by IconNMR. Data acquisition and processing were performed with TopSpin 4.1.4 (Bruker Italia Slr, Milan, Italy). Metabolic profiling 1D NMR spectra were acquired using water suppression. T2-edited Carr-Purcell-Meiboom-Gill (CPMG) experiments were acquired with d1 = 6 s; AQ = 4.92 s; FID data points = 96 k; SW = 20 ppm; ns = 32. The transmitter offset was set manually to achieve optimal suppression of the residual water signal for both experiments. FIDs were zero-filled and multiplied by an exponential weighting function corresponding to a line broadening of 0.3 Hz before Fourier transformation. Chemical shift values were referenced to the residual TSP signal (0.00 ppm).

2.2.3. Data Processing

All 1H–NMR spectra were phase-corrected using TopSpin 4.1.4 software. Then, the .zip file of all files was uploaded to the NMRProcFlow open-access web tool [18] to proceed with chemical shift calibration, baseline correction, and S/N ratio identification as well as alignment normalization and bucketing. For alignment, an interactive philosophy was performed, meaning that each interval had been chosen separately performing CluPA, as well as the least squares method [2]. All spectra were normalized using the Probabilistic Quotient Normalization (PQN) method [3], and an intelligent bucketing module was performed within the 0.8–8.5 spectra region.

2.2.4. Metabolites Screening

A series of 2D TOCSY and HSQC experiments were acquired for metabolite identification. More specifically, the above validation NMR experiments were acquired using a phase-sensitive MLEV sequence with d1 = 4 s; FID data points = 2 k (F2) and 256 (F1); SW = 20 ppm; ns = 32; and mixing time (d9) = 0.08 sec with suppression of the residual water signal. Furthermore, phase-sensitive HSQC-DEPT experiments were acquired using Echo/Antiecho-TPPI gradient selection with decoupling during acquisition (hsqcedetgpsisp2.3) with FID data points = 4 k (F2) and 288 (F1); SW = 12 ppm (F2) and 180 (F1); and ns = 160 in non-uniform sampling (NUS) acquisition mode with a NUS level of 50%.. To aid in rapid and efficient metabolite identification, Metabominer [19], an easy-to-use software tool, along with 2D TOCSY and HSQC experiments and data from the literature were used. Briefly, an automated pick-picking list from each spectrum (TOCSY and HSQC) was generated from MestreNOVA software and loaded in the Metabominer tool to (a) screen metabolites from Metabominers’ biofluids database and (b) refine the identified compounds by superimposing the spectral images to the Metabominers’ generated spectra.

2.3. Statistical Analysis

2.3.1. Demographics

As all demographic parameters did not follow the normal distribution, the Mann–Whitney U test was used for group comparisons. The statistical significance level was set at 5%. The demographics were statistically analyzed with SPSS (Statistical Package for the Social Sciences) v.26 (SPSS, Inc., Chicago, IL, USA).

2.3.2. Metabolomics

The MetaboAnalyst 5.0 was facilitated for both univariate/multivariate analyses and pathway analysis [5]. For multivariate analysis, -Pareto scaling was used, and unsupervised/supervised models (PCA/PLS-DA) were extracted at a confidence level of 95%. Feature selection of the PLS-DA model was based on variable importance projection (VIP) scores > 1.0 to reveal the variable that mostly contributes to the discrimination of the studied groups. The validation of the models was evaluated in terms of accuracy (>90.0 %), the goodness-of-fit R2 (0 ≤ R2 ≤ 1), and the predictive ability Q2 (0 ≤ Q2 ≤ 1) values. The PLS-DA results were further cross-validated by carrying out permutation tests with 1000 random permutations.
In terms of univariate analysis, non-parametric Wilcoxon t-tests were performed (p < 0.05). For the elucidation of metabolites as possible biomarkers, Receiver Operating Characteristic (ROC) curves were acquired. Finally, pathway analysis was performed by Metaboanalyst 5.0 platform.

3. Results

In total, 60 AF samples were used in the analyses; 43 were term and 17 were preterm. The mean duration of gestation was 38.5 (±0.9) weeks for the term group and 35.3 (±2.8) weeks for the preterm group (p < 0.001). Similarly, the mean birthweight was statistically different with term neonates having a mean birthweight of 3340 g (±380.77), while the preterm ones weighed 2695 g (±553.85) (p < 0.001). All parameters investigated are presented in Table 1.

3.1. NMR Analysis

A total of 27 metabolites including aliphatic and aromatic amino acids, sugars, and organic acids were assigned by the combination of 2D NMR spectroscopy (2D-HSQC, 2D TOCSY) and available online tools such as Metabominer and the literature data, as described in Section 2.2.4. The identified metabolites of the acquired AF samples are displayed in Figure 1, while the characteristic chemical shifts are represented in Supplementary Table S1.
The spectra matrix underwent intelligent bucketing to reduce the spectra data and facilitate subsequent statistical analysis. A total of 170 buckets were annotated according to the assignment procedure (Supplementary Table S2). The new matrix was subjected to both multivariate and univariate analysis.

3.2. Statistics

3.2.1. Overview of the Studied Samples

Principal Component Analysis (PCA) was conducted on the annotated NMR profile (170 annotated intelligent buckets, Supplementary Table S2) to gain an overview of the sample space. The resulting scores plot (Figure 2) revealed a pattern between the two studied groups along the first principal component, accounting for 34.3% of the metabolic variance in the studied AF samples.

3.2.2. Data Reduction Method for Unique Potential Biomarker Discovery

Digging deeper into the annotated metabolites (Supplementary Table S2), it is observed that a significant number of buckets correspond to more than one metabolite. Considering that the purpose of this study is the development of reliable potential prognostic markers for PTD, we applied a robust statistical methodology based on the annotated spectra buckets.
In detail, a gradual bucketing reduction method was implemented by including buckets that characterize unique metabolites to avoid the metabolites’ overlapping phenomenon. This procedure led to the reduction of the feature space to 94 out of 170 annotated intelligent buckets. Notably, the 94 buckets corresponded to the initial number of assigned metabolites (n = 27).
Moreover, the reduced feature space (n = 94) was subjected to biomarker analysis considering the area under the ROC curve (AUROC) and the p values promoting potential biomarkers with a good predictive ability and strong statistical significance (AUROC > 0.75 and p < 0.05). These values are in accordance with the literature [10,20].
Applying the above criteria, the process concluded with 25 intelligent buckets corresponding to 12 unique metabolites (Supplementary Table S3). Finally, to exclude multiple buckets corresponding to the same metabolite, buckets that demonstrated the higher AUROC (n = 12) were kept for further statistical analyses (Supplementary Table S4). The whole procedure is shown in Figure 3.

3.2.3. Discriminant and Pathway Analysis

Discriminant analysis was applied to define the metabolites that confirm the discrimination of the two sample categories (full term/preterm). Supervised partial least-squares discriminant analysis (PLS-DA) was employed, using the full-term/preterm classification as the response variable and the 12 potential biomarkers as the independent variables. The extracted 3D score plot manifested the separation of the studied groups across the first component (Figure 4A). The variable importance projection (VIP) plot was then retrieved to identify the variables with the greater discriminative ability (VIP > 1) among the two studied groups. Indicatively, PTB is associated with a lower abundance of dimethylglycine, glucose, and myo-inositol; while in the full-term birth, the production of succinate is fostered (Figure 4B). The model was validated by permutation test statistics (Figure 4C). The box plots of the four promoted metabolites are shown in Figure 5.
Finally, metabolites exhibiting AUROC > 0.75 (Supplementary Table S4) in AF samples were subjected to pathway analysis to relate the promoted metabolic patterns to specific pathways. The results (Supplementary Figure S1) depicted that seven metabolic pathways were significantly enriched (p < 0.05) containing at least two compounds, while two of them (alanine, aspartate, and glutamate metabolism and the citrate cycle) had the largest impact (>0.1) (Table 2).

4. Discussion

The present study explores the possibility of identifying potential concentration patterns of markers that characterize PTD. Among these 60 asymptomatic women, 43 delivered at term, while 17 delivered prematurely.
Our study supports the hypothesis that PTD is characterized by several changes in the metabolic profile of the fetus as reflected in second-trimester amniotic fluid [10]. The multivariate analysis highlighted that dimethylglycine, glucose, myo-inositol and succinate possessed the highest discriminative ability between the two studied groups. Moreover, pathway analysis revealed the most important metabolic pathways involved in our research hypothesis.
In accordance with prior research results [13], we have observed decreased glucose levels in second-trimester amniotic fluid samples obtained from women that delivered preterm. The decrease in glucose levels may be linked to increased glycolysis, potentially occurring under stressful conditions and reduced utilization of the respiratory chain pathway (as indicated by the rise in succinate levels) partly due to disrupted transplacental flow [15]. Furthermore, a meta-analysis conducted by Liu et al. (2017) supported our observations, indicating that lower glucose levels in early- or mid-trimester AF are associated with PTD [21]. Previous studies have convincingly demonstrated that decreased levels of glucose in AF are linked to the presence of microorganisms in the amniotic cavity and/or an inflammatory response [21,22,23,24]. These findings align with the notion that intraamniotic infection or inflammation (subclinical) plays a causative role in PTD [21].
Our analysis revealed that during the second trimester, women who later experienced PTD exhibited higher levels of succinate in their AF compared to women who carried their pregnancies to term. This discovery aligns with a recent study conducted by Virgiliou et al. in 2017 [14], which observed elevated fumarate levels in mid-trimester AF derived from women who delivered prematurely. Ansari et al. [22] and others [25,26,27] have proposed that succinate and fumarate, vital metabolites in both host and microbial processes, act as potent allosteric inhibitors of 2OG-dependent dioxygenases. Notably, this group of enzymes includes members of the histone demethylase family. Emerging evidence highlights a strong correlation between epigenetic events, particularly histone methylation, T-cell activation, differentiation, and commitment [22,25,26]. Hence, it is plausible that the generation, utilization, and translocation of these metabolites from the mitochondria to the cytosol, facilitating dynamic histone methylation within the nucleus, could potentially contribute to the occurrence of PTD. Given the immunomodulatory effects of succinate, as demonstrated by Al-Mushrif et al. in 2000 albeit in a slightly different setting [23], this finding needs further investigation.
Another noteworthy discovery that aligns with the existing literature [14] pertains to the reduced levels of myo-inositol found in the AF of PTD cases. Myo-inositol is a vital nutrient essential for the growth and viability of human cells. It is worth mentioning that myo-inositol has previously been identified as a significant marker in AF for women experiencing preterm labor without intraamniotic infection or inflammation [14,28,29]. Numerous assumptions may document myo-inositol’s potential role in fetal development and metabolic maturity. Myo-inositol acts as a structural basis for several secondary messengers found in eukaryotic cells, including inositol phosphates, phosphatidylinositol, and phosphatidylinositol phosphate lipids. These secondary messengers are pivotal in various cellular processes, including the regulation of intracellular calcium, gene expression, and lipid metabolism [30]. Current data suggest that a higher content of placental myo-inositol might postpone the initiation of labor by suppressing the synthesis and production of placental eicosanoids. This suppression could result from either a decrease in the availability of arachidonic acid for eicosanoid synthesis or a reduction in the activity of enzymes responsible for eicosanoid production. Moreover, myo-inositol could potentially redirect arachidonic acid metabolism towards generating anti-inflammatory eicosanoids rather than pro-inflammatory ones, thereby regulating the overall equilibrium of pro-inflammatory and anti-inflammatory factors to inhibit the onset of preterm labor. Furthermore, the administration of myo-inositol for the prevention of prematurity has been a subject of study [26].
Dimethylglycine is a derivative of the amino acid glycine. It is important to note that glycine, choline, betaine, and dimethylglycine are interconnected compounds with related metabolic pathways. While a previous metabolomic analysis of healthy pregnant individuals’ AF indicated increased levels of both choline and N-dimethylglycine during the transition from the second to the third trimester [31], there is limited information available regarding dimethylglycine levels in AF among cases of PTD. In the context of our study, dimethylglycine was found to be decreased in AF samples from women who delivered prematurely. Deciphering the role of dimethylglycine in relation to PTD is interesting since current data suggest that dimethylglycine may be engaging in different pathways that influence the regulation of the labor process itself. These effects might be related, at least in part, to its function in bolstering cellular antioxidant activity. It is well-documented that spontaneous PTD is associated with oxidative stress and disturbances in the body’s redox system, often linked with inflammation [14]. Dimethylglycine is involved in multiple metabolic pathways and can facilitate the recycling of other antioxidants, such as glutathione, which is a vital cellular antioxidant [32]. By aiding in the regeneration of other antioxidants, dimethylglycine could indirectly contribute to mitigating oxidative stress and safeguarding cells from damage caused by reactive oxygen species (ROS). Another viable scenario involves its role as a methyl donor. The transfer of a methyl group from betaine to homocysteine, catalyzed by the enzyme betaine-homocysteine methyltransferase, results in the production of methionine and dimethylglycine. Consequently, a deficiency in methyl donors may result in disturbances in metabolism and function [24].
Interestingly, as previously discussed, a recent study [16] found no significant differences in metabolite profiles associated with spontaneous PTD. This outcome may be attributed to variations in stratification and research protocols.
The strengths of this study are evident in its approach to data collection and the rigorous statistical methods applied. The collection of samples during the pre-clinical phase, when women show no symptoms, has the potential to identify individuals at high risk early on. This approach aligns with the findings of Souza et al. in 2019 [4], suggesting that it can help uncover the triggers of preterm labor. Furthermore, a meticulously statistical approach was adopted to thoroughly assess the candidate biomarkers’ diagnostic and prognostic potential. This method yielded a concise set of features that represented the most reliable and noteworthy markers, characterized by AUROC values exceeding 0.75. Regarding limitations, the study’s sample size may pose constraints on the broader applicability of its findings. Additionally, for future investigations, the integration of metabolomic techniques with advanced molecular microbiological methods could offer a deeper exploration of the intricate connections between microbiota and metabolites.

5. Conclusions

To conclude, NMR metabolomics’ analysis led to the identification of glucose, succinate, myo-inositol, and dimethylglycine as potential biomarkers of PTD. The results overall support previous data, while they underscore the dynamic nature of AF composition, reflecting the physiological changes occurring throughout pregnancy. The identified specific metabolites and their fluctuations have the potential to serve as biomarkers for monitoring pregnancy progression and identifying aberrations that may require medical attention. Furthermore, our study has highlighted the intricate interplay between maternal and fetal metabolisms within the AF environment. This intergenerational metabolic dialogue can inform our understanding of fetal development and its susceptibility to external factors.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo13111147/s1, Table S1. Table representing the characteristic chemical shifts of the identified metabolites, Table S2. Table representing the annotated buckets based on the identification procedure. The term “BX_XXXX” corresponds to ppm. Table S3. Illustration of Biomarker Analysis. In total 50 annotated buckets corresponding to unique metabolites, passed the statistical significance analysis (p < 0.05). In bold are presented the 25 annotated buckets with AUROC > 0.75 and p value < 0.05. Table S4. Buckets corresponding to the higher AUROCs and 12 unique metabolites. Figure S1: Graphical illustration of Pathway Analysis.

Author Contributions

Conceptualization, A.A., I.K., P.Z., C.K. (Charalampos Kolvatzis) and A.-M.M.; methodology, P.C., K.T., A.C., N.S.T., P.Z., A.A. and C.K. (Charalampos Kolvatzis); software P.C., K.T. and A.C.; investigation, C.K. (Charalampos Kolvatzis) and P.C.; data curation, C.K. (Charalampos Kolvatzis), P.C., K.T. and A.C.; resources, P.Z., N.S.T., A.A. and C.K. (Charalampos Kolvatzis); visualization, C.K. (Charikleia Kyrkou) and P.C.; writing—original draft preparation, C.K. (Charalampos Kolvatzis), P.C., K.T., A.C., I.T., A.-M.M. and C.K. (Charikleia Kyrkou); writing—review and editing, C.K. (Charalampos Kolvatzis), I.T., P.C., K.T., N.S.T., C.K. (Charikleia Kyrkou), A.A., I.K. and A.-M.M.; supervision, A.A., I.K., P.Z. and A.-M.M.; project administration, P.Z., A.A. and A.-M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Aristotle University of Thessaloniki Research Ethics Committee (Prot. No. 1.662/21 November 2018) for studies involving humans.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in the manuscript are available on reasonable request from the corresponding authors. The data are not publicly available due to their containing information that could compromise the privacy of research participants.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goldenberg, R.L.; Culhane, J.F.; Iams, J.D.; Romero, R. Epidemiology and causes of preterm birth. Lancet 2008, 371, 75–84. [Google Scholar] [CrossRef] [PubMed]
  2. Vu, T.N.; Laukens, K. Getting Your Peaks in Line: A Review of Alignment Methods for NMR Spectral Data. Metabolites 2013, 3, 259–276. [Google Scholar] [CrossRef] [PubMed]
  3. Wojtowicz, W.; Zabek, A.; Deja, S.; Dawiskiba, T.; Pawelka, D.; Glod, M.; Balcerzak, W.; Mlynarz, P. Serum and urine 1H NMR-based metabolomics in the diagnosis of selected thyroid diseases. Sci. Rep. 2017, 7, 9108. [Google Scholar] [CrossRef] [PubMed]
  4. Souza, R.T.; Mayrink, J.; Leite, D.F.; Costa, M.L.; Calderon, I.M.; Rocha, E.A.; Vettorazzi, J.; Feitosa, F.E.; Cecatti, J.G. Metabolomics applied to maternal and perinatal health: A review of new frontiers with a translation potential. Clinics 2019, 74, e894. [Google Scholar] [CrossRef] [PubMed]
  5. Pang, Z.; Chong, J.; Zhou, G.; de Lima Morais, D.A.; Chang, L.; Barrette, M.; Gauthier, C.; Jacques, P.-É.; Li, S.; Xia, J. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 2021, 49, W388–W396. [Google Scholar] [CrossRef] [PubMed]
  6. Gil, A.M.; Duarte, D. Biofluid metabolomics in preterm birth research. Reprod. Sci. 2018, 25, 967–977. [Google Scholar] [CrossRef] [PubMed]
  7. Ansari, A.; Bose, S.; You, Y.; Park, S.; Kim, Y. Molecular mechanism of microbiota metabolites in preterm birth: Pathological and therapeutic insights. Int. J. Mol. Sci. 2021, 22, 8145. [Google Scholar] [CrossRef]
  8. Li, X.; Huang, S.; Jiao, A.; Yang, X.; Yun, J.; Wang, Y.; Xue, X.; Chu, Y.; Liu, F.; Liu, Y. Association between ambient fine particulate matter and preterm birth or term low birth weight: An updated systematic review and meta-analysis. Environ. Pollut. 2017, 227, 596–605. [Google Scholar] [CrossRef]
  9. Michaelidou, A.-M.; Tsakoumaki, F.; Fotiou, M.; Kyrkou, C.; Athanasiadis, A.P. Metabolomics in amniotic fluid. In New Technologies and Perinatal Medicine; CRC Press: Boca Raton, FL, USA, 2019; pp. 139–148. [Google Scholar]
  10. Kolvatzis, C.; Tsakiridis, I.; Kalogiannidis, I.A.; Tsakoumaki, F.; Kyrkou, C.; Dagklis, T.; Daniilidis, A.; Michaelidou, A.-M.; Athanasiadis, A. Utilizing Amniotic Fluid Metabolomics to Monitor Fetal Well-Being: A Narrative Review of the Literature. Cureus 2023, 15, e36986. [Google Scholar] [CrossRef]
  11. Graça, G.; Diaz, S.O.; Pinto, J.; Barros, A.S.; Duarte, I.F.; Goodfellow, B.J.; Galhano, E.; Pita, C.; Almeida, M.d.C.; Carreira, I.M. Can biofluids metabolic profiling help to improve healthcare during pregnancy? Spectrosc. Int. J. 2012, 27, 515–523. [Google Scholar] [CrossRef]
  12. Graca, G.; Moreira, A.S.; Correia, A.J.V.; Goodfellow, B.J.; Barros, A.S.; Duarte, I.F.; Carreira, I.M.; Galhano, E.; Pita, C.; do Céu Almeida, M. Mid-infrared (MIR) metabolic fingerprinting of amniotic fluid: A possible avenue for early diagnosis of prenatal disorders? Anal. Chim. Acta 2013, 764, 24–31. [Google Scholar] [CrossRef] [PubMed]
  13. Graça, G.; Goodfellow, B.J.; Barros, A.S.; Diaz, S.; Duarte, I.F.; Spagou, K.; Veselkov, K.; Want, E.J.; Lindon, J.C.; Carreira, I.M. UPLC-MS metabolic profiling of second trimester amniotic fluid and maternal urine and comparison with NMR spectral profiling for the identification of pregnancy disorder biomarkers. Mol. BioSystems 2012, 8, 1243–1254. [Google Scholar] [CrossRef] [PubMed]
  14. Virgiliou, C.; Gika, H.G.; Witting, M.; Bletsou, A.A.; Athanasiadis, A.; Zafrakas, M.; Thomaidis, N.S.; Raikos, N.; Makrydimas, G.; Theodoridis, G.A. Amniotic fluid and maternal serum metabolic signatures in the second trimester associated with preterm delivery. J. Proteome Res. 2017, 16, 898–910. [Google Scholar] [CrossRef] [PubMed]
  15. Graca, G.; Duarte, I.F.; Barros, A.S.; Goodfellow, B.J.; Diaz, S.O.; Pinto, J.; Carreira, I.M.; Galhano, E.; Pita, C.; Gil, A.M. Impact of prenatal disorders on the metabolic profile of second trimester amniotic fluid: A nuclear magnetic resonance metabonomic study. J. Proteome Res. 2010, 9, 6016–6024. [Google Scholar] [CrossRef]
  16. Hallingström, M.; Barman, M.; Savolainen, O.; Viklund, F.; Kacerovsky, M.; Brunius, C.; Jacobsson, B. Metabolomic profiles of mid-trimester amniotic fluid are not associated with subsequent spontaneous preterm delivery or gestational duration at delivery. J. Matern. -Fetal Neonatal Med. 2022, 35, 2054–2062. [Google Scholar] [CrossRef]
  17. Giovannopoulou, E.; Tsakiridis, I.; Mamopoulos, A.; Kalogiannidis, I.; Papoulidis, I.; Athanasiadis, A.; Dagklis, T. Invasive prenatal diagnostic testing for aneuploidies in singleton pregnancies: A comparative review of major guidelines. Medicina 2022, 58, 1472. [Google Scholar] [CrossRef]
  18. Jacob, D.; Deborde, C.; Lefebvre, M.; Maucourt, M.; Moing, A. NMRProcFlow: A graphical and interactive tool dedicated to 1D spectra processing for NMR-based metabolomics. Metabolomics 2017, 13, 36. [Google Scholar] [CrossRef]
  19. Xia, J.; Bjorndahl, T.C.; Tang, P.; Wishart, D.S. MetaboMiner–semi-automated identification of metabolites from 2D NMR spectra of complex biofluids. BMC Bioinform. 2008, 9, 507. [Google Scholar] [CrossRef]
  20. Murgia, F.; Monni, G.; Corda, V.; Hendren, A.J.; Paci, G.; Piras, A.; Ibba, R.M.; Atzori, L. Metabolomics analysis of amniotic fluid in euploid foetuses with thickened nuchal translucency by gas chromatography-mass spectrometry. Life 2021, 11, 913. [Google Scholar] [CrossRef]
  21. Gauthier, D.W.; Meyer, W.J.; Bieniarz, A. Correlation of amniotic fluid glucose concentration and intraamniotic infection in patients with preterm labor or premature rupture of membranes. Am. J. Obstet. Gynecol. 1991, 165, 1105–1110. [Google Scholar] [CrossRef]
  22. Romero, R.; Quintero, R.; Nores, J.; Avila, C.; Hanaoka, S.; Hagay, Z.; Merchant, L.; Hobbins, J.C. Amniotic fluid white blood cell count: A rapid and simple test to diagnose microbial invasion of the amniotic cavity and predict preterm delivery. Am. J. Obstet. Gynecol. 1991, 165, 821–830. [Google Scholar] [CrossRef] [PubMed]
  23. Kirshon, B.; Rosenfeld, B.; Mari, G.; Beifort, M. Amniotic fluid glucose and intraamniotic infection. Am. J. Obstet. Gynecol. 1991, 164, 818–820. [Google Scholar] [CrossRef] [PubMed]
  24. Parthasarathy, S.; Soundararajan, P.; Sakthivelu, M.; Karuppiah, K.M.; Velusamy, P.; Gopinath, S.C.; Pachaiappan, R. The role of prognostic biomarkers and their implications in early detection of preeclampsia: A systematic review. Process Biochem. 2023, 126, 238–251. [Google Scholar] [CrossRef]
  25. Wang, M.; Wang, H.; Zheng, H.; Uhrin, D.; Dewhurst, R.J.; Roehe, R. Comparison of HPLC and NMR for quantification of the main volatile fatty acids in rumen digesta. Sci. Rep. 2021, 11, 24337. [Google Scholar] [CrossRef]
  26. Shyer, J.A.; Flavell, R.A.; Bailis, W. Metabolic signaling in T cells. Cell Res. 2020, 30, 649–659. [Google Scholar] [CrossRef]
  27. Wan, J.; Jiang, F.; Zhang, J.; Xu, Q.; Chen, D.; Yu, B.; Mao, X.; Yu, J.; Luo, Y.; He, J. Amniotic fluid metabolomics and biochemistry analysis provides novel insights into the diet-regulated foetal growth in a pig model. Sci. Rep. 2017, 7, 1–10. [Google Scholar] [CrossRef]
  28. Romero, R.; Mazaki-Tovi, S.; Vaisbuch, E.; Kusanovic, J.P.; Chaiworapongsa, T.; Gomez, R.; Nien, J.K.; Yoon, B.H.; Mazor, M.; Luo, J. Metabolomics in premature labor: A novel approach to identify patients at risk for preterm delivery. J. Matern. -Fetal Neonatal Med. 2010, 23, 1344–1359. [Google Scholar] [CrossRef]
  29. Harman, C.R. Amniotic fluid abnormalities. Semin. Perinatol. 2008, 32, 288–294. [Google Scholar] [CrossRef]
  30. Dessì, A.; Fanos, V. Myoinositol: A new marker of intrauterine growth restriction? J. Obstet. Gynaecol. 2013, 33, 776–780. [Google Scholar] [CrossRef]
  31. Orczyk-Pawilowicz, M.; Jawien, E.; Deja, S.; Hirnle, L.; Zabek, A.; Mlynarz, P. Metabolomics of human amniotic fluid and maternal plasma during normal pregnancy. PLoS ONE 2016, 11, e0152740. [Google Scholar] [CrossRef]
  32. Menon, R. Oxidative stress damage as a detrimental factor in preterm birth pathology. Front. Immunol. 2014, 5, 567. [Google Scholar] [CrossRef] [PubMed]
Figure 1. 1H NMR spectra of AF sample with annotation of the identified metabolites.
Figure 1. 1H NMR spectra of AF sample with annotation of the identified metabolites.
Metabolites 13 01147 g001
Figure 2. PCA score plot of the studied AF samples; No. of components = 2, n = 60, confidence level = 95%; red dots correspond to full-term samples, while green dots correspond to preterm samples.
Figure 2. PCA score plot of the studied AF samples; No. of components = 2, n = 60, confidence level = 95%; red dots correspond to full-term samples, while green dots correspond to preterm samples.
Metabolites 13 01147 g002
Figure 3. Overview of the buckets’ reduction methodology.
Figure 3. Overview of the buckets’ reduction methodology.
Metabolites 13 01147 g003
Figure 4. PLS-DA analysis for the AF samples of 60 volunteers. (A) Score plot of PLS-DA analysis (R2X(cum) = 0.75, Q2(cum) = 0.61, accuracy = 0.93, comp No. = 3). The green and red dots correspond to preterm and term classification, respectively; (B) VIP plot of the studied metabolites. (C) Validation of the PLS-DA analysis, by permutation test statistics, indicates that the extracted model is significantly different from a model built on random data. The permutation tests were carried out with 1000 random permutations, thus providing significance of the model at the 0.001 level.
Figure 4. PLS-DA analysis for the AF samples of 60 volunteers. (A) Score plot of PLS-DA analysis (R2X(cum) = 0.75, Q2(cum) = 0.61, accuracy = 0.93, comp No. = 3). The green and red dots correspond to preterm and term classification, respectively; (B) VIP plot of the studied metabolites. (C) Validation of the PLS-DA analysis, by permutation test statistics, indicates that the extracted model is significantly different from a model built on random data. The permutation tests were carried out with 1000 random permutations, thus providing significance of the model at the 0.001 level.
Metabolites 13 01147 g004
Figure 5. Box plot of the promoted metabolites from the discriminant analysis.
Figure 5. Box plot of the promoted metabolites from the discriminant analysis.
Metabolites 13 01147 g005
Table 1. Sample demographics and pregnancy outcome data by group (n = 60).
Table 1. Sample demographics and pregnancy outcome data by group (n = 60).
Preterm
(n = 17)
Full Term
(n = 43)
Mann–Whitney
U Test
% (n)% (n)p Value
Sex
Boy58.8% (10)51.2% (22)0.595
Girl41.2% (7)48.8% (21)
Delivery
Cesarean section17.6% (3)7% (3)0.218
Normal delivery82.4% (14)93% (40)
Mean (SD)Mean (SD)
Gestational age (weeks)35.3 (2.8)38.5 (0.9)<0.001
Amniocentesis (week)20.31 (2.46)19.68 (1.82)0.479
Age (years)36.54 (2.70)37.29 (3.63)0.548
Weight (Kg)71.82 (10.54)73.16 (9.11)0.755
Weight gain (Kg)12.9 (6.4)12.9 (5.9)0.786
Neonatal weight (g)2695.29 (553.85)3340.23 (380.77)<0.001
Neonatal weight (z-scores)−0.886 (1.06)0.350 (0.730)<0.001
Neonatal length (cm)51.2 (2.2)51.2 (1.6)0.744
SD = standard deviation.
Table 2. Results of the pathway analysis of the AF samples (in bold, the pathways of importance are depicted).
Table 2. Results of the pathway analysis of the AF samples (in bold, the pathways of importance are depicted).
No.PathwayTotalExpectedHitsRaw plog (p)Holm AdjustFDRImpact
1Aminoacyl-tRNA biosynthesis480.3440.00023.64170.01920.020.00
2Alanine, aspartate, and glutamate
metabolism
280.2030.00083.10150.06570.030.20
3Glyoxylate and dicarboxylate
metabolism
320.2330.00122.92820.09680.030.03
4Butanoate
metabolism
150.1120.00462.33970.37050.100.00
5 Citrate cycle (TCA cycle) 20 0.14 2 0.0081 2.0906 0.6494 0.14 0.12
6Glycolysis/
Gluconeogenesis
260.1820.01361.86761.00000.190.03
7Glycine, serine, and threonine
metabolism
330.2320.02141.66861.00000.260.07
8Phenylalanine,
tyrosine, and
tryptophan
biosynthesis
40.0310.02811.55111.00000.300.50
9Nitrogen
metabolism
60.0410.04191.37781.00000.350.00
10D-Glutamine and
D-glutamate
metabolism
60.0410.04191.37781.00000.350.50
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kolvatzis, C.; Christodoulou, P.; Kalogiannidis, I.; Tsiantas, K.; Tsakiridis, I.; Kyrkou, C.; Cheilari, A.; Thomaidis, N.S.; Zoumpoulakis, P.; Athanasiadis, A.; et al. Metabolomic Profiling of Second-Trimester Amniotic Fluid for Predicting Preterm Delivery: Insights from NMR Analysis. Metabolites 2023, 13, 1147. https://doi.org/10.3390/metabo13111147

AMA Style

Kolvatzis C, Christodoulou P, Kalogiannidis I, Tsiantas K, Tsakiridis I, Kyrkou C, Cheilari A, Thomaidis NS, Zoumpoulakis P, Athanasiadis A, et al. Metabolomic Profiling of Second-Trimester Amniotic Fluid for Predicting Preterm Delivery: Insights from NMR Analysis. Metabolites. 2023; 13(11):1147. https://doi.org/10.3390/metabo13111147

Chicago/Turabian Style

Kolvatzis, Charalampos, Paris Christodoulou, Ioannis Kalogiannidis, Konstantinos Tsiantas, Ioannis Tsakiridis, Charikleia Kyrkou, Antigoni Cheilari, Nikolaos S. Thomaidis, Panagiotis Zoumpoulakis, Apostolos Athanasiadis, and et al. 2023. "Metabolomic Profiling of Second-Trimester Amniotic Fluid for Predicting Preterm Delivery: Insights from NMR Analysis" Metabolites 13, no. 11: 1147. https://doi.org/10.3390/metabo13111147

APA Style

Kolvatzis, C., Christodoulou, P., Kalogiannidis, I., Tsiantas, K., Tsakiridis, I., Kyrkou, C., Cheilari, A., Thomaidis, N. S., Zoumpoulakis, P., Athanasiadis, A., & Michaelidou, A. -M. (2023). Metabolomic Profiling of Second-Trimester Amniotic Fluid for Predicting Preterm Delivery: Insights from NMR Analysis. Metabolites, 13(11), 1147. https://doi.org/10.3390/metabo13111147

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

Article Metrics

Back to TopTop