Sites of disruption within E1 and E2 genes of HPV16 and association with cervical dysplasia

Phylogenetic analysis and antigenic epitope prediction for E6 and E7 of Alpha-papillomavirus 9 in Taizhou, China

BACKGROUND

Alpha-papillomavirus 9 (α-9) is a member of the human papillomavirus (HPV) α genus, causing 75% invasive cervical cancers worldwide. The purpose of this study was to provide data for effective treatment of HPV-induced cervical lesions in Taizhou by analysing the genetic variation and antigenic epitopes of α-9 HPV E6 and E7.

METHODS

Cervical exfoliated cells were collected for HPV genotyping. Positive samples of the α-9 HPV single type were selected for E6 and E7 gene sequencing. The obtained nucleotide sequences were translated into amino acid sequences (protein primary structure) using MEGA X, and positive selection sites of the amino acid sequences were evaluated using PAML. The secondary and tertiary structures of the E6 and E7 proteins were predicted using PSIPred, SWISS-MODEL, and PyMol. Potential T/B-cell epitopes were predicted by Industrial Engineering Database (IEDB).

RESULTS

From 2012 to 2023, α-9 HPV accounted for 75.0% (7815/10423) of high-risk HPV-positive samples in Taizhou, both alone and in combination with other types. Among these, single-type-positive samples of α-9 HPV were selected, and the entire E6 and E7 genes were sequenced, including 298 HPV16, 149 HPV31, 185 HPV33, 123 HPV35, 325 HPV52, and 199 HPV58 samples. Compared with reference sequences, 34, 12, 10, 2, 17, and 17 nonsynonymous nucleotide mutations were detected in HPV16, 31, 33, 35, 52, and 58, respectively. Among all nonsynonymous nucleotide mutations, 19 positive selection sites were selected, which may have evolutionary significance in rendering α-9 HPV adaptive to its environment. Immunoinformatics predicted 57 potential linear and 59 conformational B-cell epitopes, many of which are also predicted as CTL epitopes.

CONCLUSION

The present study provides almost comprehensive data on the genetic variations, phylogenetics, positive selection sites, and antigenic epitopes of α-9 HPV E6 and E7 in Taizhou, China, which will be helpful for local HPV therapeutic vaccine development.

Prevalence and genotype distribution of HPV infections among women in Chengdu,China

BACKGROUND

Human papilloma virus (HPV) infection among female is the cause of cervical cancer and genital warts. In China, the HPV vaccination rate and the target population screening rate among females are low, and the aims of this study on the genotype distribution and prevalence of HPV infection were to provide more targeted strategies for the prevention and treatment of cervical cancer and HPV-related diseases.

METHODS

Polymerase chain reaction-reverse dot blot (PCR-RDB) was adopted for HPV genotyping test, the prevalence and 23 genotypes distribution of HPV infections among 181,705 women in Chengdu from 2013 to 2020 were analysed.

RESULTS

The overall prevalence rate of HPV infection among 181,705 cases was 23.28%, the prevalence of HR-HPV at the age group < 20 years, 60-69 years and ≥ 70 years were higher than the overall prevalence.The prevalence of HPV showed a bimodal U-shaped curve with age; the first and second peak common occurred among females < 20 years old (42.97%) and 60-69 years old (37.56%), respectively.The top five genotypes of HPV infection among females in Chengdu were HPV52/16/58/81/53. Single infection (73.26%) was the main HPV infection pattern, followed by double infection (19.17%) and multiple infection (7.57%), the infection rate of HPV showed a gradual declined as the patterns of HPV coinfections increased, low-risk and high-risk coinfection was higher in low-risk HPV infection (43.68%) and lower in high-risk HPV infection (13.59%). The prevalence of genotypes - 6 and - 81 infection was the second highest at the age group of 20 and 40-59, respectively, while the prevalence of HPV16 was the highest at the age group of ≥ 70 among 23 genotypes among the 181,705 women.

CONCLUSIONS

The prevalence of HPV infections among women in Chengdu is higher than domestic certain developed citys, among the five vaccines available, nonavalent vaccine is more suitable for Chengdu females. For young females prioritizing vaccination is essential in the current context.Double screening for HPV DNA is important in middle-aged women (30-49 years), and screening should not be lacking in older women (> 65 years). Additionally,for patients with genital warts, it is necessary to screen for high-risk HPV infection and provide appropriate management and treatment. Given the limitations of this study, future HPV research should aim to achieve full coverage of the target population, and our studies should also include cellular or pathological data of HPV-positive cases, vaccination rates, and various lifestyle details.

Long-read sequencing reveals the structural complexity of genomic integration of HPV DNA in cervical cancer cell lines

BACKGROUND

Cervical cancer (CC) causes more than 311,000 deaths annually worldwide. The integration of human papillomavirus (HPV) is a crucial genetic event that contributes to cervical carcinogenesis. Despite HPV DNA integration is known to disrupt the genomic architecture of both the host and viral genomes in CC, the complexity of this process remains largely unexplored.

