Cervical Microbiome and Cytokine Profile at Various Stages of Cervical Cancer: A Pilot Study

Abstract

Cervical cancer (CC) is caused by high-risk human papillomavirus persistence due to the immunosuppressive tumor microenvironment mediated by cytokines. Vaginal microbiota determines the presence of certain cytokines locally. We assessed the association between cervical microbiota diversity and the histopathological diagnosis of each stage of CC, and we evaluated mRNA cervical expression levels of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ across the histopathological diagnosis and specific bacterial clusters. We determined the cervical microbiota by high throughput sequencing of 16S rDNA amplicons and classified it in community state types (CST). Mean difference analyses between alpha-diversity and histopathological diagnosis were carried out, as well as a β-diversity analysis within the histological diagnosis. Cervical cytokine mRNA expression was analyzed across the CSTs and the histopathological diagnoses. We found a significant difference in microbiota's diversity in NCL-HPV negative women vs those with squamous intraepithelial lesions (SIL) and CC(p = 0.006, p = 0.036).When β-diversity was evaluated, the CC samples showed the highest variation within groups (p<0.0006) and the largest distance compared to NCL-HPV negative ones (p<0.00001). The predominant bacteria in women with normal cytology were L. crispatus and L. iners, whereas for SIL, it was Sneathia spp. and for CC, Fusobacterium spp. We found higher median cervical levels of IL-4 and TGF-β1 mRNA in the CST dominated by Fusobacterium spp. These results suggest that the cervical microbiota may be implicated in cervical cancer pathology. Further cohort studies are needed to validate these findings.

About the Authors:

Astride Audirac-Chalifour

Contributed equally to this work with: Astride Audirac-Chalifour, Kirvis Torres-Poveda

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Kirvis Torres-Poveda

Contributed equally to this work with: Astride Audirac-Chalifour, Kirvis Torres-Poveda

Current Address: Consejo Nacional de Ciencia y Tecnología de Mexico (CONACyT), Mexico City, Mexico

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Margarita Bahena-Román

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Juan Téllez-Sosa

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Jesús Martínez-Barnetche

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Bernardo Cortina-Ceballos

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Guillermina López-Estrada

Affiliation: Private Health Center for Gynecology, Cuernavaca, Morelos, Mexico

Karina Delgado-Romero

Affiliation: Centro de Atención para la Salud de la Mujer (CAPASAM) (Center for Women’s Health), Health Services of the State of Morelos, Cuernavaca, Mexico

Ana I. Burguete-García

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

David Cantú

Affiliation: Division of Clinical Research, Instituto Nacional de Cancerología (INCan), SS, Mexico City, Mexico

Alejandro García-Carrancá

Current Address: Unit of Biomedical Research in Cancer, Instituto Nacional de Cancerología (INCan), SS and Biomedical Research Institute, Universidad Nacional Autónoma de México, Mexico City, Mexico

Affiliation: Division of Basic Research, Instituto Nacional de Cancerología (INCan), SS, Mexico City, Mexico

Vicente Madrid-Marina

* E-mail: [email protected]

Affiliation: Dirección de Infecciones Crónicas y Cáncer, Centro de Investigación sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública (INSP) (Chronic Infectious Diseases and Cancer Division, Center for Research on Infectious Diseases, National Institute of Public Health), Cuernavaca, Morelos, Mexico

Introduction

Cervical cancer (CC) is the fourth most common cancer in women and the seventh overall worldwide, with an estimated 485 000 new cases and 236 000 deaths in 2013. CC caused 6.9 million disability-adjusted life-years (DALYs) in 2013 [1], and it was the second most common cause of death by cancer among Mexican women in 2011 (10.4%) [2]. CC is caused by a persistent infection with high-risk human papillomavirus (HR-HPV). However, a HR-HPV infection is considered a necessary but not a sufficient cause for CC development [3]. Mechanical factors, like vaginal douching or sexual intercourse, and biological factors, like bacterial vaginosis (VB) [4, 5], or sexually transmitted infections (STIs) [6] alter the vaginal microenvironment and have been identified as cofactors in the persistence of an HPV infection [7].

The vast majority of HR-HPV infected women never develop CC because an adequate immune response is capable of controlling the infection and preventing its progression to precancerous lesion [8]. This suggests that additional factors act in conjunction with the HPV to influence the risk of CC development. So far, the following determinant co-factors have been associated to the appearance of CC: socio-environmental factors (cultural barriers, extreme poverty, poor sanitation areas, and limited access to health care) [9], epidemiological factors (multiparity, use of oral contraceptives for more than five years, multiple sexual partners and smoking) [10, 11], and genetic factors related to the host (polymorphisms in immune response genes, which determine a deficient immune response and local immunosuppression) [12– 14]. Therefore, CC is a complex multifactorial disease, and further research is needed to improve our understanding of its etiology.

It has been recently proposed that abnormal vaginal microbiota plays an important role in the development of cervical neoplasm [15]. An epidemiologic study identified Chlamydia trachomatis as a co-factor for CC development [16]. Similarly, other studies show that some bacterial species seem to be associated with the development of other cancers, such as colon cancer; these studies also suggest that various bacterial species preferentially inhabit tumor sites [17, 18].

There are still several gaps in knowledge regarding the association between vaginal and cervical microbiomes and CC development [19]. Bacterial culture-based evidence indicates that some potential pro-oncogenic pathogens, which may be members of commensal microbiota, contribute to tumor initiation and development [20, 21].

CC is a long-standing disease, and there are stages previous to it during which the conditions of the cervical and vaginal environment are modified, including vaginal acidity and cytokine pattern that lead to a local immunosuppression state. The presence of Lactobacilli, a low vaginal pH (< 4.5) and antimicrobial peptides are part of the defense mechanisms present in the vaginal microenvironment [22]. When an imbalance of the defense system occurs, physicochemical changes arise and produce histological alterations of the vaginal mucosa and the cervical epithelium, all of which conditions a selective pressure on the microbiota [23]. In a CC microenvironment, the presence of immunosuppressive cytokines (TGF-ß1, IL-10) favors the persistenceof the HPV infection [24, 25, 26, 27]. Most studies of the female genital tract microbiome have been carried out at the vaginal level and few studies involve the cervical microbiota and cytokine profiles. A comparative genomic analysis of the vaginal microbiome has identified changes in microbiota diversity among women with genital human immunodeficiency virus (HIV) infection, with or without bacterial vaginosis (BV) [28] and in pregnant women [29]. Additionally, it has been proposed that the vaginal microbial ecosystem and the cytokine profile play a role in promoting cervical dysplasia [30], given that an abnormal vaginal microbiota has been associated with the acquisition of an HPV infection [31]. However, few studies have been carried out on the cervical microbiome as a modifier of the HPV natural history with respect to the development of cervical lesions and cervical neoplasm[32].

Considering that it has been proved that specific species of cervicovaginal microbiota modulate the inflammatory immune response in the female genital tract of healthy black South African women [33], it is possible that the cervical microbiome is involved in promoting the expression of immunosuppressive cytokines [34].We have hypothesized that the HPV infection and the development of SIL and CC are associated with changes in microbiota diversity and with cytokine expression patterns at the cervical level. Our study examined the association between cervical microbiota diversity and composition, according to an histopathological diagnosis of each stage of the natural history of CC, and the cervical expression levels of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA.

Material and Methods

Study design and population

A cross-sectional study was carried out, using samples from a biological bank built between 2008 and 2011 from newly diagnosed cases of squamous intraepithelial lesions (SIL) (n = 268 HPV-positive samples); women with a negative Papanicolaou and a normal colposcopy as controls (n = 205, 81 HPV-negative and 124 HPV-positive samples), recruited from the Women’s Health care Center in the State of Morelos (Centro de Atención para la Salud de la Mujer del Estado de Morelos) in Mexico between June 2008 and November 2011, and from cervical squamous cell carcinoma cases (n = 171 HPV-positive samples), recruited from the Gynaecology Service at the National Cancer Institute (INCan) in Mexico City, between September 2010 and December 2011. The Bioethics and Research Committees at INCan (reference number. INCan/CC/326/10CB/609) and at the National Institute of Public Health (Instituto Nacional de Salud Pública–INSP; reference number: CI814) approved the baseline study during which the biological bank was built. In addition, all participants gave their informed consent to use their biological samples for further research studies. Every subject was interviewed for lifestyle, socio-demographic and reproductive factors known to be associated with increased risk of CC [12].

