The association between obesity and weight loss after bariatric surgery on the vaginal microbiota

doi:10.21203/rs.3.rs-31292/v2

Background: Obesity and vaginal microbiome (VMB) dysbiosis are each risk factors for adverse reproductive and oncological health outcomes in women. Here we investigated the relationship between obesity, vaginal bacterial composition, local inflammation and bariatric surgery.

Methods: Vaginal bacterial composition assessed by high-throughput sequencing of bacterial 16S rRNA genes and local cytokine levels measured using a multiplexed Magnetic Luminex Screening Assay were compared between 67 obese and 42 non-obese women. We further assessed temporal changes in the microbiota and cytokines in a subset of 27 women who underwent bariatric surgery.

Results: The bacterial component of the vaginal microbiota in obese women was characterised by a lower prevalence of a Lactobacillus-dominant VMB and higher prevalence of a high diversity (Lactobacillus spp., and Gardnerella- spp. depleted) VMB, compared with non-obese subjects (p<0.001). Obese women had higher relative abundance of Dialister species (p<0.001), Anaerococcus vaginalis (p=0.021) and Prevotella timonensis (p=0.020) and decreased relative abundance of Lactobacillus crispatus (p=0.014). Local vaginal IL-1β, IL-4, IL-6, IL-8, IFNγ, MIP-1α, and TNFα levels were all higher among obese women, however only IL-1β and IL-8 correlated with VMB species diversity. In a subset of obese women undergoing bariatric surgery, there were no significant overall differences in VMB following surgery, however 75% of these women remained obese at six months. Prior to surgery there was no relationship between body mass index (BMI) and VMB structure, however post-surgery women with a Lactobacillus-dominant VMB had a significantly lower BMI than those with a high diversity VMB.

Conclusions: Obese women have a significantly different vaginal microbiota composition with increased levels of local inflammation compared to non-obese women. Bariatric surgery does not change the VMB, however, those with the greatest weight loss six-months post-surgery are most likely to have a Lactobacillus-dominant VMB.

General Microbiology

obesity

overweight

body mass index

BMI

vaginal microbiota

bariatric surgery

Obesity has become a worldwide problem, with projections of 1 in 5 women being obese and 1 in 10 morbidly obese by 2025 [1]. The gynaecological complications of obesity include menstrual disorders [2, 3], anovulation [4], polycystic ovarian syndrome [5, 6], infertility [7-9], early pregnancy loss [10-12], preterm birth [13], obstetric complications [13, 14] and increased risk for gynaecological malignancies (endometrial and ovarian) [15-18]. The mechanism by which obesity promotes several of these outcomes is not fully understood.

There is growing recognition that microbiota (i.e. the assemblage of microorganisms found at a specific environment [19]) and their functions influence disease pathophysiology [20-22]. Shifts in microbiota composition at mucosal surfaces can lead to pathobiont overgrowth and activation of innate immune responses that in turn are modified by microbial products such as short-chain fatty acids, lipids and bioamines [23-26]. The vaginal microbiota is commonly dominated by Lactobacillus spp. [27, 28], which offer protection against colonisation of pathogenic bacteria through production of lactic acid, lowering of vaginal pH, production of antimicrobial compounds and modulation of both the immunological and physical properties of cervicovaginal mucosa [29, 30]. High-throughput sequencing approaches have aided characterisation of the vaginal microbiota in health and disease states, leading to a better understanding of the factors which affect vaginal community composition including age and ethnicity [31], menstrual cycle phase [32], oestrogen levels and menopause status [33], intercourse [34], pregnancy [35, 36] and hygiene practices [37].

Bacterial vaginosis (BV), characterised by a shift from Lactobacillus spp. dominance towards high relative abundance of anaerobes, has long been associated with a range of adverse outcomes [38]. Molecular-based characterisation of vaginal bacteria composition has extended these findings to permit identification of specific taxa that modify risk of preterm birth [39-42] and preterm premature rupture of membranes [43], sexually transmitted infections [44], human papilloma virus and cervical disease [45-47].

Obesity and specific vaginal microbiome (VMB) compositions are each risk factors for adverse reproductive and oncological health outcomes in women, but there is limited evidence describing the relationship between vaginal microbiota composition and obesity and the impact that body weight may have on local inflammation, immune response and health outcomes. Brookheart and colleagues found that overweight and obese women have higher Nugent scores and a greater occurrence of BV [48]. Conversely using data from the USA National Health and Nutrition Examination Survey 2001-2004, Koumans and co-workers concluded that body mass index (BMI) was not an independent risk factor for BV[49]. Two Korean studies have examined the relationship between VMB, assessed by metataxanomic analysis, and obesity. Oh and colleagues found that obesity associates with cervical microbiota dominated by Lactobacillus iners [50]. Si and co-workers reported that discordant twin obesity associated with increased bacterial diversity and prevalence of Prevotella [51].

In this study, we investigated how obesity, defined as BMI ³30kg/m², associates with bacterial structure and composition of the vaginal microbiome compared to non-obese women and local inflammation. In a secondary analysis, we investigated the effects of bariatric surgery on vaginal bacterial diversity and local inflammation.

A total of 109 women were prospectively recruited into this study, 67 were non-obese and 42 were obese (Table 1). Non-obese women were more likely to be nulliparous (p<0.001), non-diabetic (p=0.001), using contraception (p=0.016) and Caucasian (p=0.024) compared to obese women. There was no difference in mean age (p=0.083), smoking (p=0.128) or menopause status (p=0.233) (Fisher’s exact test). No differences were found between recent antibiotic use (p=0.258) or time since last sexual intercourse prior to sample collection (p=0.070) (data not shown).

A subset of twenty-seven of the 109 recruited women were scheduled for bariatric surgery and were sampled at baseline (n=27), 3 months (n=22) and 6 months (n=19) post-surgery. Cohort characteristics are shown in Supplementary Table 1. The mean age was 48 years (SD 8.91, range 28-65 years), and the majority of women in the bariatric surgery cohort were Caucasian (18/27, 66.7%), premenopausal (17/27, 63.0%), and only one postmenopausal participant used HRT. One third of the cohort were diabetic (8/27, 29.6%), and half of the women who were diabetic required insulin. We found that 44% of the cohort were insulin resistant at baseline.

In the 16S rRNA gene sequencing analysis, a total of 3,045,614 reads were captured from 166 samples (109 baseline, 57 longitudinal samples) with an average number of reads per sample of 18,021. To avoid sequencing bias, operational taxonomic units (OTUs) were randomly sub-sampled to the lowest common read count of 1,885 with coverage of greater than 95% (Good’s coverage index) being maintained for all samples. A total of 265 taxa were detected in the study cohort microbiota after the removal of singletons and rare OTUs. Rare OTUs were defined as those present at less than 10 counts within the entire cohort, and along with singletons, were included in the final analysis to maintain consistent read counts across all samples.

Vaginal microbiota composition and local cytokine expression according to obesity, diabetes and insulin resistance status

Ward hierarchical clustering analysis of genera-level data identified three major groups, on the basis of relative bacterial abundance; i) Lactobacillus-dominant - characterised by high relative abundance of Lactobacillus spp., ii) Gardnerella-dominant - characterised by high relative abundance of Gardnerella spp. and low relative abundance of Lactobacillus species, and iii) high diversity VMB - characterised by low relative abundance of each of Lactobacillus and Gardnerella species, and increased bacterial diversity (Figure 1a). Across the whole cohort the prevalence of Lactobacillus-dominant and high diversity VMB groups were significantly different (Figure 1b). The frequency of the bacterial groups within patient groups subcategorised on the basis of obesity status, diabetes and insulin resistance status is presented in Table 2. Prevalence of the high diversity VMB was significantly greater in obese women (obese; 18/42 (42.8%) compared to non-obese; 10/67 (14.9%), p=0.002) whilst the prevalence of Lactobacillus-dominant VMB was significantly lower (Table 2). There was no difference in the prevalence of Gardnerella-dominant VMB between obese and non-obese women. There were overall no significant associations between the three major vaginal bacterial groups and diabetes (Table 2) or where diabetic status was subcategorised into obese or non-obese women (Supplementary Table 2).