RESULTS

In this study, we conducted whole-genome sequencing (WGS) at 55-65X coverage utilizing the PacBio long-read sequencing platform in SiHa and HeLa cells, followed by comprehensive analyses of the sequence data to elucidate the complexity of HPV integration. Firstly, our results demonstrated that PacBio long-read sequencing effectively identifies HPV integration breakpoints with comparable accuracy to targeted-capture Next-generation sequencing (NGS) methods. Secondly, we constructed detailed models of complex integrated genome structures that included both the HPV genome and nearby regions of the human genome by utilizing PacBio long-read WGS. Thirdly, our sequencing results revealed the occurrence of a wide variety of genome-wide structural variations (SVs) in SiHa and HeLa cells. Additionally, our analysis further revealed a potential correlation between changes in gene expression levels and SVs on chromosome 13 in the genome of SiHa cells.

CONCLUSIONS

Using PacBio long-read sequencing, we have successfully constructed complex models illustrating HPV integrated genome structures in SiHa and HeLa cells. This accomplishment serves as a compelling demonstration of the valuable capabilities of long-read sequencing in detecting and characterizing HPV genomic integration structures within human cells. Furthermore, these findings offer critical insights into the complex process of HPV16 and HPV18 integration and their potential contribution to the development of cervical cancer.

Prevalence and genotype distribution of HPV6/11/16/18 infections among 180,276 outpatient females from a Women's and Children's Central Hospital, 2015-2021, Chengdu, China

The aims of this study on human papilloma virus (HPV) 6/11/16/18 infection among females in Chengdu were to provide more targeted strategies for the prevention and treatment of cervical cancer and genital warts. In this study, the infection status of 20 genotypes was analysed by gene chip technology. The prevalence rates of HPV-6, -11, -16, and -18 infection among 180,276 cases were 0.94%, 0.57%, 3.22%, and 1.28%, respectively. The prevalence of HPV 6/11/16/18 showed a bimodal U-shaped curve with age; the first and second peak occurred among females < 20 and ≥ 60 years old, respectively. As the multiplicity of infections involving HPV6/11/16/18 increases, the infection rate decreases. The ratios of HPV16 single infection showed a yearly increase. The top five genotypes with HPV-16, -18, -6, and -11 in coinfection were HPV52/58/53/51/33, HPV 52/16/53/58/51, HPV52/16/58/51/53 and HPV16/52/58/59/18, respectively, HPV16/18/6/11 were mainly coinfected with HR-HPV. In sum, among the five vaccines available, nonavalent vaccine is more suitable for Chengdu females. For young females prioritizing vaccination is essential in the current context, while HPV screening remains an effective approach for older females. Additionally, in patients with genital warts, it is necessary to assess the presence of high-risk HPV infection and manage it appropriately in patients with genital warts.

The analysis of HPV integration sites based on nanopore sequencing and the profiling changes along the course of photodynamic therapy

OBJECTIVE

To detect the HPV genotype and integration sites in patients with high-risk HPV infection at different stages of photodynamic therapy using nanopore technology and to evaluate the treatment effect.

METHODS

Four patients with HPV infection were selected and subjected to photodynamic therapy, and cervical exfoliated cell was sampled at before treatment, after three courses of treatment and six courses of treatment, their viral abundance and insertion sites were analyzed by nanopore technology, and pathological examinations were performed before and after treatment. In this study, we developed a novel assay that combined viral sequence enrichment and Nanopore sequencing for identification of HPV genotype and integration sites at once. The assay has obvious advantages over qPCR or NGS-based methods, as it has better sensitivity after viral sequences enrichment and can generate long-reads (kb to Mb) for better detection rate of structure variations, moreover, fast turn-around time for real-time viral sequencing and analysis.

RESULTS

The pathological grade was reduced in all four patients after photodynamic therapy. Virus has been cleared in two cases after treatment, the virus amount reduced after treatment but not completely cleared in one case, and two type viruses were cleared and one type virus persisted after treatment in the last patient with multiple infection. Viral abundance and the number of integration sites were positively correlated. Gene enrichment analysis showed complete viral clearance in 1 patient and 3 patients required follow-up.

CONCLUSION

Nanopore sequencing can effectively monitor the abundance of HPV viruses and integration sites to show the presence status of viruses, and combined with the results of gene enrichment analysis, the treatment effect can be dynamically assessed.

Human Papillomavirus 16 DNA Methylation Patterns and Investigation of Integration Status in Head and Neck Cancer Cases

