Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth

Pruski, Pamela; Correia, Gonçalo D. S.; Lewis, Holly V.; Capuccini, Katia; Inglese, Paolo; Chan, Denise; Brown, Richard G.; Kindinger, Lindsay; Lee, Yun S.; Smith, Ann; Marchesi, Julian; McDonald, Julie A. K.; Cameron, Simon; Alexander-Hardiman, Kate; David, Anna L.; Stock, Sarah J.; Norman, Jane E.; Terzidou, Vasso; Teoh, T. G.; Sykes, Lynne; Bennett, Phillip R.; Takats, Zoltan; MacIntyre, David A.

doi:10.1038/s41467-021-26215-w

Download PDF

Article
Open access
Published: 13 October 2021

Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth

Nature Communications volume 12, Article number: 5967 (2021) Cite this article

10k Accesses
25 Citations
151 Altmetric
Metrics details

Subjects

Abstract

The pregnancy vaginal microbiome contributes to risk of preterm birth, the primary cause of death in children under 5 years of age. Here we describe direct on-swab metabolic profiling by Desorption Electrospray Ionization Mass Spectrometry (DESI-MS) for sample preparation-free characterisation of the cervicovaginal metabolome in two independent pregnancy cohorts (VMET, n = 160; 455 swabs; VMET II, n = 205; 573 swabs). By integrating metataxonomics and immune profiling data from matched samples, we show that specific metabolome signatures can be used to robustly predict simultaneously both the composition of the vaginal microbiome and host inflammatory status. In these patients, vaginal microbiota instability and innate immune activation, as predicted using DESI-MS, associated with preterm birth, including in women receiving cervical cerclage for preterm birth prevention. These findings highlight direct on-swab metabolic profiling by DESI-MS as an innovative approach for preterm birth risk stratification through rapid assessment of vaginal microbiota-host dynamics.

HIV-associated vaginal microbiome and inflammation predict spontaneous preterm birth in Zambia

Article Open access 20 May 2022

Joan T. Price, Bellington Vwalika, … Jeffrey S. A. Stringer

Volatile organic compound analysis, a new tool in the quest for preterm birth prediction—an observational cohort study

Article Open access 22 July 2020

Lauren Lacey, Emma Daulton, … Siobhan Quenby

Shotgun sequencing of the vaginal microbiome reveals both a species and functional potential signature of preterm birth

Article Open access 12 November 2020

Conor Feehily, David Crosby, … Paul D. Cotter

Introduction

The vaginal microbiome is a key mediator of reproductive tract pathophysiology. Unlike other mucosal surfaces such as the gut, low diversity in the vaginal microbiome and dominance by Lactobacillus species is considered a hallmark of health, particularly during reproductive years. Lactobacillus species depletion and increased microbial diversity are characteristic of bacterial vaginosis (BV)¹, and associates with both the increased risk of acquisition² and ineffective treatment of sexually transmitted infections, including HIV^3,4. High diversity, BV-associated vaginal microbiota have also been linked to HPV infection and cervical dysplasia^5,6, in vitro fertilization failure⁷, miscarriage⁸ and preterm birth^{9,10,11,12,13}. Current clinical diagnosis of vaginal infection is largely limited to subjective assessment of clinical symptoms in addition to time-consuming microscopic evaluation of vaginal swab samples, pH and, in some cases, culture^14,15,16. These approaches fail to detect many clinically relevant species and lack insight into overall community composition. Such information can be obtained using next-generation sequencing-based approaches (e.g. metagenomics and metataxonomics), but these involve complex and time-consuming sample preparation, data processing and expense that prevents their use for routine bedside testing. Further, molecular-based characterization of microbiota is unable to assess microbiota:host interactions that ultimately determine health and disease phenotypes. Therefore, there exists a need for rapid, point-of-care vaginal diagnostics to facilitate faster clinical decision making, more judicious use of antibiotics and targeted treatment strategies.

A common mechanistic pathway linking sub-optimal vaginal microbiota composition (VMC) and pathophysiology is activation of host-innate immune response and inflammation¹⁷, which can be suppressed by Lactobacillus species such as L. crispatus, through the modification of the metabolic milieu of the cervicovaginal mucosa^16,18. During pregnancy, untimely activation of inflammation in gestational tissues (e.g. cervix, foetal membranes) caused by ascending vaginal infection is thought to cause a significant proportion of infection-associated preterm birth^19,20, which remains the primary cause of death in children under 5 years of age²¹. Consistent with this, we have recently reported that cervical cerclage, a procedure used to reinforce the cervix in women at risk of preterm delivery due to cervical shortening, can induce vaginal bacterial dysbiosis, inflammatory activation and premature cervical ripening associated with increased risk of preterm birth, if performed using a braided instead of monofilament cerclage material²². In contrast, L. crispatus dominance of the vaginal niche during pregnancy has been associated with protection against preterm birth^10,12,22. Despite this, individual and ethnic variation in VMC and host response^11,23,24 as well as cost and time constraints have limited the utility of current analytical methods used for vaginal microbiota characterization to inform clinical decision-making during pregnancy. Current methods for preterm birth prediction have poor positive predictive value and fail to provide insight into underlying aetiology, which may explain low efficacy of interventions designed to prevent preterm birth²⁵.

We have recently described a method that enables rapid, objective assessment of the chemical composition of mucosal surfaces using sample preparation-free, direct on-swab desorption electrospray ionization mass spectrometry (DESI-MS)^26,27. This involves directing a pneumatically assisted electrospray of charged aqueous droplets directly onto a rotating swab, where it forms a liquid film that desorbs and ionizes molecules from the sample, which are transferred to a mass spectrometer via an atmospheric pressure ion-transfer line. With this method, metabolic profiles can be acquired rapidly (in <3 min) and directly from a swab sample without the need for laborious sample preparation or chemical extraction procedures. While lack of chromatographic separation limits assay selectivity, it provides the advantage of vastly simplifying the amount of maintenance and operator intervention required. Additionally, since the solvent stream only ablates sample from a small area of the swab, the method is virtually non-destructible supporting multi-assay use from the sampling device (e.g., cytokine/immune marker profiling or 16S rRNA sequencing)¹⁶. Collectively, these characteristics make direct on-swab DESI-MS particularly suitable for deployment in clinical point-of-care settings.

We hypothesized that vaginal microbiota:host interactions would be reflected in the metabolic milieu of the cervicovaginal mucosa during pregnancy and provide sufficient biochemical information to predict both bacterial composition and host immune response. To test this, direct on-swab DESI-MS profiling was used to characterize the cervicovaginal metabolome of two independent cohorts of women prospectively sampled throughout pregnancy (VMET, n = 160; 455 swabs; VMET II, n = 205; 573 swabs, Fig. 1A). These data were then integrated with matched metataxonomic and immuno-profiling data to identify DESI-MS metabolic signatures predictive of VMC and local inflammatory status. We then examined if these predictive models could be used to monitor changing VMC and host inflammatory responses that associate with preterm birth and clinical interventions (e.g., cervical cerclage) used during pregnancy.

**Fig. 1: Direct on-swab desorption electrospray ionization mass spectrometry (DESI-MS) metabolic profiling of cervicovaginal fluid enables robust prediction of vaginal microbiome compositions.**

Results

Baseline characteristics of the prospective study subjects

A total of 160 women recruited to the VMET study were longitudinally sampled throughout pregnancy (n = 455 swabs) (Fig. 1A). These samples were all subjected to metabolic profiling by DESI-MS and a cross-platform comparison using five complementary liquid chromatography-MS (LC-MS) assays (Fig. 1A). A second independent patient cohort (n = 205, VMET2 study) were longitudinally sampled at comparable time points (n = 573 swabs). Samples from these patients were analysed by DESI-MS and a subset were used for immuno-profiling studies. Key demographics including maternal age, BMI, gravida and ethnicity were similar between the two patient cohorts (Table 1). Preterm births rates in these high-risk populations were 18% (n = 29; VMET) and 21% (n = 44; VMET2).

Table 1 Demographics and clinical characteristics of the study cohorts.

Full size table

Prediction of VMC by direct on-swab DESI-MS

Ward’s linkage hierarchical clustering (HC) analysis with Jensen–Shannon distance metric of bacterial species (operational taxonomic counts, OTUs) was used to group vaginal samples into 11 distinct groups as evaluated with Silhouette scores (Supplementary Fig. 1) that were subsequently reduced into community state types (CSTs) consistent with previous studies (Fig. 1B)^{9,12,13,22,23,28}. Six major CSTs were identified across both cohorts: (1) CST I (L. crispatus); (2) CST II (L. gasseri); (3) CST III (L. iners); (4) CST IV (Gardnerella vaginalis); (5) CST V (L. jensenii); (6) CST VII (other Lactobacillus spp.). CST VI (Bifidobacterium breve) was identified in the VMET2 and not in the VMET cohort, due to differences in the primer sets used for amplicon generation (Supplementary Fig. 1)^29,30. Higher taxonomic classification at genera level was achieved by grouping samples from CST I, II, III, V, and VII into Lactobacillus-dominated (LDOM, VMET: n = 379; VMET2: n = 427), and samples from CST IV and VI into Lactobacillus-depleted (LDEPL, VMET: n = 49; VMET2: n = 112) categories (Supplementary Table 1).