When women who had used antibiotics within two weeks of sample collection or had sexual intercourse within 48 hours of sample collection were excluded, the prevalence of both the high diversity and Lactobacillus-dominant VMB remained significantly different in obese women in both categories (excluding intercourse, p=0.001; excluding recent antibiotic use, p=0.003) (Supplementary Table 3). When pre- and post-menopausal women were analysed separately, premenopausal obese women had a significantly higher prevalence of a high diversity VMB (p=0.001), even when oral contraceptive (OCP) users were excluded (p=0.003). There was a significant difference in the prevalence of each genus according to subcategories of normal weight, overweight and obese women (p<0.001) (Supplementary Table 3).

Consistent with an increased prevalence of high diversity VMB (Figure 2a), increased richness (number of species observed) and alpha diversity was observed in obese women (diversity, p= 0.006) (Figure 2b, 2c, Supplementary Table 4). The vaginal microbiota of obese women was characterised by a greater mean proportion of anaerobic bacterial species, specifically unclassified Dialister spp. (unclassified) (p<0.001), Anaerococcus vaginalis (p=0.021) and Prevotella timonensis (p=0.020) (Figure 2d).

To identify vaginal microbiota biomarkers specifically associated with obesity, we performed linear discriminant analysis (LDA) effect size (LefSe) modelling on the 16S rRNA gene sequence data collected from baseline samples (Supplementary Figure 1). Vaginal microbiota of obese women was enriched with members of Bacteroidales and Clostridiales, the Prevotella genus and the phylum Actinobacteria. Conversely, non-obese women were found to have enriched levels of Lactobacillales associated OTUs.

Obese women had significantly increased expression of pro-inflammatory cytokines IL-1, IL-6, IL-8, MIP-1, IFN and TNF compared to non-obese women (Figure 3a, Supplementary Table 5). The anti-inflammatory cytokine IL-4 showed increased expression among obese women. As bacterial diversity increased among obese women (depicted using non-parametric Shannon Index), the expression of IL-1 and IL-8 but not the other cytokines, increased (Figure 3b).

Metabolic and vaginal microbiota compositional changes after bariatric surgery

At 6 months post bariatric surgery, the mean body weight of participants decreased on average by 19.2%, weight loss was similar for pre- and post-menopausal women (Supplementary Table 6). This change translated into a mean reduction of BMI of 10.3 (range 3.9 to 31.3) from a pre-surgery mean of 47.4 to a mean of 35.8 six months post-surgery. Five of 19 women moved from the obese category to the overweight category. The remaining 14 women persisted within the obese category.

In the subset of obese women undergoing bariatric surgery, there were no significant overall differences in VMB following surgery (Figure 4a, Supplementary Figure 3). Neither were changes observed following bariatric surgery according to menopause status, diabetes or insulin resistance status (Table 3). Prior to surgery there was no relationship between BMI and VMB structure, however post-surgery women with Lactobacillus-dominant VMB had a significantly lower BMI than those with a high diversity VMB (Figure 4b). This difference principally applied to pre-menopausal women (Figure 4c). Local cervicovaginal cytokine levels in the bariatric surgery cohort at baseline sampling (n=27) and 6 months post-surgery (n=21) did not show any significant changes (Supplementary Figure 2).

Oestradiol levels and sex hormone binding globulin (SHBG) levels were measured in a subset of 10 bariatric surgery patients at baseline and six months post-surgery (Supplementary Table 6). There were no overall differences in serum oestradiol levels in women pre- and post-surgery. There was no relationship between oestradiol levels and magnitude of weight loss six months post-surgery. Fasting glucose and insulin serum levels were also measured and HOMA-IR (insulin resistance index) was calculated. No associations were seen between oestradiol levels and SHBG levels or HOMA-IR and VMB group.

As oestrogen levels vary during the menstrual cycle and can be affected by exogenous sources of oestrogen such as OCP and hormone replacement therapy (HRT), we performed an analysis to include only postmenopausal women that were not taking HRT (Supplementary Table 6). There was no significant difference in the mean oestradiol levels (mean at baseline sampling was 55.0pmol/L (range 37.0-82.0pmol/L) and mean 6 months post-surgery was 37.0pmol/L (range 37.0-37.0pmol/L), (p=0.371)). SHBG levels increased in the total cohort after 6 months by 39.3% (p=0.008).

In the UK, two thirds of the female population are either overweight (30%), obese (27%) or morbidly obese (4%) [52]. Obesity has been associated with a multitude of adverse health outcomes in women [2, 4, 5, 12, 18], and although the mechanisms leading to these complications of obesity in women remain unclear, the vaginal microbiota composition may be important.

In our cohort, we found three VMB groups at genus level. Approximately half of the vaginal samples were categorised as Lactobacillus-dominant VMB, whilst the remaining samples were categorised in equal proportion as either Gardnerella-dominant (with a high relative abundance of Gardnerella spp. and low relative abundance of Lactobacillus spp.), or high diversity VMB (with a low relative abundance of each of Lactobacillus and Gardnerella spp., and increased bacterial diversity). When these samples were analysed with respect to obesity status, about 70% of obese women demonstrated a Lactobacillus-dominant VMB. This proportion is consistent with other reported studies. Brotman and colleagues showed that a Lactobacillus-dominant VMB is found in 80% of premenopausal women, but only 55% of postmenopausal women [53]. Our study represents a mixture of pre- and post-menopausal women. In those who were categorised as obese, 30% of women had a Lactobacillus-dominant VMB. This prevalence is lower than that found in non-obese postmenopausal women [53]. Our study shows that obesity associates with vaginal microbiota composition, with significantly higher vaginal bacterial species diversity and increased abundance of Dialister, Prevotella and Anaerococcus among obese women. Two previous studies, both from Korean patient cohorts, have described vaginal microbiota composition in relation to obesity in a non-pregnant population. The first study reported that obese women had a greater predominance of Lactobacillus iners compared to non-obese women who were more likely to have a Lactobacillus crispatus-dominant VMB [50]. However, this study was limited to the analysis of interrelationships between cervical Lactobacillus species only with non-Lactobacillus members of the microbiota not considered. Therefore the generalisability of these findings are unclear [50]. Our findings lead us to conclude an increased prevalence of Dialister (unclassified), Prevotella timonensis and Anaerococcus vaginalis in obese women, and dominance of Lactobacillus crispatus among non-obese women (Figure 2d). The second study, by Si et al., reported that obesity was associated with increased levels of Prevotella, and reduced Lactobacillus relative abundance consistent with our data [51].

The vaginal microbiota composition is dynamic and fluctuates throughout the life cycle with relative dominance of the niche by Lactobacillus spp. mediated by oestrogen-driven vaginal epithelium thickening and glycogen deposition, which is used as a primary energy source by lactic acid producing bacteria, encouraging a Lactobacillus-dominant VMB [54]. Accordingly, prepuberty and post-menopause vaginal microbiota composition associates with reduced glycogen levels [55] and a tendency towards a high diversity vaginal microbiota [56, 57]. In our study, vaginal microbiota composition in pre- and post-menopausal women were largely consistent with the published literature. In peri- or post-menopausal women with declining ovarian function, peripheral adipose tissue becomes the major source of production of unopposed circulating oestrogen by aromatisation of adrenal androstenedione to excess endogenous oestrogen [58-60]. In this group of women, reduced systemic oestrogen levels cause a decrease in glycogen deposition resulting in a vaginal epithelia that resembles pre-puberty with a thinner mucus layer and increased incidence of high diversity VMB [61, 62]. However, oestrogen levels are not higher in premenopausal obese women compared to non-obese. Freeman and co-workers found that premenopausal obese and overweight women had significantly lower oestradiol levels compared with non-obese women, independent of age, race, or smoking [63]. Reduced circulating oestrogen concentrations is therefore a potential explanation for the lower prevalence of Lactobacillus-dominant VMB in obese women. It is also possible that due to restrictions in mobility caused by morbid obesity, female hygiene practices are affected. Local skin irritation and breakdown caused by the presence and rubbing of excess adipose tissue, together with persistent moisture may alter the local vaginal microbiota composition. The impact of women’s sexual, sanitary and hygiene practices such as douching on the vaginal microbiota is still controversial [64, 65].

