Frequent genomic structural alterations at HPV insertion sites in cervical carcinoma

Navigating the landscape of HPV-associated cancers: From epidemiology to prevention

Human Papillomavirus (HPV) is a widespread infection associated with various cancers, including cervical, oropharyngeal, anal, and genital cancers. This infection contributes to 5 % of global cancer cases annually, affecting approximately 625,600 women and 69,400 men. Cervical cancer remains the most prevalent HPV-linked cancer among females, with the highest incidence seen in low and middle-income countries (LMICs). While most HPV infections are transient, factors such as HPV variants, age, gender, and socioeconomic status influence transmission risks. HPV is categorized into high-risk (HR-HPV) and low-risk types, with strains like HPV 16 and 18 displaying distinct demographic patterns. The intricate pathogenesis of HPV involves genetic and epigenetic interactions, with HPV oncogenes (E6 and E7) and integration into host DNA playing a pivotal role in driving malignancies. Early diagnostics, utilizing HPV DNA testing with surrogate markers such as p16, and advanced molecular techniques like PCR, liquid biopsy, and NGS, significantly impact the management of HPV-induced cancers. Effectively managing HPV-related cancers demands a multidisciplinary approach, including immunotherapy, integrating current therapies, ongoing trials, and evolving treatments. Prevention via HPV vaccination and the inclusion of cervical cancer screening in national immunization programs by conventional Pap smear examination and HPV DNA testing remains fundamental.Despite the preventability of HPV-related cancers, uncertainties persist in testing, vaccination, and treatment. This review article covers epidemiology, pathogenesis, diagnostics, management, prevention strategies, challenges, and future directions. Addressing issues like vaccine hesitancy, healthcare disparities, and advancing therapies requires collaboration among researchers, healthcare providers, policymakers, and the public. Advancements in understanding the disease's molecular basis and clinical progression are crucial for early detection, proper management, and improved outcomes.

The prevalence of HPV in Chongqing, China from 2017 to 2022: a retrospective cohort study

BACKGROUND AND AIMS

Human papillomavirus (HPV) infection is widespread. The regional differences in infection rates and dominant HPV types are a significant aspect of the virus's epidemiological characteristics. In China, these differences are further influenced by the country's diverse ethnic composition, as well as varying economic and cultural conditions across regions. Additionally, factors such as population migration contribute to the exacerbation of these regional disparities in HPV infection rates.

METHODS

We conducted an analysis of the current prevalence of HPV among 108,863 patients in Chongqing from 2017 to 2022, along with an investigation into the interaction between HPV and vaginal microbiota. Additionally, we explored the relationship between HPV clearance time and disease outcomes.

RESULTS

Our findings revealed an overall HPV positive rate of 30.14% (32,813/108,863) in Chongqing, China. The top five HPV positive rates were as follows: HPV52 (7.78%), HPV16 (5.85%), HPV58 (4.41%), HPV53 (3.17%), and HPV39 (2.47%). The age-related trend of HPV infection exhibited a bimodal curve. Importantly, the positive rate of high-risk HPV (BV: 24.74%, AV: 24.38%, VVC: 18.19%, TV: 22.81%) was significantly higher than those in the control group (all P<0.05). Follow-up for 4 years, we observed an overall clearance rate of HPV at 92.20%, with a median clearance time of 13 (1-61) months.

CONCLUSION

These findings hold significant implications for the Chinese government in evaluating the prevalence of HPV in the Chongqing region and formulating effective HPV prevention strategies. Moreover, these findings could contribute to the enhancement of public health initiatives related to HPV in the region.

Identification of novel genomic hotspots and tumor-relevant genes via comprehensive analysis of HPV integration in Chinese patients of cervical cancer

Cervical cancer accounts for 10-15% of cancer-related mortality among women globally. Infection with high-risk human papillomavirus (HPV) types constitutes a significant etiological factor in the development of cervical carcinoma. The integration of HPV DNA into the host genome is considered a pivotal event in cervical carcinogenesis. Nevertheless, the precise mechanisms underlying HPV integration and its role in promoting cancer progression remain inadequately understood. Therefore, this study aims to identify potential common denominators at HPV DNA integration sites and to analyze the adjacent cellular sequences. We conducted whole-genome sequencing on 13 primary cervical cancer samples, employing the chromosomal coordinates of 537 breakpoints to assess the statistical overrepresentation of integration sites in relation to various chromatin features. Our analysis, which encompassed all chromosomes, identified several integration hotspots within the human genome, notably at 14q32.2, 10p15, and 2q37. Additionally, our findings indicated a preferential integration of HPV DNA into intragenic and gene-dense regions of human chromosomes. A substantial number of host cellular genes impacted by the integration sites were associated with cancer, including IKZF2, IL26, AHRR, and PDCD6. Furthermore, the cellular genes targeted by integration were enriched in tumor-related terms and pathways, as demonstrated by gene ontology and KEGG analysis. In conclusion, these findings enhance our understanding of HPV integration sites and provide deeper insights into the molecular mechanisms underlying the pathogenesis of cervical carcinoma.

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.

Mapping the landscape of HPV integration and characterising virus and host genome interactions in HPV-positive oropharyngeal squamous cell carcinoma

BACKGROUND

Human papillomavirus (HPV) integration into the host genome is an important factor in HPV(+)OPSCC carcinogenesis, in conjunction with HPV oncoproteins E6/E7. However, a well-studied investigation about virus-host interaction still needs to be completed. Our objective is to characterise HPV integration to investigate potential mechanisms of tumourigenesis independent of E6/E7 oncoproteins.

MATERIALS AND METHODS

High-throughput viral integration detection was performed on 109 HPV(+)OPSCC tumours with relevant clinicopathological information. Of these tumours, 38 tumours underwent targeted gene sequencing, 29 underwent whole exome sequencing and 26 underwent RNA sequencing.

RESULTS

HPV integration was detected in 94% of tumours (with a mean integration count of 337). Tumours occurring at the tonsil/oropharyngeal wall that exhibit higher PD-L1 expression demonstrated increased integration sites (p = .024). HPV exhibited a propensity for integration at genomic sites located within specific fragile sites (FRA19A) or genes associated with functional roles such as cell proliferation and differentiation (PTEN, AR), immune evasion (CD274) and glycoprotein biosynthesis process (FUT8). The viral oncogenes E2, E4, E6 and E7 tended to remain intact. HPV fragments displayed enrichment within host copy number variation (CNV) regions. However, insertions into genes related to altered homologous recombination repair were infrequent. Genes with integration had distinct expression levels. Fifty-nine genes whose expression level was affected by viral integration were identified, for example, EPHB1, which was reported to be involved in cellular protein metabolic process.

CONCLUSIONS

HPV can promote oncogenesis through recurrent integration into functional host genome regions, leading to subsequent genomic aberrations and gene expression disruption. This study characterises viral integrations and virus-host interactions, enhancing our understanding of HPV-related carcinogenesis mechanisms.

Human papilloma virus vaccines: A comprehensive narrative review

Cervical cancer is one of the most common cancers in women aged 15-44 years in the world, with more than three-quarters of cases diagnosed at a locally advanced clinical stage with minor prospects of survival. Although only a small percentage of women with Human Papilloma Virus (HPV) develop cervical cancer and most of the HPV infections are cleared subsequently at primary stage itself, but seroconversion not always guarantees that the individual is immune to HPV. The advent of the cervical carcinoma vaccine has raised the expectations that eradication of cervical carcinoma might be possible in the near future as it exhibited remarkably high efficacy against the vaccine-specific types in naive women with no serious vaccine-related adverse events. Few prophylactic HPV vaccines are currently licensed in over 100 countries. It has also been suggested that vaccinating both men and women is more beneficial than vaccinating only females. Vaccination is a cost-effective strategy to reduce the incidence of cervical cancer and mortality compared to no vaccination based on the cost of cancer treatment. Well-coordinated vaccination strategy with focus on adolescent girls and if possible, boys can lead to dramatic impact on disease reduction around the world.

Extrachromosomal Amplification of Human Papillomavirus Episomes Is a Mechanism of Cervical Carcinogenesis

SIGNIFICANCE

Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism for oncogenesis without integration.

The 3D genome and its impacts on human health and disease

