Genomic landscape of human papillomavirus-associated cancers

Impact of human papillomavirus types on uterine cervical neoplasia

Human papillomavirus (HPV) is a major cause of cervical cancer. As the natural history of HPV-associated cervical lesions is HPV genotype-dependent, it is important to understand the characteristics of these genotypes and to manage them accordingly. Among high-risk HPVs, HPV16 and 18 are particularly aggressive, together accounting for 70% of HPV genotypes detected in cervical cancer. Other than HPV16 and 18, HPV31, 33, 35, 45, 52, and 58 are also at a high risk of progression to cervical intraepithelial neoplasia (CIN)3 or higher. Recent studies have shown that the natural history of HPV16, 18, 52, and 58, which are frequently detected in Japan, depends on the HPV genotype. For example, HPV16 tends to progress in a stepwise fashion from CIN1 to CIN3, while HPV52 and 58 are more likely to persist in the CIN1 to CIN2 state. Among the high-risk HPVs, HPV18 has some peculiar characteristics different from those of other high-risk HPV types; the detection rate in precancerous lesions is much lower than those of other high-risk HPVs, and it is frequently detected in highly malignant adenocarcinoma and small cell carcinoma. Recent findings demonstrate that HPV18 may be characterized by latent infection and carcinogenesis in stem cell-like cells. In this context, this review outlines the natural history of HPV-infected cervical lesions and the characteristics of each HPV genotype.

Advancements in the Understanding of Small-Cell Neuroendocrine Cervical Cancer: Where We Stand and What Lies Ahead

Small-cell neuroendocrine cervical carcinoma (SCNCC) is a rare yet aggressive gynecological malignancy associated with dismal clinical outcomes. Its rarity has led to a limited number of retrospective studies and an absence of prospective research, posing significant challenges for evidence-based treatment approaches. As a result, most gynecologic oncology centers have limited experience with this tumor, emphasizing the urgent need for a comprehensive review and summary. This article systematically reviews the pathogenesis, immunohistochemical and molecular characteristics, prognostic factors, and clinical management of gynecologic SCNCC. We specifically focused on reviewing the distinct genomic characteristics of SCNCC identified via next-generation sequencing technologies, including loss of heterozygosity (LOH), somatic mutations, structural variations (SVs), and microRNA alterations. The identification of these actionable genomic events offers promise for discovering new molecular targets for drug development and enhancing therapeutic outcomes. Additionally, we delve deeper into key clinical challenges, such as determining the optimal treatment modality between chemoradiation and surgery for International Federation of Gynecology and Obstetrics (FIGO) stage I phase patients within a precision stratification framework, as well as the role of targeted therapy within the homologous recombination (HR) pathway, immune checkpoint inhibitors (ICIs), and prophylactic cranial irradiation (PCI) in the management of SCNCC. Finally, we anticipate the utilization of multiple SCNCC models, including cancer tissue-originated spheroid (CTOS) lines and patient-derived xenografts (PDXs), to decipher driver events and develop individualized therapeutic strategies for clinical application.

Human Papillomavirus-Associated Giant Clear Cell Acanthoma and Squamous Cell Carcinoma: A Rare Case Report and Literature Review

Clear cell acanthoma (CCA) and squamous cell carcinoma (SCC) represent distinct entities within dermatological oncology, each posing unique diagnostic and therapeutic challenges. CCA is a rare, benign epidermal growth, often not associated with human papillomavirus (HPV) infection, whereas SCC, a more aggressive form of skin cancer, has been linked to both ultraviolet (UV) exposure and HPV. Understanding the co-occurrence of these conditions in a single patient can enhance diagnostic accuracy and therapeutic outcomes. We report a 64-year-old male who underwent an operation for a verruciform lesion in the right groin, which was diagnosed as HPV-positive CCA alongside keratinised SCC. A literature search across January 2024 revealed limited evidence directly linking HPV to CCA, suggesting a need for further investigation. The speculative association between HPV and CCA warrants deeper exploration, especially considering the potential for HPV to contribute to lesion development through indirect mechanisms. The coexistence of CCA and SCC in an elderly patient presents a unique clinical scenario. This emphasises the need for vigilant diagnosis and tailored treatment strategies, highlighting the gap in understanding the pathogenesis of CCA, particularly its potential association with HPV. Further research is crucial for elucidating the complex interactions governing these conditions and for developing targeted interventions.

HPV DNA Integration at Actionable Cancer-Related Genes Loci in HPV-Associated Carcinomas

In HPV-associated carcinomas, some examples of cancer-related genes altered by viral insertion and corresponding to potential therapeutic targets have been described, but no quantitative assessment of these events, including poorly recurrent targets, has been reported to date. To document these occurrences, we built and analyzed a database comprised of 1455 cases, including HPV genotypes and tumor localizations. Host DNA sequences targeted by viral integration were classified as "non-recurrent" (one single reported case; 838 loci), "weakly recurrent" (two reported cases; 82 loci), and highly recurrent (≥3 cases; 43 loci). Whereas the overall rate of cancer-related target genes was 3.3% in the Gencode database, this rate increased to 6.5% in "non-recurrent", 11.4% in "weakly recurrent", and 40.1% in "highly recurrent" genes targeted by integration (p = 4.9 × 10-4). This rate was also significantly higher in tumors associated with high-risk HPV16/18/45 than other genotypes. Among the genes targeted by HPV insertion, 30.2% corresponded to direct or indirect druggable targets, a rate rising to 50% in "highly recurrent" targets. Using data from the literature and the DepMap 23Q4 release database, we found that genes targeted by viral insertion could be new candidates potentially involved in HPV-associated oncogenesis. A more systematic characterization of HPV/host fusion DNA sequences in HPV-associated cancers should provide a better knowledge of HPV-driven carcinogenesis and favor the development of personalize patient treatments.

Evidence of zoonotic pathogens through biophysically induced genomic variance

Zoonoses are infectious agents that are transmissible between animals and humans. Up to 60% of known infectious diseases and 75% of emergent diseases are zoonotic. Genomic variation between homeostatic populations provides a novel window into the effect of environmental pathogens on allelic distributions within the populations. Genodynamics is a biophysical approach utilizing developed metrics on biallelic single-nucleotide polymorphisms (SNPs) that can be used to quantify the adaptive influences due to pathogens. A genomic free energy that is minimized when overall population health is optimized describes the influence of environmental agents upon genomic variation. A double-blind exploration of over 100 thousand SNPs searching for smooth functional dependencies upon four zoonotic pathogens carried by four possible hosts amidst populations that live in their ancestral environments has been conducted. Exemplars that infectious agents can have significant adaptive influence on human populations are presented. One discussed SNP is likely associated with both adaptive and innate immune regulation. The adaptive response of another SNP suggests an intriguing connection between zoonoses and human cancers. The adaptive forces of the presented pathogens upon the human genome have been quantified.

Multiple HPV integration mode in the cell lines based on long-reads sequencing

Background

The integration of human papillomavirus (HPV) is closely related to the occurrence of cervical cancer. However, little is known about the complete state of HPV integration into the host genome.

Methods

In this study, three HPV-positive cell lines, HeLa, SiHa, and CaSki, were subjected to NANOPORE long-read sequencing to detect HPV integration. Analysis of viral integration patterns using independently developed software (HPV-TSD) yielded multiple complete integration patterns for the three HPV cell lines.

Results

We found distinct differences between the integration patterns of HPV18 and HPV16. Furthermore, the integration characteristics of the viruses were significantly different, even though they all belonged to HPV16 integration. The HPV integration in the CaSki cells was relatively complex. The HPV18 integration status in HeLa cells was the dominant, whereas the percentage of integrated HPV 16 in SiHa and CaSki cells was significantly lower. In addition, the virus sequences in the HeLa cells were incomplete and existed in an integrated state. We also identified a large number of tandem repeats in HPV16 and HPV18 integration. Our study not only clarified the feasibility of high-throughput long-read sequencing in the study of HPV integration, but also explored a variety of HPV integration models, and confirmed that viral integration is an important form of HPV in cell lines.

Conclusion

Elucidating HPV integration patterns will provide critical guidance for developing a detection algorithm for HPV integration, as well as the application of virus integration in clinical practice and drug research and development.