Persistent high-risk human papillomavirus (HPV) infection is a pivotal factor in the progression of cervical cancer. In recent years, an increasing interest has emerged in comprehending the influence of HPV on head and neck squamous cell carcinoma (HNSCC). Notably, it is well established that HPV-associated HNSCC show cases with distinct molecular and clinical attributes compared to HPV-negative cases. The present study delves into the epigenetic landscape of HPV16, specifically its L1 gene and untranslated region (UTR), through pyrosequencing, while the HPV16 DNA physical status was evaluated using E2/E6 ratio analysis in HPV16-positive HNSCC FFPE biopsies. Our findings reveal substantial methylation across six sites within the HPV16 L1 gene and seven sites in the UTR. Specifically, methylation percentages of two L1 CpG sites (7136, 7145) exhibit significant associations with tumor histological grade (p < 0.01), while proving concurrent methylation across multiple sites. The HPV16 DNA physical status was not correlated with the methylation of viral genome or tumor characteristics. This is the first study that examines epigenetic modifications and the HPV16 DNA physical status in Greek HNSCC patients. Our findings suggest an orchestrated epigenetic modulation among specific sites, impacting viral gene expression and intricate virus-host interactions.

Multi-omics data reveals novel impacts of human papillomavirus integration on the epigenomic and transcriptomic signatures of cervical tumorigenesis

Integration of human papilloma virus (HPV) DNA into the human genome may progressively contribute to cervical carcinogenesis. To explore how HPV integration affects gene expression by altering DNA methylation during carcinogenesis, we analyzed a multiomics dataset for cervical cancer. We obtained multiomics data by HPV-capture sequencing, RNA sequencing, and Whole Genome Bisulfite Sequencing from 50 patients with cervical cancer. We detected 985 and 485 HPV-integration sites in matched tumor and adjacent paratumor tissues. Of these, LINC00486 (n = 19), LINC02425 (n = 11), LLPH (n = 11), PROS1 (n = 5), KLF5 (n = 4), LINC00392 (n = 3), MIR205HG (n = 3) and NRG1 (n = 3) were identified as high-frequency HPV-integrated genes, including five novel recurrent genes. Patients at clinical stage II had the highest number of HPV integrations. E6 and E7 genes of HPV16 but not HPV18 showed significantly fewer breakpoints than random distribution. HPV integrations occurring in exons were associated with altered gene expression in tumor tissues but not in paratumor tissues. A list of HPV-integrated genes regulated at transcriptomic or epigenetic level was reported. We also carefully checked the candidate genes with regulation pattern correlated in both levels. HPV fragments integrated at MIR205HG mainly came from the L1 gene of HPV16. RNA expression of PROS1 was downregulated when HPV integrated in its upstream region. RNA expression of MIR205HG was elevated when HPV integrated into its enhancer. The promoter methylation levels of PROS1 and MIR205HG were all negatively correlated with their gene expressions. Further experimental validations proved that upregulation of MIR205HG could promote the proliferative and migrative abilities of cervical cancer cells. Our data provides a new atlas for epigenetic and transcriptomic regulations regarding HPV integrations in cervical cancer genome. We demonstrate that HPV integration may affect gene expression by altering methylation levels of MIR205HG and PROS1. Our study provides novel biological and clinical insights into HPV-induced cervical cancer.

Mutation Profile of HPV16 L1 and L2 Genes in Different Geographic Areas

The causal relationship between HPV and cervical cancer in association with the high prevalence of high risk HPV genotypes led to the design of HPV vaccines based on the major capsid L1 protein. In recent years, capsid protein L2 has also become a focal point in the field of vaccine research. The present review focuses on the variability of HPV16 L1 and L2 genes, emphasizing the distribution of specific amino acid changes in the epitopes of capsid proteins. Moreover, a substantial bioinformatics analysis was conducted to describe the worldwide distribution of amino acid substitutions throughout HPV16 L1, L2 proteins. Five amino acid changes (T176N, N181T; EF loop), (T266A; FG loop), (T353P, T389S; HI loop) are frequently observed in the L1 hypervariable surface loops, while two amino acid substitutions (D43E, S122P) are adjacent to L2 specific epitopes. These changes have a high prevalence in certain geographic regions. The present review suggests that the extensive analysis of the amino acid substitutions in the HPV16 L1 immunodominant loops may provide insights concerning the ability of the virus in evading host immune response in certain populations. The genetic variability of the HPV16 L1 and L2 epitopes should be extensively analyzed in a given population.

Targeted Next Generation Sequencing for Human Papillomavirus Genotyping in Cervical Liquid-Based Cytology Samples

Simple Summary

Testing for Human Papillomavirus (HPV) is currently being implemented as part of cervical cancer (CC) screening in several countries. However, infections with all but one of the HPV types classified as possibly carcinogenic cannot be detected by the assays used for CC screening today. The aim of our study was to demonstrate the use of a targeted next generation sequencing (NGS) HPV panel for CC screening—both general practitioner-collected samples and self-samples. We here show that the targeted HPV panel can detect HPV with a sensitivity and specificity similar to commercial HPV assays, one of which is used for CC screening today. However, the targeted HPV panel possess several advantages compared to the screening assays, as it enables specific detection of all relevant HPV types and can identify viral integration, variants in the HPV genome, and dominant HPV types in multi-infected cases.

