Human papillomavirus type and tobacco use as predictors of survival in early stage cervical carcinoma

Integrated analysis of cervical squamous cell carcinoma cohorts from three continents reveals conserved subtypes of prognostic significance

Human papillomavirus (HPV)-associated cervical cancer is a leading cause of cancer deaths in women. Here we present an integrated multi-omic analysis of 643 cervical squamous cell carcinomas (CSCC, the most common histological variant of cervical cancer), representing patient populations from the USA, Europe and Sub-Saharan Africa and identify two CSCC subtypes (C1 and C2) with differing prognosis. C1 and C2 tumours can be driven by either of the two most common HPV types in cervical cancer (16 and 18) and while HPV16 and HPV18 are overrepresented among C1 and C2 tumours respectively, the prognostic difference between groups is not due to HPV type. C2 tumours, which comprise approximately 20% of CSCCs across these cohorts, display distinct genomic alterations, including loss or mutation of the STK11 tumour suppressor gene, increased expression of several immune checkpoint genes and differences in the tumour immune microenvironment that may explain the shorter survival associated with this group. In conclusion, we identify two therapy-relevant CSCC subtypes that share the same defining characteristics across three geographically diverse cohorts.

Human papillomavirus (HPV) is a known cause of cervical cancer. Here, the authors perform a multi-omic analysis using published cervical squamous cell carcinoma cohorts from the USA, Europe, and SubSaharan Africa and identify two cervical squamous cell carcinoma subtypes that display prognostic differences.

Verification of the association of the cycle threshold (Ct) values from HPV testing on Cobas4800 with the histologic grades of cervical lesions using data from two population-based cervical cancer screening trials

Objective

To verify the association of high-risk human papillomavirus (hrHPV) viral load reflected by cycle threshold (Ct) values from HPV testing on Cobas4800 assay with the histologic grades of cervical lesions via analysis on the combined data from two cervical cancer screening trials and to explore the referability of Ct values in management of the abnormalities from cervical cancer primary screening.

Methods

We analyzed the data from Chinese Multi-Center Screening Trial (CHMUST) and BUJI Cervical Cancer Screening Study Project (BUJI Study). All data to be analyzed in this paper were related to provider-collected samples. One-way ANOVA was performed to compare the Ct values among different histological groups, and Kendall’s tau-b correlation was applied to examine the association between Ct values and cervical lesion grades. The stepwise incidence of CIN2+ and CIN3+ in every 100 HPV positive individuals were calculated according to the descending of the genotype specific Ct values. The highest Ct values related to CIN3+ incidence rate 4% (or 25%) were used as the cutoffs to distinguish low-Ct value cases from the high-Ct value ones.

Results

A total of 1376 women in CHUMUST and BUJI Study who were aged 30–59 and positive of hrHPV for provider-collected samples on Cobas4800 with complete data in terms of the relevant Ct values (CtV) and the histological diagnosis were included for analysis. Our data showed significant difference among different histological grades of cervical lesions in the CtV of hrHPV, HPV16-plus (positive of HPV16 only or HPV16 plus 18 and/or pooled 12-HPV), and pooled 12-HPV (P < 0.05). No significant difference was found among different grades of lesions in term of correlated CtV of HPV18-plus (positive of HPV18 only or HPV18 plus pooled 12-HPV) (P > 0.05). The CIN2+ or CIN3+ incidence per 100 positives significantly increased corresponding to the descending of the CtV of hrHPV, HPV16-plus and pooled 12-HPV. Compared with high-CtV groups (CtV > 33.2 for hrHPV, CtV > 29.6 for pooled 12-HPV), the relevant risks (RRs) of CIN2+ for hrHPV and pooled 12-HPV positive groups with low-CtV (CtV ≤ 33.2 and ≤ 29.6, respectively) were 3.2 (95%CI 2.18–4.80) and 2.3 (95%CI 1.50–3.45). Similarly, the RRs of CIN3+ for hrHPV and pooled 12-HPV positive groups with low-CtV were 6.5 (95%CI 2.83–14.80) and 2.7 (95%CI 1.15–6.39), respectively. The RRs of CIN2+ for medium- (30.3 < CtV ≤ 37.4) and low- (≤ 30.3) CtV HPV16-plus positives were 5.1 (95%CI 0.68–38.38) and 20.6 (95%CI 2.96–143.92) related to high-CtV (> 37.4) groups, and the CIN3+ incidence in low-CtV value group was nine-fold higher of that in medium-CtV ones [RRs, 9.0 (95%CI 2.89–28.10)]. In comparing with the algorithms of “HPV16-plus/18-plus + cytology ≥ ASCUS for pooled 12-HPV”, triage algorithm “HPV16-plus/18-plus + Ct value ≤ 33.2 for pooled 12-HPV” could achieve a comparable sensitivity of 93.2%.