Eukaryotic genomes are highly compacted in the cell nucleus. Two loci separated by a long linear distance can be brought into proximity in space through DNA-binding proteins and RNAs, which contributes profoundly to the regulation of gene expression. Recent technology advances have enabled the development and application of the chromosome conformation capture (3C) technique and a host of 3C-based methods that enable genome-scale investigations into changes in chromatin high-order structures during diverse physiological processes and diseases. In this review, we introduce 3C-based technologies and discuss how they can be utilized to glean insights into the impacts of three-dimensional (3D) genome organization in normal physiological and disease processes.

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.

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.

Genome-wide profiling of retroviral DNA integration and its effect on clinical pre-infusion CAR T-cell products

Background

Clinical CAR T-cell therapy using integrating vector systems represents a promising approach for the treatment of hematological malignancies. Lentiviral and γ-retroviral vectors are the most commonly used vectors in the manufacturing process. However, the integration pattern of these viral vectors and subsequent effect on CAR T-cell products is still unclear.

Methods

We used a modified viral integration sites analysis (VISA) pipeline to evaluate viral integration events around the whole genome in pre-infusion CAR T-cell products. We compared the differences of integration pattern between lentiviral and γ-retroviral products. We also explored whether the integration sites correlated with clinical outcomes.

Results

We found that γ-retroviral vectors were more likely to insert than lentiviral vectors into promoter, untranslated, and exon regions, while lentiviral vector integration sites were more likely to occur in intron and intergenic regions. Some integration events affected gene expression at the transcriptional and post-transcriptional level. Moreover, γ-retroviral vectors showed a stronger impact on the host transcriptome. Analysis of individuals with different clinical outcomes revealed genes with differential enrichment of integration events. These genes may affect biological functions by interrupting amino acid sequences and generating abnormal proteins, instead of by affecting mRNA expression. These results suggest that vector integration is associated with CAR T-cell efficacy and clinical responses.

Conclusion

We found differences in integration patterns, insertion hotspots and effects on gene expression vary between lentiviral and γ-retroviral vectors used in CAR T-cell products and established a foundation upon which we can conduct further analyses.

The Drivers, Mechanisms, and Consequences of Genome Instability in HPV-Driven Cancers

Human papillomavirus (HPV) is the causative driver of cervical cancer and a contributing risk factor of head and neck cancer and several anogenital cancers. HPV's ability to induce genome instability contributes to its oncogenicity. HPV genes can induce genome instability in several ways, including modulating the cell cycle to favour proliferation, interacting with DNA damage repair pathways to bring high-fidelity repair pathways to viral episomes and away from the host genome, inducing DNA-damaging oxidative stress, and altering the length of telomeres. In addition, the presence of a chronic viral infection can lead to immune responses that also cause genome instability of the infected tissue. The HPV genome can become integrated into the host genome during HPV-induced tumorigenesis. Viral integration requires double-stranded breaks on the DNA; therefore, regions around the integration event are prone to structural alterations and themselves are targets of genome instability. In this review, we present the mechanisms by which HPV-dependent and -independent genome instability is initiated and maintained in HPV-driven cancers, both across the genome and at regions of HPV integration.

Human Papillomavirus Integration Strictly Correlates with Global Genome Instability in Head and Neck Cancer

Human papillomavirus (HPV)-positive head and neck cancers, predominantly oropharyngeal squamous cell carcinoma (OPSCC), exhibit epidemiologic, clinical, and molecular characteristics distinct from those OPSCCs lacking HPV. We applied a combination of whole-genome sequencing and optical genome mapping to interrogate the genome structure of HPV-positive OPSCCs. We found that the virus had integrated in the host genome in two thirds of the tumors examined but resided solely extrachromosomally in the other third. Integration of the virus occurred at essentially random sites within the genome. Focal amplification of the virus and the genomic sequences surrounding it often occurred subsequent to integration, with the number of tandem repeats in the chromosome accounting for the increased copy number of the genome sequences flanking the site of integration. In all cases, viral integration correlated with pervasive genome-wide somatic alterations at sites distinct from that of viral integration and comprised multiple insertions, deletions, translocations, inversions, and point mutations. Few or no somatic mutations were present in tumors with only episomal HPV. Our data could be interpreted by positing that episomal HPV is captured in the host genome following an episode of global genome instability during tumor development. Viral integration correlated with higher grade tumors, which may be explained by the associated extensive mutation of the genome and suggests that HPV integration status may inform prognosis.

IMPLICATIONS

Our results indicate that HPV integration in head and neck cancer correlates with extensive pangenomic structural variation, which may have prognostic implications.

Strategic Significance of Low Viral Load of Human Papillomavirus in Uterine Cervical Cytology Specimens

Infection with high-risk (HR) Human Papillomavirus (HPV) is associated with the development of precancerous lesions or invasive carcinoma of the uterine cervix. Thus, the high viral load (VL) of HR-HPV DNA currently serves as a representative quantitative marker for cervical cancer. However, the clinical significance of low HPV DNA VL remains undetermined. This study aimed to evaluate the clinical association between the low HPV DNA VL and cytology/histologic diagnosis of cervical samples. We searched the electronic medical databases for the resultant analyses of HPV genotyping among patients who underwent treatment for any cervical lesion or who had undergone gynecological examinations with any positive HPV results according to the national cancer screening service between 2015 and 2016. HPV testing with genotyping and semi-quantitative VL measurement was conducted using an Anyplex^TM II H28 Detection assay (H28 assay, Seegene, Seoul, Republic of Korea). The H28 assay is a multiplex semi-quantitative real-time PCR test using the tagging of oligonucleotide cleavage and extension (TOCE) technology. The VL was semi-quantified as high (3+; positive signal before 31 PCR cycles), intermediate (2+; positive between 31 and 39 PCR cycles), or low (1+; positive after 40 PCR cycles). Out of 5940 HPV VL analyses, 356 assays (5.99%) were reported as low VL (1+) of HPV DNA. Matched cytology diagnoses were mostly negative findings (n = 347, 97.5%), except for seven cases of atypical squamous cells of undetermined significance (1.9%) and two cases of atypical glandular cells (0.6%). During the follow-up periods, abnormal cytologic diagnoses were identified, including one case of high-grade squamous intraepithelial lesion (HSIL) and two low-grade squamous intraepithelial lesions (LSILs). The matched, confirmative histologic diagnosis of HSIL cytology was compatible with chronic inflammation, wherein the two LSILs had regular check-ups. None revealed clinically concerned outcomes associated with HPV-related squamous lesions. The cytology was most likely negative for malignancy when the VL of HPV DNA was low (1+). Additional strategic monitoring and management may thus be unnecessary.

Risk Factors Associated With Human Papillomavirus Infection, Cervical Cancer, and Precancerous Lesions in Large-Scale Population Screening

Cervical cancer is the most common gynecological malignancy and screening for risk factors with early detection has been shown to reduce the mortality. In this study, we aimed to analyze the characteristics and risk factors of human papillomavirus (HPV) infection and precancerous lesions in women and provide clinical evidence for developing strategies to prevent cervical precancerous lesions and cancer in women. Furthermore, we evaluated the influencing factors for high-risk HPV infection. From April 2018 to December 2021, 10,628 women were recruited for cervical cancer screening at Liaoning Cancer Hospital, Shenyang Sujiatun District Women’s and Infants Hospital, Benxi Manchu Autonomous County People’s Hospital, and Shandong Affiliated Hospital of Qingdao University. The study participants were tested to determine if they were HPV-positive (HPV +) or underwent thinprep cytology test (TCT) for atypical squamous cells of undetermined significance (ASCUS) and above. Furthermore, colposcopies and biopsies were performed for the histopathological examination. Finally, 9991 cases were included in the statistical analysis, and the factors influencing HPV infection and those related to cervical cancer and precancerous lesions were analyzed. HPV + infection, high-grade squamous intraepithelial lesion-positive (CINII +) in cervical high-grade intraepithelial neoplasia, and early cervical cancer diagnosis rates were 12.45, 1.09, and 95.41%, respectively. The potential risk factors for HPV were education ≤ high school [odds ratio (OR) = 1.279 (1.129–1.449), P < 0.001], age at initial sexual activity ≤ 19 years [OR = 1.517 (1.080–2.129), P = 0.016], sexual partners > 1 [OR = 1.310 (1.044–1.644), P = 0.020], ASCUS and above [OR = 11.891 (10.105–13.993), P < 0.001], non-condom contraception [OR = 1.255 (1.059–1.487), P = 0.009], and HSIL and above [OR = 1.541 (1.430–1.662), P < 0.001]. Compared with women aged 56–65 and 35–45 years [OR = 0.810 (0.690–0.950), P = 0.010] the HPV infection rate was significantly lower in those aged 46–55 years [OR = 0.79 (0.683–0.915), P = 0.002]. Furthermore, ≤ high school age [OR = 1.577 (1.042–2.387), P = 0.031], not breastfeeding [OR = 1.763 (1.109–2.804), P = 0.017], ASCUS and above [OR = 42.396 (28.042–64.098), P < 0.001] were potential risk factors for cervical cancer and precancerous lesions. In women with HPV infection, ≤ high school education level, initial sexual activity at ≤ 19 years of age, number of sexual partners > 1, ASCUS and above, non-condom contraception, HSIL and above were risk factors for HPV infection. Compared with women aged 56–65 years, those aged 35–45 and 46–55 years had significantly lower HPV infection rates, and high school age and below, non-breastfeeding, and ASCUS and above were all potential risk factors for cervical cancer and precancerous lesions.