Abstract

At present, human papillomavirus (HPV) testing is replacing morphology-based cytology as the primary tool for cervical cancer screening in several countries. However, the HPV assays approved for screening lack detection for all but one of the possibly carcinogenic HPV types and do not genotype all included HPV types. This study demonstrates the use of a targeted HPV next generation sequencing (NGS) panel to detect and genotype all 25 carcinogenic, probably carcinogenic, and possibly carcinogenic HPV types as well as the low-risk types HPV6 and HPV11. The panel was validated using a cohort of 93 paired liquid-based cytology samples (general practitioner (GP)-collected cervical samples and cervico-vaginal self-samples (SS)). Overall, the targeted panel had a sensitivity (GP = 97.7%, SS = 92.1%) and specificity (GP = 98.0%, SS = 96.4%) similar to the commercial HPV assays, Cobas^® 4800 HPV DNA test (Roche) and CLART^® HPV4S assay (GENOMICA). Interestingly, of the samples that tested positive with the NGS panel, three (6.4%) of the GP-collected samples and four (9.1%) of the self-samples tested positive exclusively for HPV types only included in the NGS panel. Thus, targeted HPV sequencing has great potential to improve the HPV screening programs since, as shown here, it can identify additional HPV positive cases, cases with HPV integration, variants in the HPV genome, and which HPV type is dominant in multi-infected cases.

Genetic variability of the HPV16 early genes and LCR. Present and future perspectives

Human papillomavirus 16 (HPV16) infection is the aetiologic factor for the development of cervical dysplasia and is regarded as highly carcinogen, because it is implicated in more than 50% of cervical cancer cases, worldwide. The tumourigenic potential of HPV16 has triggered the extensive sequence analysis of viral genome in order to identify nucleotide variations and amino acid substitutions that influence viral oncogenicity and subsequently the initiation and progression of cervical cancer. Nowadays, specific mutations of HPV16 DNA have been associated with an increased risk of high-grade squamous intraepithelial lesions and invasive cervical cancer (ICC) development, including E6: Q14H, H78Y, L83V, Ε7: N29S, S63F, E2: H35Q, P219S, T310K, E5: I65V, whereas highly conserved regions of viral DNA have been extensively characterised. In addition, numerous novel HPV16 mutations are observed among the studied populations from various geographic regions, hence advocating that different HPV16 strains seem to emerge with different tumourigenic capacities. The present review focuses on the variability of the early genes and the long control region, emphasising on the association of specific mutations with the development of severe dysplasia. Finally, it evaluates whether specific regions of HPV16 DNA are able to serve as valuable biomarkers for cervical cancer risk.

A colorimetric IsoPCR for the rapid and sensitive visual detection of high-risk HPV16 in clinical samples with hydroxynaphthol blue

HPV16 infection is found in more than 50 % of cervical cancer cases worldwide, triggering the development of numerous molecular techniques for viral diagnosis. The present study focuses on the development of a colorimetric IsoPCR for HPV16 DNA detection. The methodology combines the advantages of PCR and LAMP, while the most significant aspect of the new established methodology is the visual detection of amplification products through hydroxynapthol blue dye, thus minimizing the time and labor needed. An experimental cut-off value was tested through reconstitution experiments, while the specificity was evaluated by assessing clinical samples. The analytical sensitivity of the new colorimetric IsoPCR was found to be 0.1 viral DNA copy per reaction, while the specificity was 100 % for the detection of HPV16 DNA. The assay enabled the amplification of viral DNA in cases with viral load lower than 1 copy. In conclusion, the new established colorimetric IsoPCR can be regarded as an attractive molecular tool that facilitates the specific, rapid and highly sensitive visual detection of HPV16 DNA even at the very early stages of viral infection.

Rapid and ultrasensitive detection of circulating human papillomavirus E7 cell-free DNA as a cervical cancer biomarker

Circulating cell-free DNA (cfDNA) has attracted attention as a non-invasive biomarker for diagnosing and monitoring various cancers. Given that human papillomavirus (HPV) DNA integration and overexpression of E6/E7 oncogenes are pivotal events for carcinogenesis, we sought to determine if HPV E7 cfDNA could serve as a specific biomarker for cervical cancer detection. We applied droplet digital PCR (ddPCR) to quantify HPV16/18 E7 cfDNA from the serum of patients with cervical cancer, cervical intraepithelial neoplasia, and controls. HPV16/18 E7 cfDNA was highly specific for cervical cancer, displaying 30.77% sensitivity, 100% specificity, and an area under the curve of 0.65. Furthermore, we developed a sensitive isothermal detection of HPV16/18 E7 and the PIK3CA WT reference gene based on recombinase polymerase amplification combined with a lateral flow strip (RPA-LF). The assay took less than 30 min and the detection limit was 5-10 copies. RPA-LF exhibited 100% sensitivity and 88.24% specificity towards HPV16/18 E7 cfDNA in clinical samples. The agreement between RPA-LF and ddPCR was 83.33% (κ = 0.67) for HPV16 E7 and 100% (κ = 1.0) for HPV18 E7, indicating a good correlation between both tests. Therefore, we conclude that HPV E7 cfDNA represents a potential tumor marker with excellent specificity and moderate sensitivity for minimally invasive cervical cancer monitoring. Moreover, the RPA-LF assay provides an affordable, rapid, and ultrasensitive tool for detecting HPV cfDNA in resource-limited settings.