Using linear mixed effect modelling, a total of 113 metabolite features determined in DESI-MS negative and positive ion modes were found to significantly discriminate between LDOM and LDEPL in independent analyses of both patient cohorts (Fig. 1B) (Supplementary Data 1). The vast majority of these metabolites did not change with advancing gestation (Supplementary Fig. 2). A detailed assessment of variance explained by CST, gestational age, ethnicity, BMI and maternal age for each of the metabolite biomarkers highlighted CST and ethnicity as being the major factors of variance (Supplementary Fig. 3). In addition, the proportion of total metabolome variance explained by each factor was estimated with permutational multivariate analysis of variance (PERMANOVA, Supplementary Table 2). While the R² values estimated for each covariate varied between metabolic profiling assays, their relative importance was consistent, with between-individual variability explaining the most variance (35–45.2%), followed by CST (2.2–25.2%), ethnicity (0.5–2.4%) and gestational age (0.4–1.8%). Of those, thiomalic acid and leucyl-serine were consistently higher in LDOM samples, whereas docosanoic acid and lignoceric acid were significantly higher in LDEPL samples (Fig. 1C). Random forest classifier and ROC-curve analysis were then used to assess the performance of DESI-MS to predict VMC (Fig. 2 and Supplementary Fig. 4). Robust prediction of genera-level classification (LDOM vs LDEPL) was observed across patient cohorts particularly using features derived from negative ion polarity mode (VMET/VMET2; AUC 94.1/90.6; sensitivity: 62.0/54.5; specificity: 97.8/96.4). Discrimination between the major vaginal CSTs (CSTI, III and IV) could also be readily achieved by direct on-swab DESI-MS (Fig. 2). The predictive performance of DESI-MS for all models was comparable to averaged (min/max) prediction performance of LC-MS assays (VMET/VMET2; AUC 95.76 (94.3–97.7); sensitivity: 69.74 (57.5–80.2); specificity 98.28 (97.7–98.9)) (Supplementary Table 3 and Supplementary Figs. 4 and 5).

**Fig. 2: Comparison of DESI-MS classification performance between different vaginal microbiome compositions.**

In vitro DESI-MS profiling of vaginal commensal and pathogenic bacterial isolates

We next investigated if any of the discriminatory metabolites identified in the in vivo DESI-MS analyses could also be observed in culture isolates (n = 25) of bacterial species recognized as being predominant members of major vaginal CSTs. In total, 27 of the discriminatory metabolites identified in vivo were detected by DESI-MS in swabs of culture biomasses following correction for background media concentrations (Fig. 3). Approximately half of these metabolites were detected at levels lower than that observed in media background controls, whereas the remainder were found at levels higher than media background. No clear relationship was observed between the DESI-MS detected levels of these metabolites in vitro, CST membership nor levels observed in vivo.

**Fig. 3: Detection of in vivo discriminatory metabolite features in bacterial biomasses by DESI-MS.**

Assessment of host immune response at the mucosal interface using direct on-swab DESI-MS

A panel of 22 soluble immune markers including cytokines, chemokines, immunoglobulins and members of the complement system were measured in a subset of the VMET2 cohort samples (n = 391). Random forest regression analysis was used to predict the log-transformed concentrations of each marker using DESI-MS-derived features. Robust prediction (cross-validated R² > 0.25) was observed particularly for IL-1β (CV R² = 0.51), IL-8 (CV R² = 0.37), C3b/iC3b (CV R² = 0.31), IgG3 (CV R² = 0.31), IgG2 (CV R² = 0.27) and MBL (mannose-binding lectin) (CV R² = 0.26) (Fig. 4A). Immune markers with CV R² > 0.1 (n = 9) and DESI-MS features with an R² > 0.1 for the linear regression against the immune marker (n = 23) revealed a metabolic signature strongly associated with local immune phenotype primarily characterized by altered levels of long-chain fatty acids, glycerophospholipids and ceramides (Fig. 4B and Supplementary Table 4). Both DESI-MS predicted and immunoassay-measured levels of C3b, IL-1β, IgG2 and IgG3 were elevated in Lactobacillus-depleted vaginal microbiomes indicating activation of the local innate and adaptive immune response (Fig. 4C).

**Fig. 4: Assessment of host response at the mucosal interface using direct on-swab DESI-MS profiling.**

Direct-on swab DESI-MS as a potential tool for PTB risk stratification

We next examined the relationship between pregnancy outcome, VMC and inflammatory status as predicted by DESI-MS by combining all available data from both VMET and VMET2 patient cohorts. High vaginal microbiome diversity and instability during pregnancy, as defined by shifts between Lactobacillus-dominated and Lactobacillus-depleted compositions classified with 16S rRNA gene-based metataxonomics, was associated with an increased risk of preterm birth compared to those women maintaining Lactobacillus-dominance throughout pregnancy (odds ratio (OR) 1.97, 95% confidence interval (CI) 1.03–1.66, p = 0.04) (Fig. 5A). Similarly, vaginal microbiota instability determined solely by DESI-MS was also associated with a higher risk of preterm birth although this did not reach statistical significance (OR 1.47, 95% CI 0.75–2.78, p = 0.25).

**Fig. 5: Vaginal microbiome instability and immune activation associates with preterm birth risk and poor outcomes following cervical cerclage.**

Increased levels of IL-1β, measured directly by immunoassay or predicted by DESI-MS, were associated with high-diversity VMC in both those women experiencing term birth or preterm birth. However, levels were significantly higher in preterm birth cases (preterm 59.63 pg/ml vs term 1.94 pg/ml, p = 0.006, Welch two sample t-test, Fig. 5B). Increased levels of MBL, again either measured directly by immunoassay or predicted by DESI-MS, were also observed in high-diversity VMC of women subsequently delivering preterm (term 0.042 ng/ml vs preterm 0.692 ng/ml, p = 0.007, Welch two-sample t-test, Fig. 5C).

We next tested if DESI-MS metabolic, metataxonomics and inflammatory marker profiles obtained at three different stages of gestation, could predict subsequent preterm birth. Overall, the predictive capacity of these models was poor (Supplementary Table 5). We therefore focused our subsequent analyses on preterm birth phenotypes more likely to associate with dysregulated vaginal microbiota–host interactions, as previously described²². In women at high risk of preterm birth due to cervical shortening, vaginal levels of MBL were most frequently increased in those treated with a cervical cerclage using braided suture material compared to those where cerclage had been performed with monofilament suture (10/11, 91% vs 9/21, 43%, p = 0.011, Fisher’s exact test). This was similarly detected when VMC and MBL levels were estimated using only DESI-MS analysis of vaginal swabs (7/11, 64% vs 8/21, 38%, p = 0.316; Fig. 5D). Immunoassay also indicated increased IL-1β after braided cerclage insertion, but this was less consistently detected by DESI-MS prediction (Fig. 5E). However, following treatment with braided cerclage material, DESI-MS did accurately predict low levels of IL-1β in women delivering at term compared to those who subsequently delivered preterm (Fig. 5F).

Discussion

Despite recent developments in metataxonomics and metagenomics, VMC characterization in clinical settings remains largely limited to culture and microscopy, which like molecular-based approaches, fails to capture information regarding host response. Our method, which is easily amenable to bedside point-of-care testing^31,32, addresses this limitation by leveraging information contained within the cervicovaginal mucosa to provide robust detection of VMCs and simultaneous estimation of host immune and inflammatory status.

Many of the discriminatory metabolites identified in our study have been previously reported using gas chromatography-MS or LC-MS based assays in women suffering from BV^33,34,35,36 or in HPV infection³⁷ and have biologically plausible roles in mediating vaginal health and disease. For example, increased vaginal levels of short- and long-chain fatty acids and biogenic amines are associated with activation of pro-inflammatory pathways^{36,38,39,40,41}, which contribute to reduced barrier integrity of the epithelia and consequently, increased risk of infection^{34,42,43,44,45}. Furthermore, cervicovaginal lipid species have recently been associated with tumour progression and genital inflammation defined using an aggregated score, in a study of 78 non-pregnant, HPV-negative/positive women with cervical dysplasia³⁷. Lipid changes in vaginal discharge have also been associated with vulvovaginal candidiasis⁴⁶. In our study, we profiled a larger panel of soluble immune mediators, including humoral response mediators, and identified associations with the metabolome via regression models. Despite different experimental approaches, our results provide further evidence that lipid species are a core component of the metabolic signature of inflammation and immune response in the vaginal niche. Some of the discriminatory markers used in our in vivo DESI-MS prediction models were also detected by DESI-MS in vitro swab analyses of bacterial isolate biomasses. It should be noted that the bacterial isolates analysed here were not representative of the genetic diversity of vaginal commensal and pathogenic strains; however, our analyses indicate that the metabolite signatures predictive of microbiota composition and inflammation in vivo are likely derived from both bacterial and host sources.

Some studies have estimated that around 23–30% of total cervicovaginal metabolic variation is associated with bacterial composition^36,47,48. Similar estimates obtained in our data emphasize the impact of the microbial composition on the metabolome, which is greater than the effect of gestational age, ethnicity, maternal age, or BMI. However, using the repeated measures per individual, we found that the top proportion of variance explained could be attributed to between-individual variability. The magnitude of the R² values for microbial composition was higher in the LC-MS assays compared to DESI-MS. This may represent greater coverage of relevant constituents of the metabolome that are not detected by DESI-MS or may reflect the technical variability of DESI-MS which, like other ambient ionization methods, is recognized to be larger than that of LC-MS where implementation of experimental quality control (QC), instrumental drift correction and data filtering procedures are more established⁴⁹. Despite this, our findings highlight the capacity of direct on-swab DESI-MS to rapidly capture this information from the metabolome without the need for laborious sample preparation and comparatively high per/assay time and economic costs associated with coupled chromatography-MS assays for similar applications. This lends itself to customization and automation for ease-of-use, which are important characteristics for point-of-care testing³². Our approach also offers the advantage of providing objective simultaneous assessment of microbiota composition and inflammatory state. In comparison, current ‘gold standard’ diagnostic methods for vaginal infection (e.g. Amsel criteria) are limited to a combined subjective assessment of clinical symptoms and microscopy grading of bacterial morphotypes as well as selective culture.

There is now substantial data supporting a role for the vaginal microbiome and host immune responses in shaping preterm birth in a proportion of women^10,11,12,22. However, prediction of preterm birth using DESI-MS metabolic profiles, metataxonomics or inflammatory marker data in our patient cohort was poor. This is not surprising given the fact that preterm birth is a multi-aetiological disease state that can be caused by many different factors, including non-microbial and non-immune related causes^19,20. Because of this, we focused subsequent analyses on women who receive cervical cerclage with braided suture, who we have previously shown are at increased risk of preterm birth that involves a phenotype characterized by vaginal dysbiosis and local immune activation²². Here we confirmed these findings in an independent patient population and highlight the utility of direct swab analysis by DESI-MS to rapidly detect both microbiota and inflammatory changes caused by the intervention that associate with subsequent preterm birth risk. The ability to provide such information at point-of-care would be transformative for directing clinical decision making and ultimately improving outcomes for these women and their babies.