Identify and validate circadian regulators as potential prognostic markers and immune infiltrates in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneity pathological malignant cancer with leading causes of morbidity and mortality. EGFR inhibitors, immune checkpoint inhibitors have become novel treatments. However, the mechanism still remained uncertain. Several studies have confirmed that the circadian rhythms induce multiple malignancies developing. We utilized multi-omics analysis to demonstrate the crosstalk between circadian clock genes and tumor microenvironment in HNSCC. Firstly, we performed the LASSO Cox regression analysis based on the 16 important clock genes. A 7-gene risk model was successfully established in TCGA and validated in GEO datasets. Next, CIBERSORT and ESTIMATE methods were performed to display the immune landscape of high risk and low risk groups, and the results showed that high abundance of mast cells activated, dendritic cells activated and neutrophils were positively correlated with poor OS. To further identify hub genes, Kaplan Meier plot was applied in all TCGA and GEO datasets and two hub genes (PER2, and PER3) were identified, especially PER3, which was found strongly associated with immune score, PDCD1, CD4 + and CD8 + T cells in HNSCCC. Moreover, to explore the innate mechanism of circadian-induced pathway, we constructed a circadian-related ceRNA regulatory network containing 34 lncRNAs, 3 miRNAs and 4 core circadian genes. In-vitro experiments also verified that Per2 or Per3 could suppressed the proliferation, migration and invasion of HNSC. This study unraveled the association between PER3 and prognosis in patients with HNSC and the innate mechanism remains to be elucidated.

Whole exome sequencing identifies common mutational landscape of cervix and endometrium small cell neuroendocrine carcinoma

Background

Primary small cell neuroendocrine carcinomas of the cervix and endometrium are rare gynecological malignancies with limited treatment options. This study aimed to improve the understanding of the carcinogenesis process and identify potential therapeutic targets for these two tumor types by constructing the mutational landscape at the whole exome level.

Methods

Primary tumor tissues and their matched blood samples were obtained from 10 patients with small cell cervical neuroendocrine carcinoma (NECC) and five patients with small cell endometrial neuroendocrine carcinoma (NECE). Whole exome sequencing was performed to construct the somatic mutation profiles. Mutational signature and recurrent mutated gene analysis were used to identify tumor subtypes and common carcinogenesis processes.

Results

Based on the burden of different mutational signatures, the NECCs in this work can be divided into two subtypes, including the mismatch repair deficiency like (dMMR-like) type (4/10) and the high spontaneous deamination type (6/10). Components of the PI3K/AKT signaling and RAS signaling were exclusively mutated in these two subtypes, respectively. The integration of human papillomavirus made a limited contribution to tumorigenesis in NECC (20%). The dysfunction of the mismatch repair system and microsatellite instability are the major features of NECE. PI3K/AKT, JAK/STAT signaling, and chromatin remodeling activity were the common mutated pathways in NECE. PIK3CA, WNK2, and KMT2B underwent mutations in both the dMMR-like subtype of NECC (50% - 75%) and in NECE (60% - 80%) specimens, while exhibiting infrequent mutational occurrences in publicly available data pertaining to neuroendocrine carcinomas of the lung or bladder (< 10%).

Conclusion

We identified the two subtypes of NECC with distinct mutated pathways and potential therapy targets. The dMMR-like type NECC and NECE may share a similar carcinogenesis process that include dysfunction of PI3K/AKT signaling, cell cycle, antiapoptotic processes, and chromatin remodeling activity.

Protein Profiling in Human Papillomavirus-Associated Cervical Carcinogenesis: Cornulin as a Biomarker for Disease Progression

Nearly 90% of cervical cancers are linked to human papillomavirus (HPV). Uncovering the protein signatures in each histological phase of cervical oncogenesis provides a path to biomarker discovery. The proteomes extracted from formalin-fixed paraffin-embedded tissues of the normal cervix, HPV16/18-associated squamous intraepithelial lesion (SIL), and squamous cell carcinoma (SCC) were compared using liquid chromatography-mass spectrometry (LC-MS). A total of 3597 proteins were identified, with 589, 550, and 1570 proteins unique to the normal cervix, SIL, and SCC groups, respectively, while 332 proteins overlapped between the three groups. In the transition from normal cervix to SIL, all 39 differentially expressed proteins were downregulated, while all 51 proteins discovered were upregulated in SIL to SCC. The binding process was the top molecular function, while chromatin silencing in the SIL vs. normal group, and nucleosome assembly in SCC vs. SIL groups was the top biological process. The PI3 kinase pathway appears crucial in initiating neoplastic transformation, while viral carcinogenesis and necroptosis are important for cell proliferation, migration, and metastasis in cervical cancer development. Annexin A2 and cornulin were selected for validation based on LC-MS results. The former was downregulated in the SIL vs. normal cervix and upregulated in the progression from SIL to SCC. In contrast, cornulin exhibited the highest expression in the normal cervix and lowest in SCC. Although other proteins, such as histones, collagen, and vimentin, were differentially expressed, their ubiquitous expression in most cells precluded further analysis. Immunohistochemical analysis of tissue microarrays found no significant difference in Annexin A2 expression between the groups. Conversely, cornulin exhibited the strongest expression in the normal cervix and lowest in SCC, supporting its role as a tumor suppressor and potential biomarker for disease progression.

Long-read sequencing reveals oncogenic mechanism of HPV-human fusion transcripts in cervical cancer

Integration of high-risk human papillomavirus (HPV) into the host genome is a crucial event for the development of cervical cancer, however, the underlying mechanism of HPV integration-driven carcinogenesis remains unknown. Here, we performed long-read RNA sequencing on 12 high-grade squamous intraepithelial lesions (HSIL) and cervical cancer patients, including 3 pairs of cervical cancer and corresponding para-cancerous tissue samples to investigate the full-length landscape of cross-species genome integrations. In addition to massive unannotated isoforms, transcriptional regulatory events, and gene chimerism, more importantly, we found that HPV-human fusion events were prevalent in HPV-associated cervical cancers. Combined with the genome data, we revealed the existence of a universal transcription pattern in these fusion events, whereby structurally similar fusion transcripts were generated by specific splicing in E6 and a canonical splicing donor site in E1 linking to various human splicing acceptors. Highly expressed HPV-human fusion transcripts, eg, HPV16 E6*I-E7-E1_SD880-human gene, were the key driver of cervical carcinogenesis, which could trigger overexpression of E6*I and E7, and destroy the transcription of tumor suppressor genes CMAHP, TP63 and P3H2. Finally, evidence from in vitro and in vivo experiments demonstrates that the novel read-through fusion gene mRNA, E1-CMAHP (E1C, formed by the integration of HPV58 E1 with CMAHP), existed in the fusion transcript can promote malignant transformation of cervical epithelial cells via regulating downstream oncogenes to participate in various biological processes. Taken together, we reveal a previously unknown mechanism of HPV integration-driven carcinogenesis and provide a novel target for the diagnosis and treatment of cervical cancer.

An integrative approach toward identification and analysis of therapeutic targets involved in HPV pathogenesis with a focus on carcinomas

BACKGROUND

Persistent infection of high-risk HPVs is known to cause diverse carcinomas, mainly cervical, oropharyngeal, penile, etc. However, efficient treatment is still lacking.

OBJECTIVE

Identify and analyze potential therapeutic targets involved in HPV oncogenesis and repurposing drug candidates.

METHODS

Integrative analyses were performed on the compendium of 1887 HPV infection-associated or integration-driven disrupted genes cataloged from the Open Targets Platform and HPVbase resource. Potential target genes are prioritized using STRING, Cytoscape, cytoHubba, and MCODE. Gene ontology and KEGG pathway enrichment analysis are performed. Further, TCGA cancer genomic data of CESC and HNSCC is analyzed. Moreover, regulatory networks are also deduced by employing NetworkAnalyst.

RESULTS

We have implemented a unique approach for identifying and prioritizing druggable targets and repurposing drug candidates against HPV oncogenesis. Overall, hundred key genes with 44 core targets were prioritized with transcription factors (TFs) and microRNAs (miRNAs) regulators pertinent to HPV pathogenesis. Genomic alteration profiling further substantiated our findings. Among identified druggable targets, TP53, NOTCH1, PIK3CA, EP300, CREBBP, EGFR, ERBB2, PTEN, and FN1 are frequently mutated in CESC and HNSCC. Furthermore, PIK3CA, CCND1, RFC4, KAT5, MYC, PTK2, EGFR, and ERBB2 show significant copy number gain, and FN1, CHEK1, CUL1, EZH2, NRAS, and H2AFX was marked for the substantial copy number loss in both carcinomas. Likewise, under-explored relevant regulators, i.e., TFs (HINFP, ARID3A, NFATC2, NKX3-2, EN1) and miRNAs (has-mir-98-5p, has-mir-24-3p, has-mir-192-5p, has-mir-519d-3p) is also identified.