RNA-based high-risk HPV genotyping and identification of high-risk HPV transcriptional activity in cervical tissues

Nearly all cervical cancers are initiated by a persistent infection with one of the high-risk human papillomaviruses (high-risk HPV). High-risk HPV DNA testing is highly sensitive but cannot distinguish between active, productive infections and dormant infections or merely deposited virus. A solution for this shortcoming may be the detection of transcriptional activity of viral oncogenes instead of mere presence of high-risk HPVs. In this study, fresh-frozen cervical tissues (n = 22) were subjected to high-risk HPV DNA detection using the line probe assay and to targeted RNA next-generation sequencing using single-molecule molecular inversion probes. Targeted RNA sequencing was applied for (1) RNA-based genotyping of high-risk HPV, giving information on specific HPV-subtype (2) discrimination of E2, E6, and E7 transcripts and (3) discovery of possible non-HPV cancer biomarkers. Data were analyzed using computational biology. Targeted RNA sequencing enabled reliable genotyping of high-risk HPV subtypes and allowed quantitative detection of E2, E6, and E7 viral gene expression, thereby discriminating cervical lesions from normal cervical tissues. Moreover, targeted RNA sequencing identified possible cervical cancer biomarkers other than high-risk HPV. Interestingly, targeted RNA sequencing also provided high-quality transcription profiles from cervical scrape samples, even after 1 week of dry storage or storage in Preservcyt fixative. This proof of concept study shows that targeted RNA sequencing can be used for high-risk HPV genotyping and simultaneous detection of high-risk HPV gene activity. Future studies are warranted to investigate the potential of targeted RNA sequencing for risk assessment for the development of cervical lesions, based on molecular analysis of cervical scrapes.

Human Papillomavirus Type 16 E1 Mutations Associated with Cervical Cancer in a Han Chinese Population

Human papillomavirus type 16 (HPV16) is a high-risk HPV type and a potent carcinogen. HPV E1 is one of the most highly conserved proteins and it plays a central role in initiating HPV DNA replication. In current study, we enrolled 161 HPV16-positive cervical cancer patients (case group) and 171 HPV16-positive asymptomatic individuals (control group) in a study to analyse the association between HPV16 E1 genetic mutations and cervical cancer. The samples of case group were cervical cancer tissues and the samples of control group were cervical exfoliated cells. Three variants (A4, A1-A3 and D3) were found in the case group, 68.3% of the HPV16 E1 sequences belonged to the A4 (As) sub-lineage, 29.2% belonged to the A1-A3 (EUR) sub-lineage, and 2.5% belonged to the D3 (AA1) sub-lineage. Two variants (A4 and A1-A3) occurred in the control group. The A4 (As) sub-lineage was predominant in this group as well (66.1%), followed by the A1-A3 (EUR) sub-lineage (33.9%), but the D3 (AA1) sub-lineage was not found in the control group. The distribution of the HPV16 variants between the case and control groups was significantly different (P<0.05). When the distribution of the HPV16 E1 gene mutations was compared, the distribution of twenty-seven mutations was significantly different between the case and control groups (P<0.05), and twenty-two mutations occurred only in the D3 (AA1) sub-lineage, two were found only in the A4 (As) sub-lineage, one was found in the A1-A3 (EUR) sub-lineage, two was found in both the A4 (As) and A1-A3 (EUR) sub-lineages. In the sub-lineage analysis, the differences in the T933A (A23A), T1014G (D50E) and G2160A (R432R) mutations were statistically significant between the case and control groups for the A4 (As) sub-lineage (P<0.05), and the differences in the T2232C (F456F), G2337A (M491I) and A2547G (P561P) mutations were statistically significant between the case and control groups for the A1-A3 (EUR) sub-lineage (P<0.05). In the current study, we describe specific mutations in the HPV16 E1 gene associated with cervical cancer, and our study will provide a good reference for further functional studies of the relationship between cervical cancer carcinogenesis and HPV genes.

Methylation of p16 ink4a promoter is independent of human papillomavirus DNA physical state: a comparison between cervical pre-neoplastic and neoplastic samples