In this study, DESI-MS prediction of VMCs was capable of monitoring vaginal microbiome diversity and instability longitudinally throughout pregnancy, which in the study participants was associated with increased risk of preterm birth. These results are consistent with a recent meta-analysis of metataxonomics-based studies of the vaginal microbiome in pregnancy, which reported higher variance of VMC across trimesters in women subsequently experiencing preterm birth⁵⁰. An additional strength of our approach is the ability to simultaneously capture relationships between VMC and local immune response. Innate immune activation in the vagina often accompanies disease states associated with suboptimal VMC, including BV, preterm birth, and sexually transmitted infections^51,52,53. For example, L. crispatus dominance of the vaginal niche is associated with suppressed levels of IL-1β compared to high-diversity compositions^{17,54,55,56,57,58}. Further, MBL is a recognition molecule for G. vaginalis⁵⁹ and several single-nucleotide polymorphisms in the MBL2 gene have been reported to increase the risk of BV^60,61. Our findings show that local vaginal immune responses are reflected in the cervicovaginal metabolic phenotype and can be readily detected by DESI-MS. This has clear translational benefit. We have recently reported that cervical cerclage, a procedure used to reinforce the cervix in women at risk of preterm delivery due to cervical shortening, can induce vaginal bacterial dysbiosis, inflammatory activation and premature cervical ripening associated with increased risk of preterm birth, if performed using a braided instead of monofilament cerclage material²². Consistent with these findings, DESI-MS showed capacity to detect clinically relevant inflammatory responses to cerclage insertion that associated with patient outcome. Coupled with VMC prediction, this highlights direct swab profiling by DESI-MS as an innovative platform for host–microbiota interactions during pregnancy that could potentially facilitate PTB risk stratification, optimization of preventative interventions²⁹ (e.g. targeted antibiotic treatment following preterm premature rupture of the foetal membranes and response to treatment interventions during pregnancy⁹). While our study has focused upon application of DESI-MS in pregnancy cohorts, the rapid detection of VMC and host response offered by this methodology could also inform treatment strategies in other clinical scenarios. For example, efficacy of the topical pre-exposure prophylactic, Tenofovir, is superior in preventing HIV acquisition in women with Lactobacillus-dominated VMCs compared to those dominated by G. vaginalis and BV-associated bacteria³. It is therefore conceivable that efficacy and effectiveness could improve through stratification of treatment using point-of-care metabolic profiling of vaginal swabs. By extension, our approach may also be useful for bedside monitoring of response to treatments designed to optimize VMC, such as live biotherapeutics⁶² and vaginal microbiome transplantation⁶³.

A limitation of our approach was the relatively low sensitivity of prediction. Further inspection of the 16S rRNA gene amplicon data indicated that this was largely associated with misclassification of samples harbouring mixed compositional structures, often containing G. vaginalis (Supplementary Fig. 6). This suggests that hard-clustering techniques often used for determining CSTs may under-estimate the impact of low abundance taxa on the host mucosal metabolome. Comparison of predictive performances from the LC-MS assays indicated that metabolic coverage may also impact misclassification rate with small polar or non-polar molecule-based assays marginally outperforming lipid-based assays (Supplementary Figs. 5 and 7). It is also worth noting that the analysis of VMC in this study was limited to the assessment of relative abundances of the bacterial component of microbiome and a comparatively small number of patients carrying less prevalent community compositions (e.g., L. gasseri, L. jensenii and B. breve-dominated). The detection of B. breve-dominated women in our study was limited to the VMET2 patient cohort where a mixed formulation of the 27F forward primer set was used. This primer formulation has been shown to maintain the rRNA gene ratio of key vaginal species including Lactobacillus spp. to Gardnerella spp. as well as improve detection of Bifidobacterial species, which are otherwise not detected by the ‘universal’ 27f primer often used in metataxonomics studies^30,64. Consistent with these results, a number of recent studies have also identified B. breve-dominated vaginal microbiota using metataxonomics approaches^65,66. Additional datasets that provide improved representation of more diverse compositions and rarer taxa may also eventually facilitate prediction of VMC by DESI-MS beyond major CSTs as described here, and toward species level classification. It is also important to note that recent metataxonomics studies have indicated higher false discovery rates when using relative abundance from NGS data compared to quantitative PCR⁶⁷. Future studies incorporating microbial load within a sample may further improve the sensitivity of DESI-MS predictive models, particularly those with intermediate compositions.

In conclusion, we show that direct on-swab DESI-MS permits rapid and robust assessment of vaginal microbiota:host immune interactions reflected within the cervicovaginal mucosal metabolome. While we show that this may offer an approach for preterm birth risk stratification and selective targeting of preventative treatments, we expect it to have wider application to the assessment of vaginal microbiota–host interactions in other clinical scenarios, including both pregnant and non-pregnant women.

Methods

Study subjects and sample collection

The study was conducted with approval of the NHS National Research Ethics Service (NRES) Committees London-City and East (REC 12/LO/2003) and London–Stanmore (REC 14/LO/0328), and by the North of Scotland Research Ethics Service (REC 14/NS/1078). All participants provided written informed consent prior to sampling and experiments were performed in accordance with the approved institutional guidelines. Recruitment and sampling were performed at Imperial College Healthcare NHS Trust Hospitals (Queen Charlotte’s and Chelsea and St Mary’s Hospitals), London, UK, at Chelsea & Westminster Hospital (NHS Trust, London, UK), University College London Hospital (NHS Foundation Trust, London, UK) and the Royal Infirmary of Edinburgh, Scotland, UK. Eligibility criteria were pregnant women with a singleton pregnancy, with and without risk factors for preterm birth. Exclusion criteria included women under 18 years of age, sexual intercourse within 72 h of sampling, vaginal bleeding in the preceding week, antibiotic use in the preceding 2 weeks, multiple pregnancies, HIV or hepatitis C-positive status. Detailed maternal clinical metadata were collected for all participants and birth outcome recorded following delivery. Preterm birth was defined as labour prior to 37 weeks of gestation. Cervicovaginal fluid swab samples were collected at up to five timepoints throughout pregnancy (1. 8–16, 2. 16–20, 3. 20–26, 4. 26–30, 5. 30–37 weeks) from the posterior fornix using BBL CultureSwab MaxV Liquid Amies swabs (Becton, Dickinson and Company, Oxford, UK), placed in either Amies transport media or a sterile microcentrifuge tube on ice, before long-term storage at −80 °C.

Metabolic profiling of cervicovaginal swabs using direct swab analysis by DESI-MS

All chemicals used were analytical reagent grade. HPLC grade methanol and water for DESI-MS analysis and swab sample extraction were purchased from Sigma-Aldrich (St Louis, MO). The method used for the direct analysis of vaginal swab samples by DESI MS is detailed elsewhere²⁶. Briefly, we used an LTQ-Orbitrap Discovery mass spectrometer (Thermo Scientific, Bremen, Germany) coupled with a DESI-MS source designed for direct swab analysis. Swabs were placed into a rotating holder positioned orthogonally in front of the MS inlet capillary with a swab–capillary distance of approximately 2 mm. The DESI sprayer tip was pointed to the swab centre with a tip–sample distance of 1.5–2 mm and a distance between the tip and the inlet capillary of 2 mm. The entire surface of the medical swabs was analysed by DESI-MS through clockwise rotation of the swab toward the MS capillary. The cervicovaginal mucosa was absorbed from the swab tip by gently desorbing the analytes with charged droplets of methanol/water (95:5, v/v) mixture and directed to the mass spectrometer. For each sample, 30 scan mass spectra (m/z 50–1000, R = 30,000 (FWHM)) were recorded in the negative and positive ion mode.

Metabolic profiling of cervicovaginal swabs using direct swab analysis by LC-MS analysis

Liquid extraction was performed on each swab by adding a MeOH:H₂O (1:1, v-v) solution as eluent to a final concentration of 50 mg vaginal fluid/ml. Each blank swab was extracted with 1 ml solution using a repeated sonication and vortexing cycle for 30 s each. Recovery of soluble material was achieved through centrifugation of swabs (2000 × g for 2 min) seated in a 200 µl loading tip positioned in a sterile microcentrifuge tube. Associated supernatants were pooled and centrifuged at 16,000 × g for 10 min to remove insoluble material before the resulting supernatant was divided into three aliquots. An additional extraction of lipids from the swab was performed by repeating the procedure with an isopropanol:water (4:1, v-v) solution to a final extraction of 25 mg vaginal fluid/ml. The resulting extract was pooled with one aliquot of the methanol:water extract. All extracts were evaporated using a SpeedVac for further reconstitution before analysis by LC-MS^68,69.

Reversed phase LC-MS analysis of small metabolites

Sample reconstitution was performed in 300 μl of water. In all, 250 μl aliquot of reconstituted sample material were used for individual sample preparation of 96-well plates including the addition of 25 μl of full RP 2× labelled standard mix (l-glutamine-¹³C₅; l-glutamic acid ¹³C₅; creatinine-methyl-D₃; cytidine-5,6-D₂; citric acid ¹³C₆; l-isoleucine-¹³C₆¹⁵N; l-leucine-13C₆; l-phenylalanine-¹³C₉¹⁵N; hippuric acid-D₅, benzoic acid-¹³C₆, octanoic acid-¹³C₈, l-tryptophane-¹³C₁₁¹⁵N₂). In addition, 50 μl aliquots of each sample were used for pooling and generation of QC sample. For chromatographic separation a 2 μl aliquot of extracted metabolites from each sample was injected onto a reverse-phase 150 × 2.1 mm ACQUITY 1.8-μm High Strength Silica (HSS) column (Waters Corp.) kept at 45 °C using an ACQUITY UPLC system (Waters Corp.). The mobile phase consisting of 0.1% v/v formic acid (Fisher Scientific) in water (A) and acetonitrile containing 0.1% formic acid (B, Sigma-Aldrich). Each sample was resolved for 12.65 min at a flow rate of 0.5 ml/min. The gradient consisted of 99% A and 1% B for 0.1 min, a ramp of curve 6–100% B from 0.1 to 10.70 min.