CONCLUSIONS

We have identified potential therapeutic targets, transcriptional and post-transcriptional regulators to explicate HPV pathogenesis as well as potential repurposing drug candidates. This study would aid in biomarker and drug discovery against HPV-mediated carcinoma.

Multi-region sequencing depicts intratumor heterogeneity and clonal evolution in cervical cancer

Cervical cancer is a heterogeneous malignancy mainly caused by human papillomavirus (HPV). While a few studies have revealed heterogeneity of cervical cancer in chromosome levels, the correlation between genetic heterogeneity and HPV integration in cervical cancer remains unknown. Here, we applied multi-region whole-exome sequencing and HPV integration analysis to explore intratumor heterogeneity in cervical cancer. We sequenced 20 tumor regions and 5 adjacent normal tissues from 5 cervical cancer patients, analysis based on somatic mutations and somatic copy number alterations (SCNAs) levels were performed. Variable heterogeneity was observed between the five patients with different tumor stages and HPV infection statuses. We found HPV integration has a positive effect on somatic mutation burden, but the relation to SCNAs remains unclear. Frequently mutated genes in cervical cancer were identified as trunk events, such as FBXW7, PIK3CA, FAT1 in somatic mutations and TP63, MECOM, PIK3CA, TBL1XR1 in SCNAs. New potential driver genes in cervical cancer were summarized including POU2F2, TCF7 and UBE2A. The SCNAs level has potential relation with tumor stage, and Signature 3 related to homologous recombination deficiency may be the appropriate biomarker in advanced cervical cancer. Mutation signature analysis also revealed a potential pattern that APOBEC-associated signature occurs in early stage and signatures associated with DNA damage repair arise at the later stage of cervical cancer evolution. In a conclusion, our study provides insights into the potential relationship between HPV infection and tumor heterogeneity. Those results enhanced our understanding of tumorigenesis and progression in cervical cancer.

HPV and radiosensitivity of cervical cancer: a narrative review

Background and Objective

Cervical cancer (CC), the most common gynecological malignancy, is divided into two categories: human papillomavirus-related [HPV positive (HPV+)] and non-HPV-related [HPV negative (HPV-)]. Compared with HPV- CC, HPV+ CC has better radiosensitivity and prognosis. We conducted a literature search and summarized relevant studies to explore the detailed mechanisms by which HPV+ improves the prognosis of CC compared to HPV-.

Methods

PubMed was used to search the literature on human papillomavirus, cervical cancer, and radiotherapy up to June 2022.

Key Content and Findings

Compared with HPV- CC, HPV+ CC has better radiotherapy outcomes and better prognosis. HPV improves the radiotherapy sensitivity of CC by inhibiting damaged DNA repair, increasing cell cycle arrest, reducing hypoxia, increasing cellular immune response, and other mechanisms. However, the effect of HPV on radiotherapy sensitivity of CC is not consistent and is affected by HPV type, viral load, and many other factors. Partial HPV+ CCs, due to hypoxia and other factors, are resistant to radiotherapy and have a poor prognosis. HPV- CC has poor radiotherapy sensitivity and poor prognosis. With the spread of the vaccine, HPV- CC will gradually increase, which is a cause for concern.

Conclusions

The radiosensitivity was significantly increased in patients with HPV+ CC, compared to HPV- patients. HPV improves the radiotherapy sensitivity of cervical cancer through a number of pathways. Meanwhile, the relationship between HPV and radiotherapy sensitivity is influenced by a number of factors. Some HPV+ CCs showed radiotherapy resistance, and HPV- CCs deserve further attention.

Direct Comparison of HPV16 Viral Genomic Integration, Copy Loss, and Structural Variants in Oropharyngeal and Uterine Cervical Cancers Reveal Distinct Relationships to E2 Disruption and Somatic Alteration

Squamous cell carcinoma of the oropharynx caused by HPV type 16 (HPV16+ OPSCC) is the most common HPV-associated malignancy in the USA and has many molecular differences from uterine cervical squamous cell carcinoma (UCSCC). Our understanding of HPV oncogenesis relied on studies of UCSCC revealing a consensus model reliant on HPV integration with a loss of E2. Here, we compare patterns of HPV integration in UCSCC and OPSCC by analysis of affinity capture sequencing of the HPV16 genome in 104 OPSCC and 44 UCSCC tumors. These cohorts were contemporaneously sequenced using an identical strategy. Integration was identified using discordant read pair clustering and assembly-based approaches. Viral integration sites, structural variants, and copy losses were examined. While large-scale deep losses of HPV16 genes were common in UCSCC and were associated with E2 loss, deep copy losses of the HPV16 genome were infrequent in HPV16+ OPSCC. Similarly, structural variants within HPV16 favored E2 loss in UCSCC but not OPSCC. HPV16 integration sites were non-random, with recurrent integration hot-spots identified. OPSCC tumors had many more integration sites per tumor when compared to UCSCC and had more integration sites in genomic regions with high gene density. These data show that viral integration and E2 disruption are distinct in UCSCC and OPSCC. Our findings also add to growing literature suggesting that HPV tumorigenesis in OPSCC does not follow the model developed based on UCSCC.

ERBB2 promoter demethylation and immune cell infiltration promote a poor prognosis for cancer patients

Background

Receptor tyrosine-protein kinase erbB-2 (ERBB2) expression is a critical factor for the prognosis of various cancers. ERBB2 enrichment indicates a poor prognosis in some cancer types but could be a favorable prognostic factor in others.

Methods

We analyzed DNA methylation, mRNA, protein, immune cell infiltration, and related signaling pathways using TIMER2.0, GEPIA2, STRING, and UALCAN portal datasets in tumor tissues of diverse cancer types and their matched normal tissues.

Results

ERBB2 promoter demethylation increases transcript protein amplification and promotes a poor prognosis for cancer patients. ERBB2 gain-of-function procures immune cell infiltration for tumor growth and drives away T regulatory cells, which suppress or downregulate induction and proliferation of effector T cells. The downstream signaling pathways, such as tumor proliferation, ECM-related genes, and degradation of ECM, are involved in ERBB2 gene demethylation and immune activation in cancer progression.

Conclusion

ERBB2 gene demethylation leads to a poor prognosis in cancer patients, which is strongly influenced by the composition and abundance of tumor immune cell infiltration. ERBB2 demethylation could be used in clinical practice to identify immune profiles and direct therapeutic strategies.

Differences in integration frequencies and APOBEC3 profiles of five high-risk HPV types adheres to phylogeny

Persistent infection with Human Papillomavirus (HPV) is responsible for almost all cases of cervical cancers, and HPV16 and HPV18 associated with the majority of these. These types differ in the proportion of viral minor nucleotide variants (MNVs) caused by APOBEC3 mutagenesis as well as integration frequencies. Whether these traits extend to other types remains uncertain. This study aimed to investigate and compare genomic variability and chromosomal integration in the two phylogenetically distinct Alpha-7 and Alpha-9 clades of carcinogenic HPV types. The TaME-seq protocol was employed to sequence cervical cell samples positive for HPV31, HPV33 or HPV45 and combine these with data from a previous study on HPV16 and HPV18. APOBEC3 mutation signatures were found in Alpha-9 (HPV16/31/33) but not in Alpha-7 (HPV18/45). HPV45 had significantly more MNVs compared to the other types. Alpha-7 had higher integration frequency compared to Alpha-9. An increase in integration frequency with increased diagnostic severity was found for Alpha-7. The results highlight important differences and broaden our understanding of the molecular mechanisms behind cervical cancer induced by high-risk HPV types from the Alpha-7 and Alpha-9 clades.

Integration of scRNA-Seq and TCGA RNA-Seq to Analyze the Heterogeneity of HPV+ and HPV- Cervical Cancer Immune Cells and Establish Molecular Risk Models

Background

Numerous studies support that Human papillomavirus (HPV) can cause cervical cancer. However, few studies have surveyed the heterogeneity of HPV infected or uninfected (HPV+ and HPV-) cervical cancer (CESC) patients. Integration of scRNA-seq and TCGA data to analyze the heterogeneity of HPV+ and HPV- cervical cancer patients on a single-cell level could improve understanding of the cellular mechanisms during HPV-induced cervical cancer.