BACKGROUND Epigenetic modifications in host cells, like p16 ink4a methylation, have been considered as putative complementary mechanisms for cancer development. Because only a small proportion of infected women develop cervical cancer, other factors might be involved in carcinogenesis, either independently or in association with high-risk human papillomavirus (HR-HPV) infections, including epigenetic factors. OBJECTIVES We hypothesised that p16 ink4a methylation might have a role in cancer development driven by HPV16, mainly in the presence of intact E1/E2 genes. Thus, our objectives were to assess the status of p16 ink4a methylation and the HPV16 E1/E2 integrity in samples in different stages of cervical diseases. METHODS Presence of HPV16 was determined by E6 type-specific polymerase chain reaction (PCR). Methylation status of the p16 ink4a promoter was assessed by methylation-specific PCR in 87 cervical specimens comprising 29 low-grade (LSIL), 41 high-grade (HSIL) lesions, and 17 cervical cancers (CC). Characterisation of E1 and E2 disruption (as an indirect indicator of the presence of episomal viral DNA) was performed by PCR amplifications. FINDINGS We observed a significantly increased trend (nptrend = 0.0320) in the proportion of methylated p16 ink4a in cervical samples during cancer development. Concomitant E1 and E2 disruptions were the most frequent pattern found in all groups: CC (76%), HSIL (54%), and LSIL (73%). No statistically significant differences between p16 ink4a methylation and E1/E2 integrity, in histological groups, was observed. MAIN CONCLUSIONS There was an increase in methylation of the p16 ink4a promoter from pre-neoplastic lesions to cancer. Additionally, a high frequency of E1/E2 disruptions in LSIL/HSIL suggested that viral DNA integration was an early event in cervical disease. Moreover, the methylation status was apparently independent of HPV16 integrity.

HPV16 viral characteristics in primary, recurrent and metastatic vulvar carcinoma

Vulvar carcinoma is the fourth most common gynecological malignancy. Two separate carcinogenic pathways are suggested, where one is associated with the human papillomavirus (HPV) and HPV16 the most common genotype. The aim of this study was to evaluate HPV-markers in a set of primary tumors, metastases and recurrent lesions of vulvar squamous cell carcinomas (VSCC). Ten HPV16-positive VSCC with metastatic regional lymph nodes, distant lymphoid/hematogenous metastases or local recurrent lesions were investigated for HPV genotype, HPV16 variant, HPV16 viral load, HPV16 integration and HPV16 E2BS3 and 4 methylation. In all 10 analyzed case series, the same HPV genotype (HPV16), HPV16 variant and level of viral load were detected in all lesions within a patient case. Primary tumors with a high E2/E6 ratio were found to have fewer vulvar recurrences and/or metastases after diagnosis and treatment. Also, a significantly lower viral load was evident in regional lymph nodes compared to primary tumors. The data presented strengthens the evidence for a clonal HPV-induced pathway for vulvar carcinoma.

HGT-ID: an efficient and sensitive workflow to detect human-viral insertion sites using next-generation sequencing data

Background

Transfer of genetic material from microbes or viruses into the host genome is known as horizontal gene transfer (HGT). The integration of viruses into the human genome is associated with multiple cancers, and these can now be detected using next-generation sequencing methods such as whole genome sequencing and RNA-sequencing.

Results

We designed a novel computational workflow, HGT-ID, to identify the integration of viruses into the human genome using the sequencing data. The HGT-ID workflow primarily follows a four-step procedure: i) pre-processing of unaligned reads, ii) virus detection using subtraction approach, iii) identification of virus integration site using discordant and soft-clipped reads and iv) HGT candidates prioritization through a scoring function. Annotation and visualization of the events, as well as primer design for experimental validation, are also provided in the final report. We evaluated the tool performance with the well-understood cervical cancer samples. The HGT-ID workflow accurately detected known human papillomavirus (HPV) integration sites with high sensitivity and specificity compared to previous HGT methods. We applied HGT-ID to The Cancer Genome Atlas (TCGA) whole-genome sequencing data (WGS) from liver tumor-normal pairs. Multiple hepatitis B virus (HBV) integration sites were identified in TCGA liver samples and confirmed by HGT-ID using the RNA-Seq data from the matched liver pairs. This shows the applicability of the method in both the data types and cross-validation of the HGT events in liver samples. We also processed 220 breast tumor WGS data through the workflow; however, there were no HGT events detected in those samples.

Conclusions

HGT-ID is a novel computational workflow to detect the integration of viruses in the human genome using the sequencing data. It is fast and accurate with functions such as prioritization, annotation, visualization and primer design for future validation of HGTs. The HGT-ID workflow is released under the MIT License and available at http://kalarikrlab.org/Software/HGT-ID.html.

HPV DNA methylation at the early promoter and E1/E2 integrity: A comparison between HPV16, HPV18 and HPV45 in cervical cancer

Objectives

To compare and describe type-specific characteristics of HPV16, HPV18 and HPV45 in cervical cancer with respect to 3′LCR methylation and disruption of E1/E2.

Methods

