Vaginal microbiome community state types and high-risk human papillomaviruses in cervical precancer and cancer in North-central Nigeria
BMC Cancer volume 23, Article number: 683 (2023)
High risk human papillomaviruses (HR-HPV) have a causal role in cervical oncogenesis, and HIV-mediated immune suppression allows HR-HPV to persist. We studied whether vaginal microbiome community state types (CSTs) are associated with high-grade precancer and/or invasive cervical cancer (HSIL/ICC).
This was a cross-sectional study of adult women with cervical cancer screening (CCS) at the Jos University Teaching Hospital (JUTH) in Jos, Nigeria, between January 2020 and February 2022. Cervical swabs underwent HPV genotyping (Anyplex™ II HPV28). Cervico-vaginal lavage (CVL) sample was collected for 16 S rRNA gene amplicon sequencing. We used multivariable logistic regression modelling to assess associations between CSTs and other factors associated with HSIL/ICC.
We enrolled 155 eligible participants, 151 with microbiome data for this analysis. Women were median age 52 (IQR:43–58), 47.7% HIV positive, and 58.1% with HSIL/ICC. Of the 138 with HPV data, 40.6% were negative for HPV, 10.1% had low-risk HPV, 26.8% had single HR-HPV, and 22.5% had multiple HR-HPV types. The overall prevalence of any HR-HPV type (single and multiple) was 49.3%, with a higher proportion in women with HSIL/ICC (NILM 31.6%, LSIL 46.5%, HSIL 40.8%, and 81.5% ICC; p = 0.007). Women with HIV were more likely to have HSIL/ICC (70.3% vs. 29.7% among women without HIV). In crude and multivariable analysis CST was not associated with cervical pathology (CST-III aOR = 1.13, CST-IV aOR = 1.31). However, in the presence of HR-HPV CST-III (aOR = 6.7) and CST-IV (aOR = 3.6) showed positive association with HSIL/ICC.
Vaginal microbiome CSTs were not significantly associated with HSIL/ICC. Our findings suggest however, that CST could be helpful in identifying women with HSIL/ICC and particularly those with HR-HPV. Characterization of CSTs using point-of-care molecular testing in women with HR-HPV should be studied as an approach to improve early detection and cervical cancer prevention. Future longitudinal research will improve our understanding of the temporal effect of non-optimal CST, HR-HPV, and other factors in cervical cancer development, prevention, and control.
Approximately 600,000 new cases of invasive cervical cancer (ICC) are diagnosed annually, and over 50% of these women (311,000) die within one year of diagnosis.(1) The high burden of HIV in sub-Saharan Africa and evidence that women living with HIV have 6 times greater risk of developing ICC compared to women without HIV, calls for an urgent focus on prevention and control of cervical cancer in these high-risk populations . Furthermore, epidemiologic studies have confirmed that ICC is primarily caused by human papillomaviruses (HR-HPV)  whose persistence and oncogenic potentials are enhanced by HIV-mediated decreases in cellular immunity with other co-factors playing roles in cervical carcinogenesis [3,4,5,6,7,8,9,10,11,12,13].
To improve cervical cancer prevention and control, the WHO launched the cervical cancer elimination initiative in 2018, recommending HPV vaccination, early detection of precancer through screening with a high-performing HPV test, and treatment of precancer and early invasive cancer . The HIV and HR-HPV are sexually transmissible viruses and their attribution to cervical carcinogenesis has led to experts ascribing cervical cancer as a sexually transmissible cancer,  and previous study reports have shown that cervico-vaginal microbiome play a significant role in the process [10, 16,17,18]. Additionally, natural history studies of cervical carcinogenesis have expanded our understanding of the potential mechanism by which vaginal microbiome alone or in conjunction with local pro inflammatory cytokines act either directly or mediated through HR-HPV persistence in advancing progression of cellular changes in a normal cervix (NILM), low-grade precancer (LSIL) to high-grade precancer (HSIL), and eventually invasive cervical cancer stage [19,20,21,22,23,24].
Studies of vaginal microbiome using the community state type (CST) classification system  have observed that four CSTs are dominated by bacteria from the acid producing genus Lactobacillus: CST-1 (L. crispatus), CST-II (L. gasseri), CST-III (L. iners), CST-V (L. jensenii), while the CST-IV is characterized by diverse groups of anaerobes such as bacteria from the genus Gardnerella and other taxa associated with non-optimal cervico-vaginal health . The CSTs have been shown to correlate with severity of cervical precancer and those characterized by depletion of the protective Lactobacillus crispatus found in women with high-grades precancer and invasive cervical cancer .
Most of these studies have been done in women outside of Nigeria and West Africa. A study examining persistence of low and HR-HPV in HIV positive and HIV negative women in Nigeria found a prevalence of 9% single high-risk HPV in HIV negative women with 4% persistence at 6 months visits, while the baseline prevalence of single HR-HPV was 17% in HIV positive women at baseline with 14% persistence at 6 months . The findings implied that HIV infection was associated with high baseline prevalence of both low- and high-risk HPV, with persistent infection seen in high-risk but not low-risk HPV . One study in Nigeria examined the association of vaginal microbiome and high-risk HPV and found a moderate association between prevalent HR-HPV and a low relative abundance of Lactobacillus sp, with a significant abundance of Prevotella and Leptotrichia in HIV negative women with HR-HPV infection . These findings contribute to the understanding of HR-HPV burden and persistence, and the possible role of vaginal microbiome, but did not relate their findings to the cervical pathology.
We sought to understand the predominant CSTs in Nigerian women and evaluate associations with cervical precancer and invasive cervical cancer, HIV, and HPV infection. We hypothesized that women with high-grade cervical precancer and/or invasive cervical cancer will have more of the non-optimal CSTs and high prevalence of HR-HPV types.
Study setting and design
This was a cross-sectional study of adult women attending the cervical cancer screening and gynecologic oncology unit of the department of obstetrics and gynecology, Jos Teaching Hospital in Jos, North Central Nigeria. Participants were enrolled between January 2020 and February 2022.
Study participants and data collection procedures
Screening and enrollment of eligible study participants were done at the cervical cancer screening and colposcopy clinic of the Department of Obstetrics and Gynecology of the Jos University Teaching Hospital, Jos, Nigeria. The source population were adult women who presented for cervical cancer screening, follow-up colposcopy for abnormal cytology report, or evaluation for suspected ICC. Women who were at least 30 years old, with no prior treatment for abnormal cervical precancer or invasive cancer and an intact uterus and cervix, and who were not pregnant were eligible for inclusion. Enrollment was stratified by cervical cytology. We enrolled women with a spectrum of negative intraepithelial neoplasia or malignancy (NILM), low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), and women with confirmed invasive cervical cancer (ICC). After the principal investigator explained the study purpose and procedures and obtained written informed consent, each participant had a gynecological examination in the colposcopy room. A sterile dry metal speculum was gently inserted into the vagina to visualize the cervix and the vaginal fornices. Hybrid capture 2 Digene HPV collection brushes were used to obtain samples from cervixes, and cytobrushes were used to obtained samples for cytology from each patient. Subsequently, cervico-vaginal lavage samples were collected as described below. The cervical pathological outcomes (NILM, LSIL and HSIL) were categorized and reported according to the Bethesda 2001 system , and ICC was diagnosed by histopathology. Details of colposcopic assessments have been described previous reports from Jos Nigeria . Confirmation of HIV status was done in accordance with the HIV Rapid testing algorithm in Nigeria . In brief, HIV testing was done using the national algorithm which included the use of Determine HIV-1/2 as test 1 (T1) (Abbott, California, USA), Unigold HIV-1/2 as test 2 (T2) (Trinity Biotech Plc., Ireland), and StatPak HIV-1/2 as the tie-breaker test (T3) (Chembio Diagnostic Systems, Inc., New York, USA). Individuals that were reactive with T1 and either T2 or T3 were considered HIV-positive .
Laboratory procedures for processing, preservation, and microbiome DNA extraction
Sample collection and processing
All sample collections were performed by the gynecological oncology team in the colposcopy assessment room. Cervical vaginal Lavage (CVL) sample was obtained by washing the upper vaginal walls and cervix of each of the female subjects with 10 mL of sterile Phosphate Buffer Saline applied and aspirated with a sterile disposable pipette and Drummond Pipet-Aid. The CVL sample was aseptically aspirated into 15mL sterile falcon tubes placed on ice, capped, and transported to the genomics laboratory for processing and storage within 1–4 h of collection. The 10 mL of CVL sample in the falcon tube was vortexed and aliquoted into 2mL cryovials tubes. The aliquoted 2 mL CVL sample was centrifuged at 10,000 x g for 10 min and 1 mL of the resulting supernatant was pipetted and discarded while the pellet was resuspended in 1 mL of DNA/RNA shield and stored at -800 C.
Extraction and purification of DNA
The extraction and purification of the DNA from the CVL sample were performed using Quick-DNA™ Fecal/Soil Microbe Miniprep kit obtained from Zymo Research, USA, following manufacturer directions.
Briefly, the CVL was removed from − 800 C storage and allowed to thaw before it was resuspended by vortexing. 400 µl CVL was used for DNA extraction using the Zymo quick DNA fecal/soil microbe DNA Miniprep kit. Cat No: D6010, Lot No: 207,380, 207,953 &208,889, according to the manufacturer’s protocol. The extracted DNA was quantified using a Qubit dsDNA HS Assay Kit with a Qubit 4 Fluorometer (Thermofisher, Scientific, USA), and stored at -80 °C until use.
HPV DNA extraction
As described above the CVL samples were thawed, and 500 µl was pelleted by centrifugation at 3,000 g for 10 min and re-suspended in 200 µl phosphate-buffered saline. DNA was extracted using a QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), following the manufacturer’s protocol. Each DNA extraction run contained positive and negative controls to monitor the extraction procedure and the extracted DNA was eluted in a final volume of 50 µl. The extracted DNA was quantified using a Qubit 4 Fluorometer.
Detection of HR-HPV genotypes using the real-time polymerase chain reaction (qPCR)
For each DNA sample, HPV detection and genotyping were performed using Anyplex™ II HPV28 Detection kit (Seegene, Korea). Briefly, each PCR reaction was performed in a 20-µl reaction volume consisting of 5 µl aliquot of DNA added to Anyplex™ PCR Mix (15 µl aliquot each for mixtures A and B) on a CFX96 Real-time PCR system (Bio-Rad Laboratories, Inc., Hercules, CA, USA) according to the manufacturer’s instructions. The thermal cycler conditions consisted of an initial incubation at 50 °C for 4 min, denaturation at 95 °C for 15 min, followed by 50 cycles of denaturation (30 s at 95 °C), annealing (1 min at 60 °C), and elongation (30 s at 72 °C). Cyclic-Catcher Melting Temperature Analysis (CMTA) was performed after PCR cycles 30, 40, and 50. CMTA was performed by cooling the reaction mixture to 55 °C, holding at 55 °C for 30 s, and heating from 55 to 85 °C (5 s/0.5 °C) with continuous fluorescent monitoring. The L1 gene of HPV DNA was the target of the assay, together with simultaneous targeting of a housekeeping gene (Human beta-globin) which was co-amplified as an internal control to monitor DNA purification efficiency, PCR inhibition, and cell adequacy. The results were exported and analyzed using the Seegene Viewer software provided by the manufacturer .
Microbial community characterization
A two-stage PCR protocol was used on the extracted microbiome DNA to amplify the V3-V4 variable region of bacterial 16S rRNA genes, as described previously . Briefly, gDNA was used as template for PCR amplification with the primers CS1_357wF  and CS2_806R  (ACACTGACGACATGGTTCTACACCTACGGGNGGCWGCAG and TACGGTAGCAGAGACTTGGTCTGGACTACNVGGGTWTCTAAT, respectively; linker sequences underlined). Reactions were performed in 10 µl volumes using repliQa HiFi ToughMix (QuantaBio). Cycling conditions were 2 min denaturation at 98°C, followed by 28 cycles of 98°C for 10 sec, 50°C for 1 sec, and 68°C for 1 sec. Subsequently, a second PCR amplification was performed, also in 10 microliter reactions in 96-well plates using repliQa HiFi ToughMix. Each well received a separate primer pair with a unique 10-base barcode, obtained from the Access Array Barcode Library for Illumina (Fluidigm, South San Francisco, CA; Item# 100–4876). One microliter of PCR product from the first stage amplification was used as template for the 2nd stage, without cleanup. Cycling conditions were 98°C for 2 minutes, followed by 8 cycles of 98°C for 10”, 60 °C for 1” and 68 °C for 1”. Libraries were then pooled and sequenced with a 15% phiX spike-in on an Illumina MiSeq sequencer employing V3 chemistry (2 × 300 base paired-end reads). Library preparation and sequencing were performed at the Genomics and Microbiome Core Facility (GMCF) at Rush University.
Quality control and taxonomic identification and classification into community state types (CSTs) was conducted by University of Maryland Institute for Genomic Science, as described previously . . CSTs were determined by the VALENCIA algorithm implemented by the University of Maryland Institute for Genomic Science, which uses a distance-based metric to classify each sample to a CST based on the similarity of the sample to the centroid of CSTs identified in a reference set. In brief, dominant taxa by CST are as follows: CST-I, L. crispatus; CST-II, L. gasseri; CST-III, L. iners; CST-IV, diverse; CST-V, L. jensenii. .
All data were entered in REDCap and exported as CSV file for coding and subsequent analysis on Stata version 14.1, College station, Texas, USA. The setting and use of REDCap for managing our research data has been reported in an earlier publication .
Stratification of high-risk HPV status
The Seegene AnyplexTM II HPV28 Detection kit detects 28 human papillomavirus types .
High-risk HPV (HR-HPV) types include: HPV16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 69, 73, and 82. The low-risk types HPV (LR-HPV) include: HPV6, 11, 40, 42, 43, 44, 54, 61, and 70. The HPV status were categorized as follows: category 1 as HPV negative, category 2 as low-risk risk HPV, category 3 as single high-risk HPV (any of the HR-HPV), and category 4 as multiple high risk (2 or more of the HR-HPV types). For subsequent statistical analysis we recoded HPV types into category 1 as HPV negative; category 2 as low-risk HPV positive, and category 3 as any high-risk HPV positive (single or multiple).
We compared baseline socio-demographic, sexual risk characteristics, CSTs, HPV status, and HIV status between cervical pathology categories using the Pearson’s chi squared or Fisher’s exact tests for categorical variables and ANOVA for normally distributed continuous variables. To assess the association between vaginal microbiome CSTs with grades of cervical pathology, we considered the clinical relevance of the four cervical cytology categories. We combined NILM and LSIL as one group (normal or minor grades dysplasia) while HSIL and ICC were combined as the clinically important group that requires immediate evaluation and treatment. Observations with CST-II (n = 4) were excluded (none was CST-V) from analyses as being too sparse for inference. We performed bivariable logistic regression for the primary exposure variable (CST) and other predictor variables using robust logistic regression to obtain unadjusted odd ratios, 95% confidence intervals (CI), and Wald p-values for associations between HSIL/ICC and NILM/LSIL. We used backward stepwise selection method to build final models checked by the Akaike’s Information Criteria for model selection, including age in years, HIV status, HPV status, CSTs, years of completed education, and total number of births (parity). While tobacco use has been shown to be an important co-factor in cervical cancer risk, only two participants in our sample reported ever smoking, precluding its inclusion in analysis. We then modelled the analyses stratified by high-risk HPV versus negative for HPV (the numbers for low-risk HPV were too small for stratification in the model) to assess the possible effect of CST on cervical pathology probably mediated by high-risk HPV, a known epidemiologic factor in cervical carcinogenesis [2, 38].
One hundred and fifty-five (155) participants were enrolled, 151 with microbiome data for this analysis (four excluded for having CST-II). The median age of the participants was 52 years (IQR:43–58), 47.7% were HIV positive, and 58.1% had HSIL/ICC. The HR-HPV genotyping results were available in 138 of the participants (the HPV sample of 13 participants were not processed for PCR genotyping and were excluded in all analysis with HR-HPV as an outcome). Of the 138 with HPV data, 56 (40.6%) were negative for HPV, 14 (10.1%) had low-risk HPV, 37 (26.8%) had single HR-HPV, and 31 (22.5%) had multiple HR-HPV types. Of the 31 participants with multiple HR-HPV types, 16 or 18 with any other HR-HPV type were seen in 14 (45.2%), 35 and other HR-HPV types in 4 (12.9%), while other mixed infections other than 16/18/35 were 13 (41.9%). The overall prevalence of any HR-HPV type (single and multiple) was 49.3%, with a higher proportion in women with HSIL/ICC (NILM 31.6%, LSIL 46.5%, HSIL 40.8%, and 81.5% ICC; p = 0.007).
Women with high grade cervical precancer and ICC were on average older than women with NILM and LSIL (Table 1). The distribution of CST showed that 101 participants (66.9%; 95% CI: 58.9–74.0%) had CST-IV. CST-III (L. iners dominated) was found in 33 participants (21.9%; 95% CI: 15.9–29.2), while 17 participants had CST-I VMB (11.3%; 7.1–17.4). The distribution of CSTs did not differ significantly by cervical pathology (p = 0.153, Table 1), though cervical cancer diagnosis was more common among women with CST-IV (Table 2). The association of other factors with CSTs as outcome have been summarized in Table 2. HPV status did not significantly differ by CST, though CST-I was more likely among women with HIV infection, and CST-IV was more likely with increased parity.
The distribution of both low and high-risk HPV types in the study sample is shown in Fig. 1 with HR-HPV 16 being the most common. Women with HIV were also significantly more likely to have HSIL/ICC (70.2% vs. 29.3% for women without HIV, p-value = 0.001).
The unadjusted logistic regression (Table 3) showed that age in years (OR = 1.05; 95% CI: 1.01–1.08), positive HIV status (OR = 2.67; 95% CI: 1.38–5.20), and total births (parity) > 4 (OR = 5.94; 95% CI: 1.83–19.3) were significantly associated with HSIL/ICC, while completion of 7–12 years of education (OR = 0.28; 95% CI: 0.09–0.88), or more (OR = 0.29; 95% CI: 0.10–0.85) had significantly lower likelihood for HSIL/ICC. In the final multivariable logistic regression model (Table 3), adjusting for age, HR-HPV status, HIV status, CST, parity, and years of completed education, women with HIV had over three times the odds (aOR = 3.21). Women with any HR-HPV had 2.72 times greater odds of having HSIL/ICC.
To further understand the possible associations between CSTs and HSIL/ICC, we stratified our multivariable logistic regression models by HR-HPV status (Table 4). Among women with HR-HPV detected, the effect of HIV was strengthened with aOR = 14.1 (95% CI: 2.34–85.2). It was also observed that the association between CST-III/CST-IV and HSIL/ICC were accentuated in the presence of HR-HPV, although this finding did not reach statistical significance: CST-III (aOR = 6.69; 95% CI: 0.67–66.6) and CST-IV (aOR = 3.64; 95% CI: 0.63–20.9).
This study aimed at understanding the predominant CSTs in Nigerian women and to estimate the associations with cervical precancer and invasive cervical cancer, HIV, and HPV infection found that CST- III and CST-IV vaginal microbiome (VMB) communities were more frequent among women with high-grade cervical precancer and ICC, especially in the presence of HR-HPV. In agreement with a previous report in Nigeria , we found that CST-IV was the most frequent in the study sample and did not differ significantly by HR-HPV status, although the former report did not relate the CSTs and HR-HPV to cervical pathology as in our study. Our analyses also confirmed strong associations between HIV/HR-HPV and high-grade cervical precancer or invasive cervical cancer found in several prior studies [2, 39]. Our data showed that women with HIV had 3.2 times greater odds of high-grade cervical precancer and invasive cervical cancer compared to women without HIV. Similarly, we found that the detection of any HR-HPV types increased the odds of high-grade cervical dysplasia and ICC by at least 2.7 times.
In the presence of any HR-HPV types and compared to CST-I (L. crispatus dominated), CST-III (L. iners dominated) and CST-IV (diverse, non-optimal bacterial species) vaginal microbiomes were associated with high-grade precancer or invasive cervical cancer in a large systematic review and network meta-analysis . Specifically, the meta-analysis showed that vaginal microbiota dominated by non-Lactobacilli species (CST-IV) or L. iners were associated with 3–5 times higher odds of detecting any prevalent HR-HPV and cervical precancer or invasive cervical cancer compared with a L. crispatus-dominated state . Also a study report among South African women assessing incident HR-HPV, vaginal dysbiosis and incident cervical intraepithelial neoplasia grade 2 or more (CIN 2+) found that the acquisition of HR-HPV changes the vaginal microbiome and that anaerobic dysbiosis risk seem to increase concurrently with development of CIN2 + suggesting that the role of the vaginal microbiome composition in CIN2 + development could be mediated by persistent HR-HPV infection . Our findings are therefore in agreement with the literature implying that a better understanding of the cervico-vaginal microenvironment could be a key step for HPV viral clearance . This mechanistic role of the cervico-vaginal microbiome community state types is critical and could be explored in developing therapeutic targets for cervices already infected with the HR-HPV types but yet to advance the process of cervical carcinogenesis given the high prevalence of these common sexually transmissible viral infection in the general population. Currently, the primary HPV preventive vaccines are not readily available or affordable in most LMICs making the search for therapeutic options for cervical HPV infection a promising effort.
Furthermore, the associations between CST and the two sexually transmissible viruses - HR-HPV and HIV - have public health implications for cervical cancer prevention. In 2018, the World Health Organization launched a global initiative to eliminate cervical cancer as a public health problem by adopting highly effective evidence-based interventions 2030 . These interventions includes increasing the population-level coverage for HPV vaccination of 90% of eligible young girls between 9 and 15 years; screening of 70% eligible women by 35 years and 45 years, and effective treatment of 90% women with precancer and invasive cervical cancer by 2030 . As discussed earlier, most countries in LMICs including Nigeria do not have a national program on HPV vaccination , making interventions for therapeutics, and HIV prevention, and to maintain a healthier cervico-vaginal environment effective and plausible public health approaches for cervical cancer prevention and control efforts. Indeed, some of the known epidemiologic factors for cervical cancer [15, 24, 44,45,46,47,48,49,50,51,52,53,54] such as early sexual debut, multiple sexual partners, smoking, and unsafe sexual practices are modifiable and could be controlled using public health approaches. To further support this, the Global Burden of Cancer attributable to modifiable risk factors confirmed the contribution of behavioral risk factors such as smoking, alcohol use, and unsafe sex as modifiable risk factors attributable to a huge burden of cancer in the population . Achieving population level changes in unsafe sexual practices, however, will require societal efforts in addressing the structural and system drivers of early sexual debut and multiple sex partners such as economic, legal, cultural, and religious beliefs. In the area of therapeutics, it is possible that the biological interplay of CSTs, HR-HPV, and HIV immunosuppression seen in our study could be further examined to develop point-of-care molecular testing for rapid detection and quantification of Lactobacillus and other non-Lactobacillus dominant cervico-vaginal microbiota as biomarkers for evaluation, early detection, and treatment of high-grade precancer in the population.
Our study also showed the importance of socioeconomic factors, such as more years of completed education was associated with lesser odds for high-grade cervical precancer/ICC, while age and higher number of total births (parity) were associated with greater odds for high-grade cervical precancer/ICC. Furthermore, the unadjusted analyses of age at first pregnancy showed a significant relationship favoring women whose first pregnancies were at relatively younger age (See Table 1), and these women were more likely to have a CST-IV (Table) compared to women who were older at first pregnancy. The epidemiologic findings seem to be in agreement with recent genome-wide association study which showed that women who were older at first pregnancy had a lower odds for ICC . The interplay of these sociodemographic variables in association with unfavorable CSTs and in predicting high grades precancer and/or invasive cervical cancer is an interesting finding particularly in sub-Saharan Africa where educational attainment is low, and lower educational status has been shown to be associated with higher number of births . As shown in Table 4, there is a greater odds for high-grade cervical precancer or ICC if parity is > 4 in women with negative HR-HPV. The epidemiologic nuances surrounding this relationship, and the mechanism by which repeated childbirths contributes to cervical carcinogenesis requires further investigation particularly in the absence of HR-HPV given emerging ideas that certain sub-types of cervical cancer are not caused by HPV [57, 58].
Our data revealed trends that need confirmatory longitudinal studies with larger sample size to further characterize biomarkers of the vaginal microbiome and HR-HPV detection in Nigerian women at various stages of cervical precancer to understand how they predict progression to high-grade precancer or invasive cervical cancer stages. We intuitively hypothesized that therapeutic approaches that modifies and maintains a more optimal cervico-vaginal microenvironment could reduce the acquisition or help the clearance of HR-HPV in cervices before the process of HPV-induced cervical carcinogenesis advances to high-grades precancer and invasive stages. This would be a highly translatable public health research effort and holds the prospect for reducing population level incidence of cervical precancer and/or invasive cervical cancer particularly in Nigeria and other LMICs where primary preventive HPV vaccines and cervical cancer screening coverage are still very low.
Our study findings have a few limitations, mainly in the relatively small sample size and cross-sectional design. We also lacked the data on duration for HPV infection and did not have data on HIV treatment with antiretrovirals, viral load, or immune suppression levels with CD4 count markers. We also did not have control on intravaginal practices such as douching with soap and water, and other agents on vaginal microbiome as previously reported in the literature . Nonetheless, our findings, the first of their kind from Nigeria characterizing vaginal microbiome and relating it to HR-HPV and spectrum of normal cervix, low and high-grades precancer and invasive cervical cancer pathology, suggest a potential relationship between CST, HR-HPV, and progression to high-grade precancer and/or ICC. The findings also support the well-known evidence that HR-HPV is the principal causative factor [1, 2, 60] in cervical precancer and progression to invasive cervical cancer and add to the literature in the paucity of this information from Nigeria and sub-Saharan Africa.
We conclude that non-optimal CSTs, particularly when HR-HPV is detected, may assist in identifying women with high-grade cervical precancer. Characterization of CST using point-of-care molecular testing particularly in women identified with any HR-HPV may be a useful biomarker for early detection of severe precancer and cervical cancer prevention. Longitudinal data will improve our understanding of the temporal effect of non-optimal CSTs, HR-HPV, and other factors in cervical cancer development, and may lead to translatable public health interventions towards cervical cancer prevention and elimination globally.
All the relevant data for this analysis have been presented in the body of this manuscript and the associated tables and figures. The original database of the phenotypic data is available in the institutional REDCap database of the University of Jos. The microbiome genomic raw data has been submitted to the Sequence Read Archive (SRA) database at the NCBI with Accession Number: PRJNA914231.
Stelzle D, Tanaka LF, Lee KK, Ibrahim Khalil A, Baussano I, Shah ASV, McAllister DA, Gottlieb SL, Klug SJ, Winkler AS, et al. Estimates of the global burden of cervical cancer associated with HIV. Lancet Glob Health. 2021;9(2):e161–9.
Plummer M, de Martel C, Vignat J, Ferlay J, Bray F, Franceschi S. Global burden of cancers attributable to infections in 2012: a synthetic analysis. The Lancet Global Health. 2016;4(9):e609–16.
Jonah Musa CA, Babafemi Taiwo B, Berzins O, Silas PH, Daru O, Agbaji G, Imade AS, Sagay JA, Idoko, Phyllis J. Kanki and Robert L Murphy: high-risk human papilloma virus and cervical abnormalities in HIV-infected women with normal cervical cytology. Infect Agents Cancer 2014, 9.
Mogtomo ML, Malieugoue LC, Djiepgang C, Wankam M, Moune A, Ngane AN. Incidence of cervical disease associated to HPV in human immunodeficiency infected women under highly active antiretroviral therapy. Infect Agent Cancer. 2009;4:9.
Williamson AL. The Interaction between Human Immunodeficiency Virus and Human Papillomaviruses in Heterosexuals in Africa. J Clin Med. 2015;4(4):579–92.
Ramadhani S, Chambuso SS, Salum J, Lidenge N, Mwakibete RM, Medeiros. Influence of HIV/AIDS on Cervical Cancer: a retrospective study from Tanzania. J Global Oncol 2016.
Denslow S, Rositch AF, Firnhaber C, Ting J, Smith J. Incidence and progression of cervical lesions in women with HIV: a systematic global review. Int J STD AIDS. 2014;25(3):163–77.
Menon S, Broeck DV, Rossi R, Ogbe E, Harmon S, Mabeya H. Associations between vaginal infections and potential high-risk and high-risk human papillomavirus genotypes in female sex workers in western Kenya. Clin Ther. 2016;38(12):2567–77.
Mitra A, MacIntyre D, Lee Y, Smith A, Marchesi J, Lehne B, Bhatia R, Lyons D, Paraskevaidis E, Li J, et al. Characterisation of the vaginal microbiome in cervical intraepithelial neoplasia. The Lancet. 2016;387:75.
Mitra A, MacIntyre DA, Marchesi JR, Lee YS, Bennett PR, Kyrgiou M. The vaginal microbiota, human papillomavirus infection and cervical intraepithelial neoplasia: what do we know and where are we going next? Microbiome. 2016;4(1):58.
Piyathilake CJ, Ollberding NJ, Kumar R, Macaluso M, Alvarez RD, Morrow CD. Cervical Microbiota Associated with higher Grade Cervical Intraepithelial Neoplasia in Women infected with high-risk human papillomaviruses. Cancer Prev Res (Phila). 2016;9(5):357–66.
Ghosh I, Muwonge R, Mittal S, Banerjee D, Kundu P, Mandal R, Biswas J, Basu P. Association between high risk human papillomavirus infection and co-infection with Candida spp. and Trichomonas vaginalis in women with cervical premalignant and malignant lesions. J Clin Virol. 2017;87:43–8.
Feng M, Duan R, Gao Y, Zhang H, Qiao Y, Li Q, Zhao F. Role of Epstein-Barr Virus and Human Papillomavirus Coinfection in Cervical Intraepithelial Neoplasia in Chinese Women Living with HIV. Front Cell Infect Microbiol. 2021;11:703259.
Organization WH. Global strategy to accelerate the elimination of cervical cancer as a public health problem. In. Geneva: World Health organization; 2020.
Musa J. A Sexually Transmissible Cancer: Miracle of Cure or Prevention? In: University of Jos Inaugural Lecture Series Edited by Jos Uo, vol. 94. University of Jos Printing Press: University of Jos Printing Press; 2021: 3–69.
Zeber-Lubecka N, Kulecka M, Lindner B, Krynicki R, Paziewska A, Nowakowski A, Bidzinski M, Ostrowski J. Increased diversity of a cervical microbiome associates with cervical cancer. Front Oncol. 2022;12:1005537.
Usyk M, Zolnik CP, Castle PE, Porras C, Herrero R, Gradissimo A, Gonzalez P, Safaeian M, Schiffman M, Burk RD, et al. Cervicovaginal microbiome and natural history of HPV in a longitudinal study. PLoS Pathog. 2020;16(3):e1008376.
Suehiro TT, Malaguti N, Damke E, Uchimura NS, Gimenes F, Souza RP, Sela da Silva VR, Lopes Consolaro ME. Association of human papillomavirus and bacterial vaginosis with increased risk of high-grade squamous intraepithelial cervical lesions. Int J Gynecol Cancer. 2019;29(2):242–9.
Audirac-Chalifour A, Torres-Poveda K, Bahena-Roman M, Tellez-Sosa J, Martinez-Barnetche J, Cortina-Ceballos B, Lopez-Estrada G, Delgado-Romero K, Burguete-Garcia AI, Cantu D, et al. Cervical Microbiome and Cytokine Profile at various stages of cervical Cancer: a pilot study. PLoS ONE. 2016;11(4):e0153274.
Carrero Y, Mosquera J, Callejas D, Alvarez-Mon M. In situ increased chemokine expression in human cervical intraepithelial neoplasia. Pathol Res Pract. 2015;211(4):281–5.
Dzutsev A, Goldszmid RS, Viaud S, Zitvogel L, Trinchieri G. The role of the microbiota in inflammation, carcinogenesis, and cancer therapy. Eur J Immunol. 2015;45(1):17–31.
Gradissimo A, Lam J, Attonito JD, Palefsky J, Massad LS, Xie X, Eltoum IE, Rahangdale L, Fischl MA, Anastos K et al. Methylation of high-risk human papillomavirus genomes are Associated with Cervical Precancer in HIV-positive women. Cancer Epidemiol Biomarkers Prev 2018.
Laniewski P, Barnes D, Goulder A, Cui H, Roe DJ, Chase DM, Herbst-Kralovetz MM. Linking cervicovaginal immune signatures, HPV and microbiota composition in cervical carcinogenesis in non-hispanic and hispanic women. Sci Rep. 2018;8(1):7593.
Brusselaers N, Shrestha S, van de Wijgert J, Verstraelen H. Vaginal dysbiosis and the risk of human papillomavirus and cervical cancer: systematic review and meta-analysis. Am J Obstet Gynecol. 2019;221(1):9–18e18.
Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, Karlebach S, Gorle R, Russell J, Tacket CO, et al. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4680–7.
Mitra A, MacIntyre DA, Lee YS, Smith A, Marchesi JR, Lehne B, Bhatia R, Lyons D, Paraskevaidis E, Li JV et al. Cervical intraepithelial neoplasia disease progression is associated with increased vaginal microbiome diversity. Sci Rep 2015, 5(1).
Adebamowo SN, Olawande O, Famooto A, Dareng EO, Offiong R, Adebamowo CA, Group HAAR. Persistent low-risk and high-risk human papillomavirus infections of the Uterine Cervix in HIV-Negative and HIV-Positive women. Front Public Health. 2017;5:178.
Dareng EO, Ma B, Famooto AO, Adebamowo SN, Offiong RA, Olaniyan O, Dakum PS, Wheeler CM, Fadrosh D, Yang H, et al. Prevalent high-risk HPV infection and vaginal microbiota in nigerian women. Epidemiol Infect. 2016;144(1):123–37.
Solomon D, Davey D, Kurman R, Moriarty A, O’Connor D, Prey M, Raab S, Sherman M, Wilbur D, Wright T Jr, Young N. The 2001 Bethesda System Terminology for reporting results of cervical cytology. JAMA. 2001;287:2114–9.
Patel HK, Ikpe S, Bronson M, Birhanu S, Abimiku A, Jahun I, Detorio M, Lupoli K, Yavo D, Bassey OO, et al. Performance of HIV rapid testing algorithm in Nigeria: findings from a household-based Nigeria HIV/AIDS Indicator and Impact Survey (NAIIS). PLOS Glob Public Health. 2022;2(7):e0000466.
Anyplex2 hpv28. detection assays [https://www.seegene.com/assays/anyplex2_hpv28_detection].
Naqib A, Poggi S, Wang W, Hyde M, Kunstman K, Green SJ. Making and sequencing heavily multiplexed, high-throughput 16S ribosomal RNA gene amplicon libraries using a flexible, two-stage PCR protocol. Gene expression analysis. New York, NY: Human Press; 2018: 149–69.
Herlemann DP, Labrenz M, Jurgens K, Bertilsson S, Waniek JJ, Andersson AF. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011;5(10):1571–9.
Douglas W, Fadrosh BM, Gajer P, Sengamalay N, Ott S, Rebecca M. Brotman and Jacques Ravel: an improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. Microbiome 2014:2–6.
Holm JB, Humphrys MS, Robinson CK, Settles ML, Ott S, Fu L, Yang H, Gajer P, He X, McComb E et al. Ultrahigh-Throughput Multiplexing and Sequencing of > 500-Base-Pair Amplicon Regions on the Illumina HiSeq 2500 Platform. mSystems 2019, 4(1).
France MT, Ma B, Gajer P, Brown S, Humphrys MS, Holm JB, Waetjen LE, Brotman RM, Ravel J. VALENCIA: a nearest centroid classification method for vaginal microbial communities based on composition. Microbiome. 2020;8(1):166.
Odukoya O, Nenrot D, Adelabu H, Katam N, Christian E, Holl J, Okonkwo A, Kocherginsky M, Kim KY, Akanmu S, et al. Application of the research electronic data capture (REDCap) system in a low- and middle income country- experiences, lessons, and challenges. Health Technol (Berl). 2021;11(6):1297–304.
Jan m. m. walboomers mvj, m. michele manos, f. xavier bosch, j. alain kummer, keerti v. shah, peter j. f. snijders, julian peto, chris j. l. m. meijer and nubia mun˜ oz: HUMAN PAPILLOMAVIRUS IS a NECESSARY CAUSE OF INVASIVE CERVICAL CANCER WORLDWIDE. J Pathol. 1999;189:12–9.
Liu G, Sharma M, Tan N, Barnabas RV. HIV-positive women have higher risk of human papilloma virus infection, precancerous lesions, and cervical cancer. AIDS. 2018;32(6):795–808.
Norenhag JDJ, Olovsson M, Verstraelen H, Engstrand L, Brusselaers N. The vaginal microbiota, human papillomavirus and cervical dysplasia: a systematic review and network meta-analysis. BJOG. 2020;127(2):171–80.
van de Wijgert J, Gill AC, Chikandiwa A, Verwijs MC, Kelly HA, Omar T, Delany-Moretlwe S, Segondy M, Francis S, Darby AC, et al. Human papillomavirus infection and cervical dysplasia in HIV-positive women: potential role of the vaginal microbiota. AIDS. 2020;34(1):115–25.
Cascardi E, Cazzato G, Daniele A, Silvestris E, Cormio G, Di Vagno G, Malvasi A, Loizzi V, Scacco S, Pinto V et al. Association between Cervical Microbiota and HPV: Could This Be the Key to Complete Cervical Cancer Eradication? Biology (Basel) 2022, 11(8).
Bruni LB-RL, Mena AlberoGSerranoB, GómezD M, Muñoz J, Bosch FX, de Sanjosé S, ICO/IARC Information Centre on HPV and Cancer (HPV Information Centre). Human papillomavirus and related Diseases in Nigeria. Summary Report 27 July 2017. ICO/IARC HPV Information Centre; 2017.
Bynum SA, Wigfall LT, Brandt HM, Julious CH, Glover SH, Hebert JR. Social and Structural Determinants of Cervical Health among women engaged in HIV Care. AIDS Behav. 2016;20(9):2101–9.
Dahiya N, Bachani D, Acharya AS, Sharma DN, Gupta S, Haresh KP. Socio-Demographic, Reproductive and Clinical Profile of Women diagnosed with Advanced Cervical Cancer in a Tertiary Care Institute of Delhi. J Obstet Gynaecol India. 2017;67(1):53–60.
Houlihan CF, Baisley K, Bravo IG, Kapiga S, de Sanjose S, Changalucha J, Ross DA, Hayes RJ, Watson-Jones D. Rapid acquisition of HPV around the time of sexual debut in adolescent girls in Tanzania. Int J Epidemiol. 2016;45(3):762–73.
Luhn P, Walker J, Schiffman M, Zuna RE, Dunn ST, Gold MA, Smith K, Mathews C, Allen RA, Zhang R, et al. The role of co-factors in the progression from human papillomavirus infection to cervical cancer. Gynecol Oncol. 2013;128(2):265–70.
Uzoma Ononogbu MA, Modibbo F, Lawal I, Offiong R, Olaniyan O, Dakum P. Donna Spiegelman, William Blattner and Clement Adebamowo: cervical cancer risk factors among HIV-infected nigerian women. BMC Public Health. 2013;13:582.
Yetimalar H, Kasap B, Cukurova K, Yildiz A, Keklik A, Soylu F. Cofactors in human papillomavirus infection and cervical carcinogenesis. Arch Gynecol Obstet. 2012;285(3):805–10.
Anorlu RI. Cervical cancer: the sub-saharan african perspective. Reprod Health Matters. 2008;16(32):41–9.
Daniel GO, Musa J, Akindigh TM, Shinku F, Shuaibu SI, Kwaghe B, Afolaranmi T, Okpala h, Agbaji O, Sagay A. Prevalence and predictors of precancerous cervical lesions among HIV-positive women in Jos, north‐central Nigeria. International Journal of Gynecology & Obstetrics.
Kapeu AS, Luostarinen T, Jellum E, Dillner J, Hakama M, Koskela P, Lenner P, Love A, Mahlamaki E, Thoresen S, et al. Is smoking an independent risk factor for invasive cervical cancer? A nested case-control study within nordic biobanks. Am J Epidemiol. 2009;169(4):480–8.
Kathleen McIntyre-Seltman PEC, Guido R, Schiffman M, Wheeler CM, for The ALTS Group. Smoking is a risk factor for cervical intraepithelial Neoplasia Grade 3 among Oncogenic Human Papillomavirus DNA–Positive women with equivocal or mildly abnormal cytology. Cancer Epidemiol Biomarkers Prev. 2005;14(5):1165–70.
Tran KB, Lang JJ, Compton K, Xu R, Acheson AR, Henrikson HJ, Kocarnik JM, Penberthy L, Aali A, Abbas Q, et al. The global burden of cancer attributable to risk factors, 2010–19: a systematic analysis for the global burden of Disease Study 2019. The Lancet. 2022;400(10352):563–91.
Bowden SJ, Bodinier B, Kalliala I, Zuber V, Vuckovic D, Doulgeraki T, Whitaker MD, Wielscher M, Cartwright R, Tsilidis KK, et al. Genetic variation in cervical preinvasive and invasive disease: a genome-wide association study. Lancet Oncol. 2021;22(4):548–57.
Musa J, Achenbach CJ, Evans CT, Jordan N, Daru PH, Silas O, Sagay AS, Anorlu R, Mehta SD, Wehbe F. HIV status, age at cervical Cancer screening and cervical cytology outcomes in an opportunistic screening setting in Nigeria: a 10-year Cross sectional data analysis. Infect Agents Cancer. 2019;14(1):43.
Fernandes A, Viveros-Carreno D, Hoegl J, Avila M, Pareja R. Human papillomavirus-independent cervical cancer. Int J Gynecol Cancer. 2022;32(1):1–7.
Schmeler KM, Batman SH. Human papillomavirus-independent cervical cancer: what are the implications? Int J Gynecol Cancer. 2022;32(1):8.
Daniel GA, Hu Y, Tsementzi D, Jhaney CI, Hu YJ, Yeager KA, Bai J, Dolan M, Bruner DW. Exploring the vaginal microbiome and Intravaginal Practices in Postmenopausal Women. Nurs Res. 2021;70(5):405–11.
Schiffman M, Castle PE, Jeronimo J, Rodriguez AC, Wacholder S. Human papillomavirus and cervical cancer. The Lancet. 2007;370(9590):890–907.
The contributions of Joyce Asufi, RN, BNsc, Maryam shuaibu, RN, and Patricia Asemota, RN, BNsc, of the Nursing department and the “Operation Stop” Cervical Cancer Unit of Jos University Teaching Hospital, Jos, Nigeria is hereby acknowledged. I acknowledged the contributions and dedication of the laboratory staff (Mr. Acheng Shedrach Yakubu, Mr. Rapp Chuwang Nyam, Mr. Michael Eshioramhe, Mrs. Flavia Godfrey Lamo, Mrs. Tsok Yop Gwom, and Jacob Bwalsim) of the genomics and postgraduate core facility of the College of Health Sciences, University of Jos.
I acknowledge the funding support from the Strohm’s family through the Robert J. Havey MD Institute for Global Health to provide access to cervical cancer screening for indigent women through the community cervical cancer screening project in Jos Nigeria.
The NIH/Fogarty international center provided funding support for this research and the research-protected time for writing this manuscript under the award number: K43TW011416. The NIH/NCI funded U54CA221205 provided institutional support facilities for the conduct of this research in Jos Nigeria. JM, FJM, CJN, and OAS acknowledges training support as Fogarty fellows under the NIH/FIC D43TW009575, while YZ and OAS acknowledged funding support from NIH/FIC R21TW12092. The findings and views expressed in this manuscript are those of the authors and do not necessarily represent the views of the Fogarty International Center or the National Institutes of Health.
The authors have no competing interest to declare.
Ethics approval and consent to participate
The research protocol for this study with informed consent for the enrollment procedures were reviewed and approved by the Jos University Teaching Hospital’s Health and Research Ethics Committee. The approval reference for the research protocol is: JUTH/DCS/ADM/127/XXIX/1546. All participants signed an informed consent for all enrollment procedures which were conducted in accordance with the relevant regulatory guidelines as in the Helsinki’s declaration.
Consent for publication
All participant provided informed consent for publication of research findings as summary with no identifiable personal information.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Musa, J., Maiga, M., Green, S.J. et al. Vaginal microbiome community state types and high-risk human papillomaviruses in cervical precancer and cancer in North-central Nigeria. BMC Cancer 23, 683 (2023). https://doi.org/10.1186/s12885-023-11187-5