Methods

CESC scRNA-seq data obtained from the Gene Expression Omnibus (GEO) database and the Seurat, Monocle3 package were used for scRNA-seq data analysis. The ESTIMATE package was used for single-sample gene immune score, CIBERSORT package was used to identify immune scores of cells, and the “WGCNA” package for the weighted correlation network analysis. Univariate Cox and LASSO regression were performed to establish survival and relapse signatures. KEGG and GO analyses were performed for the signature gene. Gene Expression Profiling Interactive Analysis was used for Pan-cancer analysis.

Results

In the HPV+ CESC group, CD8+ T cells and B cells were down-regulated, whereas T reg cells, CD4+ T cells, and epithelial cells were up-regulated according to scRNA-seq data. Survival analysis of TCGA-CESC revealed that increased expression of naive B cells or CD8+ T cells favors the survival probability of CESC patients. WGCNA, univariate Cox, and LASSO Cox regression established a 9-genes survival signature and a 7-gene relapse model. Pan-cancer analysis identified IKZF3, FOXP3, and JAK3 had a similar distribution and effects in HPV-associated HNSC.

Conclusion

Analysis of scRNA-seq and bulk RNA-seq of HPV+ and HPV- CESC samples revealed heterogeneity from transcriptional state to immune infiltration. Survival and relapse models were adjusted according to the heterogeneity of HPV+ and HPV- CESC immune cells to assess the prognostic risk accurately. Hub genes represent similar protection in HPV- associated HNSC while showing irrelevant to other potential HPV-related cancers.

Construction of hypoxia-immune-related prognostic model and targeted therapeutic strategies for cervical cancer

Emerging evidence indicates that hypoxia and immunity play important roles in tumorigenesis and development. However, the hypoxia-immune-related prognostic risk model has not been established in cervical cancer (CC). We aimed to construct a hypoxia-immune-based prognostic risk model with potential application in CC patient prognosis and predicting response to targeted therapy. The RNA-seq data and corresponding clinical information were retrieved from The Cancer Genome Atlas (TCGA) database. Hypoxia and immune status of CC patients were evaluated using the Consensus Clustering method and single sample gene set enrichment analysis (ssGSEA), respectively. The univariate Cox regression, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were applied to establish the prognostic risk model of CC. The chemotherapy response for six chemotherapeutic agents of each CC patient was calculated according to the Genomics of Drug Sensitivity in Cancer (GDSC). And the Connectivity Map (CMap) database was performed to screen candidate small molecule drugs. In this study, we identified 7 gene signatures (P4HA2, MSMO1, EGLN1, ZNF316, IKZF3, ISCU, MYO1B) with prognostic values. And the survival time of patients with low-risk was significantly longer than those with high-risk. Meanwhile, CC patients in the high-risk group yielded higher sensitivity to five chemotherapeutic agents. And we listed ten candidate small-molecules drugs that exhibited a high correlation with the prognosis of cervical cancer. Thus, the prognostic model can accurately predict the prognosis of patients with CC and may be helpful for the development of new hypoxia-immune prognostic markers and therapeutic strategies for CC.

Long-read sequencing unveils high-resolution HPV integration and its oncogenic progression in cervical cancer

Integration of human papillomavirus (HPV) DNA into the human genome is considered as a key event in cervical carcinogenesis. Here, we perform comprehensive characterization of large-range virus-human integration events in 16 HPV16-positive cervical tumors using the Nanopore long-read sequencing technology. Four distinct integration types characterized by the integrated HPV DNA segments are identified with Type B being particularly notable as lacking E6/E7 genes. We further demonstrate that multiple clonal integration events are involved in the use of shared breakpoints, the induction of inter-chromosomal translocations and the formation of extrachromosomal circular virus-human hybrid structures. Combined with the corresponding RNA-seq data, we highlight LINC00290, LINC02500 and LENG9 as potential driver genes in cervical cancer. Finally, we reveal the spatial relationship of HPV integration and its various structural variations as well as their functional consequences in cervical cancer. These findings provide insight into HPV integration and its oncogenic progression in cervical cancer.

The molecular mechanisms underlying cervical carcinogenesis following integration of HPV DNA into the human genome remain elusive. Here, the authors perform long-read sequencing in 16 HPV16-positive cervical tumors and identify distinct integration types, structural variations and potential driver genes.

Cervical Precancers Originate From Infected Proliferating Reserve Cells: A Comparative Histologic and Genetic Study of Thin and Thick High-grade Squamous Intraepithelial Lesions

Human papillomavirus (HPV)-induced invasive cervical squamous cell cancer (SCC) develop via high-grade squamous intraepithelial lesion (HSIL). In contrast to classic thick HSIL, thin HSIL (≤9 cell layers) are poorly documented. This study compares histology, HPV genotypes, and aberrations in 50 cancer genes of 45 thin HSIL to 45 thick HSIL, 20 pT1a SCC, and 40 ≥pT1b SCC. Thin HSIL arose from proliferating reserve cells within endocervical epithelium or immature metaplasia throughout the transformation zone after infection with high-risk HPV genotypes (36/45; 80%), and 20% non-high-risk HPV genotypes compared with 2.5% thick HSIL, pT1a SCC, and ≥pT1b SCC. Thin HSIL were multifocal proliferations with varying epithelial thickness between 1 and 2 to 9 cell layers, with occasional transitions to thick HSIL or concomitant lesions of thick HSIL. Overall, 40% thin HSIL were located distant to and most thick HSIL occurred near or at the squamocolumnar junction. Only 20% thick HSIL showed koilocytosis. All HSIL lacked somatic gene mutations, compared with 30% pT1a and 55%≥pT1b SCC. Overrepresented rare germline variants in the MET, JAK3, and FGFR3 genes occurred in all patient groups. In summary, thin and thick HSIL arose independently of somatic gene mutations. The maturation level of the squamous epithelium at the time of transforming infection determines if a thick HSIL develops directly from HPV-infected proliferating reserve cells via thin HSIL or in stratified glycogenated squamous epithelium via low-grade squamous intraepithelial lesion. These observations raise doubts about the biological relevance of separation into thin and thick HSIL. The oncogenic potential of HPV genotypes but also germline variants may influence the natural history.

Correlation of SARS‑CoV‑2 to cancer: Carcinogenic or anticancer? (Review)

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is highly infectious and pathogenic. Among patients with severe SARS-CoV-2-caused by corona virus disease 2019 (COVID-19), those complicated with malignant tumor are vulnerable to COVID-19 due to compromised immune function caused by tumor depletion, malnutrition and anti-tumor treatment. Cancer is closely related to the risk of severe illness and mortality in patients with COVID-19. SARS-CoV-2 could promote tumor progression and stimulate metabolism switching in tumor cells to initiate tumor metabolic modes with higher productivity efficiency, such as glycolysis, for facilitating the massive replication of SARS-CoV-2. However, it has been shown that infection with SARS-CoV-2 leads to a delay in tumor progression of patients with natural killer cell (NK cell) lymphoma and Hodgkin's lymphoma, while SARS-CoV-2 elicited anti-tumor immune response may exert a potential oncolytic role in lymphoma patients. The present review briefly summarized potential carcinogenicity and oncolytic characteristics of SARS-CoV-2 as well as strategies to protect patients with cancer during the COVID-19 pandemic.

Whole Genome Assembly of Human Papillomavirus by Nanopore Long-Read Sequencing

Human papillomavirus (HPV) is a causal agent for most cervical cancers. The physical status of the HPV genome in these cancers could be episomal, integrated, or both. HPV integration could serve as a biomarker for clinical diagnosis, treatment, and prognosis. Although whole-genome sequencing by next-generation sequencing (NGS) technologies, such as the Illumina sequencing platform, have been used for detecting integrated HPV genome in cervical cancer, it faces challenges of analyzing long repeats and translocated sequences. In contrast, Oxford nanopore sequencing technology can generate ultra-long reads, which could be a very useful tool for determining HPV genome sequence and its physical status in cervical cancer. As a proof of concept, in this study, we completed whole genome sequencing from a cervical cancer tissue and a CaSki cell line with Oxford Nanopore Technologies. From the cervical cancer tissue, a 7,894 bp-long HPV35 genomic sequence was assembled from 678 reads at 97-fold coverage of HPV genome, sharing 99.96% identity with the HPV sequence obtained by Sanger sequencing. A 7904 bp-long HPV16 genomic sequence was assembled from data generated from the CaSki cell line at 3857-fold coverage, sharing 99.99% identity with the reference genome (NCBI: U89348). Intriguingly, long reads generated by nanopore sequencing directly revealed chimeric cellular-viral sequences and concatemeric genomic sequences, leading to the discovery of 448 unique integration breakpoints in the CaSki cell line and 60 breakpoints in the cervical cancer sample. Taken together, nanopore sequencing is a unique tool to identify HPV sequences and would shed light on the physical status of HPV genome in its associated cancers.