For methodological convenience, we selected 32 cases for this study, non-cervical lesions (NCL: n = 10 HPV-negative; n = 10 HPV-positive), SILs (n = 4 HPV-positive) and CC (n = 8 HPV-positive). Inclusion criteria for the subsample selection related to the medical history were: patient´s recruitment on the same day of menstrual period (seven postretirement menstrual period days), the non-use of douches and no sexual activity in previous days of the sampling. Also, not having records of antibiotic or antifungal use in the last 30 days previous to sampling. Other inclusion criteria included: molecular HPV + diagnosis; DNA and RNA purity and integrity suitable for sequencing and quantifying mRNA by qRT-PCR, and being Mexican with Mexican parents and grandparents to ensure similar ancestry. The only exclusion criterion was not having sufficient reads. “We considered 1 000 as the threshold, after performing the bioinformatics sequence analysis. The studied population’s socio-demographic and reproductive-sexual characteristics are presented in Table 1. All of the study subjects were Mexican women whose age ranged from 22 to 61 years. There was a significant difference between the groups regarding the contraceptive method and the positive HPV test. CC patients reported no use of contraceptive methods more frequently than the NCL patients. Regarding the SIL cases, the most prevalent HPV genotypes were HR-HPV (non-HPV16 or 18), whereas among CC patients the most prevalent genotype was HPV 16.

[Figure omitted. See PDF.]

Table 1. Analysis of reproductive/sexual lifestyle-related risk factors in the study population (cervical lesion and cervical cancer patients).

https://doi.org/10.1371/journal.pone.0153274.t001

Characteristics of the biological samples bank

Samples were taken from women’s cervix seven days after menses withdrawal. The biological bank used for this study is composed of DNA samples and cDNA extracted from cervical epithelial scraping swabs from women diagnosed with NCL and from fresh cell biopsies from women diagnosed with SIL and CC. Genomic DNA was extracted from cervical epithelial scrapings and biopsies previously digested with proteinase K, using the Genomic DNA Purification Kit (Fermentas Life Sciences, Vilnius, Lithuania). All possible measures were taken to avoid cross-contamination of samples during DNA extraction. DNA concentration and purity were evaluated by Thermo Scientific NanoDropTM 1000 Spectrophotometer (260/280) and DNA integrity was determined by electrophoresis in agarose gels at 1%. Total RNA was isolated from cervical samples using the Trizol reagent from Invitrogen. cDNA synthesis was carried out in the presence of 200 U of M-MLV reverse transcriptase and 2.5 ug of total RNA using standard conditions. The PCR reactions were performed in a reaction volume of 25 μL containing 1 μL of cDNA, 0.2 mM of dNTPs, 15 pmol of each primer, 2.5 μL of reaction buffer and1U of Taq DNA polymerase recombinant. Primers for the human housekeeping glyceraldehyde-3-phosphate dehydrogenase GAPDH (450 pb) were used to verify cDNA integrity.

Cervical specimens were previously tested for HPV. Viral DNA fragments from the samples were amplified by PCR using consensus primers MY09/MY11 [35], LIC1/LIC2 [36], and GP5/GP6 [37], which flank the L1 region of the HPV capsid. PCR amplification of GAPDH (556pb) was used as an internal control for DNA quality. Cell lines expressing HPV-16 (SiHa) and HPV-18 (HeLa) were used as positive controls. Deionized H₂O served as a negative control. All products were analyzed by electrophoresis in acrylamide gels at 6%. Positive samples were visualized in an agarose gel at 1.5%. The DNA band obtained was extracted and purified with the MinElute Gel Extraction Kit (Qiagen, Hilden, Germany) and sequenced using the Sanger method. The HPV sequences were analyzed by BLAST. HPV was categorized according to its phylogenetic patterns into low- and high-risk types: HPV16 and HPV18. HPV status was confirmed with the Anyplex^TM II HPV HR Detection assay from Seegene^®, based on multiplex real-time PCR, TOCE and DPO primer pairs technology, according to the supplier’s instructions. [38]

Ethics statement

This study was conducted according to the principles expressed in the Declaration of Helsinki. It was approved by the Research, Ethics and Biosafety Committees at INSP (CI:1143). Written informed consent was obtained from all participants.

Analysis of cervical cytokine mRNA expression

A qRT-PCR amplification was performed in duplicate for the gene expression analysis with the following TaqMan probes: GAPDH (ID-Hs99999905_mL), IL-4 (ID-Hs00174122_mL), IL-6 (ID-Hs00174131_mL), IL-10 (ID-Hs00961622_mL), TGF-β1 (ID-Hs00961622_mL), TNF-α (ID-Hs00174128_mL) and IFN-γ(ID-Hs00174143_mL). The amplification mix was prepared by adding 100 ng of each cDNA sample to a final reaction mixture of 10 μl containing 5μl of TaqMan PCR Master Mix for expression, 0.5 μl of probe and 3.5 μl of DNase-free molecular grade water. The amplification cycles(performed on a StepOnePlus™ from Applied Biosystems, Foster City, CA, USA) were as follows: 94°C for 10 minutes, 40 cycles at 94°C for one minute, 54°C for one minute, 72°C for one minute and 30 seconds, followed by 72°C for 15 minutes. GAPDH was used to normalize the amount of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA present in each sample [39]. Peripheral blood mononuclear cells stimulated with phytohemagglutinin for 72 hours were used to determine the dynamic range curves of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ expression; all standard curves were realized in triplicate. The expression level of mRNA for each cytokin estudied was calculated using relative quantification with the comparative Ct(2-ΔCt) method, taking GAPDH as the endogenous gene. Samples were analyzed in duplicate.

High-throughput sequencing of 16S rDNA amplicons

Amplicons of ~456 bp containing V3-V4 variable regions from 16S rRNA genes were obtained for DNA libraries preparation. Primers 347F Forward 5'-GGAGGCAGCAGTRRGGAAT-3' and 803R Reverse 5'-CTACCRGGGTATCTAATCC-3', described by Nossa, were used [40]. A first, PCR was performed with the following conditions: 50 ng of template from each DNA sample were added to a final reaction mixture of 30 μl containing 1X reaction buffer, 2.5 mM MgSO₄, 1mM dNTP, 0.2 U Platinum Taq DNA Polymerase High Fidelity (Invitrogen) and 5 pmol/μL of the primers. The amplification program was the following: 94°C for 3 minutes, then 30 cycles at 94°C for 30 seconds, 55°C for 30 seconds, 68°C for 30 seconds, followed by 68°C for 5 minutes. Amplicons were visualized by electrophoresis in agarose gel at 1%. The DNA bands obtained were purified by MinElute Gel Extraction Kit (Qiagen, Hilden, Germany).

Amplicons were re-amplified and processed with the same primers linked to adapter A of the 454 sequencing protocol followed by a 6-mer multiplex identifier and by the primer 347F. The reverse primer was the same for all the reactions linked to adapter B of the 454 sequencing protocol. The reaction conditions were as follows: 0.0625 ng of each DNA sample to a final reaction mixture of 30 μl containing 1X reaction buffer, 2.5 mM MgSO₄, 1mM dNTP, 0.2 U Platinum Taq DNA Polymerase High Fidelity (Invitrogen) and 5 pmol/μL of primers 347F and 803R. Amplification cycles were the following: 94°C for 3 minutes, then 15 cycles at 94°C for 30 seconds, 55°C for 30 seconds, 68°C for 30 seconds, followed by 68°C for 5 minutes. Gel electrophoresis was performed in agarose gel at 1.5% to visualize the integrity of the amplified products. DNA fragments were purified by MinElute Gel Extraction Kit(Qiagen, Hilden, Germany) and afterwards quantified and assessed. Amplicons were pooled in eight libraries and mixed in an equimolar manner. Subsequently, the libraries were purified with Ampurebeads XP (Beckman Coulter, Inc). Emulsion PCR titration and yield calculation were then performed to determine the volume of beads to load into sequencing plates. Subsequently, massive sequencing was performed on the Genome Sequencer Titanium Roche-454 (AppliedScience) platform.

Bioinformatics analysis

The raw sequences of the V3-V4 region from the 16SrRNA gene generated for each amplicon were processed with the Roche-454 complementary software, and *.sff files were generated. High quality reads were selected with the following criteria: reads with more than five consecutive bases with scores under 20 on the Qphred scale on a 30 bases moving window were discarded. Quality control (QC) was supervised with Thomas Girke’s R script that uses a fastq Quality.R library to plot the distribution of the rating assigned to each base according to its position on the Qphred logarithmic scale [41]. Once quality files and sequences were obtained separately from the *.sff raw files, demultiplexing was conducted in QIIME to identify the sequences belonging to each sample according to their molecular identifier (MID) [42].

