Immunohistochemistry staining and 16S ribosomal RNA gene (16S rRNA gene) sequencing were performed in order to discover and distinguish bacteria in cancerous and noncancerous ovarian tissues in situ. The compositional and functional differences of the bacteria were predicted by using BugBase and Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt).
The theory of a “sterile” female upper reproductive tract has been encountering increasing opposition due to advancements in bacterial detection. However, whether ovaries contain bacteria has not yet been confirmed yet. Herein, an experiment to detect bacteria in ovarian tissues was introduced. We chose ovarian cancer patients in the cancer group and noncancerous patients in the control group. 16S rRNA gene sequencing was used to differentiate bacteria in ovarian tissues from the cancer and control groups. Furthermore, we predicted the functional composition of the identified bacteria by using BugBase and PICRUSt. This method can also be used in other viscera and tissues since many organs have been proven to harbor bacteria in recent years. The presence of bacteria in viscera and tissues may help scientists evaluate cancerous and normal tissues and may be aid in the treatment of cancer.
Recently, an increasing number of articles have been published that prove the existence of bacteria in abdominal solid viscera, such as the kidney, spleen, liver, and ovary1,2. Geller et al. found bacteria in pancreatic tumors, and these bacteria were resistant to gemcitabine, a chemotherapeutic drug2. S. Manfredo Vieira et al. concluded that Enterococcus gallinarum was portable to the lymph nodes, liver and spleen, and it could drive autoimmunity3.
Since the cervix plays a role as a defender, bacteria in the upper female reproductive tract, which contains the uterus, fallopian tubes, and ovaries, have been minimally researched. However, some new theories have been established in recent years. Bacteria may have access to the uterine cavity during the menstrual cycle due to changes in mucins4,5. Additionally, Zervomanolakis et al. confirmed that the uterus, together with the fallopian tubes, is a peristaltic pump controlled by the endocrine system of the ovaries, and this arrangement enables bacteria to enter the endometrium, fallopian tubes, and ovaries6.
The upper reproductive tract is no longer a mystery anymore thanks to the development of bacterial detection methods. Verstraelen et al. used a barcoded paired-end sequencing method to discover uterine bacteria by targeting at the V1-2 hypervariable region of the 16S RNA gene7. Fang et al. employed barcoded sequencing in patients with endometrial polyps and revealed the presence of diverse intrauterine bacteria8. Additionally, by using the 16S RNA gene, Miles et al. and Chen et al. found bacteria in the genital system of women who had undergone salpingo-oophorectomy and hysterectomy, respectively5,9.
Bacteria in tumor tissues have gained increasing attention in recent years. Banerjee et al. discovered that the microbiome signature differed between ovarian cancer patients and controls10. Anoxynatronum sibiricum was associated with tumor stage, and Methanosarcina vacuolata might be used to diagnose ovarian cancer11. In addition to ovarian cancer, other cancers, such as stomach, lung, prostate, breast, cervix, and endometrium, have been proven to be associated with bacteria12,13,14,15,16,17,18. Poore et al. proposed a new class of microbial-based oncology diagnostics, foreseeing early-stage cancer screening19. In this protocol, we investigated the differences between cancerous and normal ovarian tissues by comparing the composition and function of bacteria in these two tissues.
Immunohistochemistry staining and 16S rRNA gene sequencing were performed to confirm the presence of bacteria in the ovaries. The differences and predicted functions of the ovarian bacteria in cancerous and noncancerous ovarian tissues were studied. The results showed the existence of bacteria in ovarian tissues. Anoxynatronum sibiricum and Methanosarcina vacuolata were related to the stage and the diagnosis of ovarian cancer, respectively. Forty-six significantly different KEGG pathways that were present in both groups were compared.
This study was approved by the Medical Institutional Ethics Committee of the First Affiliated Hospital of Xi'an Jiaotong University (No. XJTUIAF2018LSK-139). Informed consent was obtained from all enrolled patients.
1. Criteria for entering the cancer group and the control group
2. Gather samples
3. Sequence the 16S rRNA gene
4. Analyze 16S rRNA gene sequencing data
5. Predict bacterial function
6. Data
Patients
A total of 16 qualified patients were included in the study. The control group included 10 women with a diagnosis of benign uterine tumor (among them, 3 patients were diagnosed with uterine myoma, and 7 patients were diagnosed with uterine adenomyosis). Meanwhile, the cancer group contained 6 women with a diagnosis of serous ovarian cancer (among them, 2 patients were diagnosed with stage II, and 2 of them were diagnosed with stage III). The following characteristics showed no differences between patients in the control group and the cancer group: age, menopausal status, parity, history of hypertension, and history of diabetes (Table 1).
Table 1: Patient statistics Please click here to download this Table.
The richness and variety of ovarian bacterial species in both groups
The presence of bacteria using immunohistochemistry staining was shown in Figure 1. The alpha diversity of the microbes was analyzed as a method to detect the richness and variety of ovarian bacterial species. The number of species observed in the ovarian cancer tissues was smaller than that of the control group, with no significant difference. The richness (represented by the Chao 1 and ACE index) and the diversity (represented by the Shannon, Simpson, and Evenness index) of the bacterial species were both not significantly different between the cancer group and the control group (Figure 2).
Figure 2: 16S rRNA gene sequencing shows differences between the cancer and control groups in bacterial richness and diversity. (A) Observed species index (P = 0.06, Mann-Whitney U test); (B) Chao 1 index (P = 0.06, Mann-Whitney U test); (C) ACE index (P = 0.06, Mann-Whitney U test); (D) Shannon index (P = 0.32, Mann-Whitney U test; E. Evenness index (P = 0.48, Mann-Whitney U test); (F) Simpson index (P = 0.46, Mann-Whitney U test). We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Description of ovarian bacteria
Deep sequencing of the V3-V4 16S rRNA gene region was performed on all samples to obtain a better understanding of the ovarian bacteria. The results showed that Proteobacteria was the most abundant phylum (67.10% in the control group and 67.20% in the cancer group), Firmicutes was the second most abundant phylum (23.77% in the control group and 23.82% in the cancer group), and the third most abundant phylum was Bacteroidetes (3.26% in the control group and 3.41% in the cancer group). When analyzing species of the control group, the main composition was consisted of Halobacteroides halobios (14.53%), followed by Gemmata obscuriglobus (11.07%) and Methyloprofundus sedimenti (10.69%). For the cancer group, Gemmata obscuriglobus was the richest in the cluster (13.89%), followed by Halobacteroides halobius (11.99%) and Methyloprofundus sedimenti (11.12%) (Figure 3).
Figure 3: Relative abundance of phyla (> 1%) and of the top 12 species in ovarian samples. (A) The relative abundance of the phyla (> 1%) in the ovaries of the patients in the control group. (B) The relative abundance of the phyla (> 1%) in the ovaries of patients with ovarian cancer. (C) The relative abundances of the 12 most abundant bacterial species in the ovaries of the control patients. (D) The relative abundances of the 12 most abundant bacterial species in the ovaries of ovarian cancer patients. We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Different compositions of ovarian bacteria between the two groups
A comparison of different bacterial communities was carried out by PCoA using PERMANOVA. The results showed that the bacteria in the control group differed from those in the cancer group, P < 0.05 (Figure 4).
Figure 4: PCoA detects clusters of communities and the relative abundances of Anoxynatronum sibiricum and Methanosarcina vacuolata. (A) Communities were clustered using PCoA. PC1 and PC2 are plotted on the x and y axes. The red block indicates a sample in the ovarian cancer group. The blue circle indicates a sample in the control group. The samples from the ovarian cancer group were separated from other samples in the control group. (B) Communities clustered using PCoA. PC1 and PC2 are plotted on the x and y axes. The red block indicates a sample in the ovarian cancer group. The blue solid circle indicates a sample from a patient with uterine myoma, and the blue hollow circle is equal to a sample of a patient with uterine adenomyosis. (C) The relative abundance of Anoxynatronum sibiricum (control group: n = 10, cancer group: n = 6, P = 0.034, Mann-Whitney U test). (D) The relative abundance of Methanosarcina vacuolata (control group: n = 10, cancer group: n = 6, P = 0.001, Mann-Whitney U test). We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Ovarian bacterial composition in cancer and control groups from different perspectives
An analysis of the ovarian bacterial composition was performed from different perspectives to further detect the differences in the identified ovarian bacteria. In Table 2, phylum, class, order, family, genus, and species levels were considered, and statistics are provided in the chart. In particular, there was an association between the relative abundance of Anoxynatronum sibiricum and the stage of the tumor, and Methanosarcina vacuolata was a specific sign when diagnosing ovarian cancer (Table 2).
Phenotypic conservation of ovarian bacteria in the two groups based on predicted functions
In the cancer group, the expression of genes related to potentially pathogenic and oxidative stress-tolerant phenotypes was increased compared with that of the control group (Wilcoxon signed-rank test, P =0.02 and P =0.002). No significant difference was found between the ovarian cancer and control groups in the following aspects: the phenotypes of aerobic, anaerobic, facultatively anaerobic, gram-positive, and gram-negative bacteria; mobile elements; and biofilm formation of the ovarian bacteria (Figure 5). Forty-six variant KEGG pathways between the bacteria in ovaries in the cancer and control groups were determined. The ovaries in the cancer group showed 26 increased pathways. Among them, the most highly related pathways were transporters. On the other hand, the bacteria in ovarian cancer tissue showed 20 reduced pathways. The most relevant functions were as follows: secretion system, unknown functions, and two-component system. The rest of the pathways are shown in Figure 6.
Figure 1: LPS immunohistochemical expression in ovaries. (A) Control group (10 x). Scale bars, 200 µm. (B) Control group (40 x). Scale bars, 50 µm. (C) Cancer group (10x). Scale bars, 200 µm. (D) Cancer group (40 x). Scale bars, 50 µm. Arrows point to LPS staining in the ovarian tissue. We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Figure 5: Predicted metagenomes analyzed by BugBase. The expression of some genes in the cancer group was increased compared with that in the control group. These genes were related to potentially pathogenic (Wilcoxon signed-rank test, P = 0.02) and oxidative stress-tolerant phenotypes of the ovaries. (Wilcoxon signed-rank test, P = 0.002). We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Figure 6: PICRUSt analysis of different KEGG pathways between the cancer and control groups. We obtained reprint permission from previous publishers. This figure has been modified from Wang et al.11. Please click here to view a larger version of this figure.
Table 2: Richness (represented by the Chao 1 and ACE index) and the diversity (represented by the Shannon, Simpson, and Evenness index) of the bacterial species Please click here to download this Table.
Ovarian cancer has a notable influence on women's fertility25. Most ovarian cancer patients are diagnosed at late stages, and the 5-year survival rate is less than 30%18. Confirmation of bacteria in the abdominal solid viscera, including the liver, pancreas and spleen, has been published. The existence of bacteria in the upper female reproductive tract occurs because the cervix is not enclosed2,3,4,5. However, whether ovaries, which are abdominal solid viscera, are sterile or not has not yet been determined. Additionally, whether bacteria in the ovaries are related to ovarian cancer is also an important question.
The significant differences in the bacteria that we found were compared between different groups. All of the procedures mentioned above were strictly germfree, including instruments, reagents, equipment, and the operation of the whole protocol. More importantly, we used ovaries from patients with benign uterine disease as the control group to counteract possible contamination. However, in this protocol, contamination cannot be avoided. Thus, since the cancer group and control group were analyzed in the same experimental environment, merely by comparing the differences between these two groups, we could obtain primary evidence about the microbiological origin of ovarian cancer.
The findings of bacteria in ovarian tissue might start a new field investigating the bacteria influencing ovarian cancer. Additionally, the unique presence and composition of bacteria in cancerous ovarian tissues might direct the carcinogenesis of ovarian cancer, and the therapeutic and prognostic targets of bacteria. Among the 46 KEGG pathways, functions related to the biosynthesis of vancomycin group antibiotics drew particular attention. This may provide further treatment options for ovarian cancer.
However, the protocol had some limitations. First, the samples could not be collected from healthy people for ethical reasons. The control group was ovaries of patients with benign uterine disease (including uterine myoma and adenomyosis). Second, the number of samples should be larger. The limited sample size of the study might hamper the accuracy of the results. Third, although the cancer group and the control group were under the same conditions, there were no negative or positive controls. In addition, contamination could not be avoided. To date, the study of the ovarian bacteria in patients with ovarian cancer is still in an early stage. A larger-scale study with more samples is needed.
The authors have nothing to disclose.
This work was supported by the Clinical Research Award of the First Affiliated Hospital of Xi'an Jiaotong University, China (XJTU1AF-2018-017, XJTU1AF-CRF-2019-002), the Major Basic Research Project of Natural Science of Shaanxi Provincial Science and Technology Department (2018JM7073, 2017ZDJC-11), the Key Research and Development Project of Shaanxi Provincial Science and Technology Department (2017ZDXM-SF-068, 2019QYPY-138), the Shaanxi Provincial Collaborative Technology Innovation Project (2017XT-026, 2018XT-002), and the Medical Research Project of Xi'an Social Development Guidance Plan (2017117SF/YX011-3). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
We thank the colleagues in the Department of Gynecology of First Affiliated Hospital of Xi'an Jiaotong University for their contributions to collecting samples.
2200 TapeStation Software | Agilgent United States |
||
AmpliSeq for Illumina Library Prep, Indexes, and Accessories | Illumina | ||
Image-pro plus 7 | Media Cybernetics | ||
Leica ASP 300S | Leica Biosystems Division of Leica Microsystems | ||
Leica EG 1150 | Leica Biosystems Division of Leica Microsystems | ||
Leica RM2235 | Leica Biosystems Division of Leica Microsystems | ||
LPS Core monoclonal antibody, clone WN1 222-5 | Hycult Biotech | ||
Mag-Bind RxnPure Plus magnetic beads | Omega Biotek | M1386-00 | |
Mag-Bind Universal Pathogen 96 Kit | Omega Biotek | M4029-01 | |
MiSeq | Illumina | SY-410-1003 | |
Silva database | Max Planck Institute for Marine Microbiology and Jacobs University | ||
the QuantiFluor dsDNA System | Promega | E2670 | |
Trimmomatic | Björn Usadel | ||
ZytoChem Plus (HRP) Anti-Rabbit (DAB) Kit | Zytomed Systems | HRP008DAB-RB |