FGF signalling facilitates cervical cancer progression

Cervical cancer is one of the most frequently diagnosed cancers in women worldwide. While cervical cancer is caused by human papillomavirus (HPV), not all females infected with HPV develop the disease, suggesting that other factors might facilitate its progression. Growing evidence supports the involvement of the fibroblast growth factor receptor (FGFR) axis in several cancers, including gynecological. However, for cervical cancer, the molecular mechanisms that underpin the disease remain poorly understood, including the role of FGFR signaling. The aim of this study was to investigate FGF(R) signaling in cervical cancer through bioinformatic analysis of cell line and patient data and through detailed expression profiling, manipulation of the FGFR axis, and downstream phenotypic analysis in cell lines (HeLa, SiHa, and CaSki). Expression (protein and mRNA) analysis demonstrated that FGFR1b/c, FGFR2b/c, FGFR4, FGF2, FGF4, and FGF7 were expressed in all three lines. Interestingly, FGFR1 and 2 localized to the nucleus, supporting that nuclear FGFRs could act as transcription factors. Importantly, 2D and 3D cell cultures demonstrated that FGFR activation can facilitate cell functions correlated with invasive disease. Collectively, this study supports an association between FGFR signaling and cervical cancer progression, laying the foundations for the development of therapeutic approaches targeting FGFR in this disease.

Exome sequencing reveals candidate mutations implicated in sinonasal carcinoma and malignant transformation of sinonasal inverted papilloma

We explored somatic mutations in dysplastic sinonasal inverted papilloma (SNIP), SNIP with concomitant sinonasal squamous cell carcinoma (SNSCC), and SNSCC without preceding SNIP. Ten SNIP and SNSCC samples were analyzed with exome sequencing and tested for human papillomavirus. The identified mutations were compared to the most frequently mutated genes in head and neck squamous cell carcinoma (HNSCC) in the COSMIC database. Exome sequencing data were also analyzed for mutations not previously linked to SNSCC. Seven of the most commonly mutated genes in HNSCC and SNSCC in COSMIC harbored mutations in our data. In addition, we identified mutations in 23 genes that are likely to contribute to SNIP and SNSCC oncogenesis.

Genome-Wide Profiling Reveals HPV Integration Pattern and Activated Carcinogenic Pathways in Penile Squamous Cell Carcinoma

Simple Summary

Penile squamous cell carcinoma (PSCC) has been regarded as an HPV-related cancer for a long time. However, the integration pattern and carcinogenic pathways of HPV in PSCC remain unclear. The results of this study provide insights into the HPV-related carcinogenic mechanism in PSCC, which may be less prone to involvement in the traditional E6/E7 carcinogenic process, and are characterized by effects on the host genome, which result in the inactivation of tumor suppressors (CADM2, etc.) and the activation of oncogenes (KLF5, etc.), thus activating oncogenic signaling pathways (MAPK, JAK/STAT, etc.). This study could enhance our understanding of HPV integration and pave the way for subsequent HPV studies in PSCC.

Abstract

Human papillomavirus (HPV) is a significant etiologic driver of penile squamous cell carcinoma (PSCC). The integration pattern of HPV and its carcinogenic mechanism in PSCC remain largely unclear. We retrospectively reviewed 108 PSCC cases who received surgery between 2008 and 2017. Using high-throughput viral integration detection, we identified 35 HPV-integrated PSCCs. Unlike cervical cancer, the HPV E2 oncogene was not prone to involvement in integration. Eleven of the 35 (31.4%) HPV-integrated PSCCs harbored intact HPV E2; these tumors had lower HPV E6 and E7 expression and higher expression of p53 and pRb proteins than those with disrupted E2 did (p < 0.001 and p = 0.024). Integration breakpoints are preferentially distributed in or near host genes, including previously reported hotspots (KLF5, etc.) and newly identified hotspots (CADM2, etc.), which are mainly involved in oncogenic signaling pathways (MAPK, JAK/STAT, etc.). Regarding the phosphorylation levels of JNK, p38 was higher in HPV-positive tumors with MAPK-associated integration than those in HPV-positive tumors with other integration and those in HPV-negative tumors. In vitro, KLF5 knockdown inhibited proliferation and invasion of PSCC cells, while silencing CADM2 promoted migration and invasion. In conclusion, this study enhances our understanding of HPV-induced carcinogenesis in PSCC, which may not only rely on the E6/E7 oncogenes, but mat also affect the expression of critical genes and thus activate oncogenic pathways.

Development of New Cancer Treatment by Identifying and Focusing the Genetic Mutations or Altered Expression in Gynecologic Cancers

With the advent of next-generation sequencing (NGS), The Cancer Genome Atlas (TCGA) research network has given gynecologic cancers molecular classifications, which impacts clinical practice more and more. New cancer treatments that identify and target pathogenic abnormalities of genes have been in rapid development. The most prominent progress in gynecologic cancers is the clinical efficacy of poly(ADP-ribose) polymerase (PARP) inhibitors, which have shown breakthrough benefits in reducing hazard ratios (HRs) (HRs between 0.2 and 0.4) of progression or death from BRCA1/2 mutated ovarian cancer. Immune checkpoint inhibition is also promising in cancers that harbor mismatch repair deficiency (dMMR)/microsatellite instability (MSI). In this review, we focus on the druggable genetic alterations in gynecologic cancers by summarizing literature findings and completed and ongoing clinical trials.

Causes and Consequences of HPV Integration in Head and Neck Squamous Cell Carcinomas: State of the Art

A constantly increasing incidence in high-risk Human Papillomaviruses (HPV)s driven head and neck squamous cell carcinomas (HNSCC)s, especially of oropharyngeal origin, is being observed. During persistent infections, viral DNA integration into the host genome may occur. Studies are examining if the physical status of the virus (episomal vs. integration) affects carcinogenesis and eventually has further-reaching consequences on disease progression and outcome. Here, we review the literature of the most recent five years focusing on the impact of HPV integration in HNSCCs, covering aspects of detection techniques used (from PCR up to NGS approaches), integration loci identified, and associations with genomic and clinical data. The consequences of HPV integration in the human genome, including the methylation status and deregulation of genes involved in cell signaling pathways, immune evasion, and response to therapy, are also summarized.

The Role of Human Papilloma Virus in Dictating Outcomes in Head and Neck Squamous Cell Carcinoma

The Human Papilloma Virus (HPV) is an oncogenic virus which is associated with the development of head and neck squamous cell carcinoma (HNSCC), predominantly within the oropharynx. Approximately 25% of oropharyngeal squamous cell carcinoma (OPSCC) cases worldwide are attributable to HPV infection, with an estimated 65% in the United States. Transmission is via exposure during sexual contact, with distinctive anatomical features of the tonsils providing this organ with a predilection for infection by HPV. No premalignant lesion is identifiable on clinical examination, thus no comparative histological features to denote the stages of carcinogenesis for HPV driven HNSCC are identifiable. This is in contrast to HPV-driven cervical carcinoma, making screening a challenge for the head and neck region. However, HPV proffers a favorable prognosis in the head and neck region, with better overall survival rates in contrast to its HPV negative counterparts. This has resulted in extensive research into de-intensifying therapies aiming to minimize the morbidity induced by standard concurrent chemo-radiotherapy without compromising efficacy. Despite the favorable prognosis, cases of recurrence and/or metastasis of HPV positive HNSCC do occur, and are linked with poor outcomes. HPV 16 is the most frequent genotype identified in HNSCC, yet there is limited research to date studying the impact of other HPV genotype with respect to overall survival. A similar situation pertains to genetic aberrations associated in those with HPV positive HNSCC who recur, with only four published studies to date. Somatic mutations in TSC2, BRIP1, NBN, TACC3, NFE2l2, STK11, HRAS, PIK3R1, TP63, and FAT1 have been identified in recurrent HPV positive OPSCC. Finding alternative therapeutic strategies for this young cohort may depend on upfront identification of HPV genotypes and mutations which are linked with worse outcomes, thus ensuring appropriate stratification of treatment regimens.