The sequences were grouped into operational taxonomic units (OTUs) using the PyNAST algorithm for sequence alignment. All sequences which had a similarity of 99% or more were considered a single OTU [43]. A representative sequence of each OTU was chosen for later taxonomic identification. Taxa were assigned with the Uclust algorithm [44], using the Green Genes database (99% release in May 2013) as reference, so that any representative sequence aligned with a 99% similarity to a reference sequence would be assigned the same species name [45].

Sequences that did not align to the reference were extracted to another file for de novo assembly and to be considered in the diversity calculus. To determine if microbiota composition was able to cluster samples by diagnosis, an unsupervised hierarchical clustering based on the Bray Curtis dissimilarity between OTU abundance of each sample (S1 Table) was performed with R packages and plotted as a heatmap. Alpha diversity was described using a phylogenetic diversity (PD) whole tree and a Shannon diversity index (H´) calculation and their rarefaction. A distance matrix and principal components (PC) were calculated with UniFrac to define beta diversity [46]. Moreover, a coordinate analysis (PCoA) was graphed with QIIME scripts and visualized in Emperor [47].

We classified microbiota composition according to community state type (CST) based on the predominant taxa found in the samples. A cervical CST is a cluster of community states (species composition and abundance of a cervical community) that are similar in terms of the kinds and relative abundances of the observed phylotypes [48]. The clustering of community states was carried out by means of a hierarchical clustering based on the Bray Curtis dissimilarity between all pairs of community states and an average linkage.

Statistical analysis

Relevant variables were compared between NCL vs SIL and NCL vs CC, using χ² and Kruskal-Wallis test for categorical and continuous variables, respectively. T-student tests were carried out to determine the mean difference between the Shannon diversity index or the phylogenetic diversity whole tree and the histopathological diagnosis. Logistic regression models were used to determine the association between the histopathological diagnosis and the Shannon diversity index and the histopathological diagnosis (NCL independently of HPV status and SIL/CC) and PD whole tree, adjusting by age, parity, contraceptive method and HPV-genotype.

We estimated alpha diversity risk as an odds ratio (OR) with 95% confidence interval (CI). The estimated mean difference of bits units (Shannon diversity index) in cervix between SIL or CC and NCL independently of HPV status was evaluated by a linear regression analysis adjusting by age, contraceptive method and HPV-genotype. We evaluated the variation of weighted UniFrac distances using a Kruskal–Wallis test to check the beta diversity within each histopathological diagnosis group. A Wilcoxon Mann-Whitney test was carried out to evaluate the statistical significance of weighted UniFrac distances between SIL/CC vs NCL-HPV negative.

Cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA was analyzed by histopathological diagnosis (NCL vs SIL and NCL vs CC) using the Wilcoxon Mann-Whitney test. Cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA was analyzed across CST clusters by means of the Kruskall Wallis test. Finally, we performed direct correlation analyses between the Shannon diversity index and the cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA. We performed all the statistical analyses using Stata statistical software, version 13.0 (StataCorp, Collage Station, TX, USA).

Results

General pattern of the cervical communities sampled

After QC, there were 311863 high-quality reads unevenly spread through the samples. To address this issue, 1000 random subsamples were obtained from the raw data of each sample, and they were used for further bioinformatics analysis. Alpha diversity rarefaction curves (Fig A in S1 File) show that after calculating H´in more than 400 reads, subsample diversity did not increase; therefore, 1000 sequences had sufficient sequencing depth for this study. On the unsupervised heatmap (Fig 1), samples from the NCL cluster together and independently of the HPV status, showed that HPV-negative women had a higher proportion of Lactobacillus crispatus (46%) and a smaller one of Lactobacillus iners (14.9%), whereas HPV-positive women had proportions of 13.3% and 2.1%, respectively. Interestingly, these proportions seemed to switch with the presence of HPV. Furthermore, L. crispatus was found in smaller proportions in SIL (14.4%) and CC (1.3%), whereas L. iners dropped to 2.1% in SIL and was not detected in CC.

[Figure omitted. See PDF.]

Fig 1. Community composition of cervical samples at the species level as determined by massively parallel sequencing on the 454 platform.

Unsupervised heatmap of the relative abundance of microbial taxa found in the cervical microbial communities of 29 subjects, based on the Bray Curtis dissimilarity metric. The species present in relative abundance of 1% in at least one sample are listed on the X axis. The first bar on the left side represents the treatment as follows: red–HPV-negative without lesion; blue–HPV-positive without lesion; orange–squamous intraepithelial cervical lesion; green–cervical cancer. CST are depicted in the second left barside; pink–CST III, dominated by Pseudomonas oleovorans; cyan–CST II dominated by L. iners; orange–CST I dominated by L. crispatus; green–CST IV dominated by Sneathia; blue–CST V dominated by G. vaginalis; yellow–CST VIII dominated by Fusobacterium spp.; red–CST VI dominated by S. agalactiae, and purple–CST VII dominated by Fusobacterium necrophorum. Sample names appear on the right side of the graph. The cladograms at the top of the species names indicate the approximate evolutionary relationships between the species. CST: Community State Type. Dx: Histopathological diagnosis.

https://doi.org/10.1371/journal.pone.0153274.g001

Other species of Lactobacillus, L. jensenii and L. vaginalis were found only in samples from women with NCL. Gardnerella vaginalis was found mainly in HPV-negative women with NCL (11.5%), and it decreased gradually across the HPV-positive (6.9%), SIL (8.1%) and CC (3.3%) groups. S. agalactiae represented more than 90% of the microbiota composition of two of the samples. Pseudomonas oleovorans was observed in a relative abundance of 19% only among HPV-positive women with NCL. Bacteria from the Fusobacteriales order were found only in the SIL and CC groups. In the SIL group, Fusobacterium spp. displayed a relative abundance of 6.3%; Sneathia spp., 26.6%; Shuttleworhia satelles, 8.7%; and Megasphaera elsdenii, 10.4%; whereas the relative abundance of the same microorganisms in the CC group was as follows: 14%, 12.9%, 0%, and 2.2%, respectively. Fusobacterium necrophorum was only observed in the CC group (14.2%). This points to the fact that microbiota diversity and composition are different among the analyzed groups. To confirm this, alpha and beta diversity were analyzed.

The cervical communities were classified into eight CSTs according to the dominant bacteria, as shown in Table 2. CST I isdominated by L. crispatus (21%), CST II by L. iners (17%), CST III by Pseudomonas oleovorans (10%), CST IV by Sneathia spp. (17%), CST V by G. vaginalis (7%), CST VI by Streptococcus agalactiae (7%), CST VII by F. necrophorum (7%), and CST VIII by Fusobacterium spp. (14%). Except for CST VI, all the samples were clustered according to the histopathological diagnosis. CST I was composed mainly of HPV-negative women with NCL; CST II and III were predominantly composed of HPV-positive women with NCL; CST IV was predominantly composed of SIL cases; CST V was composed mainly of women with NCL regardless of their HPV status; CST VIII was composed predominantly of CC cases, and CST VII included only cases of CC (Fig 2).

[Figure omitted. See PDF.]

Fig 2. Community compositions according to histopathological diagnosis groups.

Bar chart of relative abundance of species per group.

https://doi.org/10.1371/journal.pone.0153274.g002

[Figure omitted. See PDF.]

Table 2. Distribution of samples in each community state type (CST).

https://doi.org/10.1371/journal.pone.0153274.t002

Comparison of cervical microbiota diversity across CC stages

Alpha diversity.