Reversed phase LC-MS analysis for lipids

Sample reconstitution was performed in 300 μl of water/isopropanol (4:1). In all, 250 μl aliquot of reconstituted sample material were used for individual sample preparation of 96-well plates including the addition of full RP-labelled standard mix (C17:0; LPC(6:0/0:0); LPC(9:0/0:0); LPC(15:0/0:0); PC(11:0/11:0); PC(15:0/15:0); PE(15:0/15:0); PA(17:0/17:0); PG(15:0/15:0); PS(17:0/17:0); SM(d18:1/17:0); Cer(d18:1/17:0); DG(19:0/0:0/19:0); PC(23:0/23:0); TG(8:0/8:0/8:0); TG(10:0/10:0/10:0); TG(12:0/12:0/12:0); TG(14:0/14:0/14:0); TG(15:0/15:0/15:0); TG(16:0/16:0/16:0); TG(17:0/17:0/17:0); DG(18:0/20:4/0:0); DG(18:0/18:0)). In addition, 50 μl aliquots of each sample were used for pooling and generation of QC sample. For chromatographic separation a 2 μl aliquot of extracted metabolites from each sample was injected onto a Waters Acquity UPLC BEH C8, 1.7 µm, 2.1 × 100 mm column (Waters Corp.) kept at 55 °C using an ACQUITY UPLC system (Waters Corp.). The mobile phase consisting of water:isopropanol:acetonitrile (50:25:25, all high-grade LC-MS solvents from Fisher Scientific or Sigma-Aldrich) + 5 mM ammonium acetate + 0.05% acetic acid + 20 µM phosphoric acid (A) and isopropanol:acetonitrile 50:50 (Sigma-Aldrich) + 5 mM ammonium acetate + 0.05% acetic acid (B). Each sample was resolved for 13.15 min at a flow rate of 0.5 ml/min. Starting conditions were 99% A and 1% B and the gradient changed with a ramp of curve 6 as follows: decrease to 70% A and 30% B over the first 2 min, decrease to 10% A with 90% B from 2 to 11.50 min, decrease to 0.1% A with 99.9% B from 11.50 to 12.50 min, after which the solvent composition returned to starting conditions over 0.1 min until 13.15 min.

HILIC LC-MS analysis

