Insights into the mechanisms and structure of breakage-fusion-bridge cycles in cervical cancer using long-read sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cancer cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and our study applies long-read sequencing to this important chromosomal rearrangement type. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes, and only one BFB breakpoint showed chromothripsis. Five cell lines have a chromosomal region 11q BFB event, with YAP1-BIRC3-BIRC2 amplification. Indeed, YAP1 amplification is associated with a 10-year-earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that individuals with cervical cancer and YAP1-BIRC3-BIRC2 amplification, especially those of African ancestry, might benefit from targeted therapy. In summary, we uncovered valuable insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.

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.

Testing for viral DNA integration among HPV-positive women to detect cervical precancer: An observational cohort study.

OBJECTIVE: Human papillomavirus (HPV) integration is a crucial genetic step in cervical carcinogenesis. This study aimed to evaluate the performance of an HPV integration test for the triage of HPV-positive women. DESIGN: An observational cohort study. SETTING: A cervical cancer screening programme in China. POPULATION: 1393 HPV-positive women aged 25-65 years undergoing routine cervical cancer screening and HPV integration testing with 1-year follow-up. METHODS: The sensitivity, specificity, positive predictive value and negative predictive value between HPV integration and cytology were compared. MAIN OUTCOME MEASURES: Cervical intraepithelial neoplasia grade 3 or more severe (CIN3+). RESULTS: Among 1393 HPV-positive patients, 138 (9.9% [8.3-11.5%]) were HPV integration test positive compared with 537 who had abnormal cervical cytology (38.5% [36.0-41.1%]). Compared with cytology, HPV integration exhibited higher specificity (94.5% [93.3-95.8%] versus 63.8% [61.2-66.4%]) and equivalent sensitivity (70.5% [61.4-79.7%] versus 70.5% [61.4-79.7%]) for detection of CIN3+. HPV integration-negative women accounted for 90.1% (1255/1393) of the total population and had a low immediate CIN3+ risk (2.2%). At 1-year follow-up, the progression rate in the HPV integration-positive women was higher than in the HPV integration-negative women (12.0% versus 2.1%, odds ratio 5.6, 95% CI, 2.6-11.9). In 10 conservatively managed integration-negative CIN2 patients, all showed spontaneous regression and seven showed HPV clearance after 1-year follow-up. CONCLUSION: The HPV integration test may be a precise risk stratification tool for HPV-positive women and could avoid excessive use of invasive biopsies.

Preliminary study of HPV integration status on the occurrence and development of vaginal intraepithelial neoplasia.

AIM: The aim of this study was to investigate how the integration status of HPV in the vaginal epithelium affects the development of vaginal intraepithelial neoplasia (VaIN). METHODS: Twenty-four vaginal tissues were collected before applying high-throughput viral integration detection (HIVID), medical records of them were documented, including age, thin-prep cytologic test (TCT) and HPV test results, colposcopic biopsy pathology, and other clinical data, such as history of total hysterectomy for cervical lesions, whether they were infected with HPV16/18 with a follow-up span of 2 years. We summarized the distribution of HPV integration on the host chromosome and HPV type, as well as the hotspot integration gene and its role in the development of VaIN. RESULTS: In this study, 24 cases suffered from VaIN were involved. HPV integration was detected in 11 cases; furthermore, we discovered HPV 16 and 73, chromosome 1 and 2 possessed most HPV integration sites while EMBP1, CLO5A1, EHF, ELF5 as dominate hot spots. Taken clinical outcome into account, we found a significant difference between HPV integration occurrence and VaIN (p = 0.011). CONCLUSION: (1) This study found a statistical difference between HPV integration and the occurrence of VaIN; (2) HPV integration may provide a new clinical predictor for VaIN and facilitate risk assessment and stratified management of high-risk patients.

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.

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.

HPV-CCDC106 integration promotes cervical cancer progression by facilitating the high expression of CCDC106 after HPV E6 splicing.