Sustained weight loss may be brought about by significant lifestyle changes (diet and increased physical activity) or induced by bariatric surgery, that can result in improved metabolic health. It has previously been reported that bariatric surgery can reduce cancer incidence [66], improve sex hormone profiles [67], polycystic ovarian syndrome symptoms [68], spontaneous and assisted conception rates [69-72], and reduce obstetric complications such as gestational diabetes, pregnancy-induced hypertension and macrosomia [73, 74]. Reversal or reduction in obesity brought about by bariatric surgery has been shown to be accompanied by metabolic improvement and reduction of alpha-diversity of the gut microbiota within 3 months post-operatively [75-77]. The impact of bariatric surgery on the vaginal microbiome has not however been previously assessed. In our study, in the subset of obese women undergoing bariatric surgery, there were no significant overall differences in VMB following surgery, however 75% of these women remained in the obese range. Prior to surgery there was no relationship between BMI and VMB structure, however, post-surgery women with Lactobacillus-dominant VMB had a significantly lower BMI than those with a high diversity VMB. Significant weight loss following bariatric surgery therefore associates with a tendency toward an ‘optimal’ VMB. We did not, however, find any correlation between degree of weight loss and systemic oestradiol or SHBG concentrations. Additionally, SHBG increased in every woman at 6 months post-surgery which would lead to lower bioavailable oestrogen. It is therefore unlikely any that any effect of bariatric surgery upon the VMB is principally due to changes in systemic oestradiol or SHBG concentrations.

In this study we further explored how obesity may affect local inflammation and whether this is driven by increased diversity in the VMB. Previous studies have reported how a number of immune modulators in obesity (i.e. adipose tissue macrophages, cytokines and adipokines) can drive systemic inflammation [78, 79]. Other studies have shown that genital tract infection lead to changes in local immune factors that in turn have been linked to preterm birth [80] and bacterial vaginosis persistence [81]. Cervicovaginal cytokine levels (IL-1, IL-6, IL-8, TNF, MIP-1 and IFN) were increased in our obese population. This is in line with previous reports suggesting that modulation of the VMB in pregnancy increases local expression of each of these cytokines which correlates with increasing species diversity [82]. In the present study however, there was an association between increasing bacterial species diversity and expression of only two pro-inflammatory cytokines IL-1 and IL-8. This data leads us to conclude that, in obesity, factors other than the VMB act to modulate cervicovaginal inflammation. We found no significant change in cervicovaginal cytokine levels following bariatric surgery, and cytokine concentrations were not different between the different genera groups. This finding again shows that much of the local cervicovaginal inflammation associated with obesity is unrelated to the vaginal microbiota.

Previous cross-sectional studies have highlighted that vaginal microbiota composition may be affected by obesity [50, 51]. A strength of the present study is the assessment of the vaginal microbiota composition in obese pre- and postmenopausal women separately, after exclusion of OCP and HRT-users and after taking into account confounding variables and hormonal serum levels. Furthermore, this is the first study to explore the temporal changes in the vaginal microbiota in a cohort of morbidly obese women undergoing surgically induced weight loss.

Although this study is one of the largest cohorts assessing the impact of obesity on the vaginal microbiota composition in women undergoing bariatric surgery, temporal data collected at each of the three timepoints was only available for a small number of patients which limits our ability to draw meaningful conclusions for many comparisons of interest. Recruitment of much larger numbers would be practically difficult. Follow up for longer periods may show changes not revealed at 6-months, maximum weight loss is usually achieved 12-24 months after surgery [83, 84], although the majority of women who undergo bariatric surgery regain between 5-10% of their pre-operative weight by 2-3 years [84]. Furthermore, although negative controls have been included with each DNA extraction set, no amplicons were identified, and these were not subsequently sequenced. Given that vaginal samples are not considered to be of low biomass and we previously reported sequencing data of negative controls in similar contemporaneous patient cohorts [43, 85], we do not expect this to have led to spurious results. Consistent with this, we did not observe any common kit or environmental contaminants as being prominent features of any of the patient samples sequenced in this study.’

Obesity was found to be associated with higher vaginal microbiome diversity that may partly explain changes in local inflammation. Other factors beyond the VMB (i.e. endocrine) are likely to affect local inflammatory state. Surgery-induced weight loss did not change the VMB composition, although three quarters of the women remained obese six months post-surgery. Those with the greatest weight loss six-months post-surgery were more likely to have a Lactobacillus-dominant VMB. A healthier VMB following bariatric surgery-induced weight loss may create a healthier local microenvironment to promote health, but since this does not associate with a reduction in local inflammation other factors are likely to be involved and require further investigation.

Study population – inclusion and exclusion criteria

We prospectively recruited non-pregnant women attending outpatient gynaecology and bariatric surgery clinics at Imperial College NHS Healthcare Trust between 2013-2016. A subset of this population was scheduled for bariatric surgery. Women were recruited irrespective of age, menopause status, ethnicity, parity, smoking, phase in menstrual cycle and contraception use. Women who were HIV, hepatitis B or C positive, had autoimmune disorders, or had a previous hysterectomy were excluded. Ethical approval was obtained from the National Research Ethics Service Committee London – Fulham (Approval number 13/LO/0126) and the NHS West of Scotland Research Ethics Service Committee (WoSRES) (REC 14/WS/1098). All patients gave informed consent.

Sample collection and processing

Cervicovaginal secretions were collected during the clinic visit from the posterior vaginal fornix with a BBL^TMCultureSwab^TM containing liquid Amies (Becton Dickinson, Oxford, UK) using a sterile, disposable speculum, without lubricant and immediately stored at -80°C. A second transport microbiology swab (Transwab®) containing Amies gel medium was simultaneously collected for cytokine analysis. In the subset population of women planned for bariatric surgery (Gastric band, Roux-en-Y Gastric Bypass or Vertical Sleeve Gastrectomy), we collected serial vaginal swab samples on the day of surgery and at months 3 and 6 post-surgery (Supplementary Figure 4). Serial fasting blood samples were also collected on the day of surgery and at 6-month follow-up with the aim to correlate changes in vaginal microbiota to four serum markers known to be affected by surgery-induced weight loss and hyperinsulinaemia correction. Blood samples were centrifuged at 4,472 g for 10 minutes and serum collected for freezing and storage in -80°C.

A comprehensive interview and questionnaire were used to obtain all relevant gynaecological, medical and surgical history. Menopause status, type of contraception or HRT use and menstrual cycle phase (follicular or luteal) were documented. Ethnicity was self-reported as Caucasian, Asian, Black or Other.

Whole-Genomic bacterial DNA was extracted from the CultureSwab^TM using a QiAmp Mini DNA kit (Qiagen, Venlo, Netherlands) as described previously [36]. The second swab for cytokine analysis was thawed on ice and re-suspended in 350ml phosphate-buffered saline solution with protease inhibitor (5ml/ml; Sigma Aldrich). The suspension was centrifuged at 402g for 2 minutes, and supernatant collected into a new 1.5ml microcentrifuge tube. This centrifugation step was repeated to remove any remaining cellular debris. The cell-free supernatant was stored in -80°C. Negative control swabs (blank, devoid of patient sample) were processed alongside each DNA extraction set. No amplicons were observed following PCR and gel electrophoresis of the negative controls, and these were not subsequently sequenced.

Illumina MiSeq sequencing of 16S rRNA gene amplicons and data processing

The V1-V2 hypervariable regions of 16S rRNA genes were amplified by PCR using a forward and reverse fusion primer, described in detail in Supplementary Methods S1. Bacterial profiling using a MiSeq platform (Illumina, San Diego, CA, USA) was conducted at Research and Testing Laboratory (Lubbock, TX, USA). The 16S rRNA gene sequence data was analysed with bioinformatic software package Mothur [86] using the MiSeq SOP Pipeline. Sequence reads were quality checked and normalised to the lowest number of reads (n=1,855) and singleton operational taxonomic units (OTUs); samples containing fewer than 10 reads were excluded. OTU taxonomies (from Phylum to Genus) were then determined using the RDP MultiClassifier script to generate the RDP taxonomy. Taxonomy level for species of the OTUs was determined using the USEARCH algorithm with 16S rRNA gene sequences from the cultured representatives from the RDP database [87]. Rare OTUs were defined as those present at less than 10 counts within the entire cohort. Alpha and beta diversity indices were calculated from these datasets with Mothur and R statistical package using the Vegan package.