HPV16 and HPV18 type-specific APOBEC3 and integration profiles in different diagnostic categories of cervical samples

Human papillomavirus (HPV) 16 and 18 are the most predominant types in cervical cancer. Only a small fraction of HPV infections progress to cancer, indicating that additional factors and genomic events contribute to the carcinogenesis, such as minor nucleotide variation caused by APOBEC3 and chromosomal integration. We analysed intra-host minor nucleotide variants (MNVs) and integration in HPV16 and HPV18 positive cervical samples with different morphology. Samples were sequenced using an HPV whole genome sequencing protocol TaME-seq. A total of 80 HPV16 and 51 HPV18 positive samples passed the sequencing depth criteria of 300× reads, showing the following distribution: non-progressive disease (HPV16 n = 21, HPV18 n = 12); cervical intraepithelial neoplasia (CIN) grade 2 (HPV16 n = 27, HPV18 n = 9); CIN3/adenocarcinoma in situ (AIS) (HPV16 n = 27, HPV18 n = 30); cervical cancer (HPV16 n = 5). Similar numbers of MNVs in HPV16 and HPV18 samples were observed for most viral genes, with the exception of HPV18 E4 with higher numbers across clinical categories. APOBEC3 signatures were observed in HPV16 lesions, while similar mutation patterns were not detected for HPV18. The proportion of samples with integration was 13% for HPV16 and 59% for HPV18 positive samples, with a noticeable portion located within or close to cancer-related genes.

Recurrent integration of human papillomavirus genomes at transcriptional regulatory hubs

Oncogenic human papillomavirus (HPV) genomes are often integrated into host chromosomes in HPV-associated cancers. HPV genomes are integrated either as a single copy or as tandem repeats of viral DNA interspersed with, or without, host DNA. Integration occurs frequently in common fragile sites susceptible to tandem repeat formation and the flanking or interspersed host DNA often contains transcriptional enhancer elements. When co-amplified with the viral genome, these enhancers can form super-enhancer-like elements that drive high viral oncogene expression. Here we compiled highly curated datasets of HPV integration sites in cervical (CESC) and head and neck squamous cell carcinoma (HNSCC) cancers, and assessed the number of breakpoints, viral transcriptional activity, and host genome copy number at each insertion site. Tumors frequently contained multiple distinct HPV integration sites but often only one “driver” site that expressed viral RNA. As common fragile sites and active enhancer elements are cell-type-specific, we mapped these regions in cervical cell lines using FANCD2 and Brd4/H3K27ac ChIP-seq, respectively. Large enhancer clusters, or super-enhancers, were also defined using the Brd4/H3K27ac ChIP-seq dataset. HPV integration breakpoints were enriched at both FANCD2-associated fragile sites and enhancer-rich regions, and frequently showed adjacent focal DNA amplification in CESC samples. We identified recurrent integration “hotspots” that were enriched for super-enhancers, some of which function as regulatory hubs for cell-identity genes. We propose that during persistent infection, extrachromosomal HPV minichromosomes associate with these transcriptional epicenters and accidental integration could promote viral oncogene expression and carcinogenesis.

Genomic alterations caused by HPV integration in a cohort of Chinese endocervical adenocarcinomas

The association between human papillomavirus (HPV) integration and relevant genomic changes in uterine cervical adenocarcinoma is poorly understood. This study is to depict the genomic mutational landscape in a cohort of 20 patients. HPV+ and HPV− groups were defined as patients with and without HPV integration in the host genome. The genetic changes between these two groups were described and compared by whole-genome sequencing (WGS) and whole-exome sequencing (WES). WGS identified 2916 copy number variations and 743 structural variations. WES identified 6113 somatic mutations, with a mutational burden of 2.4 mutations/Mb. Six genes were predicted as driver genes: PIK3CA, KRAS, TRAPPC12, NDN, GOLGA6L4 and BAIAP3. PIK3CA, NDN, GOLGA6L4, and BAIAP3 were recognized as significantly mutated genes (SMGs). HPV was detected in 95% (19/20) of patients with cervical adenocarcinoma, 7 of whom (36.8%) had HPV integration (HPV+ group). In total, 1036 genes with somatic mutations were confirmed in the HPV+ group, while 289 genes with somatic mutations were confirmed in the group without HPV integration (HPV− group); only 2.1% were shared between the two groups. In the HPV+ group, GOLGA6L4 and BAIAP3 were confirmed as SMGs, while PIK3CA, NDN, KRAS, FUT1, and GOLGA6L64 were identified in the HPV− group. ZDHHC3, PKD1P1, and TGIF2 showed copy number amplifications after HPV integration. In addition, the HPV+ group had significantly more neoantigens. HPV integration rather than HPV infection results in different genomic changes in cervical adenocarcinoma.

Pathology of HPV-Associated Head and Neck Carcinomas: Recent Data and Perspectives for the Development of Specific Tumor Markers

A significant subset of carcinomas developed in the head and neck (H&NCs) are associated with specific human papillomaviruses (HPV) genotypes. In particular, 40–60% of oropharyngeal carcinoma cases are linked to HPV. Epidemiological studies have demonstrated that HPV oral infections are predominantly sexually transmitted and are more frequent among men (10–18%) than women (3.6–8.8%). Although there is a large diversity of HPV genotypes associated with H&NCs, HPV16 lineage represents 83% of the reported cases. The prognostic value of HPV as a biological parameter is well recognized. However, the use of HPV DNA as a diagnostic and/or predictive marker is not fully developed. Recent data reporting the physical state of the HPV genome in tumors have shown that HPV DNA integration into the tumor cell genome could lead to the alteration of cellular genes implicated in oncogenesis. Most importantly, HPV DNA corresponds to a tumor marker that can be detected in the blood of patients. Profile of the HPV DNA molecular patterns in tumor cells using New Genome Sequencing-based technologies, allows the identification of highly specific tumor markers valuable for the development of innovative diagnostic and therapeutic approaches. This review will summarize recent epidemiological data concerning HPV-associated H&NCs, the genomic characterization of these tumors, including the presence of HPV DNA in tumor cells, and will propose perspectives for developing improved care of patients with HPV-associated H&NCs, based on the use of viral sequences as personalized tumor markers and, over the longer term, as a therapeutic target.

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.

Identification of Specific Tumor Markers in Vulvar Carcinoma Through Extensive Human Papillomavirus DNA Characterization Using Next Generation Sequencing Method

OBJECTIVES

A subset of vulvar carcinomas (VC) are associated with human papillomavirus (HPV) DNA. This trait can be used to identify tumor markers for patient's follow-up. A large diversity of HPV prevalence in VC has been reported, but no data are available concerning the insertional HPV status in this tumor type. Therefore, we have used an innovative next generation sequencing (NGS)-based CaptHPV method able to provide an extensive characterization of HPV DNA in tumors.

MATERIAL AND METHODS

Tumor tissue specimens from 55 patients with VC were analyzed using p16 immunohistochemistry, in situ hybridization, polymerase chain reaction, and CaptHPV-NGS assays.