Human papillomavirus (HPV) integration and high expression of HPV oncogenes (E6 and E7) are important mechanisms for HPV carcinogenesis in cervical cancer. However, the relationship between HPV integration and HPV E6 spliced transcripts, as well as the underlying mechanisms of HPV integration in carcinogenesis after HPV E6 splicing remains unclear. We analyzed HPV-coiled-coil domain containing 106 (CCDC106) integration samples to characterize the roles of HPV integration, E6 spliceosome I (E6*I), and high CCDC106 expression in cervical carcinogenesis. We found that E6 was alternatively spliced into the E6*I transcript in HPV-CCDC016 integration samples with low p53 expression, in contrast to the role of E6*I in preventing p53 degradation in cervical cancer cells. In addition, CCDC106 was highly expressed after HPV-CCDC106 integration, and interacted with p53, resulting in p53 degradation and cervical cancer cell progression in vitro and in vivo. Importantly, when E6*I was highly expressed in cervical cancer cells, overexpression of CCDC106 independently degraded p53 and promoted cervical cancer cell progression. In this study, we explored the underlying mechanisms of HPV-CCDC106 integration in HPV carcinogenesis after HPV E6 splicing, which should provide insight into host genome dysregulation in cervical carcinogenesis.

Integration and viral oncogene expression of human papillomavirus type 16 in oropharyngeal squamous cell carcinoma and gastric cancer.

Persistent high-risk human papillomavirus (HR-HPV) infections cause cervical cancer and a fraction of head and neck cancer. To investigate whether HR-HPV infection might be also involved in the development of gastric cancer (GC), we developed a platform utilizing a rolling circle amplification (RCA)-based nested L1 polymerase chain reaction with Sanger sequencing to genotype the HPV DNA in cancer tissues of 361 GC and 89 oropharyngeal squamous cell carcinomas (OPSCC). HPV transcriptional activity was determined by E6/E7 mRNA expression and a 3' rapid amplification of cDNA ends was performed to identify HPV integration and expression of virus-host fusion transcripts. Ten of 361 GC, 2 of 89 OPSCC, and 1 of 22 normal adjacent tissues were HPV L1 DNA-positive. Five of the 10 HPV-positive GC were genotyped as HPV16 by sequencing and 1 of 2 GC with RCA/nested HPV16 E6/E7 DNA detection exhibited HPV16 E6/E7 mRNA. Two OPSCC displayed HPV16 L1 DNA and E6/E7 mRNA, of which 1 OPSCC tissue showed virus-host RNA fusion transcripts from an intron region of KIAA0825 gene. Together, our data reveal viral oncogene expression and/or integration in GC and OPSCC and a possible etiology role of HPV infections in gastric carcinogenesis.

Human papillomavirus (HPV) integration signature in cervical lesions: identification of MACROD2 gene as HPV hot spot integration site.

BACKGROUND: High-risk HPV is clearly associated with cervical cancer. Integration of HPV DNA into the host genome is considered a key event in driving cervical carcinogenesis. However, the mechanism on how HR-HPV integration influences the host genome structure has remained enigmatic. METHODS: In our study, 25 DNA samples including 11 from fresh-frozen cervical carcinomas and 14 from fresh-frozen high-grade squamous intraepithelial lesion (HSILs) were detected using the method of HPV capture combined with next generation sequencing. RESULTS: We calculated the frequency in each viral gene or region and found that breakpoints were prone to occur in L1 and L2 instead of E2 in the cervical cancer (P = 0.0004 and P = 5.15 × 10-40) and HSIL group (P = 2.1 × 10-32 and P = 7.06 × 10-13). The results revealed that HPV16 showed a strong tendency toward intronic region (P = 5.02 × 10-64) but a subtle tendency toward intergenic region (P = 0.04). The most frequent integration site was in the MACROD2 gene (introns 2, 4, 5, 6, 8 and 9), which in MACROD2 functional domain. CONCLUSION: Our results revealed that MACROD2 is HPV hot spot integration site in cervical lesions, and its deficiency alter DNA repair and sensitivity to DNA damage thought impaired PARP1 activity resulting in chromosome instability.

Human Papillomavirus Prevalence and Integration Status in Tissue Samples of Bladder Cancer in the Chinese Population.