Prognostic Markers and Driver Genes and Options for Targeted Therapy in Human-Papillomavirus-Positive Tonsillar and Base-of-Tongue Squamous Cell Carcinoma

The incidence of Human-papillomavirus-positive (HPV⁺) tonsillar and base-of-tongue squamous cell carcinoma (TSCC and BOTSCC, respectively) is increasing epidemically, but they have better prognosis than equivalent HPV-negative (HPV⁻) cancers, with roughly 80% vs. 50% 3-year disease-free survival, respectively. The majority of HPV⁺ TSCC and BOTSCC patients therefore most likely do not require the intensified chemoradiotherapy given today to head and neck cancer patients and would with de-escalated therapy avoid several severe side effects. Moreover, for those with poor prognosis, survival has not improved, so better-tailored alternatives are urgently needed. In line with refined personalized medicine, recent studies have focused on identifying predictive markers and driver cancer genes useful for better stratifying patient treatment as well as for targeted therapy. This review presents some of these endeavors and briefly describes some recent experimental progress and some clinical trials with targeted therapy.

Molecular Signatures of Gynecological Cancers: Clinicians Perspective

Large-scale molecular profiling and DNA sequencing has revolutionized cancer research. Precision medicine is a rapidly developing area in cancer care but it is not uniformly applied across different tumor types. Biomarker-based therapy is associated with improved outcomes, both in terms of progression-free survival and overall survival. Comprehensive genomic profiling (CGP) uses next-generation sequencing to analyze the complete coding sequence of hundreds of genes from a small amount of tissue. Genes included in these assays are those associated with cancer development or have diagnostic, prognostic, familial, or therapeutic implications Genomic profiling is emerging as a clinically viable tool to personalize patient's treatment. This article discusses how the insights gained through CGP can impact treatment plan in common gynecological cancers.

Human papilloma virus (HPV) integration signature in Cervical Cancer: identification of MACROD2 gene as HPV hot spot integration site

Background

Cervical cancer (CC) remains a leading cause of gynaecological cancer-related mortality with infection by human papilloma virus (HPV) being the most important risk factor. We analysed the association between different viral integration signatures, clinical parameters and outcome in pre-treated CCs.

Methods

Different integration signatures were identified using HPV double capture followed by next-generation sequencing (NGS) in 272 CC patients from the BioRAIDs study [NCT02428842]. Correlations between HPV integration signatures and clinical, biological and molecular features were assessed.

Results

Episomal HPV was much less frequent in CC as compared to anal carcinoma (p < 0.0001). We identified >300 different HPV-chromosomal junctions (inter- or intra-genic). The most frequent integration site in CC was in MACROD2 gene followed by MIPOL1/TTC6 and TP63. HPV integration signatures were not associated with histological subtype, FIGO staging, treatment or PFS. HPVs were more frequently episomal in PIK3CA mutated tumours (p = 0.023). Viral integration type was dependent on HPV genotype (p < 0.0001); HPV18 and HPV45 being always integrated. High HPV copy number was associated with longer PFS (p = 0.011).

Conclusions

This is to our knowledge the first study assessing the prognostic value of HPV integration in a prospectively annotated CC cohort, which detects a hotspot of HPV integration at MACROD2; involved in impaired PARP1 activity and chromosome instability.

DeepHPV: a deep learning model to predict human papillomavirus integration sites

Human papillomavirus (HPV) integrating into human genome is the main cause of cervical carcinogenesis. HPV integration selection preference shows strong dependence on local genomic environment. Due to this theory, it is possible to predict HPV integration sites. However, a published bioinformatic tool is not available to date. Thus, we developed an attention-based deep learning model DeepHPV to predict HPV integration sites by learning environment features automatically. In total, 3608 known HPV integration sites were applied to train the model, and 584 reviewed HPV integration sites were used as the testing dataset. DeepHPV showed an area under the receiver-operating characteristic (AUROC) of 0.6336 and an area under the precision recall (AUPR) of 0.5670. Adding RepeatMasker and TCGA Pan Cancer peaks improved the model performance to 0.8464 and 0.8501 in AUROC and 0.7985 and 0.8106 in AUPR, respectively. Next, we tested these trained models on independent database VISDB and found the model adding TCGA Pan Cancer performed better (AUROC: 0.7175, AUPR: 0.6284) than the model adding RepeatMasker peaks (AUROC: 0.6102, AUPR: 0.5577). Moreover, we introduced attention mechanism in DeepHPV and enriched the transcription factor binding sites including BHLHA15, CHR, COUP-TFII, DMRTA2, E2A, HIC1, INR, NPAS, Nr5a2, RARa, SCL, Snail1, Sox10, Sox3, Sox4, Sox6, STAT6, Tbet, Tbx5, TEAD, Tgif2, ZNF189, ZNF416 near attention intensive sites. Together, DeepHPV is a robust and explainable deep learning model, providing new insights into HPV integration preference and mechanism. Availability: DeepHPV is available as an open-source software and can be downloaded from https://github.com/JiuxingLiang/DeepHPV.git, Contact: huzheng1998@163.com, liangjiuxing@m.scnu.edu.cn, lizheyzy@163.com.

Oncogenic driver mutations predict outcome in a cohort of head and neck squamous cell carcinoma (HNSCC) patients within a clinical trial

234 diagnostic formalin-fixed paraffin-embedded (FFPE) blocks from homogeneously treated patients with locally advanced head and neck squamous cell carcinoma (HNSCC) within a multicentre phase III clinical trial were characterised. The mutational spectrum was examined by next generation sequencing in the 26 most frequent oncogenic drivers in cancer and correlated with treatment response and survival. Human papillomavirus (HPV) status was measured by p16INK4a immunohistochemistry in oropharyngeal tumours. Clinicopathological features and response to treatment were measured and compared with the sequencing results. The results indicated TP53 as the most mutated gene in locally advanced HNSCC. HPV-positive oropharyngeal tumours were less mutated than HPV-negative tumours in TP53 (p < 0.01). Mutational and HPV status influences patient survival, being mutated or HPV-negative tumours associated with poor overall survival (p < 0.05). No association was found between mutations and clinicopathological features. This study confirmed and expanded previously published genomic characterization data in HNSCC. Survival analysis showed that non-mutated HNSCC tumours associated with better prognosis and lack of mutations can be identified as an important biomarker in HNSCC. Frequent alterations in PI3K pathway in HPV-positive HNSCC could define a promising pathway for pharmacological intervention in this group of tumours.

lncRNA PART1 and MIR17HG as ΔNp63α direct targets regulate tumor progression of cervical squamous cell carcinoma

Cervical cancer (CC) remains one of the leading causes of mortality of female cancers worldwide, with more than 90% being cervical squamous cell carcinoma (CSCC). ΔNp63α is the predominant isoform expressed in cervical epithelial tissues and exerts its antitumor function in CSCC. In this study, we have identified 39 long noncoding RNAs as ΔNp63α targets in CSCC through RNA sequencing and chromatin immunoprecipitation sequencing, in which we further confirmed and focused on the two tumor‐related long noncoding RNAs, PART1 (lncPART1) and MIR17HG (lncMIR17HG). Experiments from stable overexpression/knockdown cell lines revealed that lncPART1 and lncMIR17HG regulated cell proliferation, migration, and invasion. In vivo experiments further showed that lncPART1 suppresses tumor growth in CSCC‐derived tumors. Examinations of clinical tissues indicated that the expression of lncPART1 was positively correlated with ΔNp63α expression, while lncMIR17HG was negatively correlated with ΔNp63α expression, suggesting that ΔNp63α plays a central role via regulating its direct targets in the progression of CSCC. These findings provide novel insights in targeted therapy of cervical cancers.

Accurate Detection of HPV Integration Sites in Cervical Cancer Samples Using the Nanopore MinION Sequencer Without Error Correction