RESULTS

Our analyses showed that 8 (14.5%) of 55 cases were associated with HPV 16 DNA. No other HPV genotypes were identified. The HPV genome was in a free episomal state only in one case and both episomal and integrated into the tumor cell genome in 7. There was a single insertion in 5 cases and multiple sites, scattered at different chromosomal loci in two. ISH data suggest that some of these might reflect tumor heterogeneity. Viral integration targeted cellular genes among which were TP63, CCDC148, LOC100133091, PKP1, and POLA2. Viral integration at the PKP1 locus was associated with partial gene deletion, and no PKP1 protein was detected in tumor tissue.

CONCLUSIONS

Using the NGS-based innovative capture-HPV approach, we established a cartography of HPV 16 DNA in 8 VC cases and identified novel genes targeted by integration that may be used as specific tumor markers. In addition, we established a rationale strategy for optimal characterization of HPV status in VC.

Preferential expression of a HPV genotype in invasive cervical carcinomas infected by multiple genotypes

Multiple infections by HPV genotypes are frequently detected in HPV+ cervical lesions but the interaction between each viral genotype during carcinogenesis is poorly understood. Here we carried out a comprehensive study to characterize the multiple HPV expression and integration by RNA-seq analyses of 19 invasive cervical carcinomas coinfected by several HPV genotypes. Analysis of tumor DNA by a hybridization assay indicated multiple infections ranging from 2 to 6 different HPV genotypes. RNA-seq analysis showed that a single HPV genotype was preferentially expressed. Finally, the search for HPV/human chimeric transcripts indicated integration from preferentially expressed genotypes. In conclusion, the present study indicated that, in invasive cervical carcinomas infected by multiple HPV genotypes, one HPV was preferentially expressed, supporting the hypothesis that a single HPV genotype was associated with cancer development.

3'UTR polymorphism of Thymidylate Synthase gene increased the risk of persistence of pre-neoplastic cervical lesions

BACKGROUND

Cervical cancer is caused by high-risk Human Papillomavirus (hr-HPV) infection associated with cofactors that has been analyzed as predictors of the remission or persistence of cytological abnormalities remission or persistence. These cofactors can be either environmental, epigenetic, or genetic. Polymorphism in genes of enzymes that act on one-carbon metabolism alter their activity and also may be associated with cervical carcinogenesis because they affect DNA synthesis and repair, and gene expression. Therefore, this study aimed to analyze the risk of persistence of pre-neoplastic cervical lesions according to genetic polymorphisms involved in one-carbon metabolism.

METHODS