Even though rarefaction curves for the Shannon index (Fig A in S1 File) showed a higher diversity in the CC and SIL groups than among the HPV-negative and -positive women with NCL, we found that diversity per se, which accounts only for richness and relative abundance, is not sufficient to determine the stage of developmentof CC. Nevertheless, when it comes to indices that consider phylogenetic metrics as a PD whole tree, we found a significant difference regarding phylogenetic diversity between HPV-negative NCL and SIL and between HPV-negative NCL and CC (p values: 0.006 and 0.036, respectively) (Table 3). In other words, the composition of the cervical microbiota is different between groups. The box plots in Fig 3 show the distribution of the Shannon diversity index and the PD whole tree across histopathological diagnosis groups. The Shannon diversity index shows an increasing trend but it is not significant. The PD whole tree shows a significant difference between HPV-negative NCL and SIL and between HPV-negative NCL and CC. SIL and CC groups display the most diverse cervical microbiota. We did not find an association between the histopathological diagnosis and the Shannon diversity index (NCL regardless of HPV status and SIL/CC), nor between the PD whole tree and the histopathological diagnosis according to the logistic regression analysis (Table 4). In assessing the association between the Shannon diversity index for the cervical microbiome samples and the histopathological diagnosis, the mean estimated difference of bit units between CC and NCL was 1.11 (95% CI, 0.057–2.165, (p = 0.04) (Table 5). This confirms that microbiota diversity in CC cases is higher than in the NCL group.

[Figure omitted. See PDF.]

Fig 3. Comparison of the Shannon diversity index and the PD whole tree logistic according to the histopathological diagnosis group.

The boxplots show the distribution of the H’ (bit units) and PD values across all the samples.

https://doi.org/10.1371/journal.pone.0153274.g003

[Figure omitted. See PDF.]

Table 3. Mean difference analysis between the Shannon diversity index or the phylogenetic diversity whole tree and the histopathological diagnosis.

https://doi.org/10.1371/journal.pone.0153274.t003

[Figure omitted. See PDF.]

Table 4. Association analysis between histopathological diagnosis and alpha diversity indexes.

https://doi.org/10.1371/journal.pone.0153274.t004

[Figure omitted. See PDF.]

Table 5. Estimated mean difference of bit units (Shannon diversity index) in cervix between NCL regardless of HPV status and SIL or CC.

https://doi.org/10.1371/journal.pone.0153274.t005

β-diversity.

The PCoA results consistently showed that the cervical microbiome is notably different in every stage of the natural history of CC (Fig 4A). The three PCs where plotted in 2D for pairwise comparison; PC1 accounted for 33.44% of the variance between samples; PC2 for 26.85% and PC3 for 10.38%. The presence of Fusobacterium spp., Sneathia spp. and Megasphaera spp. was related to PC1. The presence of Bifidobacteria spp. and Pseudomonas spp. was related to PC2. The absence of Pseudomonas spp., Fusobacterium spp. and the presence of bacteria from the Bifidobacteriaceae family were related to PC3, according to the calculated factor loading matrix (Table A in S1 File).

[Figure omitted. See PDF.]

Fig 4. Beta-diversity of microbial communities by histopathological diagnosis.

A. PCoA profile of the histopathological diagnosis displayed with weighted UniFrac distances. Each figure represents one sample colored according to its histopathological diagnosis. Red circles represent HPV-negative NCL samples; blue squares represent HPV-positive NCL samples; orange triangles represent SIL and green triangles represent CC. A1. Principal component (PC)-1 accounted for 33.44% of the variation in the composition of the microbiota due to the presence of Sneathia spp. and Fusobacterium spp. A2. PC-2 accounted for 26.85% of the variation inthe composition of the microbiota due to the presence of Bifidobacteria spp. and Pseudomonas spp. A3. PC-3 accounted for 10.38% of the variation in the composition of the microbiota due to the presence of Lactobacillus spp. and Streptococcus spp. B. B1. Variation of weighted UniFrac distances within each histological diagnosis group. B2. Variation of weighted UniFrac distances compared with -HPV-negativeNCL.

https://doi.org/10.1371/journal.pone.0153274.g004

When plotting PC1 vs PC2, we could see how the majority of the HPV-negative NCL samples clustered together in the third quadrant and none were in the first and forth quadrants, where CC and SIL samples were present. The plot of PC1 vs PC3 showed how all HPV-negatives NCL samples clustered together on the third quadrant, while 50% of the CC samples were only in the fourth quadrant, pointing to the presence of Sneathia spp. and Fusobacterium spp. and the absence of organisms from the Bifidobacteriaceae family. All SIL samples clustered in the first quadrant, pointing to the presence of Fusobacterium spp, Sneathia spp and Megasphaera spp. and a certain proportion of bacteria from the Bifidobacteriaceae family in this group.

When beta diversity was evaluated within each histological diagnosis group, the CC samples showed the highest variation within the groups (Fig 4B1, Kruskal–Wallis, p<0.0006) and the largest distance compared to HPV-negative NCL samples (Fig 4B2, U Mann-Whitney, p<0.00001).

Cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA according to histopathological diagnosis.

Cervical expression levels of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA from NCL vs SIL and NCL vs CC cases normalized to GAPDH mRNA are shown in Fig 5. Even though discrete differences can be observed between cervical expression levels of IL-4, IL-6, TGF-β1, TNF-α and IFN-γ mRNA normalized to GAPDH mRNA, only the median levels of IL-10 mRNA relative to GAPDH were higher in the SIL cases than in the NCL ones. The difference was statistically significant (p = 0.04).

[Figure omitted. See PDF.]

Fig 5. Cervical cytokine mRNA levels normalized to GAPDH, in NCL versus SIL and in NCL versus CC.

(A) IL-4, (B) IL-6, (C) IL-10, (D) TNFα, (E) IFN-γ, (F) TGF-β1. * p value for the Mann-Whitney test (p = 0.04).

https://doi.org/10.1371/journal.pone.0153274.g005

Evaluation across community state type clusters.

When cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA was evaluated across CST clusters, we found higher median cervical levels of IL-4 and TGF-β1 mRNA relative to GAPDH in CST VIII dominated by Fusobacterium spp., and the difference was statistically significant (p = 0.05 and p = 0.04, respectively) (Fig 6).To see if these findings could be explained mainly by the presence of Fusobacterium spp., the microbial composition of samples from CST I, CST IV and CST VIII were compared with their respective IL-4, TGF-β1 and INF-γ expression level (Fig 7). The highest level of IL-4 was found in sample CC2, whereas the highest level of TGF-β1 was found in sample CC8, which is composed of 87.95% Fusobacterium spp.

[Figure omitted. See PDF.]

Fig 6. Cervical cytokine mRNA levels normalized to GAPDH per community state type (CST).

(A) IL-4, (B) IL-6, (C) IL-10, (D) TNFα, (E) IFN-γ, (F) TGF-β1. *p value for Kruskal Wallis test (p≤0.05).

https://doi.org/10.1371/journal.pone.0153274.g006

[Figure omitted. See PDF.]

Fig 7. INF-γ, IL-4, TGF-β1 and IL-10 expression level per sample in community state type CST I, CST IV and CSTVIII and microbiome composition.

Cervical expression level of INF-γ (A), TGF-β1 (B), IL-4 (C) and IL-10 (D) normalized with GAPDH gen. (E) Microbiome composition in relative abundance.

https://doi.org/10.1371/journal.pone.0153274.g007

Finally, when we analyzed the direct correlation between the Shannon diversity index and the PD whole tree and the cervical expression of IL-4, IL-6, IL-10, TGF-β1, TNF-α and IFN-γ mRNA, we did not find any correlation.

Discussion

Several factors are implicated in the development of CC, including bacterial co-infection [49].The main findings of this study showed that the cervical microbiome is notably different in all stages of the natural history of CC. Surprisingly and in spite of the small sample size, we found higher median cervical levels of IL-4 and TGF-β1 mRNA in CST VIII, dominated by Fusobacterium spp.. As expected, we found high abundance of Lactobacillus spp. in the cervix of HPV-negative women without lesions, and L. crispatus was the most abundant species. In the cervix of HPV-infected women without lesions, the most abundant species was L. iners. Two of the CST detected in this study (CST I and II, dominated by L. crispatus and L. iners, respectively), which were described only in HPV-negative and -positive women with NCL, were similar to those reported by Ravel et al. and Smith et al. [50, 32]. Moreover, CST V, dominated by Gardnerella vaginalis, was found in a healthy Hispanic population in Smith’s study [32]. Likewise, Huang et al. reported similar CSTs in the cervix of their control group (non-pregnant Korean women)[30].

Our findings agree with the latest studies of HPV-infected women’s vaginal microbiome (Table 6). Previous reports on vaginal microbiota state that bacterial diversity among Chinese and Korean HPV-positive women is more complex than among healthy women [34, 50]. Longitudinal studies of vaginal microbiota have shown that L. crispatus is more beneficial than L. iners. [51]. L. iners has been found more often in women with vaginal dysbiosis by HIV, HPV and HSV-2 than L. crispatus [52]. Unlike L. iners, L. crispatus produces hydrogen peroxide. This could explain why L iners is found in a higher proportion than L. crispatus among HPV-infected women, since hydrogen peroxide has been recognized as critical in maintaining a healthy vaginal microbiota [53].

[Figure omitted. See PDF.]

Table 6. Comparison of microbiota composition between different studies and the present study.

https://doi.org/10.1371/journal.pone.0153274.t006

Sneathia spp., Megasphaera elsdenii and Shuttleworhia satelles were most representative in the SIL cases. Similarly to our findings, Fusobacteria, including Sneathia spp., were previously identified as a possible microbiological marker associated with the HPV-infection in a Korean twin cohort [54]. Bacteria of the genus Sneathia are arising as possible pathogens of the female reproductive tract. Sneathia spp., previously named Leptotrichia amnionii, can be part of the normal microbiota of the genitourinary tracts of men and women [55]. A correlation between colonization with Sneathia amnii (“L. amnionii”) and CC in HPV-positive subjects has been documented [56]. In our study, we found that Sneathia spp. was the most abundant species in the cervix of women with SIL and it was less abundant in patients with CC. Taken together, these results indicate that the presence of Sneathia spp. is a characteristic trait of patients with HPV-positive SIL.

Megasphaera elsdenii and Shuttleworthia satelles have not been previously reported in women with SIL. This was somewhat unexpected as Megasphaera type 1 and Sneathia amnionii are considered good predictors of BV [57], which is a classical dysbiosis associated with abnormalities in cervical smears [58]. Little is known about Shuttleworthia satelles. The only finding is that it has been isolated from the human oral cavity and from one patient with endocarditis and prosthetic valve infection [59].

Remarkably, our findings on the late stages of CC showed that Fusobacterium spp. is significantly more abundant than in the early stages (HPV-negative or -positive NCL),and F. necrophorum was only observed in the CC group. In agreement with our results, Fusobacteria and Atopobium vaginae were also predominant in the vaginal microbiota of Korean women with high risk of cervical intraepithelial neoplasia (CIN) [60].

Fusobacterium species are anaerobic Gram-negative bacteria that are part of the normal flora in the human mouth and the gut mucosa [61]. On the other hand, some species of Fusobacterium spp. have been recognized as opportunistic pathogens in inflammatory diseases in both mouth (periodontitis) and gut (bowel disease) [62]. Furthermore, Fusobacterium spp., have been linked to colorectal cancer [63, 64]. Molecular studies have confirmed these results and a signature of virulence genes from the microbiota has been identified [65–67].

Some studies have shown the molecular mechanisms by which Fusobacterium spp. induces its pathogenic effects. Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating the E-cadherin/β-catenin signaling via its FadA adhesin, thus modifying the tumor-immune microenvironment [68, 69]. The process is as follows: FadA binds to E-cadherin/β-catenin; E-cadherin is phosphorylated on the membrane and together with FadA, they are internalized, and beta-catenin is accumulated in the cytoplasm and translocated, and results in activation of the transcription factor NF-kB in the nuclei [68]. Moreover, fadA gene levels in the colorectal epithelium from colorectal carcinoma patients are higher than in normal patients [68]. Given that Fusobacterium spp. appear in high proportion in CC, it is possible that the FadA gene levels could be overexpressed in cervical cancer patients. Likewise, previous reports demonstrate an abnormal distribution of tumor-suppressor E-cadherin, which functions through beta-catenin, in different histological types of CC [70–72].

In terms of the aim of our study, finding out whether there is a relationship between cervical microbiota and a specific cytokine profile in women with different stages of CC, we found significant differences between CST VIII (conformed by SIL and CC cases) and high levels of cervical IL-4 and TGF-β1 mRNA. Several reports, including some from our laboratory, have shown high mRNA levels of TGF-β1, IL-4, and IL-10 in cervical biopsies of patients with SIL and CC [25, 34]. Similarly, high systemic levels of these cytokines have been reported in patients with CC [13]. It has been reported that oral infection by Fusobacterium enhances the systemic level of IL-4 [73]. Thus, a Fusobacterium infection could play a key role in the development of an immunosuppressive microenvironment characterized by anti-inflammatory cytokines (Th2 cytokine profile), such as IL-4 and TGF-β1, in HPV-transformed cells from the uterine cervix.

Likewise, it was demonstrated in a recent study that the microbiota regulates Th2 immunity through the induction of type 3 RORγτ+ Tregs and Th17, acting as a key element to balance the immune responses of the epithelial cells [74]. T-lymphocyte infiltration of the cervical mucosa increases according to the degree of CC malignancy and correlates with FoxP3+ lymphocyte infiltration in the later stages [75]. CD4+CD25+Foxp3+ Tregs play an important role in the pathogenesis of cervical carcinoma [76]. Therefore, the presence of some members of the cervical microbiota, such as Fusobacterium spp., may induce Th2 immunity through the RORγτ+ Treg and Th17 cells, present in the cervical epithelium [77]. The transdifferentiation of Th17 into Tregs cells was recently illustrated: It was demonstrated by changes in their transcriptional profile and the acquisition of potent regulatory capacity. Different transcriptional profiles of pre- and post-conversion Th17 cells revealed a role for canonical TGF-β signaling [78].

Another Th2 cytokine implicated in cervical immunosuppression in CC patients is IL-10 [24]. In this study, in spite of the presence of statistically significant high IL-10 mRNA levels in SIL, we did not find any association with CST VIII. The source of high cervical production of IL-10 by tumor cells, keratinocytes, macrophages and Langerhans cells in SIL and CC cases has been demonstrated by an immunohistochemical analysis [24, 79]. Several factors contribute to a high cervical production of IL-10 and TGF-β1 in HPV-transformed cells, including HPV proteins such as E2 [80] and HPV E6/E7 oncoproteins [81], which induce a transcriptional up-regulation of IL-10 and TGF-β1 expression [34].

On the other hand, cervical high levels mRNA of IFN-γ had been found in CC cases [25, 82]; in despite of this, evaluation of serum from patients with cervical lesions and from HPV-transformed cells, exert a strong immunosuppressive effect [25–27, 82]. Additionally, IFN-γ production favors the presence of iNKC in CIN patients and in HPV oncogene–driven hyperplasia and this play an important role in suppressing cellular immunity [83, 84]. This evidence supports the idea that several processes collaborate to maintain an immunosuppressive microenvironment and favor CC development.

The evidence reported in literature indicates that the expression levels of IFN-γ differ at local and systemic level in women with CC[13, 25]. Recently, we reported that women with minor allele homozygote genotype of the polymorphism -1615C>T of IFN-γ have a significant negative association with CC and these patients present lower serum levels than those observed in HPV-positive control patients [13]. Furthermore, the IFN-γ polymorphism studied extensively is +874T>A (rs2430561) and A/A genotype has been associated with low producer phenotype of IFN-γ [85]. Since mRNA IFN-γ levels are variable at cervical level in the natural history of CC [86], so that IFN-γ would be exercising two functions in the natural history of CC, in the initial phases exert an antiviral effect to control the infection and subsequent stages an immunosuppressive effect and in consequence a CC progression.

After everything that has been discussed, we arrive to the following possible mechanisms leading to cervical cancer progression, as illustrated in Fig 8. After a HR-HPV infection in a normal epithelium, the microbiome switches its composition from being dominated by L. crispatus to L. iners. With the progression to a SIL, the diversity of the microbiota increases with the appearance of Sneathia and other Fusobacterium spp. When CC appears, Fusobacterium necrophorum is present. This change in composition and diversity could be explained by the immunosuppressive microenvironment triggered by the viral infection, and it contributes to maintain a positive feed-back loop between the cytokine profile and the cervical microbiota.

[Figure omitted. See PDF.]

Fig 8. Suggested mechanism of microbiome changes during immunosuppression development.

The cervical epithelium is represented in each stage of CC as departing from a normal epithelium (left) and its longitudinal change when a HR-HPV infects it and it progresses to SIL and CC. Microbiome composition and diversity is depicted according to the main bacteria per stage. Lactobacillus are represented as purple rectangles, Pseudomonas oleovorans as red rectangles, Fusobacterium and Sneathia as shaped rods, Streptococcus agalactieas purple circles, and HPV as blue circles. After infection takes place, the microbiome changes and its diversity increases. HPV proteins E2, E6 and E7 enhance IL-10 expression and macrophages type 2 presence. The latter is also enhanced by TGFβ-1, which is in turn stimulated by the microbiota present. The diversity of the microbiota increases, through its toxins (FadA from Fusobacterium spp.), which disturb tight junctions and promote a metastasis similar to colon carcinoma.

https://doi.org/10.1371/journal.pone.0153274.g008

Conclusions

In conclusion, we suggest that some members of the cervical microbiota are possible modifiers of the cytokine profile of the cervical microenvironment during the development of SIL and CC. Currently, a Mexican cohort study is being carried out to support these findings. Further studies are needed to delve into the mechanism used by bacteria to promote cervix immunosuppression in CC. Evidence is accumulating and it points toward a major role of the microbiota in the immune system modulation of the female genital tract [31].

We consider that this opens new trends for understanding the role of Fusobacterium spp. in cervical carcinogenesis. Fusobacterium spp. either contributes by shifting Th1 immunity to Th2 or by a direct effect on the E-cadherin/β-catenin signaling pathway on cervical HPV-transformed cells. The findings reported here may allow, after validation, to build up new diagnostic strategies with microbiological markers for patient identification during the early stages of cervical cancer.

Supporting Information

[Figure omitted. See PDF.]

S1 File.

Fig A, Alpha diversity rarefaction curves Table A, Factor loading matrix for cervical microbiota composition.

https://doi.org/10.1371/journal.pone.0153274.s001

(PDF)

S1 Table. Proportions of each OTU present in samples.

https://doi.org/10.1371/journal.pone.0153274.s002

(XLSX)

Acknowledgments

Special thanks to Ken Dewar, from Genome Québec at McGill University, for his valuable feedback and advice that allowed us to strengthen our bioinformatics analysis. The authors wish to thank Juraj Lord and Adela Iglesias for their assistance with the manuscript. The authors wish to thank all CAPASAM members and the administrative personnel at INSP. This work was submitted in partial fulfillment of the requirements for the MSc. degree of A Audirac-Chalifour from the Master’s Program in Sciences at the School of Public Health in Mexico. K Torres Poveda is CONACyT Research Fellow-Instituto Nacional de Salud Pública (INSP), Cuernavaca, Morelos, Mexico.

Author Contributions

Conceived and designed the experiments: AACh KTP. Performed the experiments: AACh MBR JTS JMB BCC GLE KDR DC. Analyzed the data: AACh KTP JMB BCC AIBG AGC VMM. Contributed reagents/materials/analysis tools: VMM. Wrote the paper: AACh KTP JMB JTS AGC VMM. Participated in the design of the study, the collection of clinical materials, carried out the bioinformatics analysis of the sequenced samples and cytokine expression and the statistical analysis, and drafted the manuscript: AACh KTP. Participated in the processing of samples: MBR. Participated in the process of sequencing of samples: JTS. Participated in designing sequencing strategy and the analysis of the results: JMB. Carried out the bioinformatics analysis: BCC. Participated in the selection of study population, gynecological sampler and patient care: GLE, KDR DC. Participated in the design of the study and the statistical analysis: AIBG. Participated in the analysis of the results and drafted the manuscript: AGC VMM. Contributed with reagents and materials for experiments: VMM. Read and approved the final manuscript: AACh KTP MBR JTS JMB BCC GLE KDR AIBG DC AGC VMM.

Citation: Audirac-Chalifour A, Torres-Poveda K, Bahena-Román M, Téllez-Sosa J, Martínez-Barnetche J, Cortina-Ceballos B, et al. (2016) Cervical Microbiome and Cytokine Profile at Various Stages of Cervical Cancer: A Pilot Study. PLoS ONE 11(4): e0153274. https://doi.org/10.1371/journal.pone.0153274

References

1. Global Burden of Disease Cancer Collaboration, Fitzmaurice C, Dicker D, Pain A, Hamavid H, Moradi-Lakeh M, et al. The GlobalBurden of Cancer2013. JAMA Oncol 2015;1(4):505–27. pmid:26181261

2. INEGI 2012. Estadísticas de mortalidad. Cubos dinámicos y CONAPO 2012. Proyecciones de la población de México 2012–2050.

3. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med 2003, 348:518–27. pmid:12571259

4. Gillet E, Meys JF, Verstraelen H, Bosire C, De Sutter P, Temmerman M, et al. Bacterial vaginosis is associated with uterine cervical human papillomavirus infection: a meta-analysis. BMC Infect Dis 2011, 11:10. pmid:21223574

5. Guo YL, You K, Qiao J, Zhao YM, Geng L. Bacterial vaginosis is conducive to the persistence of HPV infection. Int J STD AIDS 2012;23(8):581–4. pmid:22930296

6. Vriend HJ, Bogaards JA, van Bergen JE, Brink AA, van den Broek IV, Hoebe CJ, et al. Incidence and persistence of carcinogenic genital human papillomavirus infections in young women with or without Chlamydia trachomatis co-infection. Cancer Med 2015 Jul 21. (Epub ahead of print)

7. Clarke MA, Rodriguez AC, Gage JC, Herrero R, Hildesheim A, Wacholder S, et al. Burk R. A large, population-based study of age-related associations between vaginal pH and human papillomavirus infection. BMC Infect Dis 2012;12:33. pmid:22316377

8. Insinga RP, Perez G, Wheeler CM, Koutsky LA, Garland SM, Leodolter S, et al. Incident cervical HPV infections in young women: transition probabilities for CIN and infection clearance. Cancer Epidemiol Biomarkers Prev 2011; 20, 287–296. pmid:21300618

9. Almonte M, Albero G, Molano M, Carcamo C, García PJ, Pérez G. Risk factors for human papillomavirus exposure, and co-factors for cervical cancer in Latin America and the Caribbean. Vaccine 2008;19;26 (Suppl 11):L16–36. pmid:18945400

10. Vaccarella S, Herrero R, Dai M, Snijders PJ, Meijer CJ, Thomas JO, et al. Reproductive factors, oral contraceptive use, and human papillomavirus infection: pooled analysis of the IARC HPV prevalence surveys. Cancer Epidemiol Biomarkers Prev 2006;15(11):2148–2153. pmid:17119039

11. Louie KS, Castellsague X, de Sanjose S, Herrero R, Meijer CJ, Shah K, et al. Smoking and passive smoking in cervical cancer risk: pooled analysis of couples from the IARC multicentric case-control studies. Cancer Epidemiol Biomarkers Prev 2011;20(7):1379–90. pmid:21610224

12. Torres-Poveda K, Burguete-García AI, Cruz M, Martínez-Nava GA, Bahena-Román M, Ortíz-Flores E, et al. The SNP at -592 of human IL-10 gene is associated with serum IL-10 levels and increased risk for human papillomavirus cervical lesion development. Infect Agent Cancer 2012;7(1):32–40. pmid:23148667

13. Torres-Poveda K, Burguete-García AI, Bahena-Román M, Méndez-Martínez R, Zurita A, López-Estrada G, et al. Risk Allelic Load in Th2 and Th3 Cytokines gene as Biomarker of Susceptibility to cervical cancer. Submitted BMC Cancer, 2015

14. Zhang X, Zhang L, Tian C, Yang L, Wang Z. Genetic variants and risk of cervical cancer: epidemiological evidence, meta-analysis and research review. BJOG 2014; 121(6):664–73. Epub 2014 Feb 19. pmid:24548744

15. Guijon F, Paraskevas M, Rand F, Heywood E, Brunham R, McNicol P. Vaginal microbial flora as a cofactor in the pathogenesis of uterine cervical intraepithelial neoplasia. Int J Gynecol Obstet 1992;37:185–191.

16. Lehtinen M, Ault KA, Lyytikainen E, Dillner J, Garland SM, Ferris DG, et al. Chlamydia trachomatis infection and risk of cervical intraepithelial neoplasia. Sex Transm Infect 2011;87(5):372–6. pmid:21471141

17. Sears CL, Garrett WS. Microbes, microbiota, and colon cancer. Cell Host Microbe 2014; 15(3):317–28. pmid:24629338

18. Dejea CM, Wick EC, Hechenbleikner EM, White JR, Mark Welch JL, Rossetti BJ, et al. Microbiota organization is a distinct feature of proximal colorectal cancers. Proc Natl Acad Sci U S A 2014; 23;111(51):18321–6. pmid:25489084

19. Chase D, Goulder A, Zenhausern F, Monk B, Herbst-Kralovetz M. The vaginal and gastrointestinal microbiomes in gynecologic cancers: A review of applications in etiology, symptoms and treatment. Gynecol Oncol 2015;138(1):190–200. pmid:25957158

20. Ito Y, Kishishita M, Yanase S. Induction of Epstein-Barr virus antigens in human lymphoblastoid P3HR-1 cells with culture fluid of Fusobacterium nucleatum. Cancer Res 1980;40(11):4329–30. pmid:6258779

21. Han YW. Fusobacterium nucleatum: a commensal-turned pathogen.CurrOpinMicrobiol2015;23:141–7.

22. Linhares IM, Summers PR, Larsen B, Giraldo PC, Witkin SS. Contemporary perspectives on vaginal pH and lactobacilli. Am J Obstet Gynecol 2011, 204:120.e1–5.

23. Kim TK, Thomas SM, Ho M, Sharma S, Reich CI, Frank JA, et al. Heterogeneity of vaginal microbial communities within individuals. J Clin Microbiol2009;47(4):1181–9. pmid:19158255

24. Bermúdez-Morales VH, Burguete AI, Gutierrez ML, Alcocer-González JM, Madrid-Marina V. Correlation between IL-10 expression and Human Papillomavirus infection in cervical cancer. A Mechanism for Immune Response Escape. Cancer Investig 2008;26:1037–43.

25. Alcocer-González JM, Berumen J, Tamez-Guerra R, Gariglio P, Hernández-Pando R, Bermúdez-Morales VH, et al. In Vivo Expression of Immunosuppressive Cytokines in Human Papillomavirus-Transformed Cancer Cells. Viral Immunology 2006;19(3):481–91. pmid:16987066

26. Giannini SL, Al-Saleh W, Piron H, Jacobs N, Doyen J, Boniver J, et al. Cytokine expression in squamous intraepithelial lesions of the uterine cervix: implications for the generation of local immunosuppression. Clin Exp Immunol. 1998;113(2):183–9. pmid:9717966

27. Sheu BC, Lin RH, Lien HC, Ho HN, Hsu SM, Huang SC. Predominant Th2/Tc2 polarity of tumor-infiltrating lymphocytes in human cervical cancer. J Immunol. 2001;167(5):2972–8. pmid:11509647

28. Spear GT, Sikaroodi M, Zariffard MR, Landay AL, French AL, Gillevet PM. Comparison of the Diversity of the Vaginal Microbiota in HIV-Infected & HIV-Uninfected women with or without Bacterial Vaginosis. Journal of Infectious Diseases 2008;198:1131–40. pmid:18717638

29. Huang YE, Wang Y, He Y, Ji Y, Wang LP, Sheng HF, et al. Homogeneity of the vaginal microbiome at the cervix, posterior fornix, and vaginal canal in pregnant Chinese women.MicrobEcol2015;69(2):407–14. Epub 2014 Sep 18.

30. Behbakht K, Friedman J, Heimler I, Aroutcheva A, Simoes J, Faro S. Role of the vaginal microbiological ecosystem and cytokine profile in the promotion of cervical dysplasia: A case-control study. Infect Dis Obstet Gynecol 2002;10:181–6. pmid:12648311

31. Gao W, Weng J, Gao Y, Chen X. Comparison of the vaginalmicrobiota diversity of women with and without human papillomavirus infection: a cross-sectional study. BMC Infect Dis 2013;13:271. pmid:23758857

32. Smith BC, McAndrew T, Chen Z, Harari A, Barris DM, Viswanathan S, et al. The cervical microbiome over 7 years and a comparison of methodologies for its characterization. PLoS One 2012;7(7):e40425. Epub 2012 Jul 9. pmid:22792313

33. Anahtar MN, Byrne EH, Doherty KE, Bowman BA, Yamamoto HS, Soumillon M et al. Cervicovaginal Bacteria Are a Major Modulator of Host Inflammatory Responses in the Female Genital Tract, Immunity 2015;42(5):965–76. pmid:25992865

34. Torres-Poveda K, Bahena-Román M, Madrid-González C, Burguete-García AI, Bermúdez-Morales VH, Peralta-Zaragoza O, et al. Role of IL-10 and TGF-β1 in local immunosuppression in HPV-associated cervical neoplasia. World J Clinical Oncology 2014; 5(4):753–63.

35. Bernard HU, Chan SY, Manos MM, Ong CK, Villa LL, Delius H, et al. Identification and assessment of known and novel human papillomaviruses by polymerase chain reaction amplification, restriction fragment length polymorphisms, nucleotide sequence, and phylogenetic algorithms. J Infect Dis 1994 70:1077–85. J Infect Dis 70:1077–85.

36. Yoshikawa H, Kawana T, Kitagawa K, Mizuno M, Yoshikura H, Iwamoto A. Detection and typing of multiple genital human papillomaviruses by DNA amplification with consensus primers. Jpn J Cancer Res 1991;82:524–31. pmid:1648051

37. de Roda Husman AM, Walboomers JM, van den Brule AJ, Meijer CJ, Snijders PJ. The use general primers GP5 and GP6 elongated at their 3´ end with adjacent highly conserved sequences improves human papillomavirus detection by PCR. J Gen Virol 1995; 76:1057–62. pmid:9049358

38. Lee D-H, Hwang NR, Lim MC, Yoo CW, Joo J, Kim J-Y, et al. Comparison of the performance of Anyplex II HPV HR, the Cobas 4800 human papillomavirus test and Hybrid Capture 2. Annals of Clinical Biochemistry: An international journal of biochemistry and laboratory medicine [Internet]. 2015 Oct 20; Available from: http://acb.sagepub.com/content/early/2016/01/23/0004563215614036.abstract

39. El-Sherif AM, Seth R, Tighe PJ, Jenkins D. Quantitative analysis of IL-10 and IFN-c mRNA levels in normal cervix and human papillomavirus type 16 associated cervical precancer. J Pathol 2001; 195:179–85. pmid:11592096

40. Nossa CW, Oberdorf WE, Yang L, Aas JA, Paster BJ, Desantis TZ, et al. Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome. World J Gastroenterol2010; 7;16(33):4135–44002E. pmid:20806429

41. Girke T. FASTQ Quality Plots, R; A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL 05-Dec-13 R Core Team. (2013)

42. 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(5):335–6. Epub 2010 Apr 11. pmid:20383131

43. Caporaso JG, Bittinger K, Bushman FD, DeSantis TZ, Andersen GL, Knight R. PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 2010; 26:266–67. pmid:19914921

44. Edgar R. C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010; 26(19), 2460–61. Epub 2010 Aug 12 pmid:20709691

45. De Santis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microb 2006; 72(7): 5069–72.

46. Lozupone C, Knight R. UniFrac: A New Phylogenetic Method for Comparing Microbial Communities. Appl Environ Microbiol2005; 71(12): 8228–35. pmid:16332807

47. Vazquez-Baeza Y, Pirrung M, Gonzalez A, Knight R. Emperor: A tool for visualizing high-throughput microbial community data. Gigascience 2013; 2(1):16. pmid:24280061

48. Gajer P, Brotman RM, Bai G, Sakamoto J, Schutte UM, Zhong X, et al. Temporal dynamics of the human vaginal microbiota. Sci Transl Med 2012; 4:132ra152.

49. Shew ML, Ermel AC, Weaver BA, Tong Y, Tu W, Kester LM, et al. Association of Chlamydia trachomatis infection with redetection of human papillomavirus after apparent clearance. J Infect Dis 2013; 208(9): 1416–21. pmid:23911713

50. Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SSK, McCulle SL, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A 2011; 108 (Suppl. 1): 4680–87. Epub 2010 Jun 3. pmid:20534435

51. van de Wijgert JH, Borgdorff H, Verhelst R, Crucitti T, Francis S, Verstraelen H, et al. The vaginal microbiota: what have we learned after a decade of molecular characterization? PLoS One 2014;9(8):e105998. pmid:25148517

52. Borgdorff H, Tsivtsivadze E, Verhelst R, Marzorati M, Jurriaans S, Ndayisaba GF, et al. Lactobacillus-dominated cervicovaginal microbiota associated with reduced HIV/STI prevalence and genital HIV viral load in African women.ISME J 2014;8(9):1781–93. Epub 2014 Mar 6. pmid:24599071

53. Dols JA, Reid G, Kort R, Schuren FH, Tempelman H, Bontekoe TR, et al. PCR-based identification of eight Lactobacillus species and 18hr-HPVgenotypes in fixed cervical samples of South African women at risk of HIV and BV. DiagnCytopathol2012;40(6):472–7. Epub 2011 Oct 24.

54. Lee JE, Lee S, Lee H, Song YM, Lee K, Han MJ, et al. Association of the Vaginal Microbiota with Human Papillomavirus Infection in a Korean Twin Cohort. PLoS One 2013; 8(5): e63514. pmid:23717441

55. Harwich MD Jr, Serrano MG, Fettweis JM, Alves JM, Reimers MA; Vaginal Microbiome Consortium (additional members), et al. Genomic sequence analysis and characterization of Sneathia amnii sp. nov.BMC Genomics 2012;13Suppl 8:S4. Epub 2012 Dec 17. pmid:23281612

56. Nawrot R, Kamieniarz K, Malinowska M, Jozefiak A, Kedzia W, Kwasniewska A, et al. The prevalence of Leptotrichia amnionii in cervical swabs of HPV positive and negative women. Eur J Gynaecol Oncol 2010; 31(4):425–28. pmid:20882886

57. Shipitsyna E, Roos A, Datcu R, Hallén A, Fredlund H, Jensen JS, et al. Composition of the vaginal microbiota in women of reproductive age—sensitive and specific molecular diagnosis of bacterial vaginosis is possible? PLoS One 2013;8(4):e60670. pmid:23585843

58. Caixeta RC, Ribeiro AA, Segatti KD, Saddi VA, Figueiredo Alves RR, Dos Santos Carneiro MA, et al. Association between the human papillomavirus, bacterial vaginosis and cervicitis and the detection of abnormalities in cervical smears from teenage girls and young women. DiagnCytopathol2015 Jul 14.

59. Shah NB, Suri RM, Melduni RM, Patel R, Kommedal O, Saebø O, et al. Shuttleworthia satelles endocarditis: evidence of non-dental human disease. J Infect 2010;60(6):491–3. Epub 2010 Feb 17. pmid:20170673

60. Oh HY, Kim BS, Seo SS, Kong JS, Lee JK, Park SY, et al. The association of uterine cervical microbiota with an increased risk for cervical intraepithelial neoplasia in Korea. Clin Microbiol Infect 2015;21(7):674.e1–9. Epub 2015 Mar 6.

61. Sagata N, Haake SK, Mannon P, Lemon KP, Waldron L, Gevers D, et al. Composition of the adult digestive tract bacterial microbiome based in seven mouth surfaces tonsil, throat and stool samples. Genome Biol 2012; 13(6):R42. pmid:22698087

62. Signat B, Roques C, Poulet P, Duffaut D. Fusobacterium nucleatum in periodontal health and disease. Curr Issues MolBiol2011;13(2):25–36.

63. 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(2):292–8. Epub 2011 Oct 18. pmid:22009990

64. Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, et al. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res 2012;22(2):299–306. Epub 2011 Oct 18. pmid:22009989

65. Ahn J, Sinha R, Pei Z, Dominianni C, Wu J, Shi J, et al. Human gut microbiome and risk for colorectal cancer. J Natl Cancer Inst 2013;105(24):1907–11. Epub 2013 Dec 6. pmid:24316595

66. Tahara T, Yamamoto E, Suzuki H, Maruyama R, Chung W, Garriga J, et al. Fusobacterium in colonic flora and molecular features of colorectal carcinoma.Cancer Res 2014;74(5):1311–8. Epub 2014 Jan 2. pmid:24385213

67. Burns MB, Lynch J, Starr TK, Knights D, Blekhman R. Virulence genes are a signature of the microbiome in the colorectal tumor microenvironment. Genome Med 2015;7(1):55. eCollection 2015. pmid:26170900

68. Rubinstein MR, Wang X, Liu W, Hao Y, Cai G, Han YW. Fusobacterium nucleatum promotes colorectal carcinogenesis by modulating E-cadherin/β-catenin signaling via its FadA adhesin. Cell Host Microbe 2013;14(2):195–206. pmid:23954158

69. Kostic AD, Chun E, Robertson L, Glickman JN, Gallini CA, Michaud M, et al. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe 2013;14(2):207–15. pmid:23954159

70. Rodríguez-Sastre MA, González-Maya L, Delgado R, Lizano M, Tsubaki G, Mohar A, et al. Abnormal distribution of E-cadherin and beta-catenin in different histologic types of cancer of the uterine cervix. Gynecol Oncol 2005;97(2):330–6. pmid:15863126

71. Branca M, Giorgi C, Ciotti M, Santini D, Di Bonito L, Costa S, et al. Down-regulation of E-cadherin is closely associated with progression of cervical intraepithelial neoplasia (CIN), but not with high-risk human papillomavirus (HPV) or disease outcome in cervical cancer. Eur J Gynaecol Oncol 2006;27(3):215–23. pmid:16800245

72. Auvinen E, Carpen O, Korpela T, Ronty M, Vaheri A, Tarkkanen J. Altered expression of ezrin, E-Cadherin and β-Catenin in cervical neoplasia. Neoplasma 2013;60(1):56–61. pmid:23067217

73. Velsko IM, Chukkapalli SS, Rivera-Kweh MF, Chen H, Zheng D, Bhattacharyya I, et al. Fusobacterium nucleatum Alters Atherosclerosis Risk Factors and Enhances Inflammatory Markers with an Atheroprotective Immune Response in ApoEnull Mice. PLoS One 2015; 10(6):e0129795. pmid:26079509

74. Ohnmacht C, Park JH, Cording S, Wing JB, Atarashi K, Obata Y, et al. The microbiota regulates type 2 immunity through RORγt+ T cells. Science 2015 Jul 9. pii: aac4263.

75. Feng Q, Wei H, Morihara J, Stern J, Yu M, Kiviat N, et al. Th2 type inflammation promotes the gradual progression of HPV-infected cervical cells to cervical carcinoma. Gynecol Oncol 2012; 127(2):412–9. pmid:22828962

76. Chen ZF, Xu Q, Ding JB, Zhang Y, Du R, Ding Y. CD4+CD25+Foxp3+ Treg and TGF-beta play important roles in pathogenesis of Uygur cervical carcinoma.Eur J Gynaecol Oncol 2012;33(5):502–7. pmid:23185797

77. Punt S, van Vliet ME, Spaans VM, de Kroon CD, Fleuren GJ, Gorter A, et al. FoxP3(+) and IL-17(+) cells are correlated with improved prognosis in cervical adenocarcinoma. Cancer Immunol Immunother 2015;64(6):745–53. pmid:25795131

78. Gagliani N, Vesely MC, Iseppon A, Brockmann L, Xu H, Palm NW, et al. Th17 cells transdifferentiate into regulatory T cells during resolution of inflammation. Nature 2015; 523(7559):221–5. Epub 2015 Apr 29. pmid:25924064

79. Prata TT, Bonin CM, Ferreira AM, Padovani CT, Fernandes CE, Machado AP, et al. Local immunosuppression induced by high viral load of human papillomavirus: characterization of cellular phenotypes producing interleukin-10 in cervical neoplastic lesions. Immunology 2015 Jun 8.

80. Bermúdez-Morales VH, Peralta-Zaragoza O, Moreno J, Alcocer-González JM, Madrid-Marina V. IL-10 expression is regulated by HPV E2 protein in cervical cancer cells. Molecular Medicine Reports 2011;4(2):369–75. pmid:21468579

81. Peralta-Zaragoza O, Bermúdez-Morales V, Gutiérrez-Xicotencatl L, Alcocer-González J, Recillas-Targa F, Madrid-Marina V. Human papilomavirus-16 E6 and E7 proteins induce activation of human TGF-1 basal promoter in epithelial cells throughout a Sp1 recognition sequence. Viral Immunology 2006;19(3):468–80.

82. Shondel SM, Helm CW, Gercel-Taylor C, Taylor DD. Differential expression of T-cell CD3-zeta chains in patients with cervical dysplasia before and after treatment. Int. J. Gynecol. Cancer.2007;17(6):1278–82. pmid:17997796

83. Hu T, Yang P, Zhu H, Chen X, Xie X, Yang M, et al. Accumulation of invariant NKT cells with increased IFN-γ production in persistent high-risk HPV-infected high-grade cervical intraepithelial neoplasia. DiagnPathol2015;10:20.

84. Gosmann C, Frazer IH, Mattarollo SR, Blumenthal A. IL-18, but not IL-12, induces production of IFN-γ in the immunosuppressive environment of HPV16 E7 transgenic hyperplastic skin. J Invest Dermatol 2014;134(10):2562–9. pmid:24756108

85. Marangon AV, Guelsin GA, Visentainer JE, Borelli SD, Watanabe MA, Consolaro ME, et al. The association of the immune response genes to human papillomavirus-related cervical disease in a Brazilian population. Biomed Res Int2013;2013:146079. pmid:23936772

86. Yang W, Song Y, Lu YL, Sun JZ, Wang HW.Increased expression of programmed death (PD)-1 and its ligand PD-L1 correlates with impaired cell-mediated immunity in high-risk human papillomavirus-related cervical intraepithelial neoplasia. Immunology 2013;139(4):513–22. pmid:23521696