During the carcinogenesis of cervical cancer, the DNA of human papillomavirus (HPV) is frequently integrated into the human genome, which might be a biomarker for the early diagnosis of cervical cancer. Although the detection sensitivity of virus infection status increased significantly through the Illumina sequencing platform, there were still disadvantages remain for further improvement, including the detection accuracy and the complex integrated genome structure identification, etc. Nanopore sequencing has been proven to be a fast yet accurate technique of detecting pathogens in clinical samples with significant longer sequencing length. However, the identification of virus integration sites, especially HPV integration sites was seldom carried out by using nanopore platform. In this study, we evaluated the feasibility of identifying HPV integration sites by nanopore sequencer. Specifically, we re-sequenced the integration sites of a previously published sample by both nanopore and Illumina sequencing. After analyzing the results, three points of conclusions were drawn: first, 13 out of 19 previously published integration sites were found from all three datasets (i.e., nanopore, Illumina, and the published data), indicating a high overlap rate and comparability among the three platforms; second, our pipeline of nanopore and Illumina data identified 66 unique integration sites compared with previous published paper with 13 of them being verified by Sanger sequencing, indicating the higher integration sites detection sensitivity of our results compared with published data; third, we established a pipeline which could be used in HPV integration site detection by nanopore sequencing data without doing error correction analysis. In summary, a new nanopore data analysis method was tested and proved to be reliable in integration sites detection compared with methods of existing Illumina data analysis pipeline with less sequencing data required. It provides a solid evidence and tool to support the potential application of nanopore in virus status identification.

Targeting DNA Double-Strand Break Repair Enhances Radiosensitivity of HPV-Positive and HPV-Negative Head and Neck Squamous Cell Carcinoma to Photons and Protons

The response of head and neck squamous cell carcinoma (HNSCC) to radiotherapy depends on human papillomavirus type 16 (HPV) status, and where improved outcome and survival is observed in HPV-positive disease. However, strategies to further radiosensitise the tumours, particularly relatively radioresistant HPV-negative HNSCC, are actively being sought. The impact of targeting the major protein kinases involved in the signaling of DNA double-strand break (DSB) repair, namely ataxia telangiectasia-mutated (ATM), ataxia telangiectasia and Rad3-related (ATR), and the catalytic subunit of DNA-dependent protein kinase (DNA-Pkcs), on the radiosensitisation of HNSCC cells was examined. The response to both conventional photon radiotherapy, but also proton beam therapy, was analysed by clonogenic assays and 3D spheroid growth. We observed that inhibition of ATM, ATR, and particularly DNA-Pkcs, caused a significant reduction in HNSCC cell survival post-irradiation with both photons and protons, with less of an impact on the most radiosensitive HPV-positive cell line. The inhibition of DNA-Pkcs and, to a lesser extent ATM, in combination with radiation was also more effective at inhibiting the growth of 3D spheroids derived from relatively radioresistant HPV-negative HNSCC. Similar effects of the inhibitors were observed comparing photon and proton irradiation, demonstrating the potential for targeting DSB repair as an effective combination treatment for HNSCC.

Management of HPV-Related Squamous Cell Carcinoma of the Head and Neck: Pitfalls and Caveat

Head and neck squamous cell carcinomas (HNSCCs) are a very heterogeneous group of malignancies arising from the upper aerodigestive tract. They show different clinical behaviors depending on their origin site and genetics. Several data support the existence of at least two genetically different types of HNSCC, one virus-related and the other alcohol and/or tobacco and oral trauma-related, which show both clinical and biological opposite features. In fact, human papillomavirus (HPV)-related HNSCCs, which are mainly located in the oropharynx, are characterized by better prognosis and response to therapies when compared to HPV-negative HNSCCs. Interestingly, virus-related HNSCC has shown a better response to conservative (nonsurgical) treatments and immunotherapy, opening questions about the possibility to perform a pretherapy assessment which could totally guide the treatment strategy. In this review, we summarize molecular differences and similarities between HPV-positive and HPV-negative HNSCC, highlighting their impact on clinical behavior and on therapeutic strategies.

Chromatin structure regulates cancer-specific alternative splicing events in primary HPV-related oropharyngeal squamous cell carcinoma

Human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV+ OPSCC) represents a unique disease entity within head and neck cancer with rising incidence. Previous work has shown that alternative splicing events (ASEs) are prevalent in HPV+ OPSCC, but further validation is needed to understand the regulation of this process and its role in these tumours. In this study, eleven ASEs (GIT2, CTNNB1, MKNK2, MRPL33, SIPA1L3, SNHG6, SYCP2, TPRG1, ZHX2, ZNF331, and ELOVL1) were selected for validation from 109 previously published candidate ASEs to elucidate the post-transcriptional mechanisms of oncogenesis in HPV+ disease. In vitro qRT-PCR confirmed differential expression of 9 of 11 ASE candidates, and in silico analysis within the TCGA cohort confirmed 8 of 11 candidates. Six ASEs (MRPL33, SIPA1L3, SNHG6, TPRG1, ZHX2, and ELOVL1) showed significant differential expression across both methods. Further evaluation of chromatin modification revealed that ASEs strongly correlated with cancer-specific distribution of acetylated lysine 27 of histone 3 (H3K27ac). Subsequent epigenetic treatment of HPV+ HNSCC cell lines (UM-SCC-047 and UPCI-SCC-090) with JQ1 not only induced downregulation of cancer-specific ASE isoforms, but also growth inhibition in both cell lines. The UPCI-SCC-090 cell line, with greater ASE expression, also showed more significant growth inhibition after JQ1 treatment. This study confirms several novel cancer-specific ASEs in HPV+OPSCC and provides evidence for the role of chromatin modifications in regulation of alternative splicing in HPV+OPSCC. This highlights the role of epigenetic changes in the oncogenesis of HPV+OPSCC, which represents a unique, unexplored target for therapeutics that can alter the global post-transcriptional landscape.

Human papillomavirus (HPV) 16 infection is not detected in rectal carcinoma

Introduction

Persistence of human papillomavirus (HPV) infections is associated with squamous cell carcinomas of different human anatomic sites. Several studies have suggested a potential role for HPV infection, particularly HPV16 genotype, in rectal cancer carcinogenesis.. The aim of this study was to assess the frequency of oncogenic HPV 16 viral DNA sequences in rectal carcinomas cases retrieved from the pathology archive of Braga Hospital, North Portuga.

Methods

TaqMan-based type-specific real-time PCR for HPV 16 was performed using primers and probe targeting HPV16 E7 region.

Results

Most of the rectal cancer patients (88.5%, n = 206 patients), were symptomatic at diagnosis. The majority of the lesions (55.3%, n = 129) presented malignancies of polypoid/vegetant phenotype. 26.8% (n = 63) had synchronic metastasis at diagnosis. 26.2% (n = 61) patients had clinical indication for neoadjuvant therapy. Most patients with rectal cancer were stage IV (19.7% patients), followed by stage IIA (19.3%) and stage I (18.5%). All cases of the present series tested negative for HPV16.

Conclusion

The total of negative tests for HPV 16 infection is a robust argument to support the assumption that HPV 16 infection, despite of previous evidences, is not involved in rectal cancer carcinogenesis and progression.

Possible Biomarkers and Therapeutic Targets for the Management of Cervical Cancer

Introduction:

Cervical cancer is the most prevalent cancer in the world due to unusual extension of cervical cell. Cervical cancer occurs due to exposure of HPV (Human papillomavirus). According to WHO, it is the 4th most ordinary cancer in women. In 2018, approx 6.6% of population was affected around the world and 570,000 new cases were reported. In low and middle-income countries, 90% of cervical cancer deaths occur.

Methods:

Despite various factors that cause cervical cancer are included exposure to HPV, dysregulation of CASPASE enzyme, elevated expression of IAPs (Inhibitor apoptotic protein), E6 and E7 gene of HPV, inhibition of p53, BAK, p16 upregulation, CDK-inactivation causing cervical cancer, role of VEGF, role of estrogen and its receptor in cervical cancer.

Results:

Cervical cancer can be screened by Pep test. There are various therapies that can be used to treat cervical cancer. As these therapies have various side effects, so the world is moving to herbal formulations to treat cervical cancer.

Conclusion:

In this study, we will discuss cervical cancer, its cause, symptoms, pathophysiology and treatments. Early screening and detection can help in reducing the overall burden of cervical cancer in the near future.

A global epidemic increase of an HPV-induced tonsil and tongue base cancer - potential benefit from a pan-gender use of HPV vaccine