Serum biomarker and local vaginal cytokine analyses

Four serum markers including oestradiol (pmol/L), insulin (mIU/L), glucose (mmol/L) and sex hormone binding globulin (SHBG, nmol/L) were quantified using ELISA at the Imperial College Healthcare NHS Trust North West London Pathology laboratory (Supplementary Methods S2).

Levels of interleukin (IL)-4, IL-6, IL-8, IL-1β, tumour necrosis factor-alpha (TNFα), interferon-gamma (IFN-γ), and macrophage inflammatory protein-1 alpha (MIP1α) in cell-free cervicovaginal secretion supernatants were determined using the Magnetic Luminex Screening Assay multiplex kit (R&D Systems, Minneapolis, MN, USA) on a MAGPIX Analyzer (Luminex® Corporation, s-Hertogenbosch, Netherlands), as per manufacturer’s instructions. Analytes were chosen based on evidence of inflammatory markers specific to adiposity [88-90].

Statistical analyses

The population was categorized into two groups of interest for the main analysis, non-obese (BMI <30kg/m²) versus obese (BMI ³30kg/m²). We performed further supplementary analyses to assess results for different obesity status subcategories and according to insulin resistance and diabetic status. The homeostatic model assessment for insulin resistance (HOMA-IR) was calculated by the following formula: the product of fasting insulin (mU/L) multiplied by fasting glucose (nmol/L) divided by 22.5 [91]. We used the 2^nd tertile value of HOMA-IR as the cut-off to determine insulin resistance status (at 2.98). Differences in categorical clinical parameters between the two main groups of interest (non-obese versus obese) were assessed using Fisher’s exact test for each of the listed characteristics; age, ethnicity, parity, smoking, menopause, menstrual cycle, use of contraception, HRT use, diabetes status and treatment, and abnormal high vaginal swab (HVS) results.

Significant differences between vaginal microbiota at genera taxonomic level were assessed using the Statistical Analysis of Metagenomic Profiles (STAMP) software package [92]. Dependent on genera hierarchical clustering analysis, Lactobacillus spp. or Gardnerella spp. abundance among selected phenotypic categories was investigated by assigning each patient sample into one of three groups (Lactobacillus-dominant, Gardnerella-dominant, or high diversity VMB). Linear discriminant analysis (LDA) effect size (LEfSe) modelling was used to identify biomarkers based on obesity status, according to relative taxonomic abundance [93].

At genus taxonomic level, prevalence of each of the three categories relating to Lactobacillus or Gardnerella presence were compared between the two phenotypic categories (obese or non-obese) using Fisher’s exact test. We performed further sub-analyses for different weight categories, as well as by the presence of diabetes and/or insulin resistance status. A sensitivity analysis assessed whether the exclusion of women that had antibiotics less than 2 weeks before sample collection or those disclosing sexual intercourse within 48 hours from sampling would affect the results. We further analysed the results for pre- and post-menopausal women separately and after exclusion of those taking oral hormonal contraception or hormone replacement therapy (HRT).

Expression levels of assayed cytokines were compared according to obesity status, prevalence of each genus group and diversity (non-parametric Shannon Index) at baseline sampling using two-way ANOVA. Where data fell outside the range of the standard curve for each analyte, either the minimal or maximal extrapolated or minimal or maximal value of the standard curve was used, where appropriate. Analyses were performed using Prism 8, p values <0.05 considered significant.

In the subset of women undergoing bariatric surgery, changes in serum markers introduced by weight loss were analysed. We further assessed the impact of surgically-induced weight loss on the prevalence of each of the three genus groups at baseline, month 3 and 6 using McNemar’s Chi square test. We analysed the results for the full cohort and separately for pre- and post-menopausal women, and according to diabetic and insulin resistance status. Transition in vaginal microbiota across genus groups correlating with weight loss from baseline sampling to 6 months post-surgery was analysed for the total bariatric cohort and pre- and postmenopausal women separately. Cytokine and serum marker expression levels were compared between baseline sampling and 6 months after surgically-induced weight loss.

Ethics approval and consent to participate

Ethical approval was obtained from the National Research Ethics Service Committee London – Fulham (Approval number 13/LO/0126) and the NHS West of Scotland Research Ethics Service Committee (WoSRES) (REC 14/WS/1098). All patients gave informed consent.

Consent for publication

Where individual patient data is included, all patients provided informed consent.

Availability of data and material

Public access to sequence data and accompanying metadata can be obtained at the European Nucleotide Archive’s (ENA) Sequence Read Archive (SRA) (accession number: PRJEB40616).

Competing interests

The authors declare no competing interests.

Funding

This work was supported by Genesis Research Trust (P55549 and Garfield Weston Foundation, Grant number P63522 to MK); Ovarian Cancer Action (Grant number PS5827 and PSA601to MG and MK); British Society of Colposcopy Cervical Pathology Jordan/Singer Award (P47773)(MK); Imperial College Healthcare Charity (P47907) (MK, AM); NIHR Academic Clinical Fellowship programme (AM); the Imperial Experimental Cancer Medicine Centre, the Cancer Research UK Imperial Centre, Imperial Healthcare NHS Trust NIHR BRC (Grant number P45272). None of the funders have had any influence over: study design, collection, analysis and interpretation of the data, in writing the report and in the decisions to submit this article for publication. The Division of Digestive Disease at Imperial College London receives financial support from the National Institute of Health Research (NIHR) Imperial Biomedical Research Centre (BRC) based at Imperial College Healthcare NHS Trust and Imperial College London. This article is independent research funded by the NIHR BRC, and the views expressed in this publication are those of the authors and not necessarily those of the NHS, NIHR, or the Department of Health.

Authors’ contributions

The study was designed by MK, DM, OR and PB. The samples and data were acquired and collated by OR, AM, SP and MK. The data was analysed by OR, DM, AM, AS, NA, YL, PB and MK. The manuscript was drafted by OR, DM, MK, PB and revised critically for important intellectual content by all authors (OR, DM, AM, YL, AS, NA, JN, SP, MG, HG, JM, PB, MK). All authors gave final approval of the version to be published and have contributed to the manuscript.

Acknowledgements

We thank all participants of the study. Our work was supported by Imperial College Healthcare Charity, Genesis Research Trust, Ovarian Cancer Action, the British Society of Colposcopy Cervical Pathology (Jordan/Singer Award), the Imperial Experimental Cancer Medicine Centre, and the Cancer Research UK Imperial Centre. The Division of Digestive Disease at Imperial College London receives financial support from the National Institute of Health Research (NIHR) Imperial Biomedical Research Centre (BRC) based at Imperial College Healthcare NHS Trust and Imperial College London.