BACKGROUND: Human papillomavirus (HPV) infection is associated with multiple types of cancer, but the evidence has not yet been fully elucidated in bladder cancer. METHODS: Frozen tissue samples collected from 146 patients aged 32 to 89 years with bladder cancer pathological diagnosis between 2015 and 2019 were analyzed. HPV genotyping and integration status determination were performed by capture-based next generation sequencing. Statistical analysis of HPV type distributions was performed according to stage, grade, sex, and age group of patients. RESULTS: Mean (SD) age of the 146 patients was 66.64 ±â€…10.06 years and 83.56% were men. Overall HPV infection rate was 28.77% (37.50% in women and 27.05% in men), with 11.90% HPV integration events. Among them, 17.12% single and 11.65% coinfections were observed. HPV18 (24.66%) was the most prevalent genotype, followed by HPV33, 16, and 39. All HPV were European lineage (A). HPV16 was more prevalent in women (P = .04). CONCLUSIONS: HPV infection may contribute to the etiology both in men and women with bladder cancer. HPV18, followed by HPV33, 16, and 39 genotypes, potentially represent the predominant oncogenic risk types for bladder carcinogenesis.

Analysis of HPV Integrations in Mexican Pre-Tumoral Cervical Lesions Reveal Centromere-Enriched Breakpoints and Abundant Unspecific HPV Regions.

Human papillomavirus (HPV) DNA integration is a crucial event in cervical carcinogenesis. However, scarce studies have focused on studying HPV integration (HPVint) in early-stage cervical lesions. Using HPV capture followed by sequencing, we investigated HPVint in pre-tumor cervical lesions. Employing a novel pipeline, we analyzed reads containing direct evidence of the integration breakpoint. We observed multiple HPV infections in most of the samples (92%) with a median integration rate of 0.06% relative to HPV mapped reads corresponding to two or more sequence breakages. Unlike cancer studies, most integrations events were unique (supported by one read), consistent with the lack of clonal selection. Congruent to other studies, we found that breakpoints could occur, practically, in any part of the viral genome. We noted that L1 had a higher frequency of rupture integration (25%). Based on host genome integration frequencies, we found previously reported integration sites in cancer for genes like FHIT, CSMD1, and LRP1B and putatively many new ones such as those exemplified in CSMD3, ROBO2, and SETD3. Similar host integrations regions and genes were observed in diverse HPV types within many genes and even equivalent integration positions in different samples and HPV types. Interestingly, we noted an enrichment of integrations in most centromeres, suggesting a possible mechanism where HPV exploits this structural machinery to facilitate integration. Supported by previous findings, overall, our analysis provides novel information and insights about HPVint.

Whole genome sequencing reveals complexity in both HPV sequences present and HPV integrations in HPV-positive oropharyngeal squamous cell carcinomas.

BACKGROUND: High risk human papillomaviruses (HPV) plays important roles in the development of cervical cancer, a number of other anogenital cancer and they are increasingly found in oropharyngeal squamous cell carcinoma (OPSCC), however there has not been comprehensive analysis about the role how these viruses play in the development of OPSCC. METHODS: To characterize the physical status of HPV within OPSCC and to determine the effect this has throughout the host genome, we have performed 30-40X whole genome sequencing (WGS) on the BGI sequencing platform on 34 OPSCCs: 28 of which were HPV positive. We then examined the sequencing data to characterize the HPV copy number and HPV physical status to determine what effect they have on both HPV and human genome structural changes. RESULTS: WGS determined the HPV copy number across the viral genome. HPV copy number ranged from 1 copy to as high as 150 copies in each individual OPSCC. Independent of HPV copy number, most tumors had either a small or a very large deletion in the viral genome. We discovered that these deletions were the result of either HPV integration into the human genome or HPV-HPV sequence junctions. WGS revealed that ~ 70% of these tumors had HPV integrations within the human genome and HPV integration occurred independent of HPV copy number. Individual HPV integrations were found to be highly disruptive resulting in structural variations and copy number changes at or around the integration sites. CONCLUSIONS: WGS reveals that there is a great complexity in both HPV sequences present and the HPV integrations events in HPV positive OPSCCs tumors. Thus HPV may be playing different roles in the development of different OPSCCs and this further challenge the HPV-driven carcinogenesis model first proposed for cervical cancer.