The methylation level of 137 cervical cancer samples (70 with HPV16, 37 with HPV18, and 30 with HPV45) of Brazilian patients was analyzed by pyrosequencing. PCR amplifications were performed to characterize E1 and E2 disruption as an episomal surrogate.

Results

The 3′LCR of HPV16 showed a higher methylation at all CpG sites (7%, 9%, 11%, 10% and 10%) than homologous HPV18 regions (4%, 5%. 6%, 9% and 5%) and HPV45 regions (7%, 7% and 5%). Presence of intact E1/E2 was associated with higher HPV16 and HPV18 methylation levels at all CpG sites (p < 0.05). Disruption of E1/E2 was more frequently found in HPV45 (97%) and HPV18 (84%) than in HPV16 DNA (30%). HPV16 disruption was more frequently found in E1 (48%) unlike HPV18, where it was found in E2 (61%). Concomitant disruption of E1/E2 was most frequent in HPV45 (72%).

Conclusions

The findings showed a higher methylation associated with intact E1/E2 for HPV16 and HPV18. The closely phylogenetic related HPV18 and HPV45 share a similar methylation level and the frequency of viral genome disruption.

Up-regulation of miR-20a by HPV16 E6 exerts growth-promoting effects by targeting PDCD6 in cervical carcinoma cells

OBJECTIVE

MicroRNAs (miRNAs/miRs) have been reported to participate in progression of multiple tumors including cervical cancer. High-risk human papillomavirus (HPV) type 16 (HPV16) is the most common and lethal HPV type, leading to exceeding 50% of cervical cancer cases. However, the relationship between miRNA and HPV-induced cervical carcinogenesis remains elusive.

RESULTS

Here, HPV16 E6 positively regulated miR-20a expression. Overexpression of miR-20a showed growth-promoting effects on C33A cells (HPV16-negative), and knockdown of miR-20a showed growth-inhibitory effects on CaSki cells (HPV16-positive). In addition, PDCD6 was identified as a target gene of miR-20a. Overexpression of PDCD6 exerted growth-inhibitory effects (opposite to miR-20a overexpression), which could be reversed by miR-20a overexpression. More importantly, activation of AKT and p38 was observed in C33A cells overexpressing miR-20a, and the growth-promoting action of miR-20a could be abated by p38 inhibition.

CONCLUSION

Up-regulation of miR-20a by HPV16 E6 exerted growth-promoting effects by targeting PDCD6 in cervical carcinoma cells. This study demonstrated miR-20a might be a potential therapeutic target in HPV16 E6 infection type of cervical cancer.

Association of an intact E2 gene with higher HPV viral load, higher viral oncogene expression, and improved clinical outcome in HPV16 positive head and neck squamous cell carcinoma

To assess the relationship of E2 gene disruption with viral gene expression and clinical outcome in human papillomavirus (HPV) positive head and neck squamous cell carcinoma, we evaluated 31 oropharyngeal and 17 non-oropharyngeal HPV16 positive carcinomas using two PCR-based methods to test for disruption of E2, followed by Sanger sequencing. Expression of HPV16 E6, E7 and E2 transcripts, along with cellular ARF and INK4A, were also assessed by RT-qPCR. Associations between E2 disruption, E2/E6/E7 expression, and clinical outcome were evaluated by Kaplan-Meier analysis for loco-regional recurrence and disease-specific survival. The majority (n = 21, 68%) of HPV16 positive oropharyngeal carcinomas had an intact E2 gene, whereas the majority of HPV16 positive non-oropharyngeal carcinomas (n = 10, 59%) had a disrupted E2 gene. Three of the oropharyngeal tumors and two of the non-oropharyngeal tumors had deletions within E2. Detection of an intact E2 gene was associated with a higher DNA viral load and increased E2/E6/E7, ARF and INK4A expression in oropharyngeal tumors. Oropharyngeal carcinomas with an intact E2 had a lower risk of loco-regional recurrence (log-rank p = 0.04) and improved disease-specific survival (p = 0.03) compared to tumors with disrupted E2. In addition, high E7 expression was associated with lower risk of loco-regional recurrence (p = 0.004) as was high E6 expression (p = 0.006). In summary, an intact E2 gene is more common in HPV16 positive oropharyngeal than non-oropharyngeal carcinomas; the presence of an intact E2 gene is associated with higher HPV viral load, higher viral oncogene expression, and improved clinical outcome compared to patients with a disrupted E2 gene in oropharyngeal cancer.

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.

HPV16 integration probably contributes to cervical oncogenesis through interrupting tumor suppressor genes and inducing chromosome instability

Background

The integration of human papilloma virus (HPV) into host genome is one of the critical steps that lead to the progression of precancerous lesion into cancer. However, the mechanisms and consequences of such integration events are poorly understood. This study aims to explore those questions by studying high risk HPV16 integration in women with cervical intraepithelial neoplasia (CIN) and cervical squamous cell carcinoma (SCC).