Conclusion

HPV viral loads reflected by Ct values for hrHPV, HPV16-plus and pooled 12-HPV from Cobas4800 HPV testing were directly associated with the severity of cervical lesions. A lower HPV genotype-specific Ct value prompted a significantly high CIN3+ risk of 4% or higher in women positive of hrHPV, HPV16-plus or pooled 12-HPV, indicating that HPV viral load reflected by Ct values on Cobas4800 may be a promising risk indicator in management of abnormalities from primary cervical cancer screening.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13027-022-00440-4.

Prognostic implication of human papillomavirus types in cervical cancer patients: a systematic review and meta-analysis

Background

To estimate the prognostic relevance of human papillomavirus (HPV) 16 and HPV 18 in patients with cervical cancer.

Method

We searched PubMed, EMBASE, American Society of Clinical Oncology (ASCO) and the European Society of Medical Oncology (ESMO), CNKI, and Wanfang databases to search primary articles illustrating the survival outcomes in cervical cancer patients with or without HPV 16/18 infection. A meta-analysis was conducted to generate a combined hazard ratio (HR) with 95% confidence intervals (CI) for progression-free survival (PFS), disease free survival (DFS) and overall survival (OS).

Results

A total of 13 studies were included. Our meta-analysis revealed that HPV 16 positive did not have any impact on OS (HR, 0.76; 95% CI = 0.37–1.54; P = 0.44). Cervical cancer patiensts infected with HPV 18 had worse OS (HR, 1.66; 95% CI = 1.28–2.17; P = 0.0001), DFS (HR, 2.10; 95% CI = 1.73–2.54; P < 0.0001) and worse PFS (HR, 2.97; 95% CI = 1.69–5.23; P = 0.00012) compared with those not infected with HPV 18. cervical cancer patiensts infected with HPV 18 had worse PFS compared with those infected with HPV 16 ((HR, 1.34; 95% CI = 1.06–1.70; P = 0.01).

Conclusion

Cervical cancer patients infected with HPV 18 had worse survival compared with cervical cancer patients with HPV 16 infection.

Analysis of Ugandan cervical carcinomas identifies human papillomavirus clade-specific epigenome and transcriptome landscapes

Cervical cancer is the most common cancer affecting sub-Saharan African women and is prevalent among HIV-positive (HIV+) individuals. No comprehensive profiling of cancer genomes, transcriptomes or epigenomes has been performed in this population thus far. We characterized 118 tumors from Ugandan patients, of whom 72 were HIV+, and performed extended mutation analysis on an additional 89 tumors. We detected human papillomavirus (HPV)-clade-specific differences in tumor DNA methylation, promoter- and enhancer-associated histone marks, gene expression and pathway dysregulation. Changes in histone modification at HPV integration events were correlated with upregulation of nearby genes and endogenous retroviruses.

Genetic variations in human papillomavirus and cervical cancer outcomes

Cervical cancer is driven by persistent infection of human papillomavirus (HPV), which is influenced by HPV type and intratypic variants, yet the impact of HPV type and intratypic variants on patient outcomes is far less understood. Here, we examined the association of cervical cancer stage and survival with HPV type, clade, lineage, and intratypic variants within the HPV E6 locus. Of 1,028 HPV-positive cases recruited through the CerGE study, 301 were in-situ and 727 were invasive cervical cancer (ICC), with an average post-diagnosis follow-up of 4.8 years. HPV sequencing was performed using tumor-isolated DNA to assign HPV type, HPV 16 lineage, clade, and intratypic variants within the HPV 16 E6 locus, of which nonsynonomous variants were functionally annotated by molecular modeling. HPV 18-related types were more prevalent in ICC compared to in-situ disease and associated with significantly worse recurrence-free survival (RFS) compared to HPV 16-related types. The HPV 16 Asian American lineage D3 and Asian lineage A4 associated more frequently with ICC than with in situ disease and women with an intratypic HPV 16 lineage B exhibited a trend toward worse RFS than those with A, C, or D lineages. Participants with intratypic E6 variants predicted to stabilize the E6-E6AP-p53 complex had worse RFS. Variants within the highly immunogenic HPV 16 E6 region (E14-I34) were enriched in ICC compared to in-situ lesions but were not associated with survival. Collectively, our results suggest that cervical cancer outcome is associated with HPV variants that affect virus-host interactions.

Prognostic implication of human papillomavirus types and species in cervical cancer patients undergoing primary treatment

High-risk human papillomavirus (HPV) types are associated with cervical cancer. It is well established that individual HPV types vary in oncogenicity, but current data on their prognostic implication remain controversial. We examined the association between HPV types/species and the survival of 236 Chinese women aged 26–87 (mean 54.4) years after receiving primary treatment for cervical cancer. Overall, 45.8% were of FIGO stage I, 41.9% stage II, and 12.3% stage III. The four most prevalent types found were HPV-16 (60.2%), HPV-18 (21.6%), HPV-52 (11.9%), and HPV-58 (9.3%). Overall, 19.5% of patients had multiple-type infections, 78.4% harboured one or more alpha-9 species, and 28.8% harboured one or more alpha-7 species. After a median follow-up of 8.0 years, 156 (66.1%) patients survived. The 3-year overall survival rate was 75.5%. Factors independently associated with a poorer 3-year overall survival were age >60 years, tumour size >4 cm, lymph node involvement and treatment with radiotherapy+/-chemotherapy. Univariate analysis showed HPV-16 single-type infection was associated with a marginally poorer disease-specific survival (71.6% vs. 87.0%, HR: 1.71, 95% CI = 1.01–2.90), whereas non-HPV-16 alpha-9 species was associated with a better disease-specific survival (90.0% vs. 76.2%, HR: 0.36, 95% CI = 0.16–0.79). However, on multivariate analysis, HPV infection status irrespective of different grouping methods, including individual types, species, single-type or co-infection, did not carry any significant prognostic significance. In conclusion, we did not observe any association between infection with a particular HPV type/species and survival. An HPV type-based stratification in treatment and follow-up plan could not be recommended.

Human papillomavirus types 16 and 18 and the prognosis of patients with stage I cervical cancer

OBJECTIVE

This study sought to evaluate the prevalence of human papillomavirus (HPV) types 16 and 18 in women with clinical stage IB cervical cancer treated by radical hysterectomy with pelvic lymphadenectomy as well as to establish a correlation between HPV type and cancer prognosis.

METHODS

A single-center cohort study was conducted with 86 patients who had undergone radical hysterectomy for stage I cervical cancer. Prognostic factors and the presence of HPV 16 and 18 were analyzed using a polymerase chain reaction assay. A univariate analysis using Kaplan-Meier curves was conducted to estimate survival.

RESULTS

The prevalence of HPV 16 in the study group was 65.3%, and the prevalence of HPV 18 was 33.3%. The prevalence of infection with both viruses was 26.9%. Overall survival at 5 years was 91% among women with HPV 18 and 96% among those without this virus type (p=0.133). Among the women with HPV 16, the overall survival was 94%, whereas this rate was 96% among those without this virus type (p=0.663). Disease-free survival was unaffected by the presence of HPV type 16 or 18.

CONCLUSION

In the present study, despite the high prevalence of HPV types 16 and 18, the presence of these virus types did not affect the prognosis of patients with stage I cervical cancer who underwent radical hysterectomy.

Does HPV type affect outcome in oropharyngeal cancer?

Background

An epidemic of human papillomavirus (HPV)-related oropharyngeal squamous cell cancer (OPSCC) has been reported worldwide largely due to oral infection with HPV type-16, which is responsible for approximately 90% of HPV-positive cases. The purpose of this study was to determine the rate of HPV-positive oropharyngeal cancer in Southwestern Ontario, Canada.

Methods

A retrospective search identified ninety-five patients diagnosed with OPSCC. Pre-treatment biopsy specimens were tested for p16 expression using immunohistochemistry and for HPV-16, HPV-18 and other high-risk subtypes, including 31,33,35,39,45,51,52,56,58,59,67,68, by real-time qPCR.

Results

Fifty-nine tumours (62%) were positive for p16 expression and fifty (53%) were positive for known high-risk HPV types. Of the latter, 45 tumors (90%) were identified as HPV-16 positive, and five tumors (10%) were positive for other high-risk HPV types (HPV-18 (2), HPV-67 (2), HPV-33 (1)). HPV status by qPCR and p16 expression were extremely tightly correlated (p < 0.001, Fishers exact test). Patients with HPV-positive tumors had improved 3-year overall (OS) and disease-free survival (DFS) compared to patients with HPV-negative tumors (90% vs 65%, p = 0.001; and 85% vs 49%, p = 0.005; respectively). HPV-16 related OPSCC presented with cervical metastases more frequently than other high-risk HPV types (p = 0.005) and poorer disease-free survival was observed, although this was not statistically significant.

Conclusion

HPV-16 infection is responsible for a significant proportion of OPSCC in Southwestern Ontario. Other high-risk subtypes are responsible for a smaller subset of OPSCC that present less frequently with cervical metastases and may have a different prognosis.

Cervical cancer survival for patients referred to a tertiary care center in Kentucky

OBJECTIVES

To identify prognostic factors influencing cervical cancer survival for patients referred to a tertiary care center in Kentucky.

METHODS

A cohort study was performed to assess predictive survival factors of cervical cancer patients referred to the University of Kentucky from January 2001 to May 2010. Eligibility criteria included those at least 18 years-old, cervical cancer history, and no prior malignancy. Descriptive statistics were compiled and univariable and multivariable Cox proportional hazard analysis were performed.

RESULTS

381 patients met entry criteria. 95% were Caucasian (N=347) and 66% (N=243) lived in Appalachian Kentucky. The following covariates showed no evidence of a statistical association with survival: race, body mass index, residence, insurance status, months between last normal cervical cytology and diagnosis, histology, tumor grade, and location of primary radiation treatment. After controlling for identified significant variables, stage of disease was a significant predictor of overall survival, with estimated relative hazards comparing stages II, III, and IV to stage I of 3.09 (95% CI: 1.30, 7.33), 18.11 (95% CI: 7.44, 44.06), and 53.03(95% CI: 18.16, 154.87), respectively. The presence of more than two comorbid risk factors and unemployment was also correlated with overall survival [HR 4.25 (95% CI: 1.00, 18.13); HR 2.64 (95% CI 1.29, 5.42), respectively].

CONCLUSIONS

Residence and location of treatment center are not an important factor in cervical cancer survival when a tertiary cancer center can oversee and coordinate care; however, comorbid risk factors influence survival and further exploration of disease comorbidity related to cervical cancer survival is warranted.

A pilot analytic study of a research-level, lower-cost human papillomavirus 16, 18, and 45 test

The analytic performance of a low-cost, research-stage DNA test for the most carcinogenic human papillomavirus (HPV) genotypes (HPV16, HPV18, and HPV45) in aggregate was evaluated among carcinogenic HPV-positive women, which might be used to decide who needs immediate colposcopy in low-resource settings ("triage test"). We found that HPV16/18/45 test agreed well with two DNA tests, a GP5+/6+ genotyping assay (Kappa = 0.77) and a quantitative PCR assay (at a cutpoint of 5000 viral copies) (Kappa = 0.87). DNA sequencing on a subset of 16 HPV16/18/45 positive and 16 HPV16/18/45 negative verified the analytic specificity of the research test. It is concluded that the HPV16/18/45 assay is a promising triage test with a minimum detection of approximately 5000 viral copies, the clinically relevant threshold.

TP53, MDM2, NQO1, and susceptibility to cervical cancer

Host genetic variability modifies the risk of cervical cancer in women infected with oncogenic human papillomavirus (HPV). Studies have reported an association of the TP53 codon 72 arginine and cervical cancer, but the results are inconsistent. We examined the association of this single nucleotide polymorphism (SNP) in women with cervical cancer and cervical intraepithelial neoplasia grade 3, using a family-based association test. We further explored SNPs in two genes that regulate p53 stability: MDM2 (SNP309) and NQO1 (SNP609, SNP465). We also examined the relationship between host genotype and tumor HPV type. We genotyped 577 patients and their biological parents and/or siblings, using PCR-RFLP or Taqman assays. HPVs were typed by sequence-based methods. The transmission/disequilibrium test was used to detect disease-susceptibility alleles. The arginine peptide of TP53 codon 72 was overtransmitted in Caucasian families (P = 0.043), and the significance of this finding was enhanced in a subgroup of women infected with HPV16- and/or HPV18-related HPVs (P = 0.026). Allele C of NQO1 SNP609 was also overtransmitted in all cases (P = 0.026). We found no association between MDM2 SNP309 or NQO1 SNP465 and cervical cancer. Our results indicate that functional polymorphisms in TP53 codon 72 and NQO1 SNP609 associate with the risk of cervical cancer especially in women infected with type 16- and/or type 18-related HPVs.

Smoking and survival among Kentucky women diagnosed with invasive cervical cancer: 1995-2005

OBJECTIVES

To determine whether smokers with cervical cancer were more likely to die from cervical cancer compared with non smokers after adjusting for confounding factors.

METHODS

A population-based survival analysis was conducted among 2661 women diagnosed with invasive cervical cancer and reported to the Kentucky Cancer Registry from 1995-2005 and linked with state vital records and the National Death Index through 12/31/2005. A standard Kaplan - Meier approach was used in this survival analysis and Cox Proportional Hazards modeling was used to estimate adjusted hazard [aHR] ratios and 95% confidence intervals [CI] for smoking and survival for all cause and cervical cancer specific cause of death.

RESULTS

Almost half of women diagnosed with cervical cancer (48.6%) were known to be current smokers based the medical record review and reporting to KCR. For another 19.4% no tobacco status was documented (missing) and 32.1% were known non smokers. After adjustment for age and stage at diagnosis, cell type, rural residence, race, insurance coverage, and treatment received, current smoker were 35% more likely to die of any cause (aHR=1.35; 95% CI=1.17-1.56) and 21% more likely to die of cervical cancer (aHR=1.21; 95% CI=1.01-1.46) compared with known non smoking cases.

CONCLUSION

These data strongly suggest that smoking reduces cervical cancer survival.

Ethnic Disparities in Cervical Cancer Survival Among Medicare Eligible Women in a Multiethnic Population

Objectives:

To determine predictors of cervical cancer survival by socioeconomic status (SES), urbanization, race/ethnicity, comorbid conditions, and treatment among elderly Medicare-eligible women whose conditions were diagnosed with cervical cancer in a multiethnic population.

Methods:

A total of 538 women with cervical cancer aged 65 years or older were identified from 1999 to 2001 from the Texas Cancer Registry and were linked with the state Medicare data and Texas Vital Records to determine survival times. All women had similar access to care through Medicare fee-for-services insurance. A composite measure of SES was created using census tract-level data as was urbanization. Treatment and comorbid conditions were available from the Medicare data. Cox proportional hazards modeling was used for all-cause and cervical cancer-specific survival analysis.

Results:

Increased age (P < 0.0001) and advanced tumor stage (P < 0.0001) were associated with poorer all-cause and cervical cancer-specific survival. Having a comorbid condition was associated with all-cause survival (P < 0.01) but not cervical cancer-specific mortality. After adjusting for confounders, women receiving some form of treatment were almost half as likely to die with cervical cancer (adjusted hazard ratio = 0.68; 95% confidence interval, 0.52-0.89). After adjustment for all confounders, Hispanic women consistently had lower all-cause and cervical cancer-specific mortality rates relative to non-Hispanic white and non-Hispanic black women.

Conclusions:

Among women with similar health care coverage, Hispanic women had consistently lower all-cause and cervical cancer-specific mortality rates than other older women whose conditions were diagnosed with this disease in Texas. The presence of comorbid conditions and treatment were important predictors of survival, yet these factors do not explain the survival advantage for Hispanic women.

CD83 gene polymorphisms increase susceptibility to human invasive cervical cancer

We previously mapped a nonrandom frequent loss of heterozygosity (LOH) region in cervical cancers to 1 Mb of 6p23. Here, we describe the identification of a novel cervical cancer susceptibility gene, CD83. The gene was identified by several complementary approaches, including a family-based association study, comparison of transcript expression in normal and cancerous tissue, and genomic sequencing of candidate. CD83 encodes an inducible glycoprotein in the immunoglobulin superfamily and is a marker for mature dendritic cells. The association study that includes 377 family trios showed that five single nucleotide polymorphisms (SNP) within 8 kb of its 3'-end showed significant allelic association that was strengthened in a subgroup of women with invasive cancers infected by high-risk human papillomavirus type 16 and 18 (rs9296925, P = 0.0193; rs853360, P = 0.0035; rs9230, P = 0.0011; rs9370729, P = 0.0012; rs750749, P = 0.0133). Investigation of CD83 uncovered three alternative transcripts in cervical tissue and cell lines, with variant 3 (lacking exons 3 and 4) being more frequent in cervical cancer than in normal cervical epithelium (P = 0.0181). Genomic sequencing on 36 paired normal and cervical tumors revealed several somatic mutations and novel SNPs in the promoter, exons, and introns of CD83. LOH was confirmed in >90% of cervical cancer specimens. Immunofluorescence colocalized CD83 protein to the Golgi apparatus and cell membrane of cervical cancer cell lines. None of seven nearby genes was differentially expressed in cervical cancer. The importance of CD83 in epithelial versus dendritic cells needs to be determined, as does its role in promoting cervical cancer.