HPV integration and cervical cancer: a failed evolutionary viral trait.

Countless efforts have been made to eradicate cervical cancer worldwide, including improving disease screening and human papillomavirus (HPV) vaccination programs. Nevertheless, cervical cancer still claims the lives of more than 300 000 women every year. Persistent infections with high-risk HPV genotypes 16 and 18 are the main cause of cancer and may result in HPV integration into the host genome. The central dogma is that HPV integration is an important step in oncogenesis, but in fact, it impedes the virus from replicating and spreading. HPV causing cervical cancer can therefore be perceived as a failed evolutionary viral trait. Here we outline the occurrence and mechanisms of HPV integration and how this process results in oncogenic transformation.

Long-read sequencing of oropharyngeal squamous cell carcinoma tumors reveal diverse patterns of high-risk Human Papillomavirus integration.

Introduction: In North America and in most European countries, Human Papillomavirus (HPV) is responsible for over 70% of oropharyngeal squamous cell carcinomas. The burden of OPSCC, in high-income countries, has been steadily increasing over the past 20 years. As a result, in the USA and in the UK, the burden of HPV-related oropharyngeal squamous cell carcinoma in men has now surpassed that of cervical cancer in women. However, the oncogenic impact of high-risk HPV integration in oropharyngeal squamous cell carcinomas hasn't been extensively studied. The present study aimed to explore the patterns of HPV integration in oropharyngeal squamous cell carcinomas and to assess the feasibility and reliability of long-read sequencing technology in detecting viral integration events in oropharyngeal head and neck cancers. Methods: A cohort of eight HPV-positive OPSCC pre-treatment patient tumors (four males and four females), were selected. All patients received a p16INK4A positive OPSCC diagnosis and were treated at the McGill University Health Centre, a quaternary center in Montreal. A minimum of 20mg of tumor tissue was used for DNA extraction. Extracted DNA was subjected to Nanopore long-read sequencing to detect and analyze for the presence of high-risk HPV sequences. PCR and Sanger sequencing experiments were performed to confirm Nanopore long-read sequencing readings. Results: Nanopore long-read sequencing showed that seven out of eight patient samples displayed either integrated or episomal high-risk HPV sequences. Out of these seven samples, four displayed verifiable integration events upon bioinformatic analysis. Integration confirmation experiments were designed for all four samples using PCR-based methods. Sanger sequencing was also performed. Four distinct HPV integration patterns were identified: concatemer chromosomal integration in a single chromosome, bi-chromosomal concatemer integration, single chromosome complete integration and bi-chromosomal complete integration. HPV concatemer integration also proved more common than full HPV integration events. Conclusion and relevance: Long-read sequencing technologies can be effectively used to assess HPV integration patterns in OPSCC tumors. Clinically, more research should be conducted on the prognostication value of high-risk HPV integration in OPSCC tumors using long-read sequencing technology.

Effects of CpG sites methylation modification of HPV16 integration essential gene on the proliferation of cervical cancer cells.

PURPOSE: The mechanism of methylation of HPV CpG sites in the occurrence and prognosis of cervical carcinogenesis remains unclear. We investigated the effects of demethylation of the CpG sites of E2 and E6, essential genes of HPV16 integration, on cervical cancer cell expression, integration, and proliferation. MATERIALS AND METHODS: HPV16-positive (Caski) cells were treated with different concentrations of the demethylation compound 5-aza-dc (0, 5, 10, 20 µmol/l) in vitro. After the intervention, the methylation statuses of HPV16 E2 and E6 were detected by TBS, the expression levels of E2 and E6 mRNA and protein were detected by real-time PCR and western blot, cell proliferation activity was detected by CCK8, and cell cycle and apoptosis were determined by FCM. GraphPad Prism version 8.4.2 and R version 4.2.3 were used for relevant data analyses. RESULTS: The methylation levels of HPV16 E2 and E6 CpG sites decreased gradually with increasing 5-aza-dc intervention concentrations. With decreasing E2 and E6 methylation rates, E2 expression increased, the E2/E6 ratio increased, E6 expression decreased, and the growth inhibition rate of Caski cells increased. E2 and E6 expression were negatively and positively correlated with their degrees of methylation respectively, while the E2/E6 mRNA to protein ratio was negatively correlated with the methylation degrees of E2 and E6. CONCLUSION: Demethylation can be used as a prospective treatment to affect HPV expression and persistent infection, providing a new theoretical basis for the clinical treatment of viral infections.