In 2007, human papillomavirus (HPV) type 16 was finally recognized as a risk factor, besides smoking and alcohol, for oropharyngeal squamous cell carcinoma (OPSCC), including tonsillar squamous cell carcinoma (TSCC), by the International Agency for Research against Cancer. Just before, in 2006, the Food and Drug Administration had approved Gardasil, the first vaccine against HPV16, 18, 6 and 11, for preventive vaccination women against cervical cancer. Concurrently, some Western countries, where smoking was decreasing, disclosed an epidemic increase in the incidence of OPSCC, especially of TSCC and base of tongue cancer (BOTSCC), together accounting for 80-90% of all OPSCCs, and mainly affecting men. The epidemic was later revealed to be due to a rise in HPV-positive cases, and scientists in the field suggested HPV vaccination also of boys. Globally, there are roughly 96 000 incident OPSCC cases/year of which 20-24% are caused by HPV, thereby accounting for around 22 000 OPSCC cases annually. Of these cases, 80-90% are due to HPV16 infection and would be prevented with the presently registered HPV vaccines. In Western countries, such as Sweden (with almost 400 TSCC and BOTSCC cases per year) and the United States, HPV prevalence in OPSCC is higher and around 70%. HPV vaccination of girls has been initiated in many countries, and the vaccines have been efficient and their side effects limited. HPV vaccination of boys has, however, been the exception, but should definitely not be delayed any further. It would benefit both girls and boys directly, and result in better and more robust herd immunity. Today, we have the possibility to eliminate several high-risk HPV types in the younger generations and avoid more than 600 000 cancer cases annually worldwide, and this possibility should be embraced by offering global pan-gender HPV vaccination.

Whole-exome and RNA sequencing reveal novel insights into the pathogenesis of HPV associated cervical cancer

BACKGROUND

Worldwide, cervical cancer is the fouth leading cause of deaths in gynecological oncology. Although the causes of cervical cancer have been extensively investigated, understanding of its exact pathogenesis remains incomplete.

OBJECTIVE

This study aimed to identify alterations of genome and transcriptome of HPV associated cervical cancer pathogenesis using multi-omics approaches.

METHODS

Cervical cancer and matched adjacent non-tumor specimens of one HPV16+ and two HPV- patients were sampled for whole-exome sequencing (WES) and RNA sequencing to characterize DNA mutations and gene expression profiles. WES and Affymetrix SNP 6.0 arrays data were analyzed from 6 HPV- and 93 HPV16+ cervical cancer patients in the cancer genome atlas (TCGA) database, as an independent validation group.

RESULTS

WES identified 64 somatic mutation genes in tumors of 3 patients. HPV16+ tumor got fewer somatic mutated genes than HPV- tumors, which was validated by TCGA results. In this study, somatic mutated profile, CNV and gene expression heat map presented that HPV16+ tumors was distinct with HPV- tumors. The most significant altered pathways and GO terms were both related with cell cycle. Integrated analysis of multi-omics showed positive correlation between gene expression level and copy numbers.

CONCLUSIONS

The results of this study provided novel insights into the pathogenesis of HPV associated cervical cancer.

Integrative analysis of DNA methylation and gene expression identified cervical cancer-specific diagnostic biomarkers

Cervical cancer is the leading cause of death among women with cancer worldwide. Here, we performed an integrative analysis of Illumina HumanMethylation450K and RNA-seq data from TCGA to identify cervical cancer-specific DNA methylation markers. We first identified differentially methylated and expressed genes and examined the correlation between DNA methylation and gene expression. The DNA methylation profiles of 12 types of cancers, including cervical cancer, were used to generate a candidate set, and machine-learning techniques were adopted to define the final cervical cancer-specific markers in the candidate set. Then, we assessed the protein levels of marker genes by immunohistochemistry by using tissue arrays containing 93 human cervical squamous cell carcinoma samples and cancer-adjacent normal tissues. Promoter methylation was negatively correlated with the local regulation of gene expression. In the distant regulation of gene expression, the methylation of hypermethylated genes was more likely to be negatively correlated with gene expression, while the methylation of hypomethylated genes was more likely to be positively correlated with gene expression. Moreover, we identified four cervical cancer-specific methylation markers, cg07211381 (RAB3C), cg12205729 (GABRA2), cg20708961 (ZNF257), and cg26490054 (SLC5A8), with 96.2% sensitivity and 95.2% specificity by using the tenfold cross-validation of TCGA data. The four markers could distinguish tumors from normal tissues with a 94.2, 100, 100, and 100% AUC in four independent validation sets from the GEO database. Overall, our study demonstrates the potential use of methylation markers in cervical cancer diagnosis and may boost the development of new epigenetic therapies.

Immune Signature-Based Subtypes of Cervical Squamous Cell Carcinoma Tightly Associated with Human Papillomavirus Type 16 Expression, Molecular Features, and Clinical Outcome

Substantial heterogeneity exists within cervical cancer that is generally infected by human papillomavirus (HPV). However, the most common histological subtype of cervical cancer, cervical squamous cell carcinoma (CSCC), is poorly characterized regarding the association between its heterogeneity and HPV oncoprotein expression. We filtered out 138 CSCC samples with infection of HPV16 only as the first step; then we compressed HPV16 E6/E7 expression as HPV_pca and correlated HPV_pca with the immunological profiling of CSCC based on supervised clustering to discover subtypes and to characterize the differences between subgroups in terms of the HPV_pca level, pathway activity, epigenetic dysregulation, somatic mutation frequencies, and likelihood of responding to chemo/immunotherapies. Supervised clustering of immune signatures revealed two HPV16 subtypes (namely, HPV16-IMM and HPV16-KRT) that correlated with HPV_pca and clinical outcomes. HPV16-KRT is characterized by elevated expression of genes in keratinization, biological oxidation, and Wnt signaling, whereas HPV16-IMM has a strong immune response and mesenchymal features. HPV16-IMM exhibited much more epigenetic silencing and significant mutation at FBXW7, while MUC4 and PIK3CA were mutated frequently for HPV16-KRT. We also imputed that HPV16-IMM is much more sensitive to chemo/immunotherapy than is HPV16-KRT. Our characterization tightly links the expression of HPV16 E6/E7 with biological and clinical outcomes of CSCC, providing valuable molecular-level information that points to decoding heterogeneity. Together, these results shed light on stratifications of CSCC infected by HPV16 and shall help to guide personalized management and treatment of patients.

High olive oil diets enhance cervical tumour growth in mice: transcriptome analysis for potential candidate genes and pathways

BACKGROUND

Numerous epidemiologic studies have found a close association between obesity and cancer. Dietary fat is a fundamental contributor to obesity and is a risk factor for cancer. Thus far, the impact of dietary olive oil on cancer development remains inconclusive, and little is known about its underlying mechanisms.

METHODS

Nude mouse xenograft models were used to examine the effects of high olive oil diet feeding on cervical cancer (CC) development and progression. Cell proliferation, migration and invasion were observed by the methods of EdU incorporation, Wound healing and Transwell assay, separately. RNA-sequencing technology and comprehensive bioinformatics analyses were used to elucidate the molecular processes regulated by dietary fat. Differentially expressed genes (DEGs) were identified and were functionally analyzed by Gene Ontology (GO), Kyoto Enrichment of Genes and Genomes (KEGG). Then, protein-protein interaction (PPI) network and sub-PPI network analyses were conducted using the STRING database and Cytoscape software.

RESULTS

A high olive oil diet aggravated tumourigenesis in an experimental xenograft model of CC. Oleic acid, the main ingredient of olive oil, promoted cell growth and migration in vitro. Transcriptome sequencing analysis of xenograft tumour tissues was then performed to elucidate the regulation of molecular events regulated by dietary fat. Dietary olive oil induced 648 DEGs, comprising 155 up-regulated DEGs and 493 down-regulated DEGs. GO and pathway enrichment analysis revealed that some of the DEGs including EGR1 and FOXN2 were involved in the transcription regulation and others, including TGFB2 and COL4A3 in cell proliferation. The 15 most strongly associated DEGs were selected from the PPI network and hub genes including JUN, TIMP3, OAS1, OASL and EGR1 were confirmed by real-time quantitative PCR analysis.

CONCLUSIONS

Our study suggests that a high olive oil diet aggravates CC progression in vivo and in vitro. We provide clues to build a potential link between dietary fat and cancerogenesis and identify areas requiring further investigation.