- Open Access
Potential contribution of the uterine microbiome in the development of endometrial cancer
Genome Medicinevolume 8, Article number: 122 (2016)
Endometrial cancer studies have led to a number of well-defined but mechanistically unconnected genetic and environmental risk factors. One of the emerging modulators between environmental triggers and genetic expression is the microbiome. We set out to inquire about the composition of the uterine microbiome and its putative role in endometrial cancer.
We undertook a study of the microbiome in samples taken from different locations along the female reproductive tract in patients with endometrial cancer (n = 17), patients with endometrial hyperplasia (endometrial cancer precursor, n = 4), and patients afflicted with benign uterine conditions (n = 10). Vaginal, cervical, Fallopian, ovarian, peritoneal, and urine samples were collected aseptically both in the operating room and the pathology laboratory. DNA extraction was followed by amplification and high-throughput next generation sequencing (MiSeq) of the 16S rDNA V3-V5 region to identify the microbiota present. Microbiota data were summarized using both α-diversity to reflect species richness and evenness within bacterial populations and β-diversity to reflect the shared diversity between bacterial populations. Statistical significance was determined through the use of multiple testing, including the generalized mixed-effects model.
The microbiome sequencing (16S rDNA V3-V5 region) revealed that the microbiomes of all organs (vagina, cervix, Fallopian tubes, and ovaries) are significantly correlated (p < 0.001) and that there is a structural microbiome shift in the cancer and hyperplasia cases, distinguishable from the benign cases (p = 0.01). Several taxa were found to be significantly enriched in samples belonging to the endometrial cancer cohort: Firmicutes (Anaerostipes, ph2, Dialister, Peptoniphilus, 1–68, Ruminococcus, and Anaerotruncus), Spirochaetes (Treponema), Actinobacteria (Atopobium), Bacteroidetes (Bacteroides and Porphyromonas), and Proteobacteria (Arthrospira). Of particular relevance, the simultaneous presence of Atopobium vaginae and an uncultured representative of the Porphyromonas sp. (99 % match to P. somerae) were found to be associated with disease status, especially if combined with a high vaginal pH (>4.5).
Our results suggest that the detection of A. vaginae and the identified Porphyromonas sp. in the gynecologic tract combined with a high vaginal pH is statistically associated with the presence of endometrial cancer. Given the documented association of the identified microorganisms with other pathologies, these findings raise the possibility of a microbiome role in the manifestation, etiology, or progression of endometrial cancer that should be further investigated.
The causative or triggering agents for endometrial cancer remain elusive despite continued research along the PI3K/PTEN/mTOR/HIF axis in type I  and the p53 tumor-suppressor system in type II endometrial cancer . Host genetics explain only 20 % of endometrial cancer incidence through microsatellite instability (MSI)  or abnormalities in aerobic glycolysis . The efforts to identify the cause of the remaining 80 % of cases have led to studies of a number of environmental and host factors including hormones , obesity , and diabetes . However, these alone do not address the question of tumorigenic mechanism. There is a need to examine potential causative agents, studies of which bring the promise of developing targeted prevention strategies.
Here, we explore a major source of environmental influence on the uterine microenvironment—the microbiome. Microbial influence on the etiology and progression of cancer has already been well established for Helicobacter pylori and gastric cancer . Recent high-throughput sequencing assays have revealed associations between colorectal cancer and infection with Fusobacteria  and Porphyromonas  that are suggestive of a broader microbiome role in cancerous processes. Like the two examples above, endometrial cancer also often arises from a pro-inflammatory profile . We sought to explore the potential microbial triggers for inflammation and tumorigenesis through examination of the uterine microbiome in participants with endometrial cancer.
The microbial partners along the female reproductive tract have been long known to play an important role in health and disease along the woman’s reproductive tract. Lactic acid producing microbes have a strong role in determining the microbial community membership of the vaginal microbiome and have been shown to protect against infection . Gynecologic pathogens associated with bacterial vaginosis, such as Atopobium vaginae and Gardnerella vaginalis have been associated with obstetric complications, such as preterm labor . However, few studies have directly probed the microbes within the uterine environment and how these microbes could influence cancer within the endometrial lining. Given the inflammatory profile in endometrial cancer manifestation, we hypothesized that there is a microbiome component in the malignancy and that its signature in patients diagnosed with the disease is distinguishable from that of patients without malignancy.
We report the results from 31 participants enrolled at the Gynecologic Division, Mayo Clinic, Rochester, MN under an IRB approval protocol (12–004445). The inclusion criteria were the following: 18 years of age or older; women undergoing hysterectomy by any standard surgical approach; undergoing hysterectomy for benign disease, hyperplasia, or any stage of endometrial cancer. Patients with any of the following criteria were excluded from our study: women who were pregnant or nursing; had taken antibiotics within two weeks preceding surgery; surgeon using morcellation during the hysterectomy procedure, due to the size of the uterus or for any other reason. Upon enrollment the participants were requested to fill out an optional questionnaire about sexual and reproductive health and history. The metadata from the questionnaires was stored at REDCap . Cancer participants were also requested to provide a stool sample for the search for putative endometrial cancer signatures.
Vaginal and cervical samples
All participants were requested not to douche with betadine on the day of surgery or the day immediately preceding it. All the vaginal and cervical swabs and scrapes were collected by the surgeon (with guidance on site by the research team) immediately after the administration of anesthesia and immediately preceding the standard pre-surgical betadine douche. Both the vaginal and cervical swabs were performed with three sterile Dacron swabs each and placed in a sterile tube with 1 mL of Tris-EDTA (TE) buffer kept on dry ice until storage at –80 °C. One of the vaginal swabs was used for immediate on-site vaginal pH measurement with a Hydrion measuring pH tape. The scrapes were performed using sterilized (autoclaved at 121 °C for 20 min) pap smear spatulas and placed in sterile tubes with TE buffer kept in dry ice until storage at –80 °C.
Uterine, Fallopian, and ovarian samples
Once removed, the uterus, Fallopian tubes, and ovaries were handed by the surgeon to the instrumentalist nurse who placed them inside a sterile transport bag and into a closed sterile container. The research team then transported the container to the pathology lab (within the same clean area) where the organs were handed to a pathologists’ assistant (PA) to be processed under sterile conditions. The grossing station where the specimen was processed was sterilized by the research team, including all the tools needed by the PA for handling. The PA used surgical gloves and mask when handling the specimen. The PA performed a bilateral cut of the uterus and splayed it. The research team advanced to the collection of the uterine swabs (Dacron) and scrapes (sterilized pap smear spatulas) and documentation (by placement of push pins in sampled locations and digital photograph). The PA then proceeded to the aseptic collection of samples needed for the diagnosis and, once complete, the research team collected the uterine, Fallopian, and ovarian biopsies (approximately 4 mm of tissue was collected per biopsy by the use of a pair of sterile tweezers, scalpel, and surgical ruler). Each collected sample was placed in a sterile tube with 1 mL of TE buffer and kept on dry ice until storage at –80 °C. A petri dish with Lysogeny broth (LB) was kept open on the grossing station during sample collection to detect any possible airborne contamination of the specimen. The LB was swabbed and the swab was stored in a tube with 1 mL of TE and kept on dry ice until storage along with all the other samples.
Once thawed, the swab and scrape samples were vortexed to bring the collected material into solution. The biopsy samples were macerated by the use of sterile pestles. The swab and scrape samples were centrifuged for 10 min at 10,000 g to collect the bacterial cells and the supernatant was discarded. All genomic DNA extractions were performed by using the MoBio PowerSoil Kit (MoBio Laboratories, Inc., Carlsbad, CA, USA) as described by the manufacturer; however, instead of vortexing, an MP FastPrep (MP Biomedicals, Solon, OH, USA) was used instead, for 60 s at 6.0 m/s, to obtain a more effective and rapid lysis of the cells. After extraction the DNA content was measured using High Sensitivity Qubit (Life Technologies Corporation, Carlsbad, CA, USA). The V3-V5 region of the 16S rDNA was then amplified through a polymerase chain reaction (PCR) as follows: 25 μL of Kapa HiFi (Kapa Biosystems, Woburn, MA, USA), 1.5 μL (10 uM) forward primer, 1.5 μL (10 uM) reverse primer, 50 ng of DNA with the remaining volume being added by molecular grade water (up to a final volume of 50 μL per reaction). The forward primer was the universal primer 357 F (5’GTCCTACGGGAGGCAGCAG3’) with the added construct on the 5’ end of the 5’ Illumina Adapter (5’AATGATACGGCGACCACCGAGATCTACAC3’) + Forward Primer Pad (5’TATGGTAATT3’) to a total sequence: 5’AATGATACGGCGACCACCGAGATCTACACTATGGTAATTGTCCTACGGGAGGCAGCAG3’ and the universal bacterial reverse primer was 926R (5’CCGTCAATTCMTTTRAGT3’) with an added construct on the 5’ end of the reverse complement of 3’ Illumina adapter (5’CAAGCAGAAGACGGCATACGAGATGCCGCATTCGAT3’) + Barcode (12 base pairs) to a total sequence: 5’CAAGCAGAAGACGGCATACGAGATGCCGCATTCGATXXXXXXXXXXXXCCGTCAATTCMTTTRAGT3’. The barcode introduced in the reverse primer construct was unique to each sample, functioning as a genetic ID for sequencing. The PCR cycle was the following: 95 °C for 3 min, 98 °C for 20 s, 70 °C for 15 s, 72 °C for 15 s, cycle repeated 34 times, and 72 °C for 5 min. The products of the amplification were verified by a TapeStation D1K Tape (2200 TapeStation Instrument, Agilent Technologies, Santa Clara, CA, USA) to be free of contamination and to contain the expected amplification size, approximately 700 base pairs. If the amplification was unsuccessful, the parameters of the reaction or cycle were adjusted in repeated attempts. In some cases (mostly biopsy samples) the amplification was not successful even after repeated attempts. The reduced number of microorganisms present in the upper reproductive tract is likely to justify this result and attests for the success of the sterile collection of the samples. In samples that failed 16S rDNA amplification, NEBNext Microbiome DNA Enrichment Kit (New England Biolabs Inc., Ipswitch, MA, USA) was used to separate the microbiome from the human DNA to increase the odds of a successful amplification from samples naturally enriched with human DNA (mostly tissue samples). Controls of both the DNA extraction and Microbiome Enrichment processes were performed and are shown in Supplement 5. Upon verification the PCR products were purified using Agencourt AMPure (Beckman Coulter, Brea, CA, USA). After purification the concentrations were measured using Qubit High Sensitivity. The 16S rDNA sequencing was performed by the MGF (Medical Genome facility at Mayo Clinic, Rochester) using a high-throughput next-generation Illumina MiSeq (San Diego, CA, USA) sequencing platform.
Sequence reads were aligned with our own custom multiple alignment tool known as the Illinois-Mayo Taxon Operations for RNA Dataset Organization (IM-TORNADO) that merges paired end reads into a single multiple alignment and obtains taxa calls . IM-TORNADO then clusters sequences into operational taxonomic units (OTUs) using AbundantOTU+ .
A total of 16,366,472 sequence reads (17,657–828,181 reads per sample) were obtained (mean of 199,591 ± 190,153 reads) after quality control. Further processing for visualization was performed using QIIME  and METAGENassist .
α-diversity and β-diversity analysis
To compare the microbiota composition between cohorts, we summarized the data using both α-diversity and β-diversity. α-diversity reflects species richness and evenness within bacterial populations. Two α-diversity metrics, the observed OTU number and the Shannon index, were investigated. Rarefaction curves were used to compare the α-diversity measures. The observed OTU number reflects species richness, whereas the Shannon index measures both species richness and evenness. β-diversity reflects the shared diversity between bacterial communities in terms of ecological distance between samples; different distance metrics provide distinctive views of community structure. Two β-diversity measures (unweighted and weighted UniFrac distances) were calculated using the OTU table and a phylogenetic tree (“GUniFrac” function in the R package GUniFrac) . The unweighted UniFrac reflects differences in community membership (i.e. the presence or absence of an OTU), whereas the weighted UniFrac captures this information and also differences in abundance. Rarefaction was performed on the OTU table before calculating the distances.
To assess the association with α-diversity, we fitted a linear mixed effects model (LME) to the α-diversity metrics with a random intercept for each subject (“lme” function in R package “nlme”), adjusting for covariates if necessary. Wald test was used to assess the significance. To assess the association with β-diversity measures, we used a variant of PERMANOVA procedure (“adonis” function in the R “vegan” package), which is a multivariate analysis of variance based on distance matrices and permutation . To retain the within-subject correlation, we used a block-permutation scheme, where samples from the same participant were assigned a different subject ID. Significance was assessed by 1000 permutations and the covariate was adjusted if necessary. Ordination plots were generated using non-metric multidimensional scaling (NMDS) as implemented in R (“metaMDS” function in the R “vegan” package).
To test for the correlation between organs, we used a permutation test based on Bray-Curtis distance with the test statistic calculated as the distance between the organs from different participants minus the distance between the organs from the same participant. We next permuted each participant for the same organ type using the same block-permutation scheme as above. The p value was calculated as the percentage of permutations that produce a test statistic more extreme than what is observed. To identify the taxa shared by both organs, we used a taxon-specific Euclidean distance, defined based on the presence and absence of a given taxon, and applied the same permutation test. To test whether the distance from cohort 1 to cohort 2 is greater than the distance from cohort 1 to cohort 3, we used a permutation test with the test statistic as the difference between these two distances and block-permutation was used for assessing the significance.
Differential abundance analysis
We conducted differential abundance analysis at phylum, family, and genus levels and filtered rare taxa with prevalence less than 20 % to reduce the number of the tests. We fit a generalized mixed-effects model to the taxa count data using the PQL method, assuming a random intercept for each participant to account for within-subject correlation (“glmmPQL” in R “MASS” package). We fitted an overdispersed Poisson to the counts if the zero proportion is less than 25 % and an overdispersed Binomial model (presence/absence) otherwise. For the overdispersed Poisson model, we included the log of library size as an offset to account for variable sequencing depth. In the overdispersed Binomial model, the log of library size was included as a covariate to account for potential dependence of occurrence probability with sequencing depth. We used the winsorized data (97 % upper quantile) to reduce the potential impact of outliers upon the parameter estimates. To improve power to detect differential taxa, which show consistent change in both the uterus and lower tract microbiome, we pooled the uterus and lower tract data and included the sampling site (uterus/lower tract) as a covariate in the model. The same analyses were also repeated for both datasets separately to confirm the source of the identified signals using pooled data. Statistical significance was assessed based on the Wald test. False discovery rate (FDR) control (B-H procedure, “p.adjust” in standard R packages) was used for correcting for multiple testing, and FDR-adjusted p values or q values will be reported. All statistical analyses were performed in R 3.0.2 (R Development Core Team, Vienna, Austria). The receiver operating characteristic (ROC) curve and area under the curve (AUC) were generated using the median of the replicates with the software generated by Johns Hopkins. (http://www.rad.jhmi.edu/jeng/javarad/roc/).
A total of 31 Caucasian patients undergoing hysterectomy were included in this study. Of those, ten women were diagnosed with a benign gynecologic condition (control cohort), four women were diagnosed with endometrial hyperplasia (cancer precursor, hyperplasia cohort), and 17 women were diagnosed with endometrial cancer (cancer cohort). All diagnoses were made based on the final surgical pathology following hysterectomy. Healthy, asymptomatic women were not included in our study because hysterectomies (surgical removal of the uterus) are not performed on healthy individuals. The inclusion of this population in our study would mandate a different collection protocol to assess the uterine environment and involve the inclusion of multiple confounding variables that could influence the microbiome data. Nevertheless, the inclusion of a variety of benign uterine conditions in our control group provides an assessment of the microbiome that is specifically associated with a cancerous condition and not simply the result of a diseased state. Patients diagnosed with endometrial cancer were significantly older, predominantly postmenopausal, and hypertensive (Table 1).
In order to characterize the microbiome of the patients we collected vaginal and cervical samples (lower genital tract) in the operating room and endometrial, Fallopian, and ovarian samples in the pathology laboratory (collection details are provided in the “Methods” section). The deep-sequencing of the V3-V5 16S rDNA region of all 238 collected samples resulted in the identification of 3545 OTUs. The endometrial microbiome was dominated by Shigella and Barnesiella, with Staphylococcus, Blautia, and Parabacteroides particularly relevant in the benign cohort and Bacteroides and Faecalibacterium more relevant in the endometrial cancer cohort (Fig. 1). The uterine microbiome results are consistent with the very limited number of studies that have assessed the human microbiome composition through culture-based methods, where Escherichia, Streptococcus, Staphylococcus, and Enterococcus were found to be the predominant taxa in women with chronic endometritis and dysfunctional bleeding . The very recent 16S rDNA assessment of the uterine microbiome via transcervical collection is also consistent with Bacteroides being a dominant uterine taxa . In the lower genital tract (vagina and cervix), Prevotella and Lactobacillus were the dominant taxa, with Stenotrophomonas and Shigella more characteristic in the benign cohort and Porphyromonas more common in the endometrial cancer cohort (Fig. 2). These results are also consistent with the pre-menopausal and post-menopausal profiles reported by others [23, 24], with the exception of Stenotrophomonas. Because our benign population is not gynecologically healthy, but instead presented with a variety of conditions (pelvic pain, abnormal bleeding, fibroids, and prolapse), it is possible that Stenotrophomonas may be more prominent in this patient population than in an asymptomatic group of participants. Although it is also possible that this could be the result of contamination, we did not find this taxon to be prominent in our controls (Additional file 1). It is therefore unlikely that this is the case. In the Fallopian tubes, Shigella and Bacteroides were the most dominant taxa, with Staphylococcus, Lactobacillus, Barnesiella, and Pseudomonas commonly appearing in the benign cohort (Fig. 3). In the ovaries, Stenotrophomas, Xanthomonas, and Lactobacillus dominated the benign cohort, while Bacteroides dominated the endometrial cancer cohort (Fig. 4). There is no current literature on the human microbiome composition of Fallopian tubes or ovaries.
Organ microbiome correlation
We first started by assessing if the microbiomes between the different organs were correlated. For instance, whether the vaginal microbiome of a given patient resembled the uterine microbiome of that particular patient more than the uterine microbiome of any other patient. The results showed a very significant correlation between all organs based on a distance-based permutation test (See “Methods” and Table 2). The correlation was also significant, though to a lesser degree, for the stool samples when compared to all organs. The correlation structure held for both benign and cancer cohorts (Additional file 2). Genus level analysis revealed several genera that were significantly shared between the lower genital tract and uterus (Additional file 3). These results are indicative of an overall host specific microbiome effect (host selection effect) and/or transfer of microbiomes across the different organs (microbial movement across organs). The correlation between organs also suggests a potential gain in statistical power by a combined analysis. We thus performed both combined (uterus + lower genital tract) and separate analyses when assessing the microbiota between different disease states.
Overall microbiome structure difference between benign, hyperplasia, and endometrial cancer
We first compared the overall microbiota structure between disease states by investigating the α-diversity and β-diversity. The α-diversity (number of observed OTUs and Shannon index) in the cancer cohort was significantly higher than in the benign cohort (p = 0.003 and 0.01 for the two α-diversity metrics, LME) and the difference was much stronger in uterus (p = 0.03 and 0.01, Fig. 5) than in the lower genital tract (p = 0.17 and 0.31, Additional file 4). The endometrial α-diversity of the hyperplasia cohort was similar to the cancer cohort and was also significantly higher than the benign cohort (p = 0.07 and 0.04, Fig. 5). β-diversity analysis revealed a significant difference in the overall microbiota structure between the three cohorts (p = 0.01, unweighted UniFrac, PERMANOVA, Fig. 6). Consistent with the α-diversity analysis, the difference was mainly observed in the uterus (p = 0.05 and 0.11 for uterus and lower genital tract, unweighted UniFrac). We next conducted pairwise comparisons using the endometrial samples. The endometrial microbiome of both endometrial cancer and hyperplasia cohorts displayed some level of difference from the benign cohort (p = 0.09 and 0.07, unweighted UniFrac). In contrast, the hyperplasia cohort was not distinguishable from the endometrial cancer cohort (p = 0.23, unweighted UniFrac) (Fig. 6). Comparison of the distance between the benign and hyperplasia cohort to the distance between cancer and hyperplasia cohort reveals that hyperplasia is closer to the cancer cohort (p = 0.05, unweighted UniFrac, permutation test; Additional file 5). Interestingly, the distance between the benign and hyperplasia cohort is also significantly larger than that between the benign and cancer cohort (p = 0.05, unweighted UniFrac, Additional file 5). Because endometrial hyperplasia can be a clinical precursor to endometrial cancer, and the uterine microbiome of the four patients diagnosed with endometrial hyperplasia is distinct from the benign cohort and presents some but not complete clustering with an endometrial cancer subgroup, we removed these patients from the primary analysis. This allowed us to compare the benign and endometrial cancer cohorts without the impact of the hyperplasia cases. These were later introduced in a secondary analysis.
The dataset also contains Fallopian and ovarian samples. We therefore tested the microbiota difference between the benign and cancer cohorts for these two organs. Interestingly, we identified a significant difference for the ovaries (p = 0.003, unweighted UniFrac, Additional file 6) suggesting a microbiome connection between the ovarian microniche and endometrial cancer presence/absence.
Endometrial cancer microbiome signature
After the overall microbiome assessment, we performed taxa analysis to determine whether the benign and endometrial cancer cohort displayed differential microbiota. We first performed a combined analysis pooling the samples from both the uterus and lower genital tract. At the genus level there were 12 taxa significantly enriched in the endometrial cancer cohort (Table 3 and Additional file 7, q < 0.10). When we further inquired at a finer level (OTU), we found eight OTUs significantly associated with endometrial cancer (Table 4, q < 0.05). OTU 8 (Atopobium sp.) and OTU 9 (Porphyromonas sp.) became of particular relevance since they were pervasive across samples recovered from endometrial cancer patients and largely absent from the samples recovered from patients in the benign cohort. The Atopobium V3-V5 16S rDNA signature matches (100 %) that of Atopobium vaginae, a well-known vaginal pathogen . The Porphyromonas signature is a close match (99 % sequence identity) to Porphyromonas somerae (Fig. 7), a described pathogen recovered from soft tissue and bone infections . Separate analyses of endometrial and lower genital tract samples revealed a high concordance of the identified genera from the pooled analysis, indicating that both uterine and lower genital tract microbiota may be associated with cancer diagnosis (Table 3).
Vaginal pH and endometrial cancer
Vaginal pH was significantly correlated with an endometrial cancer diagnosis (p = 0.0053), with endometrial cancer patients typically displaying a high vaginal pH (>4.5). However, the vaginal pH is known to raise in approximately 95 % of postmenopausal women  due to physiological and microbiological changes . Therefore, the correlation between endometrial cancer and high vaginal pH could not be detangled from age effects alone. Nevertheless, we were able to determine that the microbiome pH effects were independent of the microbiome disease effects in the uterus since the vaginal pH level was not significantly correlated with the uterine microbiome (p = 0.22 and 0.29, unweighted and weighted UniFrac, PERMANOVA), indicating that they can be used as distinct factors.
Lower tract microbiome association with endometrial cancer
In the lower genital tract, the association of Atopobium vaginae and the identified Porphyromonas sp. with a diagnosis of endometrial cancer has a sensitivity of 73–93 %, and specificity of 67–90 % (Fig. 8). The sensitivity is improved if the vaginal pH is factored in, although specificity is decreased (Table 5; sensitivity – 100 %, specificity – 60 %).
Endometrial hyperplasia microbiome
We had four patients with a final diagnosis of endometrial hyperplasia, which is a known endometrial cancer precursor, in particular in the case of complex hyperplasia with atypia. Three of our patients had simple hyperplasia without atypia (H07, H08, and H63) and one had complex hyperplasia with atypia (H72). Interestingly, the Atopobium vaginae and the Porphyromonas sp. presence/absence profile of the vaginal microbiome of these four patients more closely resembled a benign microbiome signature (Table 5), while the uterine microbiome signature of two of them (H63 and H72) were closer to an endometrial cancer signature.
Snapshots of progression
The correlation and variation between the microbiomes recovered is illustrated in the snapshots, which demonstrate the variable microbiome landscape within and between patients (Fig. 9). We were able to successfully amplify bacterial DNA from 94 % of the lower genital tract samples (vaginal/cervical), 87 % of uterine samples, 50 % of the Fallopian, 61 % of ovarian, 29 % of urine, and 17 % of peritoneal or ascites samples. This progression is likely representative of the bacterial burden in the different body sites.
Here we present a pilot high-throughput microbiome assessment of the female reproductive tract of patients diagnosed with a variety of benign uterine conditions warranting a hysterectomy (abnormal bleeding, fibroids, uterine prolapse, and pelvic pain), endometrial hyperplasia (with and without atypia), and an endometrial cancer diagnosis (endometrioid, mucinous, serous, squamous, and carcinosarcoma). The dominant taxa in the vaginal and cervical (lower tract) microbiome were Prevotella and Lactobacillus, which is consistent with current vaginal microbiome literature . The dominant bacteria in the uterine microbiome were Shigella and Barnesiella, which is also consistent with the existent culture-based literature of the uterine microbial composition in women with endometritis and abnormal bleeding . However, in contrast with our data, these authors reported low concordance between the vaginal and uterine bacteria. Despite the low number of patients in our study (n = 31), the microbiome correlation between all organs (vagina/cervix, uterus, Fallopian tubes, ovaries) was very significant, including the stool recovered from each patient. We interpret these results to indicate that there is a strong selective host effect on the microbiome and/or that there is movement or transfer of microorganisms across the different body sites. Since the samples were not collected simultaneously or in the same location or by the same personnel, we do not believe this correlation to be the result of an artifact. In addition, the stool samples were collected by the patient and processed separately, further supporting the veracity of the correlation. The microbiome organ correlations were equally significant for benign and endometrial cancer patients. Because we do not have healthy asymptomatic patients in this study we cannot assess whether this correlation is generally present or if it may be indicative of a diseased status in itself. We were unable to amplify a positive bacterial signal in about 40–50 % of all Fallopian and ovarian tissue biopsies collected. We conjecture that these results indicate a very low bacterial load (below detection) rather than a truly sterile environment. Further advancements in our ability to detect and amplify bacterial DNA from tissue samples should improve this success rate.
Our results indicate that endometrial hyperplasia can be distinguished from a benign uterine condition based on its microbiome structure, suggesting either a microbiome role in the early phases of cellular transformation or a notorious response to physiologic or chemical gradient shifts within the host’s cellular microenvironment. This structural differentiation is not apparent between hyperplasia and endometrial cancer patients or between benign and endometrial cancer patients. This could be indicative of a transient microbial ecological disturbance that is later normalized to a new equilibrium state, overall closer to the starting point.
Our results suggest that the detection of A. vaginae and the identified Porphyromonas sp. in the gynecologic tract is associated with the presence of endometrial cancer, especially if combined with a high vaginal pH (>4.5). Though cause and effect cannot be discerned from association alone, examination of the hyperplasia cases shows that although these microbes are absent from the lower tract, they are present in half the cases in the uterus, supporting an early disease role for these microbes. A. vaginae has been increasingly recognized as a prominent gynecologic and obstetric pathogen, being positively associated with Nugent scores and bacterial vaginosis , intrauterine infections , and other invasive infections of the female genital tract . While we did not anticipate finding Atopobium to be associated with endometrial cancer, we provided all patients with a comprehensive questionnaire about present and past gynecologic and obstetric diseases or conditions. One question directly asked if the patient had a current or past diagnosis of bacterial vaginosis. Only one patient (H72), a patient with hyperplasia, declared having had bacterial vaginosis in the past (more than six months away from present time). Through verification of the medical records, which is allowed by our institutional IRB, no additional patients had a previous diagnosis of bacterial vaginosis in their medical record. While bacterial vaginosis is a condition difficult to diagnose and not always medically reported, we believe we did all we could reasonably do to assess the prevalence of this disease in our study population. Given that only one patient indicated a past occurrence of bacterial vaginosis, we do not have evidence that this disease impacted our results significantly or differently among our benign and study cohort. While the specific Porphyromonas sp. has yet to be characterized in the literature, the association of members of the Porphyromonas genus with cancers has been recently verified. Porphyromonas gingivalis has been found to be an accurate biomarker for risk of death due to orodigestive cancer, independently of periodontal disease , and the carriage of Porphyromonas species has also been found to be associated with colorectal cancer participants . Given the correlation of these two microorganisms with the disease along with their association with other pathologies, it is possible that they are involved in the etiology or aggravation of conditions leading up to the development of endometrial cancer. Based on the documented association of A. vaginae  with bacterial vaginosis, it is possible that this microbe causes a chronic inflammatory profile that eventually leads to local immune dysregulation and facilitates intracellular infection by Porphyromonas species. Porphyromonas species have been found intracellularly  and it is conceivable that the microorganism we found is capable of disrupting normal cell regulatory functions that may eventually lead to a carcinogenic trigger. The trigger could then be reinforced by the anoxic microenvironment fomented by these microorganisms. We believe this to be a hypothesis worthy of mechanistic investigation.
Our study is limited by a small sample size, due to the technical difficulty of collecting these specimens in real time in the operating room and in the pathology laboratory within a sterile field while guaranteeing that the diagnosis and diagnostic time is not impact by our research procedures. While sample size is always a valid concern, based on the microbiome-based power calculation, the current sample size is powered to detect a relatively large overall effect: 90 % power for an ω2 = 0.04, unweighted UniFrac, an effect size similar to that of antibiotics . Even at this small sample size, we were still able to identify significant microbiome differences between disease states and identify differential abundant taxa after multiple testing correction. Though the increased power may be partially due to the inclusion of multiple replicates for each participant and the use of sophisticated statistical models to address the correlation between replicates and thus reduce the sampling error, these significant results nevertheless indicate that there is a large difference between benign and cancer states.
We investigated potential sources of confounding in this comparative study. The study involved slightly different methods of sampling the microbiota including variation in bacterial DNA enrichment (used to separate bacterial DNA from human DNA in tissue samples), collection type (swab versus scrape versus biopsy), and sampling position (posterior versus superior). Based on marginal PERMANOVA tests, we were able to detect significant effects of bacterial DNA enrichment method and collection type in profiling the microbiota (p < 0.001, unweighted UniFrac) while the sampling position was not significant (p = 0.28 and 0.67, unweighted and weighted UniFrac). However, these technical variables were not true confounders in this comparative study due to roughly equal proportions of different sampling methods in both benign and cancer cohorts (p > 0.3, Fisher’s exact test). In fact, if these technical variables were adjusted in the model, we achieved a similar level of statistical significance in testing the microbiota difference between the two cohorts (data not shown). We thus ruled out the potential confounding effects of these technical variables. Among the demographical and clinical variables, age, BMI, vaginal pH level, menopausal status, and history of hypertension were potential confounders, which had different distributions in the benign and cancer cohorts (Table 1). Marginal PERMANOVA tests on the uterus samples revealed that these variables had less significant effects on the endometrial microbiota than the cohort effect (Additional file 8), indicating that the observed difference could not be completely explained by these potential confounders. However, a larger sample study may be needed to disentangle these confounding effects with confidence. In order to specifically address concerns that the observed findings could be impacted by the age differential between our benign and cancer cohorts, we performed a subset analysis where we extracted patients in the age range of 48–60 years with both a cancer and benign diagnosis. We repeated the analysis with this age-matched subset and we still observed the same trend for both A. vaginae and Porphyromonas, though less significant due to the reduction in the sample size (6 cancer versus 5 benign, Additional file 9). Although age differential is a direct reflection of the patient populations, enrollment targeted efforts will be made in future studies to decrease this gap.
Lastly, while our enrollment exclusion and inclusion criteria did not specify nor exclude any ethnic or racial description, our study population is entirely Caucasian. Although this is a representative reflection of the composition of our patient population at Mayo Clinic, Rochester, MN, it is not representative of the country’s demographics. In future studies we will seek an ethnically diverse patient population to investigate whether our results extend to other populations.
Extending this study to a larger number of patients will allow for the verification of the findings and increase the statistical power. Culturing the identified Porphyromonas sp. and investigating its effects on endometrial cells and their immunological pro-inflammatory profile response, especially in the co-presence of A. vaginae, is warranted. Because of the modifiable nature of the microbiome, these findings also hold promise to endometrial cancer prevention.
We found a distinct microbiome signature in patients with endometrial cancer and hyperplasia. We have shown that in our study population the detection of A. vaginae and the identified Porphyromonas sp. in the gynecologic tract is associated with the presence of endometrial cancer, especially if combined with a high vaginal pH (>4.5). These findings provide important insights into the etiology or manifestation of the disease with broad implications for biomarker development in the early detection of and screening for endometrial cancer.
Area under the curve
Body mass index
False discovery rate
Illinois-Mayo Taxon Operations for RNA Dataset Organization
Linear mixed effects model
Mammalian target of rapamycin
Operational taxonomic unit
Polymerase chain reaction
Phosphatase and tensin homolog
Receiver operating characteristic
King A, Selak M, Gottlieb E. Succinate dehydrogenase and fumarate hydratase: linking mitochondrial dysfunction and cancer. Oncogene. 2006;25:4675–82.
Bansal N, Yendluri V, Wenham RM. The molecular biology of endometrial cancers and the implications for pathogenesis, classification, and targeted therapies. Cancer Control. 2009;16:8.
Hampel H, Frankel W, Panescu J, Lockman J, Sotamaa K, Fix D, et al. Screening for Lynch syndrome (hereditary nonpolyposis colorectal cancer) among endometrial cancer patients. Cancer Res. 2006;66:7810–7.
Zanssen S, Schon EA. Mitochondrial DNA mutations in cancer. PLoS Med. 2005;2:e401.
Collaborators MWS. Endometrial cancer and hormone-replacement therapy in the Million Women Study. Lancet. 2005;365:1543–51.
Fader AN, Arriba LN, Frasure HE, von Gruenigen VE. Endometrial cancer and obesity: epidemiology, biomarkers, prevention and survivorship. Gynecol Oncol. 2009;114:121–7.
Friberg E, Orsini N, Mantzoros C, Wolk A. Diabetes mellitus and risk of endometrial cancer: a meta-analysis. Diabetologia. 2007;50:1365–74.
Nomura A, Stemmermann GN, Chyou P-H, Kato I, Perez-Perez GI, Blaser MJ. Helicobacter pylori infection and gastric carcinoma among Japanese Americans in Hawaii. N Engl J Med. 1991;325:1132–6.
Kostic AD, Gevers D, Pedamallu CS, Michaud M, Duke F, Earl AM, et al. Genomic analysis identifies association of Fusobacterium with colorectal carcinoma. Genome Res. 2012;22:292–8.
Ahn J, Sinha R, Pei Z, Dominianni C, Wu J, Shi J, et al. Human gut microbiome and risk of colorectal cancer. J Natl Cancer Inst. 2013;105:1907–11.
Dossus L, Rinaldi S, Becker S, Lukanova A, Tjonneland A, Olsen A, et al. Obesity, inflammatory markers, and endometrial cancer risk: a prospective case–control study. Endocr Relat Cancer. 2010;17:1007–19.
Braundmeier AG, Lenz KM, Inman KS, Chia N, Jeraldo P, Walther-António MR, et al. Individualized medicine and the microbiome in reproductive tract. Front Physiol. 2015;6:97.
Menard JP, Mazouni C, Salem-Cherif I, Fenollar F, Raoult D, Boubli L, et al. High vaginal concentrations of Atopobium vaginae and Gardnerella vaginalis in women undergoing preterm labor. Obstet Gynecol. 2010;115:134–40.
Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG. Research electronic data capture (REDCap)—a metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform. 2009;42:377–81.
Sipos M, Jeraldo P, Chia N, Qu A, Dhillon AS, Konkel ME, et al. Robust computational analysis of rRNA hypervariable tag datasets. PLoS One. 2010;5:e15220.
Ye Y. Identification and quantification of abundant species from pyrosequences of 16S rRNA by consensus alignment. Proceedings (IEEE Int Conf Bioinformatics Biomed). 2010;2011:153–7.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Arndt D1, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, Sinelnikov I, Budwill K, Nesbø CL, Wishart DS. METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012;40:88–95.
Chen J, Bittinger K, Charlson ES, Hoffmann C, Lewis J, Wu GD, Collman RG, Bushman FD, Li H. Associating microbiome composition with environmental covariates using generalized UniFrac distances. Bioinformatics. 2012;28:2106–13.
McArdle BH, Anderson MJ. Fitting multivariate models to community data: a comment on distancebased redundancy analysis. Ecology. 2001;82(1):290–7.
Cicinelli E, De Ziegler D, Nicoletti R, Colafiglio G, Saliani N, Resta L, et al. Chronic endometritis: correlation among hysteroscopic, histologic, and bacteriologic findings in a prospective trial with 2190 consecutive office hysteroscopies. Fertil Steril. 2008;89:677–84.
Verstraelen H, Vilchez-Vargas R, Desimpel F, Jauregui R, Vankeirsbilck N, Weyers S, et al. Characterisation of the human uterine microbiome in non-pregnant women through deep sequencing of the 16S rRNA gene. PeerJ. 2015;4, e1602.
Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci. 2011;108:4680–7.
Brotman RM, Shardell MD, Gajer P, Fadrosh D, Chang K, Silver MI, et al. Association between the vaginal microbiota, menopause status, and signs of vulvovaginal atrophy. Menopause (New York, NY). 2014;21:450–8.
Marconi C, Cruciani F, Vitali B, Donders GG. Correlation of atopobium vaginae amount with bacterial vaginosis markers. J Lower Genit Tract Dis. 2012;16:127–32.
Summanen PH, Durmaz B, Väisänen M-L, Liu C, Molitoris D, Eerola E, et al. Porphyromonas somerae sp. nov., a pathogen isolated from humans and distinct from Porphyromonas levii. J Clin Microbiol. 2005;43:4455–9.
Freedman M. Vaginal pH, estrogen and genital atrophy. Menopause Manag. 2008;17:9–13.
Farage M, Maibach H. Lifetime changes in the vulva and vagina. Arch Gynecol Obstet. 2006;273:195–202.
Knoester M, Lashley LE, Wessels E, Oepkes D, Kuijper EJ. First report of Atopobium vaginae bacteremia with fetal loss after chorionic villus sampling. J Clin Microbiol. 2011;49:1684–6.
Chan JF, Lau SK, Curreem SO, To KK, Leung SS, Cheng VC, et al. First report of spontaneous intrapartum Atopobium vaginae bacteremia. J Clin Microbiol. 2012;50:2525–8.
Ahn J, Segers S, Hayes RB. Periodontal disease, Porphyromonas gingivalis serum antibody levels and orodigestive cancer mortality. Carcinogenesis. 2012;33:1055–8.
Onderdonk AB, Delaney ML, Fichorova RN. The human microbiome during bacterial vaginosis. Clin Microbiol Rev. 2016;29:223–38.
Rudney JD, Chen R, Sedgewick GJ. Intracellular Actinobacillus actinomycetemcomitans and Porphyromonas gingivalis in buccal epithelial cells collected from human subjects. Infect Immun. 2001;69:2700–7.
Kelly BJ, Gross R, Bittinger K, Sherrill-Mix S, Lewis JD, Collman RG, et al. Power and sample-size estimation for microbiome studies using pairwise distances and PERMANOVA. Bioinformatics. 2015;31:2461–8.
MWA, HN, and NC acknowledge support from the Center for Individualized Medicine, Microbiome Program at Mayo Clinic. MWA acknowledges support from the Microbial Community Metabolism and Health Program (MCMHP) Career Development Award Education/Disease. AM and NC acknowledge support from the John W. Anderson Foundation. FM acknowledges support from the Division of Obstetrics and Gynecology, University of Cagliari, Italy.
Availability of data and materials
The dataset of sequence files supporting the conclusions of this article are available in SRA repository (SRP064295).
MWA, DC, HN, AM, and NC initiated and designed the study. MWA, FM, AH, TD, and EHC acquired the data. JC performed the statistical analysis. MWA, JC, AH, HN, AM, and NC analyzed and interpreted the data. MWA, JC, FM, AH, GK, HN, AM, and NC wrote or revised the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Consent was obtained from each participant included in the study; however, images and data included in the study have been de-identified.
Ethics approval and consent to participate
Participants were consented under IRB #12-004445, which was approved and reviewed by the Mayo Clinic Institutional Review Board. All participants provided written consent.
Sequencing results from controls. A total of 14 controls were performed, with five of them not retrieving any sequence reads. sV1, sV2 – V-shaped scrape with Microbiome Enrichment Nebnext kit; sP1, sP2 – Pap scrape with Microbiome Enrichment Nebnext kit; sT1, sT2 – TE buffer with Microbiome Enrichment Nebnext kit; sS1 – Swab with Microbiome Enrichment Nebnext kit; sV1B, sV2B – V-shaped scrape. (XLSX 26 kb)
Organ correlation with benign and endometrial cancer patients. (XLSX 43 kb)
Genera correlation between lower tract and uterus. (XLSX 51 kb)
α-diversity for the lower tract. A Observed OTU number. B Shannon index. (PDF 384 kb)
Unweighted UniFrac distance between benign, hyperplasia, and endometrial cancer cohorts. (PDF 8 kb)
Unweighted UniFrac distance for ovary between benign and endometrial cancer cohort. (PDF 9 kb)
Genera significantly enriched in the endometrial cancer cohort. (PDF 11 kb)
Confounder analysis. (XLSX 37 kb)
Cohort age differential subset analysis. (XLS 33 kb)