Our sample consisted of 106 women, divided into two groups - Remission (n = 60), i.e., with the presence of pre-neoplastic lesions at first meeting (T₁) and normal cytology after 6 months of follow-up (T₂), and Persistence (n = 46), i.e., with the presence of pre-neoplastic lesions at T₁ and T₂. We obtained cervical samples for cytological analysis (T₁ and T₂), HPV detection (T₁), and evaluation of polymorphism C667T of Methylenetetrahydrofolate Reductase (MTHFR C677T), A2756G of Methionine Synthase (MS A2756G), A66G of Methionine Synthase Reductase (MTRR A66G), double or triple 28 bp tandem repeat in 5'-untranslated enhanced region of Thymidylate Synthase (TSER), and 6 bp deletion at nucleotide 1494 in TS 3'-untranslated region (TS3'UTR). To analyze all genetic polymorphisms simultaneously, we calculated the Genetic Risk Score (GRS).

RESULTS

We observed no differences between the Remission and Persistence groups regarding the GRS. Also, there were no differences in the genotypic and allelic distribution of MTHFR C677T and MS A2756G polymorphisms. However, the risk of persistence was higher among women with the heterozygote genotype - ins/del [OR (IC95%): 3.22 (1.19-8.69), p = 0.021], or the polymorphic genotype - del/del [OR (IC95%): 6.50 (1.71-24.70), p = 0.006] of TS3'UTR.

CONCLUSIONS

The presence of the TS3'UTR polymorphism increased the risk of persistence of cervical abnormalities. This genetic variant could be a potential marker of cervical carcinogenesis and therefore assist the follow-up of women with persistent pre-neoplastic cervical lesions.

Insertional oncogenesis by HPV70 revealed by multiple genomic analyses in a clinically HPV-negative cervical cancer

Cervical carcinogenesis, the second leading cause of cancer death in women worldwide, is caused by multiple types of human papillomaviruses (HPVs). To investigate a possible role for HPV in a cervical carcinoma that was HPV-negative by PCR testing, we performed HPV DNA hybridization capture plus massively parallel sequencing. This detected a subgenomic, URR-E6-E7-E1 segment of HPV70 DNA, a type not generally associated with cervical cancer, inserted in an intron of the B-cell lymphoma/leukemia 11B (BCL11B) gene in the human genome. Long range DNA sequencing confirmed the virus and flanking BCL11B DNA structures including both insertion junctions. Global transcriptomic analysis detected multiple, alternatively spliced, HPV70-BCL11B, fusion transcripts with fused open reading frames. The insertion and fusion transcripts were present in an intraepithelial precursor phase of tumorigenesis. These results suggest oncogenicity of HPV70, identify novel BCL11B variants with potential oncogenic implications, and underscore the advantages of thorough genomic analyses to elucidate insights into HPV-associated tumorigenesis.

The association of integration patterns of human papilloma virus and single nucleotide polymorphisms on immune- or DNA repair-related genes in cervical cancer patients

The present study investigated the association between single nucleotide polymorphisms (SNPs) in immune- or DNA repair-related genes and the integration pattern of human papillomavirus (HPV), a promising prognostic marker in cervical cancer. The HPV integration patterns of cervical cancer patients were determined by polymerase chain reaction and in situ hybridization, and categorized as episomal (group A), single-copy or multi-copy tandem repetition integrated (group B), and undetectable HPV types (group C). After sample and SNP quality control, 166,505 SNPs in 161 samples (38, 111, and 12 patients in groups A, B, and C, respectively) were examined. None of the SNPs reached genome-wide significance, and several candidate SNPs for future study were selected, including rs10999435 on chromosome 10q22, rs1322054 on chromosome 9q32-33, and rs10902171 on chromosome 11p15. Luciferase assay identified rs1322054 as the primary functional variant to regulate gene expression in immune cell. Further studies are needed to determine the genetic background of different integration patterns of HPV in cervical cancer patients.

Lack of association between methylenetetrahydrofolate reductase C677T polymorphism, HPV infection and cervical intraepithelial neoplasia in Brazilian women

BACKGROUND

Cervical cancer has high prevalence and mortality rates in worldwide female population. Persistent infection by high-risk Human Papillomavirus (hr-HPV) is the main cause of this cancer. However, many environmental, genetical, and epigenetical cofactors can modulate viral infection and cervical carcinogenesis. Methylenetetrahydrofolate reductase (MTHFR) C677T polymorphism is a genetic factor that has been associated with many pathologies, including cancer. Nevertheless, studies with cervical cancer presented controversial results, and varied according to ethnicity. Thus, the aim of this study was to determine association between MTHFR C677T polymorphism, Human Papillomavirus (HPV) infection and cervical cancer.

METHODS

A case-control study was performed with 150 histological cervical samples. Case group were divided in Cervical Intraepithelial Neoplasia (CIN) grade I (n = 30), CIN II (n = 30), CIN III (n = 30), and Squamous Cervical Carcinoma (SCC) (n = 30). Control group was composed by 30 samples without lesion, presenting cervicitis. HPV detection was performed by conventional Polymerase Chain Reaction (PCR) with SPF primers set, and by real-time PCR specific for HPV 16 and hr-HPV. MTHFR C677T polymorphism was analyzed by PCR followed by Restriction Fragment Length Polymorphism (RFLP).

RESULTS

Frequency of MTHFR CC genotype was 72.7% (n = 109), CT 23.3% (n = 35) and TT 4.0% (n = 6). Polymorphic T allele frequency was 15.7%. No statistically significant association was observed between MTHFR C677T polymorphism and presence of pre-neoplastic or neoplastic cervical lesions. Similar frequencies of T allele was observed in control (23.3%) and cases (13.3%) groups (p = 0.174). In addition, there was no statistically significant association between MTHFR C677T polymorphism and viral infection, even considering hr-HPV or HPV 16 positivity.

CONCLUSION

MTHFR C677T polymorphism was not associated with cervical cancer and HPV infection.

HPV DNA integration site as proof of the origin of ovarian metastasis from endocervical adenocarcinoma: three case reports

Background

Most endocervical adenocarcinomas are human papillomavirus (HPV)-related cancers associated with p16 immunostaining. Ovarian metastasis from cervical cancer is a rare phenomenon, the mechanism of dissemination remains unclear. The diagnosis of metastasis may be difficult to establish when the ovarian neoplasm presents features consistent with primary tumor. Immunohistochemical expression of p16 in ovarian tumors can guide the diagnosis of metastasis from HPV-related cervical cancer, but p16 positivity is nonspecific. Identical HPV genotype in the paired endocervical and ovarian tumors is a better marker for cervical origin, which may also be confirmed by identical HPV integration site.

Case presentation

Two women presented with HPV18 cervical adenocarcinoma. No signs of disease were visible on MRI after treatment. After several years of follow-up, mucinous ovarian tumors were discovered in both patients. Molecular analyses showed that the ovarian lesions were HPV18-positive; indicating a primary cervical origin. A third woman was diagnosed with grade 1 ovarian endometrioid carcinoma with no peritoneal carcinomatosis. Final histological examination and HPV genotyping revealed HPV18-related in situ endometrioid adenocarcinoma in the endocervix and HPV18-related invasive endometrioid adenocarcinoma in the endometrium and both ovaries. Additional molecular analyses performed in two patients identified the same HPV integration sites in both the ovarian and cervical tumors, confirming that the ovarian mass was a metastasis from the cervical adenocarcinoma.

Conclusion

We report three new cases of ovarian neoplasia in which the diagnosis of metastasis from cervical cancer was supported by the same HPV genotype and the same integration site in the paired cervical and ovarian tumors. To our knowledge, this is the first report of molecular evidence of the cervical origin of an ovarian metastasis. HPV screening should be performed in ovarian tumors for all patients with history of cervical neoplasia.

Characterization of six Merkel cell polyomavirus-positive Merkel cell carcinoma cell lines: Integration pattern suggest that large T antigen truncating events occur before or during integration

Merkel cell carcinoma (MCC), an aggressive neuroendocrine skin tumor, is a polyomavirus-induced human cancer. To study the causal relationship of MCC carcinogenesis with the integrated Merkel cell polyomavirus (MCPyV) in detail, well-characterized MCC cell lines are needed. Consequently, in the current study, we established and characterized six MCPyV-positive MCC cell lines. Microarray-based comparative genomic hybridization revealed a stable genome carrying only a limited number of chromosomal gains and deletions. All cell lines expressed MCC markers Keratin-20 and neuron-specific enolase as well as truncated MCPyV-encoded large T antigen (LT). For five cell lines, we were able to identify the MCPyV-integration sites in introns of different genes. The LT-truncating stop codon mutations and integration sites were affirmed in the respective clinical patient samples. Inverse PCR suggested that three of the cell lines contained MCPyV genomes as concatemers. This notion was confirmed for the two cell lines with known integration sites. Importantly, our observation of distinct stop codon mutations in cell lines with concatemeric MCPyV integration indicates that these LT-truncating mutations occur before integration. In summary, we provide the detailed characterization of six MCPyV-positive MCC cell lines, which are likely to serve as valuable tools in future MCC research.

HPV insertional pattern as a personalized tumor marker for the optimized tumor diagnosis and follow-up of patients with HPV-associated carcinomas: a case report

BACKGROUND

In clinical oncology, only a few applications have been developed using HPV as a personalized tumor marker, a lack most probably related to the limited information obtained by the classical Polymerase Chain Reaction (PCR) approach. To overcome this limitation, we have recently developed the capture-based Next-Generation Sequencing (NGS) "CaptHPV" assay, designed to provide an extensive and comprehensive molecular characterization of HPV DNA sequences associated with neoplasias, ie the sequence of the viral genome (245 genotypes), its physical state, viral load, integration site and genomic alterations at integration locus. These data correspond to highly specific tumor markers that can be used to improve diagnosis and patient's follow-up.

CASE PRESENTATION

We report here a case that is a straightforward and practical illustration of the power of the CaptHPV method. A patient developed successively a carcinoma of the anal canal and of the tongue. The two tumors were squamous cell carcinoma, found associated with HPV16 using PCR. In order to document a possible metastasis to the tongue from the anal cancer, we performed CaptHPV analysis on the two tumors. The analysis of the anal carcinoma found 55 viral/human hybrid reads allowing the identification of the HPV16 DNA integration in the 4q25 chromosomal band locus with a 178,808 bp deletion in the cell genome. Molecular analysis of the tongue tumor disclosed 6110 reads of HPV16, with a viral pattern strictly identical to that of the anal tumor. A total of 131 hybrid reads between HPV16 and the cell genome were found, corresponding exactly to the same locus of integration of viral DNA at the 4q25 site. The 178,808 bp genomic deletion was also found in the lingual tumor. The exact identity of HPV insertional signatures in the two tumors, demonstrates unambiguously that the tongue tumor derived from the anal cancer whereas neither histological immunophenotyping nor classical viral analysis using PCR could allow a definitive diagnosis.

CONCLUSION

Our observation indicates that the establishment of a detailed cartography of HPV DNA sequences in a tumor specimen provides crucial information for the design of specific biomarkers that can be used for diagnostic, prognostic or predictive purposes.

Genome-wide profiling of human papillomavirus DNA integration in liquid-based cytology specimens from a Gabonese female population using HPV capture technology

Human papillomavirus (HPV) is recognised as the cause of precancerous and cancerous cervical lesions. Furthermore, in high-grade lesions, HPV is frequently integrated in the host cell genome and associated with the partial or complete loss of the E1 and E2 genes, which regulate the activity of viral oncoproteins E6 and E7. In this study, using a double-capture system followed by high-throughput sequencing, we determined the HPV integration status present in liquid-based cervical smears in an urban Gabonese population. The main inclusion criteria were based on cytological grade and the detection of the HPV16 genotype using molecular assays. The rate of HPV integration in the host genome varied with cytological grade: 85.7% (6/7), 71.4% (5/7), 66.7% (2/3) 60% (3/5) and 30.8% (4/13) for carcinomas, HSIL, ASCH, LSIL and ASCUS, respectively. For high cytological grades (carcinomas and HSIL), genotypes HPV16 and 18 represented 92.9% of the samples (13/14). The integrated form of HPV16 genotype was mainly found in high-grade lesions in 71.4% of samples regardless of cytological grade. Minority genotypes (HPV33, 51, 58 and 59) were found in LSIL samples, except HPV59, which was identified in one HSIL sample. Among all the HPV genotypes identified after double capture, 10 genotypes (HPV30, 35, 39, 44, 45, 53, 56, 59, 74 and 82) were detected only in episomal form. Our study revealed that the degree of HPV integration varies with cervical cytological grade. The integration event might be a potential clinical prognostic biomarker for the prediction of the progression of neoplastic lesions.

TaME-seq: An efficient sequencing approach for characterisation of HPV genomic variability and chromosomal integration

HPV genomic variability and chromosomal integration are important in the HPV-induced carcinogenic process. To uncover these genomic events in an HPV infection, we have developed an innovative and cost-effective sequencing approach named TaME-seq (tagmentation-assisted multiplex PCR enrichment sequencing). TaME-seq combines tagmentation and multiplex PCR enrichment for simultaneous analysis of HPV variation and chromosomal integration, and it can also be adapted to other viruses. For method validation, cell lines (n = 4), plasmids (n = 3), and HPV16, 18, 31, 33 and 45 positive clinical samples (n = 21) were analysed. Our results showed deep HPV genome-wide sequencing coverage. Chromosomal integration breakpoints and large deletions were identified in HPV positive cell lines and in one clinical sample. HPV genomic variability was observed in all samples allowing identification of low frequency variants. In contrast to other approaches, TaME-seq proved to be highly efficient in HPV target enrichment, leading to reduced sequencing costs. Comprehensive studies on HPV intra-host variability generated during a persistent infection will improve our understanding of viral carcinogenesis. Efficient identification of both HPV variability and integration sites will be important for the study of HPV evolution and adaptability and may be an important tool for use in cervical cancer diagnostics.

Association of Base Excision Repair Gene hOGG1 Ser326Cys Polymorphism with Susceptibility to Cervical Squamous Cell Carcinoma and High-Risk Human Papilloma Virus Infection in a Chinese Population

AIM

This study investigated the association of the human 8-oxoguanine glycosylase 1 (hOGG1) Ser326Cys polymorphism with risk of cervical squamous cell carcinoma (CSCC) and high-risk human papilloma virus (HR-HPV) infection.

BACKGROUND

The hOGG1 Ser326Cys polymorphism is reported to be correlated with the risk of several cancers. However, there are reports that have found no significant differences in the frequency of the hOGG1 Ser326Cys between cervical carcinoma patients and controls.

METHODS

hOGG1 Ser326Cys was genotyped through modified allele mismatch amplification polymerase chain reaction in 1200 healthy controls, 400 cervical intraepithelial neoplasia (CIN) grade III cases, and 400 CSCC cases.

RESULTS

The homozygous genotype of hOGG1 Cys326Cys (GG) was associated with increased risk of CIN III (odds ratio, OR = 1.81 [1.31-2.49], p < 0.001) and CSCC (OR = 3.05 [2.2-4.20], p < 0.001). The G allele or G carrier (GG + CG) genotype was a highly-significant risk factor for CSCC (OR = 1.49 [1.14-1.97], p = 0.004). In the HR-HPV-positive group, the homozygous genotype of hOGG1 GG was associated with increased risk of CSCC (OR = 3.66 [2.02-6.62], p < 0.001) and risk of CIN III (OR = 1.82 [1.08-3.06], p = 0.024). The proportion of G allele carriers was significantly increased in CIN III (51.9%, [322/620], OR = 1.33 [1.03-1.72], p = 0.028) and CSCC (62.1% [221/356], OR = 2.02 [1.51-2.71], p < 0.001). The GG and GC genotypes were consistently identified as significant risk factors for CSCC (OR = 1.73 [1.06-2.83], p = 0.029) in the HR-HPV infected group. We further observed enrichment of the hOGG1 Ser326Cys polymorphism in the CIN III (p = 0.021) and CSCC (p < 0.001) stratified by age at first intercourse, with more significant enrichment (p = 0.036) in the HR-HPV infection group.

CONCLUSIONS

Our findings support associations of the hOGG1 Ser326Cys polymorphism with CSCC carcinogenesis and susceptibility to HR-HPV infection. The hOGG1 Ser326Cys polymorphism may serve as a potential genetic biomarker of susceptibility to cervical cancer and HR-HPV infection.

Characterization of HPV integration, viral gene expression and E6E7 alternative transcripts by RNA-Seq: A descriptive study in invasive cervical cancer

Scarce data are available on the expression of papillomavirus genome and the frequency of alternatively spliced E6E7 mRNAs in invasive cervical cancer. We carried out a comprehensive characterization of HPV expression by RNA-Seq analysis in 22 invasive cervical cancer with HPV16 or HPV18, characterizing the presence of integrated/episomal viral DNA, the integration sites in human genome and the proportion of alternative splicing products of E6 and E7 genes. The expression patterns suggested the presence of episomal and/or integrated viral DNA, with integration detected in most tumors, frequently occurring within human genes in HPV18+ and in intergenic regions in HPV16+ tumors. Alternative splicing of E6E7 transcripts showed E6*I as the most frequent isoform for both viral types, followed by E6*II and E6/E7 (unspliced) transcripts in HPV16+, and by E6/E7 in HPV18+ tumors. Previously described E6*VI and E6*V transcript isoforms for HPV16, and E6*X for HPV18, were rare or not detected.

The precision prevention and therapy of HPV-related cervical cancer: new concepts and clinical implications

Cervical cancer is the third most common cancer in women worldwide, with concepts and knowledge about its prevention and treatment evolving rapidly. Human papillomavirus (HPV) has been identified as a major factor that leads to cervical cancer, although HPV infection alone cannot cause the disease. In fact, HPV-driven cancer is a small probability event because most infections are transient and could be cleared spontaneously by host immune system. With persistent HPV infection, decades are required for progression to cervical cancer. Therefore, this long time window provides golden opportunity for clinical intervention, and the fundament here is to elucidate the carcinogenic pattern and applicable targets during HPV-host interaction. In this review, we discuss the key factors that contribute to the persistence of HPV and cervical carcinogenesis, emerging new concepts and technologies for cancer interventions, and more urgently, how these concepts and technologies might lead to clinical precision medicine which could provide prediction, prevention, and early treatment for patients.

Nuclear factor I X is a recurrent target for HPV16 insertions in anal carcinomas

Anal carcinomas (AC) are associated with human papillomavirus (HPV) DNA sequences, but little is known about the physical state of the viral genome in carcinoma cells. To define the integration status and gene(s) targeted by viral insertions in AC, tumor DNAs extracted from 35 tumor specimen samples in patients with HPV16-associated invasive carcinoma were analyzed using the detection of integrated papillomavirus sequences-PCR approach. The genomic status at integration sites was assessed using comparative genomic hybridization-array assay and gene expression using reverse transcription quantitative PCR (RT-qPCR). HPV16 DNA was found integrated in 25/35 (71%) cases and the integration locus could be determined at the molecular level in 19 cases (29 total integration loci). HPV DNA was inserted on different chromosomes, but 5 cases harbored viral sequences at 19p13.2, within the nuclear factor I X (NFIX) locus. Viral DNA mapped between the most distal and the two proximal alternatively expressed exons of this gene in three cases (CA21, CA04, and CA35) and upstream of this gene (663 kb and 2.3 Mb) in the others. CGH arrays showed genomic gains/amplifications at the NFIX region, associated with HPV within the gene and RT-qPCR, revealed NFIX mRNA overexpression. Other genes targeted by integration were IL20RB, RPS6KA2, MSRA1, PIP5K1B, SLX4IP, CECR1, BCAR3, ATF6, CSNK1G1, APBA2, AGK, ILF3, PVT1, TRMT1, RAD51B, FASN, CCDC57, DSG3, and ZNF563. We identified recurrent targeting of NFIX by HPV16 insertion in anal carcinomas, supporting a role for this gene in oncogenesis, as reported for non-HPV tumors.

Patterns and mechanisms of structural variations in human cancer

Next-generation sequencing technology has enabled the comprehensive detection of genomic alterations in human somatic cells, including point mutations, chromosomal rearrangements, and structural variations (SVs). Using sophisticated bioinformatics algorithms, unbiased catalogs of SVs are emerging from thousands of human cancer genomes for the first time. Via careful examination of SV breakpoints at single-nucleotide resolution as well as local DNA copy number changes, diverse patterns of genomic rearrangements are being revealed. These "SV signatures" provide deep insight into the mutational processes that have shaped genome changes in human somatic cells. This review summarizes the characteristics of recently identified complex SVs, including chromothripsis, chromoplexy, microhomology-mediated breakage-induced replication (MMBIR), and others, to provide a holistic snapshot of the current knowledge on genomic rearrangements in somatic cells.

Next generation sequencing and its applications in HPV-associated cancers

Approximately 18% of all human cancers have a viral etiology, and human papillomavirus (HPV) has been identified as one of the most prevalent viruses that plays causative role in nearly all cervical cancers and, in addition, in subset of head and neck, anal, penile and vulvar cancers. The recent introduction of next generation sequencing (NGS) and other omics approaches have resulted in comprehensive knowledge on the pathogenesis of HPV-driven tumors. Specifically, these approaches have provided detailed information on genomic HPV integration sites, disrupted genes and pathways, and common and distinct genetic and epigenetic alterations in different human HPV-associated cancers. This review focuses on HPV integration sites, its concomitantly disrupted genes and pathways and its functional consequences in both cervical and head and neck cancers. Integration of NGS data with other omics and clinical data is crucial to better understand the pathophysiology of each individual malignancy and, based on this, to select targets and to design effective personalized treatment options.

Synthetic siRNA targeting human papillomavirus 16 E6: a perspective on in vitro nanotherapeutic approaches

High-risk human papillomaviruses infect skin and mucosa, causing approximately 5% of cancers worldwide. In the search for targeted nanotherapeutic approaches, siRNAs against the viral E6 transcript have been molecules of interest but have not yet seen successful translation into the clinic. By reviewing the past approximately 15 years of in vitro literature, we identify the need for siRNA validation protocols which concurrently evaluate ranges of key treatment parameters as well as characterize downstream process restoration in a methodical, quantitative manner and demonstrate their implementation using our own data. We also reflect on the future need for more appropriate cell culture models to represent patient lesions as well as the application of personalized approaches to identify optimal treatment strategies.

Identification of reliable biomarkers of human papillomavirus 16 methylation in cervical lesions based on integration status using high-resolution melting analysis

BACKGROUND

The dynamic methylation of human papillomavirus (HPV) 16 DNA is thought to be associated with the progression of cervical lesions. Previous studies that did not consider the physical status of HPV 16 may have incorrectly mapped HPV 16 methylomes. In order to identify reliable biomarkers for squamous cervical cancer (SCC), we comprehensively evaluated the methylation of HPV 16 depending on the integration incidence of each sample.

METHODS

Based on the integration status of 115 HPV 16-infected patients (50 SCC, 30 high-grade squamous intraepithelial lesion [HSIL], and 35 low-grade squamous intraepithelial lesion [LSIL]) and HPV 16-infected Caski cell lines by PCR detection of integrated papillomavirus sequences, we designed a series of primers that would not be influenced by breakpoints for a high-resolution melting (HRM) PCR method to detect the genome methylation.

RESULTS

A few regions with recurrent interruptions were identified in E1, E2/E4, L1, and L2 despite scattering of breakpoints throughout all eight genes of HPV 16. Frequent integration sites often occurred concomitantly with methylated CpG sites. The HRM PCR method showed 100% agreement with pyrosequencing when 3% was set as the cutoff value. A panel of CpG sites such as nt5606, nt5609, nt5615, and nt5378 can be combined in reweighing calculations to distinguish SCC from HSIL and LSIL patients which have high sensitivity and specificity (88% and 92.31%, respectively).

CONCLUSIONS

Our research shows that combination of CpG sites nt5606, nt5609, nt5615, and nt5378 can be used as potential diagnosis biomarkers for SCC, and the HRM PCR method is suitable for clinical methylation analysis.

HPV integration hijacks and multimerizes a cellular enhancer to generate a viral-cellular super-enhancer that drives high viral oncogene expression

Integration of human papillomavirus (HPV) genomes into cellular chromatin is common in HPV-associated cancers. Integration is random, and each site is unique depending on how and where the virus integrates. We recently showed that tandemly integrated HPV16 could result in the formation of a super-enhancer-like element that drives transcription of the viral oncogenes. Here, we characterize the chromatin landscape and genomic architecture of this integration locus to elucidate the mechanisms that promoted de novo super-enhancer formation. Using next-generation sequencing and molecular combing/fiber-FISH, we show that ~26 copies of HPV16 are integrated into an intergenic region of chromosome 2p23.2, interspersed with 25 kb of amplified, flanking cellular DNA. This interspersed, co-amplified viral-host pattern is frequent in HPV-associated cancers and here we designate it as Type III integration. An abundant viral-cellular fusion transcript encoding the viral E6/E7 oncogenes is expressed from the integration locus and the chromatin encompassing both the viral enhancer and a region in the adjacent amplified cellular sequences is strongly enriched in the super-enhancer markers H3K27ac and Brd4. Notably, the peak in the amplified cellular sequence corresponds to an epithelial-cell-type specific enhancer. Thus, HPV16 integration generated a super-enhancer-like element composed of tandem interspersed copies of the viral upstream regulatory region and a cellular enhancer, to drive high levels of oncogene expression.

HPV Integration in Head and Neck Squamous Cell Carcinomas: Cause and Consequence

Human papillomaviruses (HPVs) are a necessary cause of anogenital squamous cell carcinomas (SCC) and a subgroup of head and neck SCC, i.e., those originating in the oropharynx. The key events in high-risk HPV (HRHPV)-associated neoplastic progression include persistent infection, deregulated expression of virus early genes in basal epithelial cells, local immune suppression and the accumulation of chromosomal alterations. Evidence for these events particularly comes from studies of uterine cervical carcinogenesis; primary premalignant HRHPV-positive lesions of the head and neck mucosa are seldomly detected. Integration of virus DNA into host chromosomes is considered an important driver of carcinogenesis and observed in 40 up to 90 % of uterine cervical SCC (UCSCC) and oropharyngeal SCC (OPSCC), dependent on the integration detection method used and HRHPV type. In OPSCC, > 90 % HPV-positive tumors are infected with HPV16. Ten up to 60 % of HPV-positive tumors thus contain extrachromosomal (episomal) virus. In this chapter, causes and consequences of HPV integration are summarized from the literature, with special focus on the site of HPV integration in the cellular genome, and its effect on expression of viral oncogenes (particularly E6 and E7), on human (tumor) gene expression and on deregulation of cell proliferation, apoptosis and cell signaling pathways. Also data on DNA methylation, viral load and clinical outcome in relation to HPV integration are provided.

Molecular approaches for HPV genotyping and HPV-DNA physical status

Persistent infection with high-risk human papillomavirus (HPV) genotypes is the leading cause of cervical cancer development. To this end several studies have focused on designing molecular assays for HPV genotyping, which are considered as the gold standard for the early diagnosis of HPV infection. Moreover, the tendency of HPV DNA to be integrated into the host chromosome is a determining event for cervical oncogenesis. Thus, the establishment of molecular techniques was promoted in order to investigate the physical status of the HPV DNA and the locus of viral insertion into the host chromosome. The molecular approaches that have been developed recently facilitate the collection of a wide spectrum of valuable information specific to each individual patient and therefore can significantly contribute to the establishment of a personalised prognosis, diagnosis and treatment of HPV-positive patients. The present review focuses on state of the art molecular assays for HPV detection and genotyping for intra-lesion analyses, it examines molecular approaches for the determination of HPV-DNA physical status and it discusses the criteria for selecting the most appropriate regions of viral DNA to be incorporated in HPV genotyping and in the determination of HPV-DNA physical status.

Comprehensive mapping of the human papillomavirus (HPV) DNA integration sites in cervical carcinomas by HPV capture technology

Integration of human papillomavirus (HPV) DNA into the host genome can be a driver mutation in cervical carcinoma. Identification of HPV integration at base resolution has been a longstanding technical challenge, largely due to sensitivity masking by HPV in episomes or concatenated forms. The aim was to enhance the understanding of the precise localization of HPV integration sites using an innovative strategy. Using HPV capture technology combined with next generation sequencing, HPV prevalence and the exact integration sites of the HPV DNA in 47 primary cervical cancer samples and 2 cell lines were investigated. A total of 117 unique HPV integration sites were identified, including HPV16 (n = 101), HPV18 (n = 7), and HPV58 (n = 9). We observed that the HPV16 integration sites were broadly located across the whole viral genome. In addition, either single or multiple integration events could occur frequently for HPV16, ranging from 1 to 19 per sample. The viral integration sites were distributed across almost all the chromosomes, except chromosome 22. All the cervical cancer cases harboring more than four HPV16 integration sites showed clinical diagnosis of stage III carcinoma. A significant enrichment of overlapping nucleotides shared between the human genome and HPV genome at integration breakpoints was observed, indicating that it may play an important role in the HPV integration process. The results expand on knowledge from previous findings on HPV16 and HPV18 integration sites and allow a better understanding of the molecular basis of the pathogenesis of cervical carcinoma.

Functional variants of human papillomavirus type 16 demonstrate host genome integration and transcriptional alterations corresponding to their unique cancer epidemiology

BACKGROUND

Human papillomaviruses (HPVs) are a worldwide burden as they are a widespread group of tumour viruses in humans. Having a tropism for mucosal tissues, high-risk HPVs are detected in nearly all cervical cancers. HPV16 is the most common high-risk type but not all women infected with high-risk HPV develop a malignant tumour. Likely relevant, HPV genomes are polymorphic and some HPV16 single nucleotide polymorphisms (SNPs) are under evolutionary constraint instigating variable oncogenicity and immunogenicity in the infected host.

RESULTS

To investigate the tumourigenicity of two common HPV16 variants, we used our recently developed, three-dimensional organotypic model reminiscent of the natural HPV infectious cycle and conducted various "omics" and bioinformatics approaches. Based on epidemiological studies we chose to examine the HPV16 Asian-American (AA) and HPV16 European Prototype (EP) variants. They differ by three non-synonymous SNPs in the transforming and virus-encoded E6 oncogene where AAE6 is classified as a high- and EPE6 as a low-risk variant. Remarkably, the high-risk AAE6 variant genome integrated into the host DNA, while the low-risk EPE6 variant genome remained episomal as evidenced by highly sensitive Capt-HPV sequencing. RNA-seq experiments showed that the truncated form of AAE6, integrated in chromosome 5q32, produced a local gene over-expression and a large variety of viral-human fusion transcripts, including long distance spliced transcripts. In addition, differential enrichment of host cell pathways was observed between both HPV16 E6 variant-containing epithelia. Finally, in the high-risk variant, we detected a molecular signature of host chromosomal instability, a common property of cancer cells.

CONCLUSIONS

We show how naturally occurring SNPs in the HPV16 E6 oncogene cause significant changes in the outcome of HPV infections and subsequent viral and host transcriptome alterations prone to drive carcinogenesis. Host genome instability is closely linked to viral integration into the host genome of HPV-infected cells, which is a key phenomenon for malignant cellular transformation and the reason for uncontrolled E6 oncogene expression. In particular, the finding of variant-specific integration potential represents a new paradigm in HPV variant biology.

Human papillomavirus (HPV) genome status & cervical cancer outcome--A retrospective study

Background & objectives:

Persistent infections with high-risk (HR) human papillomaviruses such as HPV 16, 18, 31, 33 and 45 have been identified as the major aetiological factor for cervical cancer. The clinical outcome of the disease is often determined by viral factors such as viral load, physical status and oncogene expression. The aim of the present study was to evaluate the impact of such factors on clinical outcome in HPV16 positive, locally advanced cervical cancer cases.

Methods:

One hundred and thirty two pretreatment cervical tumour biopsies were selected from patients undergoing radiotherapy alone (n=63) or concomitant chemo-radiation (n=69). All the samples were positive for HPV 16. Quantitative real time-PCR was carried out to determine viral load and oncogene expression. Physical status of the virus was determined for all the samples by the ratio of E2_copies/E7_copies; while in 73 cases, the status was reanalyzed by more sensitive APOT (amplification of papillomavirus oncogene transcripts) assay. Univariate analysis of recurrence free survival was carried out using Kaplan-Meier method and for multivariate analysis the Cox proportional hazard model was used.

Results:

The median viral load was 19.4 (IQR, 1.9- 69.3), with viral integration observed in 86 per cent cases by combination of the two methodologies. Both univariate and multivariate analyses identified viral physical status as a good predictor of clinical outcome following radiation treatment, with episomal form being associated with increased recurrence free survival.

Interpretation & conclusions:

The present study results showed that viral physical status might act as an important prognostic factor in cervical cancer.

Genomic characterization of viral integration sites in HPV-related cancers

Persistent infection with carcinogenic human papillomaviruses (HPV) causes the majority of anogenital cancers and a subset of head and neck cancers. The HPV genome is frequently found integrated into the host genome of invasive cancers. The mechanisms of how it may promote disease progression are not well understood. Thoroughly characterizing integration events can provide insights into HPV carcinogenesis. Individual studies have reported limited number of integration sites in cell lines and human samples. We performed a systematic review of published integration sites in HPV-related cancers and conducted a pooled analysis to formally test for integration hotspots and genomic features enriched in integration events using data from the Encyclopedia of DNA Elements (ENCODE). Over 1,500 integration sites were reported in the literature, of which 90.8% (N = 1,407) were in human tissues. We found 10 cytobands enriched for integration events, three previously reported ones (3q28, 8q24.21 and 13q22.1) and seven additional ones (2q22.3, 3p14.2, 8q24.22, 14q24.1, 17p11.1, 17q23.1 and 17q23.2). Cervical infections with HPV18 were more likely to have breakpoints in 8q24.21 (p = 7.68 × 10(-4) ) than those with HPV16. Overall, integration sites were more likely to be in gene regions than expected by chance (p = 6.93 × 10(-9) ). They were also significantly closer to CpG regions, fragile sites, transcriptionally active regions and enhancers. Few integration events occurred within 50 Kb of known cervical cancer driver genes. This suggests that HPV integrates in accessible regions of the genome, preferentially genes and enhancers, which may affect the expression of target genes.

Circulating human papillomavirus DNA detected using droplet digital PCR in the serum of patients diagnosed with early stage human papillomavirus-associated invasive carcinoma

Specific human papillomavirus genotypes are associated with most ano‐genital carcinomas and a large subset of oro‐pharyngeal carcinomas. Human papillomavirus DNA is thus a tumour marker that can be detected in the blood of patients for clinical monitoring. However, data concerning circulating human papillomavirus DNA in cervical cancer patients has provided little clinical value, due to insufficient sensitivity of the assays used for the detection of small sized tumours. Here we took advantage of the sensitive droplet digital PCR method to identify circulating human papillomavirus DNA in patients with human papillomavirus‐associated carcinomas.

A series of 70 serum specimens, taken at the time of diagnosis, between 2002 and 2013, were retrospectively analyzed in patients with human papillomavirus‐16 or human papillomavirus‐18‐associated carcinomas, composed of 47 cases from the uterine cervix, 15 from the anal canal and 8 from the oro‐pharynx. As negative controls, 18 serum samples from women with human papillomavirus‐16‐associated high‐grade cervical intraepithelial neoplasia were also analyzed. Serum samples were stored at −80°C (27 cases) or at −20°C (43 cases). DNA was isolated from 200 µl of serum or plasma and droplet digital PCR was performed using human papillomavirus‐16 E7 and human papillomavirus‐18 E7 specific primers.

Circulating human papillomavirus DNA was detected in 61/70 (87%) serum samples from patients with carcinoma and in no serum from patients with cervical intraepithelial neoplasia. The positivity rate increased to 93% when using only serum stored at −80°C. Importantly, the two patients with microinvasive carcinomas in this series were positive. Quantitative evaluation showed that circulating viral DNA levels in cervical cancer patients were related to the clinical stage and tumour size, ranging from 55 ± 85 copies/ml (stage I) to 1774 ± 3676 copies/ml (stage IV).

Circulating human papillomavirus DNA is present in patients with human papillomavirus‐associated invasive cancers even at sub‐clinical stages and its level is related to tumour dynamics. Droplet digital PCR is a promising method for circulating human papillomavirus DNA detection and quantification. No positivity was found in patients with human papillomavirus‐associated high grade cervical intraepithelial neoplasia.

Mechanistic signatures of HPV insertions in cervical carcinomas

To identify new personal biomarkers for the improved diagnosis, prognosis and biological follow-up of human papillomavirus (HPV)-associated carcinomas, we developed a generic and comprehensive Capture-HPV method followed by Next Generation Sequencing (NGS). Starting from biopsies or circulating DNA samples, this Capture-NGS approach rapidly identifies the HPV genotype, HPV status (integrated, episomal or absence), the viral-host DNA junctions and the associated genome rearrangements. This analysis of 72 cervical carcinomas identified five HPV signatures. The first two signatures contain two hybrid chromosomal-HPV junctions whose orientations are co-linear (2J-COL) or non-linear (2J-NL), revealing two modes of viral integration associated with chromosomal deletion or amplification events, respectively. The third and fourth signatures exhibit 3-12 hybrid junctions, either clustered in one locus (MJ-CL) or scattered at distinct loci (MJ-SC) while the fifth signature consists of episomal HPV genomes (EPI). Cross analyses between the HPV signatures and the clinical and virological data reveal unexpected biased representation with respect to the HPV genotype, patient age and disease outcome, suggesting functional relevance(s) of this new classification. Overall, our findings establish a facile and comprehensive rational approach for the molecular detection of any HPV-associated carcinoma and definitive personalised sequence information to develop sensitive and specific biomarkers for each patient.