Methods

Specifically, HPV integration status of 13 HPV16-infected patients were investigated by ligation-mediated PCR (DIPS-PCR) followed by DNA sequencing.

Results

In total, 8 HPV16 integration sites were identified inside or around genes associated with cancer development. In particular, the well-studied tumor suppressor genes SCAI was found to be integrated by HPV16, which would likely disrupt its expression and therefore facilitate the migration of tumor. On top of that, we observed several cases of chromosome translocation events coincide with HPV integration, which suggests the existence of chromosome instability. Additionally, short overlapping sequences were observed between viral derived and host derived fragments in viral-cellular junctions, indicating that integration was mediated by micro homology-mediated DNA repair pathway.

Conclusions

Overall, our study suggests a model in which HPV16 might contribute to oncogenesis not only by disrupting tumor suppressor genes, but also by inducing chromosome instability.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-016-0454-4) contains supplementary material, which is available to authorized users.

Analysis of E2 gene integrity in HPV16 and HPV58 viruses isolated from women with cervical pathology

Integration of human papillomavirus (HPV) DNA into human cells accompanied by the disruption of the viral genome has been described as a prerequisite for cancer development. This study aimed to investigate E2 gene integrity of HPV16 and HPV58 viruses isolated from infected women with cervical lesions. Forty-two HPV16- and 31 HPV58-positive samples were analysed. E2 integrity was assumed when all fragments covering the E2 gene were amplified with specific polymerase chain reaction primers. Overall, in 59% of the samples, at least one fragment was not amplified in HPV16- (57%) and HPV58-positive samples (61%). Samples from high-grade squamous intraepithelial lesions had the highest frequency of E2 gene disruptions (73%), followed by samples from low-grade squamous intraepithelial lesions (63%) and, finally, samples from invasive cervical cancer (35%). Association between the integrity status of the E2 gene, and lesion grade was assessed by the chi-squared test applied to the combined set of viruses (p = 0.6555) or to populations of the same virus type (HPV58, p = 0.3101; HPV16, p = 0.3024). In conclusion, in this study, no association was found between the presence of E2 gene disruptions and the grade of cervical lesions caused by HPV16 and HPV58.

E1 detection as prognosticator in human papillomavirus-positive head and neck cancers

PURPOSE

HPV-related locally advanced head and neck cancers (LA-HNCs) show a good prognosis. This study aimed to investigate the HPV prevalence in LA-HNCs and compare the prognostic value of E1, E6 and L1 genomic viral fragments and p16, individually and in combination, in order to find the best prognosticator in terms of overall survival (OS) and progression-free survival (PFS).

PATIENTS AND METHODS

HPV16 was searched in 255 LA-HNC formalin-fixed paraffin-embedded tumor tissues, 89 oropharyngeal cancers (OPCs), and 166 non-OPCs by DNA-PCR with 3 primer pairs. p16 was analyzed by immunohistochemistry in 235 patients.

RESULTS

The prevalence of positive samples decreased constantly from E6 to L1 and E1 in both OPCs and non-OPCs. Each LA-HNC patient highlighted variable positivity for each fragment. OPCs showed a higher prevalence of positive samples compared to non-OPCs.Positive coexistence of all the fragments was more common in OPCs (31.5%) than non-OPCs (4.2%), and E1 detection was always associated with E6 and L1. E1-positive OPCs showed improved OS (p = 0.012) and PFS (p = 0.036), while L1- or E6-positive ones did not. p16-positive patients were more prevalent in the OPC (29.8%) than the non-OPC group (7.3%) (p&lt;0.0001) and its prognostic value was not superior to that of E1. However, the multivariate Cox analysis which included E1, L1, E6 status and p16 expression did not show a significant p value.

CONCLUSIONS

Though HPV16 positivity measured by DNA-PCR was higher for L1 and E6, they performed weakly as prognosticators; E1 might become a strong prognostic marker for OS and PFS in OPCs.

Identification of rearranged sequences of HPV16 DNA in precancerous and cervical cancer cases

Integration of HPV16 DNA into the host chromosome is considered to be a crucial step towards genomic instability and cervical cancer development. Aim of the present study was to investigate the presence of HPV16 rearranged intra-viral sequences in HPV16-positive normal, precancerous and cervical cancer samples using the method of Restriction Site-PCR (RS-PCR). Sequence analysis of HPV16 integrants revealed for the first time in clinical samples two distinct rearranged intra-viral sequences, concerning the conjunction of E2 and L1 genes and the conjunction of E1 and L1 genes with inverted orientation. Furthermore mapping analysis of the E1 and E2 genes in cervical samples with rearranged intra-viral sequences of HPV16 genome was conducted in order to determine the integrity of viral genes. The identification of intra-viral rearrangements provides valuable information regarding the HPV16 DNA integration, and may be a significant biomarker for the presence of chromosomal instability and DNA damages in clinical samples.