Trends in adult body-mass index in 200 countries from 1975 to 2014: a pooled analysis of 1698 population-based measurement studies with 19.2 million participants. Lancet 2016, 387(10026):1377-1396.
Lash MM, Armstrong A: Impact of obesity on women's health. Fertil Steril 2009, 91(5):1712-1716.
Ko KM, Han K, Chung YJ, Yoon KH, Park YG, Lee SH: Association between Body Weight Changes and Menstrual Irregularity: The Korea National Health and Nutrition Examination Survey 2010 to 2012. Endocrinology and metabolism (Seoul, Korea) 2017, 32(2):248-256.
Brewer CJ, Balen AH: The adverse effects of obesity on conception and implantation. Reproduction (Cambridge, England) 2010, 140(3):347-364.
Lim SS, Norman RJ, Davies MJ, Moran LJ: The effect of obesity on polycystic ovary syndrome: a systematic review and meta-analysis. Obes Rev 2013, 14(2):95-109.
Tang T, Glanville J, Hayden CJ, White D, Barth JH, Balen AH: Combined lifestyle modification and metformin in obese patients with polycystic ovary syndrome. A randomized, placebo-controlled, double-blind multicentre study. Hum Reprod 2006, 21(1):80-89.
Rich-Edwards JW, Goldman MB, Willett WC, Hunter DJ, Stampfer MJ, Colditz GA, Manson JE: Adolescent body mass index and infertility caused by ovulatory disorder. Am J Obstet Gynecol 1994, 171(1):171-177.
Maheshwari A, Stofberg L, Bhattacharya S: Effect of overweight and obesity on assisted reproductive technology--a systematic review. Hum Reprod Update 2007, 13(5):433-444.
Balen AH, Anderson RA: Impact of obesity on female reproductive health: British Fertility Society, Policy and Practice Guidelines. Hum Fertil (Camb) 2007, 10(4):195-206.
Lashen H, Fear K, Sturdee DW: Obesity is associated with increased risk of first trimester and recurrent miscarriage: matched case-control study. Hum Reprod 2004, 19(7):1644-1646.
Metwally M, Ong KJ, Ledger WL, Li TC: Does high body mass index increase the risk of miscarriage after spontaneous and assisted conception? A meta-analysis of the evidence. Fertil Steril 2008, 90(3):714-726.
Boots C, Stephenson MD: Does obesity increase the risk of miscarriage in spontaneous conception: a systematic review. Semin Reprod Med 2011, 29(6):507-513.
Marchi J, Berg M, Dencker A, Olander EK, Begley C: Risks associated with obesity in pregnancy, for the mother and baby: a systematic review of reviews. Obes Rev 2015, 16(8):621-638.
Sebire NJ, Jolly M, Harris JP, Wadsworth J, Joffe M, Beard RW, Regan L, Robinson S: Maternal obesity and pregnancy outcome: a study of 287,213 pregnancies in London. Int J Obes Relat Metab Disord 2001, 25(8):1175-1182.
Kyrgiou M, Kalliala I, Markozannes G, Gunter MJ, Paraskevaidis E, Gabra H, Martin-Hirsch P, Tsilidis KK: Adiposity and cancer at major anatomical sites: umbrella review of the literature. BMJ (Clinical research ed) 2017, 356:j477.
Poorolajal J, Jenabi E: The association between BMI and cervical cancer risk: a meta-analysis. Eur J Cancer Prev 2016, 25(3):232-238.
Aune D, Navarro Rosenblatt DA, Chan DS, Abar L, Vingeliene S, Vieira AR, Greenwood DC, Norat T: Anthropometric factors and ovarian cancer risk: a systematic review and nonlinear dose-response meta-analysis of prospective studies. Int J Cancer 2015, 136(8):1888-1898.
Raglan O, Kalliala I, Markozannes G, Cividini S, Gunter MJ, Nautiyal J, Gabra H, Paraskevaidis E, Martin-Hirsch P, Tsilidis KK et al: Risk factors for endometrial cancer: An umbrella review of the literature. Int J Cancer 2019, 145(7):1719-1730.
Marchesi JR, Ravel J: The vocabulary of microbiome research: a proposal. Microbiome 2015, 3:31.
Kostic AD, Chun E, Robertson L, Glickman JN, Gallini CA, Michaud M, Clancy TE, Chung DC, Lochhead P, Hold GL et al: Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe 2013, 14(2):207-215.
Marchesi JR, Dutilh BE, Hall N, Peters WH, Roelofs R, Boleij A, Tjalsma H: Towards the human colorectal cancer microbiome. PLoS One 2011, 6(5):e20447.
Alekseyenko AV, Perez-Perez GI, De Souza A, Strober B, Gao Z, Bihan M, Li K, Methe BA, Blaser MJ: Community differentiation of the cutaneous microbiota in psoriasis. Microbiome 2013, 1(1):31.
Rose WA, 2nd, McGowin CL, Spagnuolo RA, Eaves-Pyles TD, Popov VL, Pyles RB: Commensal bacteria modulate innate immune responses of vaginal epithelial cell multilayer cultures. PLoS One 2012, 7(3):e32728.
Bercik P, Denou E, Collins J, Jackson W, Lu J, Jury J, Deng Y, Blennerhassett P, Macri J, McCoy KD et al: The intestinal microbiota affect central levels of brain-derived neurotropic factor and behavior in mice. Gastroenterology 2011, 141(2):599-609, 609 e591-593.
Ohland CL, Kish L, Bell H, Thiesen A, Hotte N, Pankiv E, Madsen KL: Effects of Lactobacillus helveticus on murine behavior are dependent on diet and genotype and correlate with alterations in the gut microbiome. Psychoneuroendocrinology 2013, 38(9):1738-1747.
McDonald D, Hornig M, Lozupone C, Debelius J, Gilbert JA, Knight R: Towards large-cohort comparative studies to define the factors influencing the gut microbial community structure of ASD patients. Microb Ecol Health Dis 2015, 26:26555.
Romero R, Hassan SS, Gajer P, Tarca AL, Fadrosh DW, Nikita L, Galuppi M, Lamont RF, Chaemsaithong P, Miranda J et al: Correction: The composition and stability of the vaginal microbiota of normal pregnant women is different from that of non-pregnant women. Microbiome 2014, 2(1):10.
Rogosa M, Sharpe ME: Species differentiation of human vaginal lactobacilli. J Gen Microbiol 1960, 23:197-201.
Champer M, Wong AM, Champer J, Brito IL, Messer PW, Hou JY, Wright JD: The role of the vaginal microbiome in gynaecological cancer. BJOG 2017.
Pruski P, Lewis HV, Lee YS, Marchesi JR, Bennett PR, Takats Z, MacIntyre DA: Assessment of microbiota:host interactions at the vaginal mucosa interface. Methods 2018, 149:74-84.
Borgdorff H, van der Veer C, van Houdt R, Alberts CJ, de Vries HJ, Bruisten SM, Snijder MB, Prins M, Geerlings SE, Schim van der Loeff MF et al: The association between ethnicity and vaginal microbiota composition in Amsterdam, the Netherlands. PLoS One 2017, 12(7):e0181135.
Gajer P, Brotman RM, Bai G, Sakamoto J, Schutte UM, Zhong X, Koenig SS, Fu L, Ma ZS, Zhou X et al: Temporal dynamics of the human vaginal microbiota. Science translational medicine 2012, 4(132):132ra152.
Brotman RM, Shardell MD, Gajer P, Fadrosh D, Chang K, Silver MI, Viscidi RP, Burke AE, Ravel J, Gravitt PE: Association between the vaginal microbiota, menopause status, and signs of vulvovaginal atrophy. Menopause 2014, 21(5):450-458.
Shannon B, Yi TJ, Perusini S, Gajer P, Ma B, Humphrys MS, Thomas-Pavanel J, Chieza L, Janakiram P, Saunders M et al: Association of HPV infection and clearance with cervicovaginal immunology and the vaginal microbiota. Mucosal Immunol 2017.
Aagaard K, Riehle K, Ma J, Segata N, Mistretta TA, Coarfa C, Raza S, Rosenbaum S, Van den Veyver I, Milosavljevic A et al: A metagenomic approach to characterization of the vaginal microbiome signature in pregnancy. PLoS One 2012, 7(6):e36466.
MacIntyre DA, Chandiramani M, Lee YS, Kindinger L, Smith A, Angelopoulos N, Lehne B, Arulkumaran S, Brown R, Teoh TG et al: The vaginal microbiome during pregnancy and the postpartum period in a European population. Sci Rep 2015, 5:8988.
Brotman RM, Klebanoff MA, Nansel TR, Andrews WW, Schwebke JR, Zhang J, Yu KF, Zenilman JM, Scharfstein DO: A longitudinal study of vaginal douching and bacterial vaginosis--a marginal structural modeling analysis. Am J Epidemiol 2008, 168(2):188-196.
Murphy K, Mitchell CM: The Interplay of Host Immunity, Environment and the Risk of Bacterial Vaginosis and Associated Reproductive Health Outcomes. J Infect Dis 2016, 214 Suppl 1:S29-35.
Kindinger LM, Bennett PR, Lee YS, Marchesi JR, Smith A, Cacciatore S, Holmes E, Nicholson JK, Teoh TG, MacIntyre DA: The interaction between vaginal microbiota, cervical length, and vaginal progesterone treatment for preterm birth risk. Microbiome 2017, 5(1):6.
Callahan BJ, DiGiulio DB, Goltsman DSA, Sun CL, Costello EK, Jeganathan P, Biggio JR, Wong RJ, Druzin ML, Shaw GM 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 U S A 2017, 114(37):9966-9971.
Fettweis JM, Serrano MG, Brooks JP, Edwards DJ, Girerd PH, Parikh HI, Huang B, Arodz TJ, Edupuganti L, Glascock AL et al: The vaginal microbiome and preterm birth. Nat Med 2019, 25(6):1012-1021.
Elovitz MA, Gajer P, Riis V, Brown AG, Humphrys MS, Holm JB, Ravel J: Cervicovaginal microbiota and local immune response modulate the risk of spontaneous preterm delivery. Nat Commun 2019, 10(1):1305.
Brown RG, Marchesi JR, Lee YS, Smith A, Lehne B, Kindinger LM, Terzidou V, Holmes E, Nicholson JK, Bennett PR et al: Vaginal dysbiosis increases risk of preterm fetal membrane rupture, neonatal sepsis and is exacerbated by erythromycin. BMC Med 2018, 16(1):9.
Martin HL, Richardson BA, Nyange PM, Lavreys L, Hillier SL, Chohan B, Mandaliya K, Ndinya-Achola JO, Bwayo J, Kreiss J: Vaginal lactobacilli, microbial flora, and risk of human immunodeficiency virus type 1 and sexually transmitted disease acquisition. J Infect Dis 1999, 180(6):1863-1868.
Mitra A, MacIntyre DA, Marchesi JR, Lee YS, Bennett PR, Kyrgiou M: The vaginal microbiota, human papillomavirus infection and cervical intraepithelial neoplasia: what do we know and where are we going next? Microbiome 2016, 4(1):58.
Kyrgiou M, Mitra A, Moscicki AB: Does the vaginal microbiota play a role in the development of cervical cancer? Transl Res 2017, 179:168-182.
Mitra A, MacIntyre DA, Ntritsos G, Smith A, Tsilidis KK, Marchesi JR, Bennett PR, Moscicki A-B, Kyrgiou M: The vaginal microbiota associates with the regression of untreated cervical intraepithelial neoplasia 2 lesions. Nature Communications 2020, 11(1):1999.
Brookheart RT, Lewis WG, Peipert JF, Lewis AL, Allsworth JE: Association between obesity and bacterial vaginosis as assessed by Nugent score. Am J Obstet Gynecol 2019.
Koumans EH, Sternberg M, Bruce C, McQuillan G, Kendrick J, Sutton M, Markowitz LE: The prevalence of bacterial vaginosis in the United States, 2001-2004; associations with symptoms, sexual behaviors, and reproductive health. Sex Transm Dis 2007, 34(11):864-869.
Oh HY, Seo SS, Kong JS, Lee JK, Kim MK: Association between Obesity and Cervical Microflora Dominated by Lactobacillus iners in Korean Women. Journal of clinical microbiology 2015, 53(10):3304-3309.
Si J, You HJ, Yu J, Sung J, Ko G: Prevotella as a Hub for Vaginal Microbiota under the Influence of Host Genetics and Their Association with Obesity. Cell host & microbe 2017, 21(1):97-105.
Statistics on Obesity, Physical Activity and Diet - England, 2018. NHS Digital [https://digital.nhs.uk/data-and-information/publications/statistical/statistics-on-obesity-physical-activity-and-diet/statistics-on-obesity-physical-activity-and-diet-england-2018#key-facts]
Brotman RM, Shardell MD, Gajer P, Fadrosh D, Chang K, Silver MI, Viscidi RP, Burke AE, Ravel J, Gravitt PE: Association between the vaginal microbiota, menopause status, and signs of vulvovaginal atrophy. Menopause 2018, 25(11):1321-1330.
Spear GT, French AL, Gilbert D, Zariffard MR, Mirmonsef P, Sullivan TH, Spear WW, Landay A, Micci S, Lee BH et al: Human alpha-amylase present in lower-genital-tract mucosal fluid processes glycogen to support vaginal colonization by Lactobacillus. J Infect Dis 2014, 210(7):1019-1028.
Farage M, Maibach H: Lifetime changes in the vulva and vagina. Arch Gynecol Obstet 2006, 273(4):195-202.
Hammerschlag MR, Alpert S, Rosner I, Thurston P, Semine D, McComb D, McCormack WM: Microbiology of the vagina in children: normal and potentially pathogenic organisms. Pediatrics 1978, 62(1):57-62.
Hammerschlag MR, Alpert S, Onderdonk AB, Thurston P, Drude E, McCormack WM, Bartlett JG: Anaerobic microflora of the vagina in children. Am J Obstet Gynecol 1978, 131(8):853-856.
Kaaks R, Lukanova A, Kurzer MS: Obesity, endogenous hormones, and endometrial cancer risk: a synthetic review. Cancer Epidemiol Biomarkers Prev 2002, 11(12):1531-1543.
Key TJ, Allen NE, Verkasalo PK, Banks E: Energy balance and cancer: the role of sex hormones. Proc Nutr Soc 2001, 60(1):81-89.
Zhao H, Zhou L, Shangguan AJ, Bulun SE: Aromatase expression and regulation in breast and endometrial cancer. J Mol Endocrinol 2016, 57(1):R19-33.
Muhleisen AL, Herbst-Kralovetz MM: Menopause and the vaginal microbiome. Maturitas 2016, 91:42-50.
Gliniewicz K, Schneider GM, Ridenhour BJ, Williams CJ, Song Y, Farage MA, Miller K, Forney LJ: Comparison of the Vaginal Microbiomes of Premenopausal and Postmenopausal Women. Front Microbiol 2019, 10:193.
Freeman EW, Sammel MD, Lin H, Gracia CR: Obesity and reproductive hormone levels in the transition to menopause. Menopause 2010, 17(4):718-726.
Noyes N, Cho KC, Ravel J, Forney LJ, Abdo Z: Associations between sexual habits, menstrual hygiene practices, demographics and the vaginal microbiome as revealed by Bayesian network analysis. PLoS One 2018, 13(1):e0191625.
Hassan S, Chatwani A, Brovender H, Zane R, Valaoras T, Sobel JD: Douching for perceived vaginal odor with no infectious cause of vaginitis: a randomized controlled trial. J Low Genit Tract Dis 2011, 15(2):128-133.
MacKintosh ML, Derbyshire AE, McVey RJ, Bolton J, Nickkho-Amiry M, Higgins CL, Kamieniorz M, Pemberton PW, Kirmani BH, Ahmed B et al: The impact of obesity and bariatric surgery on circulating and tissue biomarkers of endometrial cancer risk. Int J Cancer 2019, 144(3):641-650.
Legro RS, Dodson WC, Gnatuk CL, Estes SJ, Kunselman AR, Meadows JW, Kesner JS, Krieg EF, Jr., Rogers AM, Haluck RS et al: Effects of gastric bypass surgery on female reproductive function. J Clin Endocrinol Metab 2012, 97(12):4540-4548.
Jamal M, Gunay Y, Capper A, Eid A, Heitshusen D, Samuel I: Roux-en-Y gastric bypass ameliorates polycystic ovary syndrome and dramatically improves conception rates: a 9-year analysis. Surg Obes Relat Dis 2012, 8(4):440-444.
Maggard MA, Yermilov I, Li Z, Maglione M, Newberry S, Suttorp M, Hilton L, Santry HP, Morton JM, Livingston EH et al: Pregnancy and fertility following bariatric surgery: a systematic review. JAMA 2008, 300(19):2286-2296.
Tan O, Carr BR: The impact of bariatric surgery on obesity-related infertility and in vitro fertilization outcomes. Semin Reprod Med 2012, 30(6):517-528.
Moran LJ, Norman RJ: The effect of bariatric surgery on female reproductive function. J Clin Endocrinol Metab 2012, 97(12):4352-4354.
Merhi ZO: Impact of bariatric surgery on female reproduction. Fertil Steril 2009, 92(5):1501-1508.
Guelinckx I, Devlieger R, Vansant G: Reproductive outcome after bariatric surgery: a critical review. Hum Reprod Update 2009, 15(2):189-201.
Uzoma A, Keriakos R: Pregnancy management following bariatric surgery. J Obstet Gynaecol 2013, 33(2):109-114.
Liu R, Hong J, Xu X, Feng Q, Zhang D, Gu Y, Shi J, Zhao S, Liu W, Wang X et al: Gut microbiome and serum metabolome alterations in obesity and after weight-loss intervention. Nat Med 2017, 23(7):859-868.
Li JV, Ashrafian H, Bueter M, Kinross J, Sands C, le Roux CW, Bloom SR, Darzi A, Athanasiou T, Marchesi JR et al: Metabolic surgery profoundly influences gut microbial-host metabolic cross-talk. Gut 2011, 60(9):1214-1223.
Li JV, Reshat R, Wu Q, Ashrafian H, Bueter M, le Roux CW, Darzi A, Athanasiou T, Marchesi JR, Nicholson JK et al: Experimental bariatric surgery in rats generates a cytotoxic chemical environment in the gut contents. Front Microbiol 2011, 2:183.
Ouchi N, Parker JL, Lugus JJ, Walsh K: Adipokines in inflammation and metabolic disease. Nat Rev Immunol 2011, 11(2):85-97.
Hotamisligil GS: Inflammation, metaflammation and immunometabolic disorders. Nature 2017, 542(7640):177-185.
Christiaens I, Zaragoza DB, Guilbert L, Robertson SA, Mitchell BF, Olson DM: Inflammatory processes in preterm and term parturition. J Reprod Immunol 2008, 79(1):50-57.
Thurman AR, Kimble T, Herold B, Mesquita PM, Fichorova RN, Dawood HY, Fashemi T, Chandra N, Rabe L, Cunningham TD et al: Bacterial Vaginosis and Subclinical Markers of Genital Tract Inflammation and Mucosal Immunity. AIDS Res Hum Retroviruses 2015, 31(11):1139-1152.
Kindinger LM, MacIntyre DA, Lee YS, Marchesi JR, Smith A, McDonald JA, Terzidou V, Cook JR, Lees C, Israfil-Bayli F et al: Relationship between vaginal microbial dysbiosis, inflammation, and pregnancy outcomes in cervical cerclage. Sci Transl Med 2016, 8(350):350ra102.
Maciejewski ML, Arterburn DE, Van Scoyoc L, Smith VA, Yancy WS, Jr., Weidenbacher HJ, Livingston EH, Olsen MK: Bariatric Surgery and Long-term Durability of Weight Loss. JAMA surgery 2016, 151(11):1046-1055.
Sjostrom L, Narbro K, Sjostrom CD, Karason K, Larsson B, Wedel H, Lystig T, Sullivan M, Bouchard C, Carlsson B et al: Effects of bariatric surgery on mortality in Swedish obese subjects. N Engl J Med 2007, 357(8):741-752.
Brown RG, Al-Memar M, Marchesi JR, Lee YS, Smith A, Chan D, Lewis H, Kindinger L, Terzidou V, Bourne T et al: Establishment of vaginal microbiota composition in early pregnancy and its association with subsequent preterm prelabor rupture of the fetal membranes. Transl Res 2019, 207:30-43.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD: Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol 2013, 79(17):5112-5120.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26(19):2460-2461.
Himbert C, Delphan M, Scherer D, Bowers LW, Hursting S, Ulrich CM: Signals from the Adipose Microenvironment and the Obesity-Cancer Link-A Systematic Review. Cancer Prev Res (Phila) 2017, 10(9):494-506.
Gilbert CA, Slingerland JM: Cytokines, obesity, and cancer: new insights on mechanisms linking obesity to cancer risk and progression. Annu Rev Med 2013, 64:45-57.
Schmidt FM, Weschenfelder J, Sander C, Minkwitz J, Thormann J, Chittka T, Mergl R, Kirkby KC, Fasshauer M, Stumvoll M et al: Inflammatory cytokines in general and central obesity and modulating effects of physical activity. PLoS One 2015, 10(3):e0121971.
Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC: Homeostasis model assessment: insulin resistance and beta-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 1985, 28(7):412-419.
Parks DH, Beiko RG: Identifying biologically relevant differences between metagenomic communities. Bioinformatics (Oxford, England) 2010, 26(6):715-721.
Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C: Metagenomic biomarker discovery and explanation. Genome Biol 2011, 12(6):R60.

Table 1. Patient characteristics of the whole cohort at baseline sample collection (n=109) according to obesity status.

Characteristics	Non-obese (n=67)	Obese (n=42)	Total (n=109)	P value^a
Age (years)				0.083
Mean (SD, range)	44 (11.79, 20-75)	46 (11.26, 28-72)	44 (11.79, 20-75)
Ethnicity, n/N (%)				0.024
Caucasian	48/67 (71.6)	22/42 (52.4)	70/109 (64.2)
Asian	3/67 (4.5)	5/42 (11.9)	8/109 (7.3)
Black	12/67 (17.9)	15/42 (35.7)	27/109 (24.8)
Other	4/67 (6.0)	0/42 (0)	4/109 (3.7)
Parity, n/N (%)				<0.001
Nulliparous	44/67 (65.7)	7/42 (16.7)	51/109 (46.8)
Parous	23/67 (34.3)	35/42 (83.3)	58/109 (53.2)
Smoking status, n/N (%)				0.128
Current smoker	11/67 (16.4)	2/42 (4.8)	13/109 (11.9)
Non-smoker	56/67 (83.6)	40/42 (95.2)	96/109 (88.1)
HVS results n/N (%)				0.062
Normal	46/67 (68.7)	36/42 (85.7)	82/109 (75.2)
Abnormal	15/67 (22.4)	6/42 (14.3)	21/109 (19.3)
Unknown	6/67 (8.9)	0/42 (0)	6/109 (5.5)
Abnormal HVS results, n/N (%)				0.076
Bacterial vaginosis	3/15 (20.0)	1/6 (16.7)	4/21 (19.0)
E.coli	1/15 (6.7)	0/6 (0)	1/21 (4.8)
S.aureus	1/15 (6.7)	0/6 (0)	1/21 (4.8)
Group B streptococcus (S.agalactiae)	2/15 (13.3)	3/6 (50.0)	5/21 (23.8)
Yeast	0/15 (0)	1/6 (16.7)	1/21 (4.8)
Mixed coliforms	0/15 (0)	1/6 (16.6)	1/21 (4.8)
Unknown	8/15 (53.3)	0/6 (0)	8/21 (38.0)
Menopause status, n/N (%)				0.233
Premenopausal	50/67 (74.6)	26/42 (61.9)	76/109 (70.0)
Postmenopausal	17/67 (25.4)	16/42 (38.1)	33/109 (30.0)
Phase of menstrual cycle (PrMP), n/N (%)				0.115
Luteal	24/50 (48.0)	7/26 (26.9)	31/76 (40.8)
Follicular	15/50 (30.0)	8/26 (30.8)	23/76 (30.3)
Ovulation	0/50 (0)	1/26 (3.8)	1/76 (1.3)
Unknown	11/50 (22.0)	10/26 (38.5)	21/76 (27.6)
Use of contraception (PrMP), n/N (%)				0.016
Nil	27/50 (54.0)	20/26 (77.0)	47/76 (61.9)
Condoms	5/50 (10.0)	0/26 (0)	5/76 (6.6)
COCP	12/50 (24.0)	1/26 (3.8)	13/76 (17.1)
POP	1/50 (2.0)	0/26 (0)	1/76 (1.3)
Copper IUD	0/50 (0)	0/26 (0)	0/76 (0)
Mirena IUS	3/50 (6.0)	5/26 (19.2)	8/76 (10.5)
Vaginal ring	0/50 (0)	0/26 (0)	0/76 (0)
Contraceptive implant	2/50 (4.0)	0/26 (0)	2/76 (2.6)
Contraceptive injection	0/50 (0)	0/26 (0)	0/76 (0)
Use of HRT (PoMP), n/N (%)				0.175
Yes	5/17 (29.4)	1/16 (6.2)	6/33 (18.2)
No	12/17 (70.6)	15/16 (93.8)	27/33 (81.8)
Diabetes status, n/N (%)				0.001
Non-diabetic	65/67 (97.0)	32/42 (76.2)	97/109 (89.0)
Diabetic	2/67 (3.0)	10/42 (23.8)	12/109 (11.0)
Diabetic treatment, n/N (%)				0.212
Diet control only	1/2 (50.0)	3/10 (30.0)	4/12 (33.3)
Metformin alone	0/2 (0)	3/10 (30.0)	3/12 (25.0)
Metformin combined 2nd diabetic medication	0/2 (0)	0/10 (0)	0/12 (0)
Other oral antiglycaemic medication	1/2 (50.0)	0/10 (0)	1/12 (8.4)
Insulin (alone or with oral medication)	0/2 (0)	4/10 (40.0)	4/12 (33.3)
HOMA-IR^b, n/N (%)				0.136
Insulin resistant	2/67 (3.0)	5/42 (11.9)	7/109 (6.4)
Non-insulin resistant	5/67 (7.5)	5/42 (11.9)	10/109 (9.2)
Unknown insulin resistance status	60/67 (89.5)	32/42 (76.2)	92/109 (84.4)

^aCalculated using Fisher’s exact test. A p-value of less than 0.05 demonstrates a significant difference in the distribution of the demographic of interest (e.g. ethnicity), according to obesity status.

^b HOMA-IR was calculated according to the formula: fasting insulin (mU/L) multiplied by fasting glucose (nmol/L)/22.5. The 2nd tertile was used as the cut-off to determine insulin resistance status. Insulin resistance cut-off value: 2.98

BMI, body mass index; COCP, combined oral contraceptive pill; E.coli, Escherichia coli; HOMA-IR, homeostatic model of assessment - insulin resistance; HVS, high vaginal swab; IUD, intrauterine device; IUS, intrauterine system; PoMP, postmenopausal; POP, progesterone-only pill; PrMP, premenopausal; S.aureus, staphylococcus aureus; SD, standard deviation.

Table 2. Prevalence of genus group present according to obesity status, diabetes and insulin resistance status for the whole cohort (n=109) at baseline sampling.

	*Lactobacillus* dominant n/N (%)	*Gardnerella* dominant n/N (%)	High diversity n/N (%)	Total n/N (%)
Obesity status
Non-obese (BMI <30.0kg/m²)	41/67 (61.2)	16/67 (23.9)	10/67 (14.9)	67/67 (100)
Obese (≥30.0kg/m²)	12/42 (28.6)	12/42 (28.6)	18/42 (42.8)	42/42 (100)
Total	53/109 (48.6)	28/109 (25.7)	28/109 (25.7)	109/109 (100)
P value^a				0.002
Diabetes status
Non-diabetic	49/97 (50.5)	25/97 (25.8)	23/97 (23.7)	97/97 (100)
Diabetic	4/12 (33.3)	3/12 (25.0)	5/12 (41.7)	12/12 (100)
Total	53/109 (48.6)	28/109 (25.7)	28/109 (25.7)	109/109 (100)
P value				0.451
Insulin resistance status^b
Non-insulin resistant	4/10 (40.0)	2/10 (20.0)	4/10 (40.0)	10/10 (100)
Insulin resistant	2/7 (28.6)	2/7 (28.6)	3/7 (42.9)	7/7 (100)
Total	6/17 (35.3)	4/17 (23.5)	7/17 (41.2)	17/17 (100)
P value				1.000

^aA p-value of less than 0.05 demonstrates a significant difference in the distribution of genus group present according to obesity, diabetic and insulin resistance status. Fisher’s exact test employed as small numbers were present for each group (e.g. diabetic).

^b Where concomitant fasting serum samples were available, fasting glucose and fasting insulin levels were identified. Using these values, we were able to calculate the HOMA-IR, according to the formula: the product of fasting insulin (mU/L) multiplied by fasting glucose (nmol/L) divided by 22.5. The 2nd tertile was used as the cut-off to determine insulin resistance status. Insulin resistance cut-off value: 2.98.

BMI, body mass index

Table 3. Prevalence of genus group present over serial sample collection timepoints in the total bariatric cohort and according to menopausal status, diabetes status and insulin resistance status.

	*Lactobacillus* dominant n/N (%)	*Gardnerella* dominant n/N (%)	High diversity n/N (%)	TOTAL n/N (%)
Total bariatric surgery cohort
Baseline (day of surgery) (n=27)	10/27 (37.0)	5/27 (18.6)	12/27 (44.4)	27/27 (100)
3 months post-surgery (n=22)	7/22 (31.8)	6/22 (27.3)	9/22 (40.9)	22/22 (100)
6 months post-surgery (n=19)	10/19 (52.6)	2/19 (10.5)	7/19 (36.9)	19/19 (100)
Total (n=84)	31/84 (36.9)	16/84 (19.1)	37/84 (44.0)	84/84 (100)
p value^a				0.506
Pre-menopausal women only
Baseline (day of surgery) (n=17)	6/17 (35.3)	4/17 (23.5)	7/17 (41.2)	17/17 (100)
3 months post-surgery (n=15)	5/15 (33.3)	4/15 (26.7)	6/15 (40.0)	15/15 (100)
6 months post-surgery (n=14)	8/14 (57.1)	2/14 (14.3)	4/14 (28.6)	14/14 (100)
Total (n=58)	22/58 (37.9)	12/58 (20.7)	24/58 (41.4)	58/58 (100)
p value				0.572
Post-menopausal women only
Baseline (day of surgery) (n=10)	4/10 (40.0)	1/10 (10.0)	5/10 (50.0)	10/10 (100)
3 months post-surgery (n=7)	2/7 (28.6)	2/7 (28.6)	3/7 (42.8)	7/7 (100)
6 months post-surgery (n=5)	2/5 (40.0)	0/5 (0)	3/5 (60.0)	5/5 (100)
Total (n=26)	9/26 (34.6)	4/26 (15.4)	13/26 (50.0)	26/26 (100)
p value				0.999
Non-diabetic women only
Baseline (day of surgery) (n=19)	6/19 (31.6)	3/19 (15.8)	10/19 (52.6)	19/19 (100)
3 months post-surgery (n=16)	6/16 (37.5)	4/16 (25.0)	6/16 (37.5)	16/16 (100)
6 months post-surgery (n=14)	8/14 (57.1)	1/14 (7.2)	5/14 (35.7)	14/14 (100)
Total (n= 58)	21/58 (36.2)	9/58 (15.5)	28/58 (48.3)	58/58 (100)
p value				0.572
Diabetic women only
Baseline (day of surgery) (n=8)	4/8 (50.0)	2/8 (25.0)	2/8 (25.0)	8/8 (100)
3 months post-surgery (n=6)	1/6 (16.7)	2/6 (33.3)	3/6 (50.0)	6/6 (100)
6 months post-surgery (n=5)	2/5 (40.0)	1/5 (20.0)	2/5 (40.0)	5/5 (100)
Total (n=26)	10/26 (38.5)	7/26 (26.9)	9/26 (34.6)	26/26 (100)
p value**				n/a
Women without insulin resistance^c
Baseline (day of surgery) (n=4)	1/4 (25.0)	1/4 (25.0)	2/4 (50.0)	4/4 (100)
3 months post-surgery (n=2)	0/2 (0)	0/2 (0)	2/2 (100)	2/2 (100)
6 months post-surgery (n=7)	4/7 (57.1)	1/7 (14.3)	2/7 (28.6)	7/7 (100)
Total (n=13)	5/13 (38.5)	2/13 (15.4)	6/13 (46.1)	13/13 (100)
p value				n/a
Women with insulin resistance
Baseline (day of surgery) (n=5)	2/5 (40.0)	1/5 (20.0)	2/5 (40.0)	5/5 (100)
3 months post-surgery (n=1)	0/1 (0)	0/1 (0)	1/1 (100)	1/1 (100)
6 months post-surgery (n=1)	1/1 (100)	0/1 (0)	0/1 (0)	1/1 (100)
Total (n=7)	3/7 (42.9)	1/7 (14.3)	3/7 (42.9)	7/7 (100)
p value^b				n/a

^a A p-value of less than 0.05 demonstrates a significant difference between the proportion (%) of genus group present (Lactobacillus-dominant, Gardnerella-dominant or high diversity VMB) over serial timepoints 0 to 6 months (McNemar’s Chi-square test).

^b Sample size too small to compute p-value using McNemar’s Chi-square test.

^c Where concomitant fasting serum samples were available, fasting glucose and fasting insulin levels were identified. Using these values, we were able to calculate the HOMA-IR, according to the formula: the product of fasting insulin (mU/L) multiplied by fasting glucose (nmol/L) divided by 22.5. The 2nd tertile was used as the cut-off to determine insulin resistance status. Insulin resistance cut-off value: 2.98.

SupplementaryMaterialsclean.docx

The association between obesity and weight loss after bariatric surgery on the vaginal microbiota

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Results

Discussion

Conclusions

Materials And Methods

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 2