Multi-omics characterization of silent and productive HPV integration in cervical cancer.

Cervical cancer (CC) that is caused by high-risk human papillomavirus (HPV) remains a significant public health problem worldwide. HPV integration sites can be silent or actively transcribed, leading to the production of viral-host fusion transcripts. Herein, we demonstrate that only productive HPV integration sites were nonrandomly distributed across both viral and host genomes, suggesting that productive integration sites are under selection and likely to contribute to CC pathophysiology. Furthermore, using large-scale, multi-omics (clinical, genomic, transcriptional, proteomic, phosphoproteomic, and single-cell) data, we demonstrate that tumors with productive HPV integration are associated with higher E6/E7 proteins and enhanced tumor aggressiveness and immunoevasion. Importantly, productive HPV integration increases from carcinoma in situ to advanced disease. This study improves our understanding of the functional consequences of HPV fusion transcripts on the biology and pathophysiology of HPV-driven CCs, suggesting that productive HPV integration should be evaluated as an indicator of high risk for progression to aggressive cancers.

Characteristics of HPV integration in cervical adenocarcinoma and squamous carcinoma.

PURPOSE: HPV integration usually occurs in HPV-related cancer, and is the main cause of cancer. But the carcinogenic mechanism of HPV integration is unclear. The study aims to provide a theoretical basis for understanding the pathogenesis of cervical adenocarcinoma (AC) and cervical squamous carcinoma (SCC). METHODS: We used HPV capture sequencing to obtain HPV integration sites in AC and SCC, and analyzed cytobands, distribution of genetic and genomic elements, identified integration hotspot genes, clinicopathological parameters, breakpoints of HPV16 and performed pathway analysis. Then we conducted immunohistochemical (IHC) assay to preliminarily verify the expression of most frequently integrated genes in AC, STARD3 and ERBB2. RESULTS: The results revealed that the most frequently observed integrated cytoband was 17q12 in AC and 21p11.2 in SCC, respectively. The breakpoints in both AC and SCC were more tended to occur within gene regions, compared to intergenetic regions. Compared to SCC samples, AC samples had a higher prevalence of genomic elements. In AC, HPV integration has no significantly difference with clinicopathological parameters, but in SCC integration correlated with differentiation (P < 0.05). Breakpoints of HPV in SCC located in LCR more frequently compared to AC, which destroyed the activation of promoter p97. Hotspot genes of HPV integration were STARD3 and ERBB2 in AC, and RNA45S rDNA and MIR3648-1 in SCC, respectively. Meanwhile, we preliminarily proved that the expression of STARD3 and ERBB2, the most frequently integrated genes, would increase after integration. CONCLUSION: These results suggested that HPV may utilize the powerful hosts' promoters to express viral oncogenes and overexpression of viral oncogenes plays a significant role in the carcinogenesis of SCC. In AC, HPV integration may affect hosts' oncogenes, and the dysregulation of oncogenes may primarily contribute to progression of AC.

MIR-376A-2-5P as a potential prognostic marker for advanced penile squamous cell carcinomas trough HPV-dependent pathways.