Sample reconstitution was performed in 200 μl of water:acetonitrile mixture (1:3.6, v-v). 115 μl aliquot of reconstituted sample material were used for individual sample preparation of 96-well plates including the addition of 2.5 μl of 48× labelled IS standard mix ((uracil-2-¹³C¹⁵N₂; N-benzoyl-d₅-glycine, adenosine-2-D₁, adenine-2-D₁, taurine-¹⁵N, l-tryptophan-D₅, l-phenylalanine-¹³C₉¹⁵N, creatine-(methyl-D₃)-monohydrate, l-arginine-¹³C₆-hydrochloride). In addition, 20 μl aliquots of each sample were used for pooling and generation of QC sample. HILIC chromatography analysis was performed using a Waters Acquity UPLC BEH HILIC (1.7 μm, 2.1 × 150 mm) column (Waters Corporation, Milford, MA, USA) kept at 40 °C. The mobile phases consisted of 0.1% formic acid and 20 mM ammonium formate in water (A, Fisher Scientific), and 0.1% formic acid in acetonitrile (B, Sigma-Aldrich) at a flow rate of 0.5 ml/min. Starting conditions were 5% A and 95% B and the gradient changed with a ramp of curve 6 as follows: increase to 20% A and 80% B over the first 5.5 min; increase to 50% A with 50% B from 5.5 to 7.1 min, after which the solvent composition returned to starting conditions over 0.1 min until 12.65 min. The injection volume was 2 μl of extracted metabolites for negative ion mode analysis.

LC-MS instrumental operation and analysis

The column eluent was introduced directly into the mass spectrometer by electrospray. MS was performed on a Waters Xevo G2-S QTOF mass spectrometer (Waters Ltd, Elstree, UK) operating in either negative-ion (ESI−) or positive-ion (ESI+) electrospray ionization mode with a capillary voltage of 1 kV for RP and 1.5 kV for HILIC, and a sampling cone voltage of 20 V. The desolvation gas flow was set to 1000 l/h and the desolvation temperature was set to 600 °C. The cone gas flow was 150 l/h, and the source temperature was 120 °C. Accurate mass was maintained by introduction of LockSpray interface of Leucine Enkephalin (m/z 556.2771 in ESI+, m/z 554.2615 in ESI−) at a concentration of 200 pg/μl in 50% aqueous acetonitrile with a scan time of 0.07 s over 4 scans, after each interval of 60 s. Data were collected in centroid mode from 50 to 1200m/z in MS scanning mode. To ensure system suitability and stability, a study reference (SR) QC sample was prepared by combining equal aliquots of all the samples and injected at regular intervals throughout the analytical run. This SR sample was also used to condition the column (30 injections) prior to the analysis of both the ESI+ and ESI mode batches. Blank samples (i.e. injection of the reconstitution solvent) were also run to check the presence of artefact or contaminant peaks. A dilution series of the SR sample was also acquired at the beginning and end of injection sequence. Data-dependent acquisition (DDA) and MS^E analysis of the SR sample were performed for structural elucidation.

Mass spectral data processing

DESI-MS raw mass spectral data were converted from.raw files into.mzXML format via the ProteoWizard msConverterGUI⁷⁰. The first 30 spectra per sample were averaged into a single spectrum and saved as.csv file using the MALDIquant R package (v1.19.3). The detailed parameter settings and algorithm used in the R package MALDIquant⁷¹ are as follows: Peak detection was performed using ‘detectPeaks’, with SNR set to 3 and half window size to 10, and using the median absolute deviation (MAD) method. Peaks were aligned with a half window size set to 20, tolerance set to 0.002, signal-to-noise-ratio (SNR) set to 2 using the warping method LOWESS. The reference peaks were created by calling the function ‘referencePeaks’ from MALDIquant with method ‘strict’, ‘minFrequency’ set to 0.01 and peak binning tolerance 0.002.

Raw LC-MS data files were first converted to the.mzML open format and the Proteowizard msconvert..mzML files were processed in R (v4.0.3) using the XCMS (v3.10.2) package⁷². Peaks were detected with the centWave algorithm (ppm = 25, mzdiff = 0.001, prefilter = c(4, 1000), noise = 100, snthresh = 5), and grouped with the ‘density’ method (bw = 3, mzwid = 0.0007, minFrac = 0.4). The centWave peakwidth parameter was set depending on the chromatographic method (c(1.5, 14) for HILIC-LC-MS (−) c(3, 12) in RP-LC-MS Lipid (+/−) and c(1.5, 5) in the RP-LC-MS Metabolite (+/−). Non-detected peaks were filled with fillPeaks, using the default arguments, and no retention time correction was done. The XCMS datasets were further filtered and the feature intensities drift corrected in Python (v3.8.5) using the nPYc-Toolbox (v1.2.4)⁷³. Drift correction was performed using LOWESS trend-line model fitted on the SR samples. Features where the coefficient of variation measured on repeated injections of the SR sample were larger than 30% or whose Pearson correlation with dilution was less than 0.7 (measured using the SR dilution series) were excluded from the data matrix.

Metabolite identification

Target m/z features were putatively annotated using online databases including HMDB, Metlin, MMCD, KEGG and Lipidmaps with a 5 ppm tolerance for each compound. To account for potential imprecision in m/z measurement due to peak binning and data processing artefacts, raw m/z value were confirmed by inspection of the.raw data. Metabolite annotation was performed by searching the measured m/z ratios against METLIN (http://metlin.scripps.edu), Lipidmaps (http://www.lipidmaps.org) and the HMDB (http://www.hmdb.ca)^74,75,76. Further structural elucidation was performed using MS/MS experiments via collision-induced dissociation on the LTQ-Discovery MS instrument (Thermo Scientific), and with DDA of the precursor ion by the quadrupole of the Xevo G2 XS Q-TOF mass spectrometer (Waters corporation). For the annotation of metabolites, the MS/MS spectra were matched against spectral libraries from HMDB, NIST and METLIN that were compiled with either authentic standards or theoretical assignment. Identification of metabolites where MS/MS reference spectra were not available were annotated using chemical fragmentation rules. For the structural assignment of glycerophospholipid, fragments of the polar head group or the fatty acyl chains were investigated to confirm the annotation proposed by the databases and discriminate isomers. MSI levels of compound identification were further reported as suggested by the Metabolomics Standards Initiative (MSI)⁷⁷.

DNA extraction and sequencing of 16S rRNA amplicons

Extraction of bacterial DNA was performed as previously described⁷⁸. For the VMET cohort, V1–V3 hypervariable regions of bacterial 16S rRNA genes were amplified using a forward and reverse fusion primer. The forward primer was constructed with the Illumina i5 adapter (5′-AATGATACGGCGACCACCGAGATCTACAC-3′), an 8-base pair (bp) bar code, a primer pad (forward, 5′-TATGGTAATT-3′), and the 28F primer (5′-GAGTTTGATCNTGGCTCAG-3′) (64). The reverse fusion primer consisted of the Illumina i7 adapter (5′-CAAGCAGAAGACGGCATACGAGAT-3′), an 8-bp bar code, a primer pad (reverse, 5′-AGTCAGTCAG-3′), and the 519R primer (5′-GTNTTACNGCGGCKGCTG-3′). For the VMET2 cohort, the V1–V2 hyper variable regions were amplified with the forward primer set (28f-YM) consisting of a mixture of the following primers mixed at a 4:1:1:1 ratio: 28F-Borrellia GAGTTTGATCCTGGCTTAG; 28F-Chlorflex GAATTTGATCTTGGTTCAG; 28F-Bifido GGGTTCGATTCTGGCTCAG; 28F GAGTTTGATCNTGGCTCAG. The reverse primer consisted of 388R TGCTGCCTCCCGTAGGAGT³⁰. Sequencing was performed at RTL Genomics (Lubbock, TX, USA) using an Illumina MiSeq platform (Illumina Inc.). Resulting sequence data were analysed using the MiSeq SOP Pipeline of the Mothur package⁷⁹. The Silva bacterial database (www.arb-silva.de/) was used for sequence alignment and classification was performed using the RDP (Ribosomal Database Project) database reference sequence files⁸⁰. Determination of operational taxonomic unit taxonomies (phylum to genus) and species-level taxonomies was performed with USEARCH with 16S rRNA gene sequences from the cultured representatives from the RDP database⁸¹. Species-level taxonomies were complemented using information from the STIRRUPS database⁸².

Counts from all OTUs assigned to the same species were summed to generate matrices of total counts per species. Species with less than 50 total counts were excluded. Community state types were then assigned to each sample with HC, using Ward-linkage and the Jensen–Shannon distance. The number of clusters which maximized the mean silhouette score were selected. The relative abundance of each species in a cluster was inspected and the clusters matched to CST reported in previous publications, by comparing the relative abundances of the top five OTUs obtained for each cluster. Samples where the most abundant species was a Lactobacillus spp. other than L. crispatus, L. iners, L. gasseri or L. jensenii, where manually assigned to a separate cluster (CST VII). HC analyses and heatmaps were performed using Python (v3.8.5), ‘scipy’⁸³ (v1.5.2), the ‘matptlotlib’⁸⁴ (v3.3.2) and ‘seaborn’⁸⁵ (v0.11.0) libraries.

Immune/inflammatory profiling of cervicovaginal samples

Swabs were thawed on ice and re-suspended in 350 ml of phosphate-buffered saline solution containing protease inhibitor (5 µL/ml; Sigma-Aldrich) before being centrifuged at 7000 × g for 10 min to remove cellular debris. Customized multiplex assays (R&D Systems) were used together with a Bio-Plex 200 system (Bio-Rad Laboratories Ltd) to quantify levels of interleukin-1β (IL-1β), IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IL-18, interferon-γ, GM-CSF and tumour necrosis factor-α, with a 10-fold dilution of 1L-8 performed using Calibrator Diluent RD6-52 prior to assaying. The Human Complement Magnetic Bead panel 1 (Milliplex®) was used for measurement of C5, C5a and MBL, while panel 2 was used for measurement of C3b. Levels of immunglobulins IgA, IgG1, IgG2, IgG3, IgG4 and IgM were assessed using the ProcartaPlex Human Antibody Isotyping Panel 7plex assay (Thermo Fisher Scientific).

DESI-MS profiling of bacterial isolates

Bacterial isolates analysed in this study were derived from the DSMZ-German Collection of Microorganism and Cell Culture GmbH (DSM number, Supplementary Table 6) or clinical samples received by the Imperial College NHS Healthcare Trust Diagnostic Microbiology laboratory (NHS number, Supplementary Table 6) at Charing Cross Hospital, London, after the completion of standard identification workflows. Identification of isolates from clinical samples was performed using a MALDI Biotyper instrument (Bruker, UK)⁸⁶. Isolates were stored on beads and in glycerol broth at −80 °C. The bacterial culture library consisted of Atopobium vaginae, Bifidobacterium breve, Gardnerella vaginalis, Lactobacillus crispatus (5× isolates), Lactobacillus gasseri (5× isolates), Lactobacillus iners, Lactobacillus gasseri (5× isolates), Lactobacillus jensenii (5× isolates), Prevotella amni, Prevotella disiens, Prevotella timonensis, Staphylococcus aureus, Streptococcus agalactiae and Streptococcus anginosus. Microorganisms were grown from either fresh cultures or beads with each isolate cultured 10 times using optimal culture conditions (Supplementary Table 5).

Bacterial isolates were sampled directly from solid agar plates using a medical rayon swab before being transferred into a sterile microcentrifuge tube and stored at −80 °C. Three ‘blank’ swab samples of each solid agar media were also collected. Direct swab analysis using DESI-MS was performed in both positive and negative ion mode. Spectra were processed and assembled into a single data matrix. For each isolate, the intensities per m/z feature of biological replicates were compared against their corresponding culture media background samples using the Wilcoxon–Mann–Whitney test. The list of significant features was matched to the set of identified features from the analysis of human cervicovaginal swab samples. The mean fold change for matched peaks (tolerance < 5 ppm) was estimated as the ratio of the mean intensity within bacterial biomass samples to the mean intensity in the background culture media samples. The heatmap in Supplementary Fig 6 was generated in Python (v3.8.5), with the ‘matptlotlib’ (v3.3.2) and ‘seaborn’ (v0.11.0) libraries.

Statistical analysis of metabolomic profiling data

The linear mixed-effects modelling analyses were performed in R using the ‘lme4’⁸⁷ package (v1.1.25). For all metabolic profile variables in each data matrix, a linear mixed-effects model with the following ‘lme4’ formula was fitted: Metabolite ~ GestationalAge + CST + MaternalAge + BMI + Ethnicity + (GestationalAge||SubjectID). This model structure accounts for the repeated measures by using a random intercept per pregnancy and a random slope per pregnancy. No random effect term was used to model correlation between the random slopes and random intercepts. Gestational age, CST, age, BMI and ethnicity were modelled as fixed effects. All models were fitted using restricted maximum likelihood (‘REML = TRUE’ in lmer function call). The ‘emmeans’⁸⁸ package (v1.5.2.1) was used to perform contrast coding and obtain effect size estimates and P values for contrasts and trends of interest. For the comparisons between different CST’s and LDOM vs LDEPL, the mean levels and contrasts obtained from the model were estimated assuming maternal age = 30 years, BMI = 23 and gestational age = 20 weeks, and averaged over all levels of ethnicity. For the gestational age trend estimates were further averaged across all CSTs. Detailed tables with effect size estimates and P values for these analyses and contrasts are available via a GitHub repository (see ‘Code availability’). The comparisons between Lactobacillus dominant vs depleted was encoded as a comparison between the grand mean of CST I, II III, V and VII levels vs the mean of CST IV. The P values for each contrast tested were calculated using the Kenward–Roger approximation, implemented in the ‘pbkrtest’ package (v0.4.8.6)⁸⁹. The P values for all contrasts within a single metabolic signal were not corrected for multiple testing. Instead, for each of the main contrasts of interest (early gestation vs late gestation, CST I vs III, CST I vs IV, CST III vs IV and LDOM vs LDEPL), all P values obtained across all metabolic variables in an assay were pooled and FDR corrected together as a signature, using Benjamini–Hochberg false discovery rate correction and selecting a 5% FDR cut-off. The heatmaps in Fig. 1B were created in Python (v3.8.5), ‘matptlotlib’ (v3.3.2) and ‘seaborn’ (v0.11.0). Only features which were statistically significant after false discovery rate correction and replicated in both VMET and VMET2 datasets were plotted. To account for differences in m/z calibration between datasets, a feature was only considered to be replicated between datasets if it was possible to find a marker with an m/z error of less than 5 part-per-million in the final FDR corrected signature from the other study. The linear mixed effect model, semi-partial r² and conditional R² measures were calculated using the R packages ‘r2glmm’ (v0.1.2) and ‘MuMIn’ (v1.43.17), respectively. PERMANOVA analyses⁹⁰ were performed with the ‘vegan’ package (v2.5.7), using the adonis2 function, with method = ‘euclidean’, permutations = 999, and by adding terms sequentially following their order in the model formula log(MetaboliteMatrix + 1) ~ CST + GestationalAge + Ethnicity + BMI + Age + SubjectID.

For the immune-metabolite association analysis, a series of linear models with the formula lm(log(ImmuneMarker + 1) ~ Metabolite) were fitted in R (v4.0.3) for all immune-marker cytokine pair combinations, and the corresponding F-test P values calculated. Benjamini–Hochberg false discovery rate correction was used independently for each immune marker analysis, with a 5% FDR cut-off. The heatmaps in Fig. 2A, B were generated using Python (v3.8.5), ‘matptlotlib’ (v3.3.2) and ‘seaborn’ (v0.11.0), using the results of the linear model analysis and the random forest regression models (detailed in the next section).

Prediction of CST, preterm birth and immune markers from the metabolomic data

The CST type was predicted from the metabolomic profiles using random forest classifiers. For the Lactobacillus dominant vs depleted comparison, samples which were assigned to CST I, II, III, V or VII were assigned to the Lactobacillus dominant class, and samples with CST IV or VI label were assigned to Lactobacillus depleted. Random forest classifiers were trained to predict the CST or Lactobacillus depletion status using the metabolic profile variables, without using any other clinical or demographical variables, including gestational age. An assessment of the impact on model performance of using gestational age as covariate was also performed (see the analyses under ‘CST Typing by DESI-MS’ in the GitHub code repository).

Prediction of preterm birth was also performed using random forest classifiers. The centred-log-ratio transformation was applied to the 16S data matrices with the R package ‘propr’ (v4.2.6)⁹¹. The gestational period was divided in three gestational windows (1st timepoint: 0–14 weeks; 2nd timepoint: 14–24 weeks; 3rd timepoint 24–40 weeks) and a random forest classifier models fitted at each window using only the 16S data, immune markers, or DESI-MS profiles as predictors. Only one sample per patient (selected as the closest to the gestational age window mid-point) was used in each model. For the immune marker prediction, immune marker concentrations were initially log-transformed after addition of a constant offset log(x + 1). A random forest regressor was trained to predict a single log-transformed immune marker level. All random forest models were trained in R (v4.0.3), using the ‘randomForest’ (v4.6-14)⁹² package in combination with the ‘caret’ package (v6.0-86)⁹³, for model fitting, performance metric calculation and cross-validation. RF classifiers were trained using constant default parameters: number of trees ‘ntree’ = 1000, and the number of variables ‘mtry = number of metabolic variables/3’. For the random forest regression models, ‘mtry = sqrt(number of metabolic variables)’. A repeated (n = 15) stratified fivefold cross-validation procedure was used for all models. This approach ensures that the train and test sets contain a similar proportion of samples from each class.

The performance of the RF models was assessed by calculating the R² using only data left out from each CV round test sets (performed by default by ‘caret’). For the classifiers, ROC curves, precision–recall curves, accuracy, sensitivity, specificity, positive-predictive value and negative predictive values were calculated and plotted with the R packages ‘caret’, ‘pROC’ (v1.16.2)⁹⁴, ‘plotROC’ (v2.2.1)⁹⁵, ‘precrec’ (v0.11.2)⁹⁶ and ‘ggplot2’ (v3.3.2)⁹⁷. For the regression analysis, the R² was calculated with ‘caret’ and predictions plotted with ‘ggplot2’.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

The metabolic profiling data generated in this study have been deposited and made publicly available in the MetaboLights database under study identifier MTBLS717. The sequence data for the study are publicly available through the European Nucleotide Archive [https://www.ebi.ac.uk/ena] under accession numbers PRJEB 11895, 12577 and (https://www.ebi.ac.uk/ena/browser/view/PRJEB41427). Relevant clinical and patient metadata are publicly available in the GitHub repository at https://www.github.com/gscorreia89/desims-cst-analysis/, digital resource identifier: https://doi.org/10.5281/zenodo.5513501⁹⁸.

Code availability

Key analysis code and processed datasets required to reproduce the statistical analyses presented in this study are deposited and publicly available in the GitHub repository at https://www.github.com/gscorreia89/desims-cst-analysis/, digital resource identifier: https://doi.org/10.5281/zenodo.5513501⁹⁸.

References

Ravel, J. et al. Daily temporal dynamics of vaginal microbiota before, during and after episodes of bacterial vaginosis. Microbiome 1, 29 (2013).
Article PubMed PubMed Central Google Scholar
Borgdorff, H. et al. Lactobacillus-dominated cervicovaginal microbiota associated with reduced HIV/STI prevalence and genital HIV viral load in African women. ISME J. 8, 1781–1793 (2014).
Article CAS PubMed PubMed Central Google Scholar
Klatt, N. R. et al. Vaginal bacteria modify HIV tenofovir microbicide efficacy in African women. Science 356, 938–945 (2017).
Article ADS CAS PubMed Google Scholar
van de Wijgert, J. The vaginal microbiome and sexually transmitted infections are interlinked: consequences for treatment and prevention. PLoS Med. 14, e1002478 (2017).
Article PubMed PubMed Central Google Scholar
Norenhag, J. et al. The vaginal microbiota, human papillomavirus and cervical dysplasia: a systematic review and network meta-analysis. BJOG 127, 171–180 (2020).
Article CAS PubMed Google Scholar
Mitra, A. et al. The vaginal microbiota associates with the regression of untreated cervical intraepithelial neoplasia 2 lesions. Nat. Commun. 11, 1999 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Koedooder, R. et al. The vaginal microbiome as a predictor for outcome of in vitro fertilization with or without intracytoplasmic sperm injection: a prospective study. Hum. Reprod. 34, 1042–1054 (2019).
Article CAS PubMed Google Scholar
Al-Memar, M. et al. The association between vaginal bacterial composition and miscarriage: a nested case-control study. BJOG 127, 264–274 (2020).
Article CAS PubMed Google Scholar
Brown, R. G. et al. Vaginal dysbiosis increases risk of preterm fetal membrane rupture, neonatal sepsis and is exacerbated by erythromycin. BMC Med. 16, 9 (2018).
Article PubMed PubMed Central CAS Google Scholar
Callahan, B. J. et al. Replication and refinement of a vaginal microbial signature of preterm birth in two racially distinct cohorts of US women. Proc. Natl Acad. Sci. USA 114, 9966–9971 (2017).
Article CAS PubMed PubMed Central Google Scholar
Elovitz, M. A. et al. Cervicovaginal microbiota and local immune response modulate the risk of spontaneous preterm delivery. Nat. Commun. 10, 1305 (2019).
Article ADS PubMed PubMed Central CAS Google Scholar
Fettweis, J. M. et al. The vaginal microbiome and preterm birth. Nat. Med. 25, 1012–1021 (2019).
Article CAS PubMed PubMed Central Google Scholar
Kindinger, L. M. et al. The interaction between vaginal microbiota, cervical length, and vaginal progesterone treatment for preterm birth risk. Microbiome 5, 6 (2017).
Article PubMed PubMed Central Google Scholar
Amsel, R. et al. Nonspecific vaginitis. Diagnostic criteria and microbial and epidemiologic associations. Am. J. Med. 74, 14–22 (1983).
Article CAS PubMed Google Scholar
Nugent, R. P., Krohn, M. A. & Hillier, S. L. Reliability of diagnosing bacterial vaginosis is improved by a standardized method of gram stain interpretation. J. Clin. Microbiol. 29, 297–301 (1991).
Article CAS PubMed PubMed Central Google Scholar
Pruski, P. et al. Assessment of microbiota:host interactions at the vaginal mucosa interface. Methods 149, 74–84 (2018).
Article CAS PubMed Google Scholar
Anahtar, M. N. et al. Cervicovaginal bacteria are a major modulator of host inflammatory responses in the female genital tract. Immunity 42, 965–976 (2015).
Article CAS PubMed PubMed Central Google Scholar
Hearps, A. C. et al. Vaginal lactic acid elicits an anti-inflammatory response from human cervicovaginal epithelial cells and inhibits production of pro-inflammatory mediators associated with HIV acquisition. Mucosal Immunol. 10, 1480–1490 (2017).
Article CAS PubMed Google Scholar
Goldenberg, R. L., Culhane, J. F., Iams, J. D. & Romero, R. Epidemiology and causes of preterm birth. Lancet 371, 75–84 (2008).
Article PubMed PubMed Central Google Scholar
Romero, R., Dey, S. K. & Fisher, S. J. Preterm labor: one syndrome, many causes. Science 345, 760–765 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Liu, L. et al. Global, regional, and national causes of under-5 mortality in 2000-15: an updated systematic analysis with implications for the sustainable development goals. Lancet 388, 3027–3035 (2016).
Article PubMed PubMed Central Google Scholar
Kindinger, L. M. et al. Relationship between vaginal microbial dysbiosis, inflammation, and pregnancy outcomes in cervical cerclage. Sci. Transl. Med. 8, 350ra102 (2016).
Article PubMed Google Scholar
Ravel, J. et al. Vaginal microbiome of reproductive-age women. Proc. Natl Acad. Sci. USA 108(Suppl 1), 4680–4687 (2011).
Article ADS CAS PubMed Google Scholar
Serrano, M. G. et al. Racioethnic diversity in the dynamics of the vaginal microbiome during pregnancy. Nat. Med. 25, 1001–1011 (2019).
Article CAS PubMed PubMed Central Google Scholar
Georgiou, H. M., Di Quinzio, M. K., Permezel, M. & Brennecke, S. P. Predicting preterm labour: current status and future prospects. Dis. Markers 2015, 435014 (2015).
Article PubMed PubMed Central Google Scholar
Pruski, P. et al. Medical swab analysis using desorption electrospray ionization mass spectrometry: a noninvasive approach for mucosal diagnostics. Anal. Chem. 89, 1540–1550 (2017).
Article CAS PubMed Google Scholar
Takats, Z., Wiseman, J. M., Gologan, B. & Cooks, R. G. Mass spectrometry sampling under ambient conditions with desorption electrospray ionization. Science 306, 471–473 (2004).
Article ADS CAS PubMed Google Scholar
DiGiulio, D. B. et al. Temporal and spatial variation of the human microbiota during pregnancy. Proc. Natl Acad. Sci. USA 112, 11060–11065 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Brown, R. G. et al. Prospective observational study of vaginal microbiota pre- and post-rescue cervical cerclage. BJOG 126, 916–925 (2019).
Article CAS PubMed PubMed Central Google Scholar
Frank, J. A. et al. Critical evaluation of two primers commonly used for amplification of bacterial 16 S rRNA genes. Appl Environ. Microbiol. 74, 2461–2470 (2008).
Article ADS CAS PubMed PubMed Central Google Scholar
Ma, X. & Ouyang, Z. Ambient ionization and miniature mass spectrometry system for chemical and biological analysis. Trends Anal. Chem. 85, 10–19 (2016).
Article CAS Google Scholar
Ferreira, C. R. et al. Ambient ionization mass spectrometry for point-of-care diagnostics and other clinical measurements. Clin. Chem. 62, 99–110 (2016).
Article CAS PubMed Google Scholar
McMillan, A. et al. A multi-platform metabolomics approach identifies highly specific biomarkers of bacterial diversity in the vagina of pregnant and non-pregnant women. Sci. Rep. 5, 14174 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Srinivasan, S. et al. Metabolic signatures of bacterial vaginosis. MBio 6, e00204-15 (2015).
Article PubMed PubMed Central CAS Google Scholar
Vitali, B. et al. Vaginal microbiome and metabolome highlight specific signatures of bacterial vaginosis. Eur. J. Clin. Microbiol. Infect. Dis. 34, 2367–2376 (2015).
Article CAS PubMed Google Scholar
Nelson, T. M. et al. Vaginal biogenic amines: biomarkers of bacterial vaginosis or precursors to vaginal dysbiosis? Front. Physiol. 6, 253 (2015).
Article PubMed PubMed Central Google Scholar
Ilhan, Z. E. et al. Deciphering the complex interplay between microbiota, HPV, inflammation and cancer through cervicovaginal metabolic profiling. EBioMedicine 44, 675–690 (2019).
Article PubMed PubMed Central Google Scholar
Al-Mushrif, S., Eley, A. & Jones, B. M. Inhibition of chemotaxis by organic acids from anaerobes may prevent a purulent response in bacterial vaginosis. J. Med. Microbiol. 49, 1023–1030 (2000).
Article CAS PubMed Google Scholar
Aldunate, M. et al. Vaginal concentrations of lactic acid potently inactivate HIV. J. Antimicrob. Chemother. 68, 2015–2025 (2013).
Article CAS PubMed PubMed Central Google Scholar
Mirmonsef, P. et al. Short-chain fatty acids induce pro-inflammatory cytokine production alone and in combination with toll-like receptor ligands. Am. J. Reprod. Immunol. 67, 391–400 (2012).
Article CAS PubMed Google Scholar
O’Hanlon, D. E., Moench, T. R. & Cone, R. A. In vaginal fluid, bacteria associated with bacterial vaginosis can be suppressed with lactic acid but not hydrogen peroxide. BMC Infect. Dis. 11, 200 (2011).
Article PubMed PubMed Central CAS Google Scholar
Wolrath, H., Boren, H., Hallen, A. & Forsum, U. Tirimethylamine content in vaginal secretion and its relation to bacterial vaginosis. Apmis 110, 819–824 (2002).
Article CAS PubMed Google Scholar
Yeoman, C. J. et al. A multi-omic systems-based approach reveals metabolic markers of bacterial vaginosis and insight into the disease. PLoS ONE 8, e56111 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Nazli, A. et al. Exposure to HIV-1 directly impairs mucosal epithelial barrier integrity allowing microbial translocation. PLoS Pathog. 6, e1000852 (2010).
Article PubMed PubMed Central CAS Google Scholar
Rutkowsky, J. M. et al. Acylcarnitines activate proinflammatory signaling pathways. Am. J. Physiol. Endocrinol. Metab. 306, E1378–1387 (2014).
Article CAS PubMed PubMed Central Google Scholar
Sanches, J. M. et al. Vaginal lipidomics of women with vulvovaginal candidiasis and cytolytic vaginosis: a non-targeted LC-MS pilot study. PLoS ONE 13, e0202401 (2018).
Article PubMed PubMed Central CAS Google Scholar
Borgogna, J. C. et al. Authors’ reply: The vaginal metabolome and microbiota of cervical HPV-positive and HPV-negative women: a cross-sectional analysis. BJOG 127, 773–774 (2020).
Article PubMed Google Scholar
Borgogna, J. C. et al. The vaginal metabolome and microbiota of cervical HPV-positive and HPV-negative women: a cross-sectional analysis. BJOG 127, 182–192 (2020).
Article CAS PubMed Google Scholar
Broadhurst, D. et al. Guidelines and considerations for the use of system suitability and quality control samples in mass spectrometry assays applied in untargeted clinical metabolomic studies. Metabolomics 14, 72 (2018).
Article PubMed PubMed Central CAS Google Scholar
Kosti, I., Lyalina, S., Pollard, K. S., Butte, A. J. & Sirota, M. Meta-analysis of vaginal microbiome data provides new insights into preterm birth. Front. Microbiol. 11, 476 (2020).
Article PubMed PubMed Central Google Scholar
Witkin, S. S., Linhares, I. M. & Giraldo, P. Bacterial flora of the female genital tract: function and immune regulation. Best. Pract. Res. Clin. Obstet. Gynaecol. 21, 347–354 (2007).
Article PubMed Google Scholar
Jespers, V. et al. A longitudinal analysis of the vaginal microbiota and vaginal immune mediators in women from sub-Saharan Africa. Sci. Rep. 7, 11974 (2017).
Article ADS PubMed PubMed Central CAS Google Scholar
Campisciano, G., Zanotta, N., Licastro, D., De Seta, F. & Comar, M. In vivo microbiome and associated immune markers: New insights into the pathogenesis of vaginal dysbiosis. Sci. Rep. 8, 2307 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar
Hedges, S. R., Barrientes, F., Desmond, R. A. & Schwebke, J. R. Local and systemic cytokine levels in relation to changes in vaginal flora. J. Infect. Dis. 193, 556–562 (2006).
Article CAS PubMed Google Scholar
Mattsby-Baltzer, I., Platz-Christensen, J. J., Hosseini, N. & Rosen, P. IL-1 beta, IL-6, TNF alpha, fetal fibronectin, and endotoxin in the lower genital tract of pregnant women with bacterial vaginosis. Acta Obstet. Gynecol. Scand. 77, 701–706 (1998).
CAS PubMed Google Scholar
Platz-Christensen, J. J., Mattsby-Baltzer, I., Thomsen, P. & Wiqvist, N. Endotoxin and interleukin-1 alpha in the cervical mucus and vaginal fluid of pregnant women with bacterial vaginosis. Am. J. Obstet. Gynecol. 169, 1161–1166 (1993).
Article CAS PubMed Google Scholar
Spandorfer, S. D., Neuer, A., Giraldo, P. C., Rosenwaks, Z. & Witkin, S. S. Relationship of abnormal vaginal flora, proinflammatory cytokines and idiopathic infertility in women undergoing IVF. J. Reprod. Med. 46, 806–810 (2001).
CAS PubMed Google Scholar
Wennerholm, U. B. et al. Interleukin-1 alpha, interleukin-6 and interleukin-8 in cervico/vaginal secretion for screening of preterm birth in twin gestation. Acta Obstet. Gyn Scan 77, 508–514 (1998).
CAS Google Scholar
Pellis, V. et al. Mannose binding lectin and C3 act as recognition molecules for infectious agents in the vagina. Clin. Exp. Immunol. 139, 120–126 (2005).
Article CAS PubMed PubMed Central Google Scholar
Giraldo, P. C. et al. Mannose-binding lectin gene polymorphism, vulvovaginal candidiasis, and bacterial vaginosis. Obstet. Gynecol. 109, 1123–1128 (2007).
Article CAS PubMed Google Scholar
Kalia, N., Singh, J., Sharma, S. & Kaur, M. Impact of SNPs interplay across the locus of MBL2, between MBL and Dectin-1 gene, on women’s risk of developing recurrent vulvovaginal infections. Cell Biosci. 9, 35 (2019).
Article PubMed PubMed Central Google Scholar
Cohen, C. R. et al. Randomized trial of Lactin-V to prevent recurrence of bacterial vaginosis. N. Engl. J. Med. 382, 1906–1915 (2020).
Article CAS PubMed PubMed Central Google Scholar
Lev-Sagie, A. et al. Vaginal microbiome transplantation in women with intractable bacterial vaginosis. Nat. Med. 25, 1500–1504 (2019).
Article CAS PubMed Google Scholar
Walker, A. W. et al. 16S rRNA gene-based profiling of the human infant gut microbiota is strongly influenced by sample processing and PCR primer choice. Microbiome 3, 26 (2015).
Article PubMed PubMed Central Google Scholar
France, M. T. et al. VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition. Microbiome 8, 166 (2020).
Article PubMed PubMed Central Google Scholar
Lee, S. et al. Community state types of vaginal microbiota and four types of abnormal vaginal microbiota in pregnant Korean women. Front Public Health 8, 507024 (2020).
Article PubMed PubMed Central Google Scholar
Tettamanti Boshier, F. A. et al. Complementing 16S rRNA gene amplicon sequencing with total bacterial load to infer absolute species concentrations in the vaginal microbiome. mSystems 5, https://doi.org/10.1128/mSystems.00777-19 (2020).
Izzi-Engbeaya, C. et al. The effects of kisspeptin on beta-cell function, serum metabolites and appetite in humans. Diabetes Obes. Metab. 20, 2800–2810 (2018).
Article CAS PubMed PubMed Central Google Scholar
Lewis, M. R. et al. Development and application of ultra-performance liquid chromatography-TOF MS for precision large scale urinary metabolic phenotyping. Anal. Chem. 88, 9004–9013 (2016).
Article CAS PubMed Google Scholar
Kessner, D., Chambers, M., Burke, R., Agus, D. & Mallick, P. ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics 24, 2534–2536 (2008).
Article CAS PubMed PubMed Central Google Scholar
Gibb, S. & Strimmer, K. MALDIquant: a versatile R package for the analysis of mass spectrometry data. Bioinformatics 28, 2270–2271 (2012).
Article CAS PubMed Google Scholar
Smith, C. A., Want, E. J., O’Maille, G., Abagyan, R. & Siuzdak, G. XCMS: processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 78, 779–787 (2006).
Article CAS PubMed Google Scholar
Sands, C. J. et al. The nPYc-Toolbox, a Python module for the pre-processing, quality-control and analysis of metabolic profiling datasets. Bioinformatics 35, 5359–5360 (2019).
Article CAS PubMed PubMed Central Google Scholar
Smith, C. A. et al. METLIN: a metabolite mass spectral database. Ther. Drug Monit. 27, 747–751 (2005).
Article CAS PubMed Google Scholar
Fahy, E. A comprehensive classification system for lipids (vol 46, pg 839, 2005). J. Lipid Res. 51, 1618–1618 (2010).
CAS Google Scholar
Wishart, D. S. et al. HMDB 3.0—the human metabolome database in 2013. Nucleic Acids Res. 41, D801–D807 (2013).
Article CAS PubMed Google Scholar
Sumner, L. W. et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics 3, 211–221 (2007).
Article CAS PubMed PubMed Central Google Scholar
MacIntyre, D. A. et al. The vaginal microbiome during pregnancy and the postpartum period in a European population. Sci. Rep. 5, 8988 (2015).
Article CAS PubMed PubMed Central Google Scholar
Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541 (2009).
Article ADS CAS PubMed PubMed Central Google Scholar
Wang, Q., Garrity, G. M., Tiedje, J. M. & Cole, J. R. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261 (2007).
Article ADS CAS PubMed PubMed Central Google Scholar
Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010).
Article CAS PubMed Google Scholar
Fettweis, J. M. et al. Species-level classification of the vaginal microbiome. BMC Genomics 13(Suppl 8), S17 (2012).
Article PubMed PubMed Central Google Scholar
Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272 (2020).
Article CAS PubMed PubMed Central Google Scholar
Hunter, J. D. Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90–95 (2007).
Article Google Scholar
Waskom, M. et al. mwaskom/seaborn (v0.11.0). Zenodo https://doi.org/10.5281/zenodo.4019146 (2020).
Schulthess, B. et al. Identification of Gram-positive cocci by use of matrix-assisted laser desorption ionization-time of flight mass spectrometry: comparison of different preparation methods and implementation of a practical algorithm for routine diagnostics. J. Clin. Microbiol. 51, 1834–1840 (2013).
Article CAS PubMed PubMed Central Google Scholar
Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48 (2015).
Article Google Scholar
Lenth, R. et al. Estimated marginal means, aka least-squares means (v1.5.2-1). https://cran.r-project.org/web/packages/emmeans (2020).
Halekoh, U. & Højsgaard, S. A Kenward-Roger approximation and parametric bootstrap methods for tests in linear mixed models the R Package pbkrtest. J. Stat. Softw. 59, (2012).
McArdle, B. H. & Anderson, M. J. Fitting multivariate models to community data: a comment on distance-based redundancy analysis. Ecology 82, 290–297 (2001).
Article Google Scholar
Quinn, T. P., Richardson, M. F., Lovell, D. & Crowley, T. M. propr: an R-package for identifying proportionally abundant features using compositional data analysis. Sci. Rep. 7, 16252 (2017).
Article ADS PubMed PubMed Central CAS Google Scholar
Liaw, A. & Wiener, M. randomForest: Breiman and Cutler’s Random Forests for classification and regression. https://cran.r-project.org/web/packages/randomForest (2018).
Kuhn, M. caret: classification and regression training (v6.0-86). https://cran.r-project.org/web/packages/caret (2020).
Robin, X. et al. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics 12, 77 (2011).
Article PubMed PubMed Central Google Scholar
Sachs, M. C. plotROC: a tool for plotting roc curves. J. Stat. Softw. 79, 2 (2017).
Saito, T. & Rehmsmeier, M. Precrec: fast and accurate precision-recall and ROC curve calculations in R. Bioinformatics 33, 145–147 (2017).
Article CAS PubMed Google Scholar
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).
Correia, G. D., Pruski, P. & MacIntyre, D. A. Gscorreia89/desims-cst-analysis: supporting data analysis repository for the manuscript “Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth” (v1.0). Zenodo https://doi.org/10.5281/zenodo.5513501 (2021).

Download references

Acknowledgements

We would like to thank all patients who have participated in this study and members of the Women’s Health Research Centre who facilitated and coordinated study recruitment and sample collection. We also acknowledge Lili Herendi, Verena Horneffer van der Sluis and Maria Gomez for coordinating and running LC-MS based assays at the Clinical Phenotyping Center (CPC) at Imperial College London, UK, and the NIHR Imperial Biomedical Research Centre (BRC), which provided infrastructure support. This work was supported by the European Research Council (Grant No. 617896), the Imperial College NIHR Biomedical Research Centre (Grant No. 58434), the Medical Research Council (Grant No. MR/L009226/1), NIHR Clinical Lectureship Scheme, Genesis Research Trust and the March of Dimes European Preterm Birth Research Centre at Imperial College London. A.L.D. is supported by the NIHR University College London Hospitals Biomedical Research Centre. The views expressed in this publication are those of the authors and not necessarily those of the NHS, NIHR or Department of Health and Social Care.

Author information

These authors contributed equally: Pamela Pruski, Gonçalo D. S. Correia.

Authors and Affiliations

Division of Systems Medicine, Department of Metabolism, Digestion and Reproduction, Faculty of Medicine Imperial College London, London, UK
Pamela Pruski, Gonçalo D. S. Correia, Paolo Inglese, Julian Marchesi, Simon Cameron, Kate Alexander-Hardiman & Zoltan Takats
National Phenome Centre, Imperial College London, London, UK
Gonçalo D. S. Correia, Paolo Inglese & Zoltan Takats
March of Dimes Prematurity Research Centre at Imperial College London, London, UK
Holly V. Lewis, Katia Capuccini, Denise Chan, Yun S. Lee, Julian Marchesi, Vasso Terzidou, Lynne Sykes, Phillip R. Bennett, Zoltan Takats & David A. MacIntyre
Imperial College Parturition Research Group, Institute of Reproductive and Developmental Biology, Department of Metabolism, Digestion and Reproduction, Imperial College London, London, UK
Holly V. Lewis, Katia Capuccini, Denise Chan, Richard G. Brown, Lindsay Kindinger, Yun S. Lee, Vasso Terzidou, Lynne Sykes, Phillip R. Bennett & David A. MacIntyre
Queen Charlotte’s & Chelsea Hospital, Imperial College London, London, UK
Holly V. Lewis, Denise Chan, Richard G. Brown, Lynne Sykes & Phillip R. Bennett
Elizabeth Garrett Anderson Institute for Women’s Health, University College London, London, UK
Lindsay Kindinger & Anna L. David
Faculty of Health and Applied Sciences, University West of England, Bristol, UK
Ann Smith
MRC Centre for Molecular Bacteriology and Infection, Imperial College London, London, UK
Julie A. K. McDonald
School of Biological Sciences, Institute for Global Food Security, Queen’s University Belfast, Belfast, UK
Simon Cameron
MRC Centre for Reproductive Health, University of Edinburgh, Edinburgh, UK
Sarah J. Stock & Jane E. Norman
Faculty of Health Sciences, University of Bristol, Bristol, UK
Jane E. Norman
Chelsea & Westminster Hospital, NHS Trust, London, UK
Vasso Terzidou
St Mary’s Hospital, Imperial College Healthcare NHS Trust, London, UK
T. G. Teoh
Tommy’s National Centre for Miscarriage Research, Imperial College London, London, UK
Phillip R. Bennett & David A. MacIntyre

Authors

Pamela Pruski
View author publications
You can also search for this author in PubMed Google Scholar
Gonçalo D. S. Correia
View author publications
You can also search for this author in PubMed Google Scholar
Holly V. Lewis
View author publications
You can also search for this author in PubMed Google Scholar
Katia Capuccini
View author publications
You can also search for this author in PubMed Google Scholar
Paolo Inglese
View author publications
You can also search for this author in PubMed Google Scholar
Denise Chan
View author publications
You can also search for this author in PubMed Google Scholar
Richard G. Brown
View author publications
You can also search for this author in PubMed Google Scholar
Lindsay Kindinger
View author publications
You can also search for this author in PubMed Google Scholar
Yun S. Lee
View author publications
You can also search for this author in PubMed Google Scholar
Ann Smith
View author publications
You can also search for this author in PubMed Google Scholar
Julian Marchesi
View author publications
You can also search for this author in PubMed Google Scholar
Julie A. K. McDonald
View author publications
You can also search for this author in PubMed Google Scholar
Simon Cameron
View author publications
You can also search for this author in PubMed Google Scholar
Kate Alexander-Hardiman
View author publications
You can also search for this author in PubMed Google Scholar
Anna L. David
View author publications
You can also search for this author in PubMed Google Scholar
Sarah J. Stock
View author publications
You can also search for this author in PubMed Google Scholar
Jane E. Norman
View author publications
You can also search for this author in PubMed Google Scholar
Vasso Terzidou
View author publications
You can also search for this author in PubMed Google Scholar
T. G. Teoh
View author publications
You can also search for this author in PubMed Google Scholar
Lynne Sykes
View author publications
You can also search for this author in PubMed Google Scholar
Phillip R. Bennett
View author publications
You can also search for this author in PubMed Google Scholar
Zoltan Takats
View author publications
You can also search for this author in PubMed Google Scholar
David A. MacIntyre
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

P.P., G.D.S.C., L.S., P.R.B., Z.T. and D.A.M. conceived and designed the study. Clinical sampling and coordination of metadata was performed by H.V.L., D.C., R.G.B., L.K., A.L.D., S.J.S., L.S., V.T. and P.R.B. Bacterial isolates were obtained from S.C., K.A-H., J.A.K.M.D. and J.M. Experiments were performed by P.P., H.V.L., K.C., D.C., R.B., L.K., Y.S.L., J.A.K.M.D. and K.H. Data processing and analysis was performed G.D.S.C., P.P., P.I., A.S., Z.T. and D.A.M. P.P., G.D.S.C. and D.A.M. prepared all figures and tables and wrote the first draft of the manuscript. All authors critically reviewed, read and approved the final manuscript.

Corresponding authors

Correspondence to Zoltan Takats or David A. MacIntyre.

Ethics declarations

Competing interests

P.P., P.R.B., Z.T. and D.A.M. hold patents for the use of rapid evaporative ionization mass spectrometry and desorption electrospray ionization mass spectrometry analysis of swabs and biopsy samples (US10026599B2, EP3265817B1). P.P., G.C., P.R.B., Z.T. and D.A.M. have filed a provisional patent for the use of rapid evaporative ionization mass spectrometry and desorption electrospray ionization mass spectrometry analysis of swabs for prediction of vaginal microbiota composition and inflammatory status (GB2110293.4). The remaining authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

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

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Dataset 1

Reporting Summary

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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Pruski, P., Correia, G.D.S., Lewis, H.V. et al. Direct on-swab metabolic profiling of vaginal microbiome host interactions during pregnancy and preterm birth. Nat Commun 12, 5967 (2021). https://doi.org/10.1038/s41467-021-26215-w

Download citation

Received: 14 March 2021
Accepted: 23 September 2021
Published: 13 October 2021
DOI: https://doi.org/10.1038/s41467-021-26215-w

This article is cited by

Effect of transvaginal Lactobacillus supplementation on reversing lower genital tract dysbiosis and improving perinatal outcomes in PCOS patients after IVF-FET: a study protocol for a multicenter randomized controlled trial
- Huixi Chen
- Yaoyao Tu
- He-Feng Huang
Trials (2023)
Microbial signatures in amniotic fluid at preterm birth and association with bronchopulmonary dysplasia
- Birte Staude
- Silvia Gschwendtner
- Harald Ehrhardt
Respiratory Research (2023)
Small molecule metabolites: discovery of biomarkers and therapeutic targets
- Shi Qiu
- Ying Cai
- Aihua Zhang
Signal Transduction and Targeted Therapy (2023)
Preterm birth is associated with xenobiotics and predicted by the vaginal metabolome
- William F. Kindschuh
- Federico Baldini
- Tal Korem
Nature Microbiology (2023)
N-glycosylation of cervicovaginal fluid reflects microbial community, immune activity, and pregnancy status
- Gang Wu
- Paola Grassi
- Stuart M. Haslam
Scientific Reports (2022)

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.