Penile cancer (PeCa) is a rare tumor, generally associated with socioeconomic conditions in low-income countries. Hence, a delay in diagnosis and treatment leads in more advanced tumors, to higher comorbidity, and mortality. Human papillomavirus (HPV) infection has been identified as one of the major risk factors for PeCa. In addition, viral integration sites have been related to copy number alterations, impacting miRNAs/mRNA interactions and, consequently, the molecular pathways related to them. Nonetheless, studies on differentially expressed miRNAs (miRDEs) in PeCa are still scarce, especially in PeCa associated with high-risk HPV (hrHPV). To investigate the role of these gene regulators in PeCa progression, 827 miRNAs (Nanostring Technologies™, Seattle, WA, USA) were evaluated in 22 hrHPV-associated penile squamous cell carcinomas and five non-tumor penile tissues. For functions of miRNAs/target genes and relationship with HPV we conducted an integrated analysis by Diana Tools, KEGG, HPVbase, and InterSPPI-HVPPI platforms. We found that 25 miRNAs of the most differentially expressed impact 43 top molecular pathways, of which the fatty acid biosynthesis pathway, prions, miRNAs in cancer and hippo signaling (P<1.0-325, for each) were the most statistically significant. Notably, 23 out of 25 are located at HPV integration sites (HPVis). MiR-1206, miR-376b-3p and miR-495-3p were downregulated and associated with perineural invasion. In addition, a comparison between advanced and early diseases revealed 143 miRDEs. ROC analysis of a single (miR-376a-2-5p), paired (miR-376a-2-5p, miR-551b-3p) or combination of five miRDEs (miR-99a-5p, miR-150-5p, miR-155-5p, let-7c-5p, miR-342-3p) showed robust discriminatory power (AUC = 0.9; P = 0.0114, for each). Strikingly, miR-376a-2-5p exhibited the highest values of sensitivity and specificity, with 100% and 83.3%, respectively, indicating this miRNA as a potential prognostic marker in hrHPV-penile carcinogenesis.

Cis-regulatory effect of HPV integration is constrained by host chromatin architecture in cervical cancers.

Human papillomavirus (HPV) infections are the primary drivers of cervical cancers, and often HPV DNA gets integrated into the host genome. Although the oncogenic impact of HPV encoded genes is relatively well known, the cis-regulatory effect of integrated HPV DNA on host chromatin structure and gene regulation remains less understood. We investigated genome-wide patterns of HPV integrations and associated host gene expression changes in the context of host chromatin states and topologically associating domains (TADs). HPV integrations were significantly enriched in active chromatin regions and depleted in inactive ones. Interestingly, regardless of chromatin state, genomic regions flanking HPV integrations showed transcriptional upregulation. Nevertheless, upregulation (both local and long-range) was mostly confined to TADs with integration, but not affecting adjacent TADs. Few TADs showed recurrent integrations associated with overexpression of oncogenes within them (e.g. MYC, PVT1, TP63 and ERBB2) regardless of proximity. Hi-C and 4C-seq analyses in cervical cancer cell line (HeLa) demonstrated chromatin looping interactions between integrated HPV and MYC/PVT1 regions (~ 500 kb apart), leading to allele-specific overexpression. Based on these, we propose HPV integrations can trigger multimodal oncogenic activation to promote cancer progression.

Insights into the Mechanisms and Structure of Breakage-Fusion-Bridge Cycles in Cervical Cancer using Long-Read Sequencing.

Cervical cancer is caused by human papillomavirus (HPV) infection, has few approved targeted therapeutics, and is the most common cause of cancer death in low-resource countries. We characterized 19 cervical and four head and neck cell lines using long-read DNA and RNA sequencing and identified the HPV types, HPV integration sites, chromosomal alterations, and cancer driver mutations. Structural variation analysis revealed telomeric deletions associated with DNA inversions resulting from breakage-fusion-bridge (BFB) cycles. BFB is a common mechanism of chromosomal alterations in cancer, and this is one of the first analyses of these events using long-read sequencing. Analysis of the inversion sites revealed staggered ends consistent with exonuclease digestion of the DNA after breakage. Some BFB events are complex, involving inter- or intra-chromosomal insertions or rearrangements. None of the BFB breakpoints had telomere sequences added to resolve the dicentric chromosomes and only one BFB breakpoint showed chromothripsis. Five cell lines have a Chr11q BFB event, with YAP1/BIRC2/BIRC3 gene amplification. Indeed, YAP1 amplification is associated with a 10-year earlier age of diagnosis of cervical cancer and is three times more common in African American women. This suggests that cervical cancer patients with YAP1/BIRC2/BIRC3-amplification, especially those of African American ancestry, might benefit from targeted therapy. In summary, we uncovered new insights into the mechanisms and consequences of BFB cycles in cervical cancer using long-read sequencing.