NGS-based approach to determine the presence of HPV and their sites of integration in human cancer genome

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.

Integrated mutational landscape analysis of poorly differentiated high-grade neuroendocrine carcinoma of the uterine cervix

High-grade neuroendocrine cervical cancers (NETc) are exceedingly rare, highly aggressive tumors. We analyzed 64 NETc tumor samples by whole-exome sequencing (WES). Human papillomavirus DNA was detected in 65.6% (42/64) of the tumors. Recurrent mutations were identified in PIK3CA, KMT2D/MLL2, K-RAS, ARID1A, NOTCH2, and RPL10. The top mutated genes included RB1, ARID1A, PTEN, KMT2D/MLL2, and WDFY3, a gene not yet implicated in NETc. Somatic CNV analysis identified two copy number gains (3q27.1 and 19q13.12) and five copy number losses (1p36.21/5q31.3/6p22.2/9q21.11/11p15.5). Also, gene fusions affecting the ACLY-CRHR1 and PVT1-MYC genes were identified in one of the eight samples subjected to RNA sequencing. To resolve evolutionary history, multiregion WES in NETc admixed with adenocarcinoma cells was performed (i.e., mixed-NETc). Phylogenetic analysis of mixed-NETc demonstrated that adenocarcinoma and neuroendocrine elements derive from a common precursor with mutations typical of adenocarcinomas. Over one-third (22/64) of NETc demonstrated a mutator phenotype of C > T at CpG consistent with deficiencies in MBD4, a member of the base excision repair (BER) pathway. Mutations in the PI3K/AMPK pathways were identified in 49/64 samples. We used two patient-derived-xenografts (PDX) (i.e., NET19 and NET21) to evaluate the activity of pan-HER (afatinib), PIK3CA (copanlisib), and ATR (elimusertib) inhibitors, alone and in combination. PDXs harboring alterations in the ERBB2/PI3K/AKT/mTOR/ATR pathway were sensitive to afatinib, copanlisib, and elimusertib (P < 0.001 vs. controls). However, combinations of copanlisib/afatinib and copanlisib/elimusertib were significantly more effective in controlling NETc tumor growth. These findings define the genetic landscape of NETc and suggest that a large subset of these highly lethal malignancies might benefit from existing targeted therapies.

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.

Viral oncogenes, viruses, and cancer: a third-generation sequencing perspective on viral integration into the human genome

The link between viruses and cancer has intrigued scientists for decades. Certain viruses have been shown to be vital in the development of various cancers by integrating viral DNA into the host genome and activating viral oncogenes. These viruses include the Human Papillomavirus (HPV), Hepatitis B and C Viruses (HBV and HCV), Epstein-Barr Virus (EBV), and Human T-Cell Leukemia Virus (HTLV-1), which are all linked to the development of a myriad of human cancers. Third-generation sequencing technologies have revolutionized our ability to study viral integration events at unprecedented resolution in recent years. They offer long sequencing capabilities along with the ability to map viral integration sites, assess host gene expression, and track clonal evolution in cancer cells. Recently, researchers have been exploring the application of Oxford Nanopore Technologies (ONT) nanopore sequencing and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) sequencing in cancer research. As viral integration is crucial to the development of cancer via viruses, third-generation sequencing would provide a novel approach to studying the relationship interlinking viral oncogenes, viruses, and cancer. This review article explores the molecular mechanisms underlying viral oncogenesis, the role of viruses in cancer development, and the impact of third-generation sequencing on our understanding of viral integration into the human genome.

Genome sequencing reveals molecular subgroups in oral epithelial dysplasia

This study aimed to analyze the molecular characteristics of oral epithelial dysplasia (OED), highlighting the pathways and variants of genes that are frequently mutated in oral squamous cell carcinoma (OSCC) and other cancers. Ten archival OED cases were retrieved for retrospective clinicopathological analysis and exome sequencing. Comparative genomic analysis was performed between high-grade dysplasia (HGD) and low-grade dysplasia (LGD), focusing on 57 well-known cancer genes, of which 10 were previously described as the most mutated in OSCC. HGD cases had significantly more variants; however, a similar mutational landscape to OSCC was observed in both groups. CASP8+FAT1/HRAS, TP53, and miscellaneous molecular signatures were also present. FAT1 is the gene that is most affected by pathogenic variants. Hierarchical divisive clustering showed division between the two groups: "HGD-like cluster" with 4HGD and 2LGD and "LGD-like cluster" with 4 LGD. MLL4 pathogenic variants were exclusively in the "LGD-like cluster". TP53 was affected in one case of HGD; however, its pathway was usually altered. We describe new insights into the genetic basis of epithelial malignant transformation by genomic analysis, highlighting those associated with FAT1 and TP53. Some LGDs presented a similar mutational landscape to HGD after cluster analysis. Perhaps molecular alterations have not yet been reflected in histomorphology. The relative risk of malignant transformation in this molecular subgroup should be addressed in future studies.

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

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

Detection of Microorganisms in Body Fluid Samples

Next-generation sequencing (NGS) has been widely applied to the identification of microbiome in body fluids. The methodology of 16S rRNA amplicon sequencing is simple, fast, and cost-effective. It overcomes the problem that some microorganisms cannot be isolated or cultured. Low abundant bacteria can also be amplified and sequenced, but the resolution of classification can hardly reach species or sub-species level; moreover, this methodology is mainly used to identify bacterial populations, and other microorganisms like viruses or fungi cannot be sequenced. On the other hand, the microbiome profiling obtained by shotgun metagenomic sequencing is more comprehensive with better resolution, and more accurate classification can be expected due to higher coverage of genomic sequences from microorganisms. By combining the capture-based method with metagenomic sequencing, we can further enrich and detect low abundant microorganisms and identify the viral integration sites in host gDNA at once.

HPV genotyping by L1 amplicon sequencing of archived invasive cervical cancer samples: a pilot study

BACKGROUND

Human papillomavirus (HPV) is the primary cause of invasive cervical cancer (ICC). The prevalence of various HPV genotypes, ranging from oncogenically low- to high-risk, may be influenced by geographic and demographic factors, which could have critical implications for the screening and prevention of HPV infection and ICC incidence. However, many technical factors may influence the identification of high-risk genotypes associated with ICC in different populations.

METHODS

We used high-throughput sequencing of a single amplicon within the HPV L1 gene to assess the influence of patient age, race/ethnicity, histological subtype, sample type, collection date, experimental factors, and computational parameters on the prevalence of HPV genotypes detected in archived DNA (n = 34), frozen tissue (n = 44), and formalin-fixed paraffin-embedded (FFPE) tissue (n = 57) samples collected in the Los Angeles metropolitan area.

RESULTS

We found that the percentage of off-target human reads and the concentration of DNA amplified from each sample varied by HPV genotype and by archive type. After accounting for the percentage of human reads and excluding samples with especially low levels of amplified DNA, the HPV prevalence was 95% across all ICC samples: HPV16 was the most common genotype (in 56% of all ICC samples), followed by HPV18 (in 21%). Depending upon the genotyping parameters, the prevalence of HPV58 varied up to twofold in our cohort. In archived DNA and frozen tissue samples, we detected previously established differences in HPV16 and HPV18 frequencies based on histological subtype, but we could not reproduce those findings using our FFPE samples.

CONCLUSIONS

In this pilot study, we demonstrate that sample collection, preparation, and analysis methods can influence the detection of certain HPV genotypes and must be carefully considered when drawing any biological conclusions based on HPV genotyping data from ICC samples.

Fusobacterium nucleatum is associated with inflammation and poor survival in early-stage HPV-negative tongue cancer

Persistent pathogen infection is a known cause of malignancy, although with sparse systematic evaluation across tumor types. We present a comprehensive landscape of 1060 infectious pathogens across 239 whole exomes and 1168 transcriptomes of breast, lung, gallbladder, cervical, colorectal, and head and neck tumors. We identify known cancer-associated pathogens consistent with the literature. In addition, we identify a significant prevalence of Fusobacterium in head and neck tumors, comparable to colorectal tumors. The Fusobacterium-high subgroup of head and neck tumors occurs mutually exclusive to human papillomavirus, and is characterized by overexpression of miRNAs associated with inflammation, elevated innate immune cell fraction and nodal metastases. We validate the association of Fusobacterium with the inflammatory markers IL1B, IL6 and IL8, miRNAs hsa-mir-451a, hsa-mir-675 and hsa-mir-486-1, and MMP10 in the tongue tumor samples. A higher burden of Fusobacterium is also associated with poor survival, nodal metastases and extracapsular spread in tongue tumors defining a distinct subgroup of head and neck cancer.

Influence of tumor mutational burden, inflammatory gene expression profile, and PD-L1 expression on response to pembrolizumab in head and neck squamous cell carcinoma

Background

To characterize genomic determinants of response to pembrolizumab in recurrent/metastatic (R/M) head and neck squamous cell carcinoma (HNSCC) in the KEYNOTE-012 study.

Methods

Associations between biomarkers (tumor mutational burden (TMB), neoantigen load (NL), 18-gene T-cell-inflamed gene expression profile (Tcell_infGEP), and PD-L1 combined positive score (CPS)) and clinical outcomes with pembrolizumab were assessed in patients with R/M HNSCC (n=192). Tumor human papillomavirus (HPV) status was also evaluated with the use of p16 immunohistochemistry and whole exome sequencing (WES; HPV⁺, mapping >20 HPV reads) in pretreatment tumor samples (n=106).

Results

TMB, clonality-weighted TMB, and Tcell_infGEP were significantly associated with objective response (p=0.0276, p=0.0201, and p=0.006, respectively), and a positive trend was observed between NL and PD-L1 CPS and clinical response (p=0.0550 and p=0.0682, respectively). No correlation was observed between TMB and Tcell_infGEP (Spearman ρ=–0.026) or TMB and PD-L1 (Spearman ρ=0.009); a correlation was observed between Tcell_infGEP and PD-L1 (Spearman ρ=0.511). HPV status by WES and p16 immunohistochemistry showed concordance (84% ҡ=0.573) among patients whose HPV results were available using both methods.

Conclusions

TMB and inflammatory biomarkers (Tcell_infGEP and PD-L1) may represent distinct and complementary biomarkers predicting response to anti-programmed death 1 therapies in HNSCC; further study of these relationships in randomized clinical trials is needed.

Trial registration number

NCT01848834.

The Prevalence of High- and Low-Risk Types of HPV in Patients with Squamous Cell Carcinoma of the Head and Neck, Patients with Chronic Tonsillitis, and Healthy Individuals Living in Poland

Human papillomavirus (HPV) is a virus with the potential to infect human epithelial cells and an etiological agent of many types of cancer, including head and neck cancer. The aim of the study was to determine the prevalence of HPV infection in patients with head and neck squamous cell carcinoma (HNSCC), patients with chronic tonsillitis, and healthy individuals, and to establish high- and low-risk HPV genotypes in these groups. The objectives also comprised the delineation of the relationship between the infection with high- or low-risk HPV subtypes and clinicopathological and demographic characteristics of the study groups. This study was composed of 76 patients diagnosed with HNSCC, 71 patients with chronic tonsillitis, and 168 cases without either of these conditions (the control group). HPV detection and identification of subtypes were performed on isolated DNA using a test which allowed detection of 33 common high-risk and low-risk HPV subtypes. The prevalence of HPV infection was 42.1%, 25.4%, and 37.5% in HNSCC, chronic tonsillitis, and control groups, respectively. HPV 16 was the most prevalent genotype in all groups and the non-oncogenic HPV 43/44 was frequent in HNSCC patients. This analysis provides insight into the prevalence of oral oncogenic and non-oncogenic HPVs in patients with head and neck cancer, patients with chronic tonsillitis and healthy individuals, and leads to the conclusion that further investigations are warranted to examine a larger cohort of patients focusing on high- and low-risk HPV genotypes. Efforts should be focused on screening and prevention strategies, and therefore, it is important to introduce tools for effective detection of HPV genotypes. Furthermore, given the role of vaccines against oral HPV infection, our observations lead to the suggestion that HPV vaccination should be of considerable importance in public health strategies.

A NGS-based Blood Test For the Diagnosis of Invasive HPV-associated Carcinomas with Extensive Viral Genomic Characterization

PURPOSE

Use of circulating tumor DNA (ctDNA) for diagnosis is limited regarding the low number of target molecules in early-stage tumors. Human papillomavirus (HPV)-associated carcinomas represent a privileged model using circulating viral DNA (ctHPV DNA) as a tumor marker. However, the plurality of HPV genotypes represents a challenge. The next-generation sequencing (NGS)-based CaptHPV approach is able to characterize any HPV DNA sequence. To assess the ability of this method to establish the diagnosis of HPV-associated cancer via a blood sample, we analyzed ctHPV DNA in HPV-positive or HPV-negative carcinomas.

EXPERIMENTAL DESIGN

Patients (135) from France and Senegal with carcinoma developed in the uterine cervix (74), oropharynx (25), oral cavity (19), anus (12), and vulva (5) were prospectively registered. Matched tumor tissue and blood samples (10 mL) were taken before treatment and independently analyzed using the CaptHPV method.

RESULTS

HPV prevalence in tumors was 60.0% (81/135; 15 different genotypes). Viral analysis of plasmas compared with tumors was available for 134 patients. In the group of 80 patients with HPV-positive tumors, 77 were also positive in plasma (sensitivity 95.0%); in the group of 54 patients with HPV-negative tumors, one was positive in plasma (specificity 98.1%). In most cases, the complete HPV pattern observed in tumors could be established from the analysis of ctHPV DNA.

CONCLUSIONS

In patients with carcinoma associated with any HPV genotype, a complete viral genome characterization can be obtained via the analysis of a standard blood sample. This should favor the development of noninvasive diagnostic tests providing the identification of personalized tumor markers. See related commentary by Rostami et al., p. 5158.

Genomic analyses of high-grade neuroendocrine gynecological malignancies reveal a unique mutational landscape and therapeutic vulnerabilities

High‐grade neuroendocrine carcinoma of gynecologic origin (NEC‐GYN) is a highly aggressive cancer that often affects young women. The clinical management of NEC‐GYN is typically extrapolated from its counterpart, small cell carcinoma of the lung (SCLC), but, unfortunately, available therapies have limited benefit. In our NEC‐GYN cohort, median progression‐free survival (PFS) and overall survival (OS) were 1 and 12 months, respectively, indicating the highly lethal nature of this cancer. Our comprehensive genomic analyses unveiled that NEC‐GYN harbors a higher mutational burden with distinct mutational landscapes from SCLC. We identified 14 cancer driver genes, including the most frequently altered KMT2C (100%), KNL1 (100%), NCOR2 (100%), and CCDC6 (93%) genes. Transcriptomic analysis identified several novel gene fusions; astonishingly, the MALAT1 lincRNA gene was found in ˜ 20% of all fusion events in NEC‐GYN. Furthermore, NEC‐GYN exhibited a highly immunosuppressive state, intact RB1 expression, and was uniquely enriched with the YAP1^high molecular subtype. Our study identifies several potential therapeutic targets and suggests an urgent need to re‐evaluate the treatment options for NEC‐GYN.

Humans with inherited T cell CD28 deficiency are susceptible to skin papillomaviruses but are otherwise healthy

We study a patient with the human papilloma virus (HPV)-2-driven "tree-man" phenotype and two relatives with unusually severe HPV4-driven warts. The giant horns form an HPV-2-driven multifocal benign epithelial tumor overexpressing viral oncogenes in the epidermis basal layer. The patients are unexpectedly homozygous for a private CD28 variant. They have no detectable CD28 on their T cells, with the exception of a small contingent of revertant memory CD4+ T cells. T cell development is barely affected, and T cells respond to CD3 and CD2, but not CD28, costimulation. Although the patients do not display HPV-2- and HPV-4-reactive CD4+ T cells in vitro, they make antibodies specific for both viruses in vivo. CD28-deficient mice are susceptible to cutaneous infections with the mouse papillomavirus MmuPV1. The control of HPV-2 and HPV-4 in keratinocytes is dependent on the T cell CD28 co-activation pathway. Surprisingly, human CD28-dependent T cell responses are largely redundant for protective immunity.

Human Papillomavirus Detection by Whole-Genome Next-Generation Sequencing: Importance of Validation and Quality Assurance Procedures

Next-generation sequencing (NGS) yields powerful opportunities for studying human papillomavirus (HPV) genomics for applications in epidemiology, public health, and clinical diagnostics. HPV genotypes, variants, and point mutations can be investigated in clinical materials and described in previously unprecedented detail. However, both the NGS laboratory analysis and bioinformatical approach require numerous steps and checks to ensure robust interpretation of results. Here, we provide a step-by-step review of recommendations for validation and quality assurance procedures of each step in the typical NGS workflow, with a focus on whole-genome sequencing approaches. The use of directed pilots and protocols to ensure optimization of sequencing data yield, followed by curated bioinformatical procedures, is particularly emphasized. Finally, the storage and sharing of data sets are discussed. The development of international standards for quality assurance should be a goal for the HPV NGS community, similar to what has been developed for other areas of sequencing efforts including microbiology and molecular pathology. We thus propose that it is time for NGS to be included in the global efforts on quality assurance and improvement of HPV-based testing and diagnostics.

Whole-Exome Sequencing Analysis of Oral Squamous Cell Carcinoma Delineated by Tobacco Usage Habits

Oral squamous cell carcinoma (OSCC) is a common cancer of the oral cavity in India. Cigarette smoking and chewing tobacco are known risk factors associated with OSCC. However, genomic alterations in OSCC with varied tobacco consumption history are not well-characterized. In this study, we carried out whole-exome sequencing to characterize the mutational landscape of OSCC tumors from subjects with different tobacco consumption habits. We identified several frequently mutated genes, including TP53, NOTCH1, CASP8, RYR2, LRP2, CDKN2A, and ATM. TP53 and HRAS exhibited mutually exclusive mutation patterns. We identified recurrent amplifications in the 1q31, 7q35, 14q11, 22q11, and 22q13 regions and observed amplification of EGFR in 25% of samples with tobacco consumption history. We observed genomic alterations in several genes associated with PTK6 signaling. We observed alterations in clinically actionable targets including ERBB4, HRAS, EGFR, NOTCH1, NOTCH4, and NOTCH3. We observed enrichment of signature 29 in 40% of OSCC samples from tobacco chewers. Signature 15 associated with defective DNA mismatch repair was enriched in 80% of OSCC samples. NOTCH1 was mutated in 36% of samples and harbored truncating as well as missense variants. We observed copy number alterations in 67% of OSCC samples. Several genes associated with non-receptor tyrosine kinase signaling were affected in OSCC. These molecules can serve as potential candidates for therapeutic targeting in OSCC.

Circulating HPV DNA in the Management of Oropharyngeal and Cervical Cancers: Current Knowledge and Future Perspectives

Human papillomaviruses (HPVs) are associated with invasive malignancies, including almost 100% of cervical cancers (CECs), and 35–70% of oropharyngeal cancers (OPCs). HPV infection leads to clinical implications in related tumors by determining better prognosis and predicting treatment response, especially in OPC. Currently, specific and minimally invasive tests allow for detecting HPV-related cancer at an early phase, informing more appropriately therapeutical decisions, and allowing for timely disease monitoring. A blood-based biomarker detectable in liquid biopsy represents an ideal candidate, and the use of circulating HPV DNA (ct-DNA) itself could offer the highest specificity for such a scope. Circulating HPV DNA is detectable in the greatest part of patients affected by HPV-related cancers, and studies have demonstrated its potential usefulness for CEC and OPC clinical management. Unfortunately, when using conventional polymerase chain reaction (PCR), the detection rate of serum HPV DNA is low. Innovative techniques such as droplet-based digital PCR and next generation sequencing are becoming increasingly available for the purpose of boosting HPV ct-DNA detection rate. We herein review and critically discuss the most recent and representative literature, concerning the role of HPV ctDNA in OPC and CEC in the light of new technologies that could improve the potential of this biomarker in fulfilling many of the unmet needs in the clinical management of OPC and CEC patients.

Integrative sequencing discovers an ATF1-motif enriched molecular signature that differentiates hyalinizing clear cell carcinoma from mucoepidemoid carcinoma

OBJECTIVES

Salivary gland tumors are comprised of a diverse group of malignancies with widely varying prognoses. These cancers can be difficult to differentiate, especially in cases with limited potential for immunohistochemistry (IHC)-based characterization. Here, we sought to define the molecular profile of a rare salivary gland cancer called hyalinizing clear cell carcinoma (HCCC), and identify a molecular gene signature capable of distinguishing between HCCC and the histopathologically similar disease, mucoepidermoid carcinoma (MEC).

MATERIALS AND METHODS

We performed the first integrated full characterization of five independent HCCC cases.

RESULTS

We discovered insulin-like growth factor alterations and aberrant IGF2 and/or IGF1R expression in HCCC tumors, suggesting a potential dependence on this pathway. Further, we identified a 354 gene signature that differentiated HCCC from MEC, and was significantly enriched for genes with an ATF1 binding motif in their promoters, supporting a transcriptional pathogenic mechanism of the characteristic EWSR1-ATF1 fusion found in these tumors. Of the differentially expressed genes, IGF1R, SGK1 and SGK3 were found to be elevated in the HCCCs relative to MECs. Finally, analysis of immune checkpoints and subsequent IHC demonstrated that CXCR4 protein was elevated in several of the HCCC cases.

CONCLUSION

Collectively, our data identify an ATF1-motif enriched gene signature that may have clinical utility for molecular differentiation of HCCCs from other salivary gland tumors and discover potential actionable alterations that may benefit the clinical care of recurrent HCCC patients.

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.

An integrated approach to determine the abundance, mutation rate and phylogeny of the SARS-CoV-2 genome

The analysis of the SARS-CoV-2 genome datasets has significantly advanced our understanding of the biology and genomic adaptability of the virus. However, the plurality of advanced sequencing datasets-such as short and long reads-presents a formidable computational challenge to uniformly perform quantitative, variant or phylogenetic analysis, thus limiting its application in public health laboratories engaged in studying epidemic outbreaks. We present a computational tool, Infectious Pathogen Detector (IPD), to perform integrated analysis of diverse genomic datasets, with a customized analytical module for the SARS-CoV-2 virus. The IPD pipeline quantitates individual occurrences of 1060 pathogens and performs mutation and phylogenetic analysis from heterogeneous sequencing datasets. Using IPD, we demonstrate a varying burden (5.055-999655.7 fragments per million) of SARS-CoV-2 transcripts across 1500 short- and long-read sequencing SARS-CoV-2 datasets and identify 4634 SARS-CoV-2 variants (~3.05 variants per sample), including 449 novel variants, across the genome with distinct hotspot mutations in the ORF1ab and S genes along with their phylogenetic relationships establishing the utility of IPD in tracing the genome isolates from the genomic data (as accessed on 11 June 2020). The IPD predicts the occurrence and dynamics of variability among infectious pathogens-with a potential for direct utility in the COVID-19 pandemic and beyond to help automate the sequencing-based pathogen analysis and in responding to public health threats, efficaciously. A graphical user interface (GUI)-enabled desktop application is freely available for download for the academic users at http://www.actrec.gov.in/pi-webpages/AmitDutt/IPD/IPD.html and for web-based processing at http://ipd.actrec.gov.in/ipdweb/ to generate an automated report without any prior computational know-how.

Different Methods in HPV Genotyping of Anogenital and Oropharyngeal Lesions: Comparison between VisionArray® Technology, Next Generation Sequencing, and Hybrid Capture Assay

(1) Background: Human papillomaviruses (HPVs) are known to be related to the development of about 5% of all human cancers. The clinical relevance of HPV infection has been deeply investigated in carcinomas of the oropharyngeal area, uterine cervix, and anogenital area. To date, several different methods have been used for detecting HPV infection. The aim of the present study was to compare three different methods for the diagnosis of the presence of the HPV genome. (2) Methods: A total of 50 samples were analyzed. Twenty-five of them were tested using both next generation sequencing (NGS) and VisionArray® technology, the other 25 were tested using Hybrid Capture (HC) II assay and VisionArray® technology. (3) Results: A substantial agreement was obtained using NGS and VisionArray® (κ = 0.802), as well as between HC II and VisionArray® (κ = 0.606). In both analyses, the concordance increased if only high risk HPVs I(HR-HPVs) were considered as “positive”. (4) Conclusions: Our data highlighted the importance of technical choice in HPV characterization, which should be guided by the clinical aims, costs, starting material, and turnaround time for results.

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.

Genomic landscape analyses in cervical carcinoma and consequences for treatment

Where we are on the road to 'tailor-made' precision medicine for drug-resistant cervical carcinoma? We explored studies about analyses of viral and human genomes, epigenomes and transcriptomes, DNA mutation analyses, their importance in detecting HPV sequences, mechanisms of drug resistance to established and targeted therapies with small molecule or therapeutic antibodies, to radiosensitivity and to chemoradiotherapy. The value of repurposing of old drugs initially approved for other disease indications and now considered for cervix cancer therapy is also discussed. The microbiome influences drug response and survival too. HPV genomic integration sites were less significant. Nomograms (Lee et al., 2013) even outperformed FIGO staging regarding prediction of five-year overall survival times. We conclude that there are still many loose threads to be followed up, before coherent conclusions for individualized therapy of drug-resistant cervical carcinoma can be drawn.

Molecular Profiling Reveals Common and Specific Development Processes in Different Types of Gynecologic Cancers

BACKGROUND

Gynecologic cancers have become a major threat to women's health. The molecular biology of gynecologic cancers is not as well understood as that of breast cancer, and precision targeting is still new. Although viewed collectively as a group of cancers within the female reproductive system, they are more often studied separately. A comprehensive within-group comparison on molecular profiles is lacking.

METHODS

We conducted a whole-exome sequencing study of cervical/endometrial/ovarian cancer samples from 209 Chinese patients. We combined our data with genomic and transcriptomic data from relevant TCGA cohorts to identify and verify common/exclusive molecular changes in cervical/endometrial/ovarian cancer.

RESULTS

We identified shared molecular features including a COSMIC signature of deficient mismatch repair (dMMR), four recurrent copy-number variation (CNV) events, and extensive alterations in PI3K-Akt-mTOR signaling and cilium component genes; we also identified transcription factors and pathways that are exclusively altered in cervical/endometrial/ovarian cancer. The functions of the commonly/exclusively altered genomic circuits suggest (1) a common reprogramming process during early tumor initiation, which involves PI3K activation, defects in mismatch repair and cilium organization, as well as disruption in interferon signaling and immune recognition; (2) a cell-type specific program at late-stage tumor development that eventually lead to tumor proliferation and migration.

CONCLUSION

This study describes, from a molecular point of view, how similar and how different gynecologic cancers are, and it provides a hypothesis about the causes of the observed similarities and differences.

Merkel Cell Polyomavirus in Merkel Cell Carcinoma: Integration Sites and Involvement of the KMT2D Tumor Suppressor Gene

Merkel cell carcinoma (MCC) is an uncommon, lethal cancer of the skin caused by either Merkel cell polyomavirus (MCPyV) or UV-linked mutations. MCPyV is found integrated into MCC tumor genomes, accompanied by truncation mutations that render the MCPyV large T antigen replication incompetent. We used the open access HPV Detector/Cancer-virus Detector tool to determine MCPyV integration sites in whole-exome sequencing data from five MCC cases, thereby adding to the limited published MCPyV integration site junction data. We also systematically reviewed published data on integration for MCPyV in the human genome, presenting a collation of 123 MCC cases and their linked chromosomal sites. We confirmed that there were no highly recurrent specific sites of integration. We found that chromosome 5 was most frequently involved in MCPyV integration and that integration sites were significantly enriched for genes with binding sites for oncogenic transcription factors such as LEF1 and ZEB1, suggesting the possibility of increased open chromatin in these gene sets. Additionally, in one case we found, for the first time, integration involving the tumor suppressor gene KMT2D, adding to previous reports of rare MCPyV integration into host tumor suppressor genes in MCC.

HPV-EM: an accurate HPV detection and genotyping EM algorithm

Accurate HPV genotyping is crucial in facilitating epidemiology studies, vaccine trials, and HPV-related cancer research. Contemporary HPV genotyping assays only detect < 25% of all known HPV genotypes and are not accurate for low-risk or mixed HPV genotypes. Current genomic HPV genotyping algorithms use a simple read-alignment and filtering strategy that has difficulty handling repeats and homology sequences. Therefore, we have developed an optimized expectation–maximization algorithm, designated HPV-EM, to address the ambiguities caused by repetitive sequencing reads. HPV-EM achieved 97–100% accuracy when benchmarked using cell line data and TCGA cervical cancer data. We also validated HPV-EM using DNA tiling data on an institutional cervical cancer cohort (96.5% accuracy). Using HPV-EM, we demonstrated HPV genotypic differences in recurrence and patient outcomes in cervical and head and neck cancers.

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

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

Simultaneous detection and comprehensive analysis of HPV and microbiome status of a cervical liquid-based cytology sample using Nanopore MinION sequencing

Human papillomavirus (HPV) is a major pathogen that causes cervical cancer and many other related diseases. HPV infection related cervical microbiome could be an induce factor of cervical cancer. However, it is uncommon to find a single test on the market that can simultaneously provide information on both HPV and the microbiome. Herein, a novel method was developed in this study to simultaneously detect HPV infection and microbiota composition promptly and accurately. It provides a new and simple way to detect vaginal pathogen situation and also provide valuable information for clinical diagnose. This approach combined multiplex PCR, which targeted both HPV16 E6E7 and full-length 16S rRNA, and Nanopore sequencing to generate enough information to understand the vagina condition of patients. One HPV positive liquid-based cytology (LBC) sample was sequenced and analyzed. After comparing with Illumina sequencing, the results from Nanopore showed a similar microbiome composition. An instant sequencing evaluation showed that 15 min sequencing is enough to identify the top 10 most abundant bacteria. Moreover, two HPV integration sites were identified and verified by Sanger sequencing. This approach has many potential applications in pathogen detection and can potentially aid in providing a more rapid clinical diagnosis.

HPViewer: sensitive and specific genotyping of human papillomavirus in metagenomic DNA

Motivation

Shotgun DNA sequencing provides sensitive detection of all 182 HPV types in tissue and body fluid. However, existing computational methods either produce false positives misidentifying HPV types due to shared sequences among HPV, human and prokaryotes, or produce false negative since they identify HPV by assembled contigs requiring large abundant of HPV reads.

Results

We designed HPViewer with two custom HPV reference databases masking simple repeats and homology sequences respectively and one homology distance matrix to hybridize these two databases. It directly identified HPV from short DNA reads rather than assembled contigs. Using 100 100 simulated samples, we revealed that HPViewer was robust for samples containing either high or low number of HPV reads. Using 12 shotgun sequencing samples from respiratory papillomatosis, HPViewer was equal to VirusTAP, and Vipie and better than HPVDetector with the respect to specificity and was the most sensitive method in the detection of HPV types 6 and 11. We demonstrated that contigs-based approaches had disadvantages of detection of HPV. In 1573 sets of metagenomic data from 18 human body sites, HPViewer identified 104 types of HPV in a body-site associated pattern and 89 types of HPV co-occurring in one sample with other types of HPV. We demonstrated HPViewer was sensitive and specific for HPV detection in metagenomic data.

Availability and implementation

HPViewer can be accessed at https://github.com/yuhanH/HPViewer/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Establishment, molecular and biological characterization of HCB-514: a novel human cervical cancer cell line

Cervical cancer is the fourth most common cancer in women. Although cure rates are high for early stage disease, clinical outcomes for advanced, metastatic, or recurrent disease remain poor. To change this panorama, a deeper understanding of cervical cancer biology and novel study models are needed. Immortalized human cancer cell lines such as HeLa constitute crucial scientific tools, but there are few other cervical cancer cell lines available, limiting our understanding of a disease known for its molecular heterogeneity. This study aimed to establish novel cervical cancer cell lines derived from Brazilian patients. We successfully established one (HCB-514) out of 35 cervical tumors biopsied. We confirmed the phenotype of HCB-514 by verifying its’ epithelial and tumor origin through cytokeratins, EpCAM and p16 staining. It was also HPV-16 positive. Whole-exome sequencing (WES) showed relevant somatic mutations in several genes including BRCA2, TGFBR1 and IRX2. A copy number variation (CNV) analysis by nanostring and WES revealed amplification of genes mainly related to kinases proteins involved in proliferation, migration and cell differentiation, such as EGFR, PIK3CA, and MAPK7. Overexpression of EGFR was confirmed by phospho RTK-array and validated by western blot analysis. Furthermore, the HCB-514 cell line was sensitive to cisplatin. In summary, this novel Brazilian cervical cancer cell line exhibits relevant key molecular features and constitutes a new biological model for pre-clinical studies.

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.

Copy number variation: A prognostic marker for young patients with squamous cell carcinoma of the oral tongue

Background

The incidence of squamous cell carcinoma of the oral tongue (SCCOT) is increasing in people under age 40. There is an urgent need to identify prognostic markers that help identify young SCCOT patients with poor prognosis in order to select these for individualized treatment.

Materials and methods

To identify genetic markers that can serve as prognostic markers for young SCCOT patients, we first investigated four young (≤40 years) and five elderly patients (≥50 years) using global RNA sequencing and whole‐exome sequencing. Next, we combined our data with data on SCCOT from the cancer genome atlas (TCGA), giving a total of 16 young and 104 elderly, to explore the correlations between genomic variations and clinical outcomes.

Results

In agreement with previous studies, we found that SCCOT from young and elderly patients was transcriptomically and also genomically similar with no significant differences regarding cancer driver genes, germline predisposition genes, or the burden of somatic single nucleotide variations (SNVs). However, a disparate copy number variation (CNV) was found in young patients with distinct clinical outcome. Combined with data from TCGA, we found that the overall survival was significantly better in young patients with low‐CNV (n = 5) compared to high‐CNV (n = 11) burden (P = 0.044).

Conclusions

Copy number variation burden is a useful single prognostic marker for SCCOT from young, but not elderly, patients. CNV burden thus holds promise to form an important contribution when selecting suitable treatment protocols for young patients with SCCOT.

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.

Single-cell whole-genome sequencing identifies human papillomavirus integration in cervical tumour cells prior to and following radiotherapy

Single-cell sequencing technology is a promising systematic and comprehensive approach to delineate clonal associations between cells. The present study collected 13 and 12 cervical cells from fresh tumour tissue prior to and following radiotherapy, respectively, from a 46-year-old female patient with exogenous-type cervical carcinoma. Next, single-cell whole-genome sequencing analysis was performed on each cell. Examination revealed that normal cells could be clearly distinguished from tumour cells among the 25 cells. Tumour cells prior to and following radiotherapy almost represented two independent clones, with the main subpopulation prior to radiotherapy being killed and the minor subpopulation prior to radiotherapy becoming the main subpopulation following radiotherapy. A human papillomavirus (HPV) integration site was detected in POU class 5 homeobox 1B (POU5F1B) in tumour cells following radiotherapy, which has been reported to be a frequent HPV integration site in cervical carcinoma. These results indicate that tumour cells with HPV integration in POU5F1B survive radiotherapy, and that tumour cells prior to and following radiotherapy exhibit distinct characteristics.

ViFi: accurate detection of viral integration and mRNA fusion reveals indiscriminate and unregulated transcription in proximal genomic regions in cervical cancer

The integration of viral sequences into the host genome is an important driver of tumorigenesis in many viral mediated cancers, notably cervical cancer and hepatocellular carcinoma. We present ViFi, a computational method that combines phylogenetic methods with reference-based read mapping to detect viral integrations. In contrast with read-based reference mapping approaches, ViFi is faster, and shows high precision and sensitivity on both simulated and biological data, even when the integrated virus is a novel strain or highly mutated. We applied ViFi to matched genomic and mRNA data from 68 cervical cancer samples from TCGA and found high concordance between the two. Surprisingly, viral integration resulted in a dramatic transcriptional upregulation in all proximal elements, including LINEs and LTRs that are not normally transcribed. This upregulation is highly correlated with the presence of a viral gene fused with a downstream human element. Moreover, genomic rearrangements suggest the formation of apparent circular extrachromosomal (ecDNA) human-viral structures. Our results suggest the presence of apparent small circular fusion viral/human ecDNA, which correlates with indiscriminate and unregulated expression of proximal genomic elements, potentially contributing to the pathogenesis of HPV-associated cervical cancers. ViFi is available at https://github.com/namphuon/ViFi.

Utility of high-throughput DNA sequencing in the study of the human papillomaviruses

The Papillomaviridae family is probably the most diverse group of viruses that affect vertebrates. The study of the relationship between infection by certain types of human papillomavirus (HPV) and the development of neoplastic epithelial lesions is of particular interest because of the high prevalence of HPV-related carcinomas in populations of developing countries. To understand the mechanisms of infection and their association with different clinical manifestations, molecular tools play an important role in the description of new types of HPV, the characterization of effector properties of the viral factors, the specific diagnosis and monitoring of HPV types, and the alteration patterns at genetic level in the host. Technological advances in the field of DNA sequencing have led to the development of different next-generation sequencing systems, allowing obtaining a large amount of data and broadening the applications to study viral diseases. In this review, we summarize the main approaches and their perspectives where the use of massively parallel sequencing has been proved as a useful tool in the research of the HPV infection.

Genomic characterization of tobacco/nut chewing HPV-negative early stage tongue tumors identify MMP10 asa candidate to predict metastases

OBJECTIVES

Nodal metastases status among early stage tongue squamous cell cancer patients plays a decisive role in the choice of treatment, wherein about 70% patients can be spared from surgery with an accurate prediction of negative pathological lymph node status. This underscores an unmet need for prognostic biomarkers to stratify the patients who are likely to develop metastases.

MATERIALS AND METHODS

We performed high throughput sequencing of fifty four samples derived from HPV negative early stage tongue cancer patients habitual of chewing betel nuts, areca nuts, lime or tobacco using whole exome (n=47) and transcriptome (n=17) sequencing that were analyzed using in-house computational tools. Additionally, gene expression meta-analyses were carried out for 253 tongue cancer samples. The candidate genes were validated using qPCR and immuno-histochemical analysis in an extended set of 50 early primary tongue cancer samples.

RESULTS AND CONCLUSION

Somatic analysis revealed a classical tobacco mutational signature C:G>A:T transversion in 53% patients that were mutated in TP53, NOTCH1, CDKN2A, HRAS, USP6, PIK3CA, CASP8, FAT1, APC, and JAK1. Similarly, significant gains at genomic locus 11q13.3 (CCND1, FGF19, ORAOV1, FADD), 5p15.33 (SHANK2, MMP16, TERT), and 8q24.3 (BOP1); and, losses at 5q22.2 (APC), 6q25.3 (GTF2H2) and 5q13.2 (SMN1) were observed in these samples. Furthermore, an integrated gene-expression analysis of 253 tongue tumors suggested an upregulation of metastases-related pathways and over-expression of MMP10 in 48% tumors that may be crucial to predict nodal metastases in early tongue cancer patients. In overall, we present the first descriptive portrait of somatic alterations underlying the genome of tobacco/nut chewing HPV-negative early tongue cancer, and identify MMP10 asa potential prognostic biomarker to stratify those likely to develop metastases.

An Algorithm to Evaluate Suspected Lung Metastases in Patients with HPV-Associated Oropharyngeal Cancer

Distinguishing between distantly metastatic and metachronous lung primary carcinoma is challenging for patients with a history of head and neck cancer. There are implications for registry data, prognosis and related counseling, and management options, including eligibility for precision oncology trials. Patients with human papillomavirus (HPV)-associated oropharyngeal squamous cell carcinoma who were treated under a uniform clinical protocol and achieved a complete response were identified in a single-institution prospective head and neck cancer epidemiology database (n = 205). Fifteen patients presented with pulmonary nodule(s) after completion of therapy. We describe our algorithm for the evaluation of these patients, including histopathology, p16 immunohistochemistry, and HPV in situ hybridization.

Comprehensive Genomic Profiling of Metastatic Squamous Cell Carcinoma of the Anal Canal

Squamous cell carcinoma of the anal canal (SCCA) is a rare gastrointestinal malignancy with an increasing annual incidence globally. The majority of cases are linked to prior infection with the human papillomavirus (HPV). For patients with metastatic SCCA, no consensus standard treatment exists. Identification of relevant targeted agents as novel therapeutic approaches for metastatic SCCA has been limited by a lack of comprehensive molecular profiling. We performed whole-exome sequencing on tumor-normal pairs from 24 patients with metastatic SCCA. Tumor tissue from 17 additional patients was analyzed using a 263-gene panel as a validation cohort. Gene expression profiling was performed on available frozen tissue to assess for differential expression patterns. Based on these findings, patient-derived xenograft (PDX) models of SCCA were generated to test targeted therapies against PI3K and EGFR. Despite a low mutation burden, mutations in PIK3CA, MLL2, and MLL3 were among the most commonly mutated genes. An association between TP53 mutations and HPV-negative SCCA tumors was observed. Gene expression analysis suggested distinct tumor subpopulations harboring PIK3CA mutations and for which HPV had integrated into the host genome. In vivo studies demonstrated improvement with anti-EGFR treatment. Gene mutation frequencies, tumor mutation burden, and gene expression patterns for metastatic SCCA appear similar to other HPV-associated malignancies.Implications: This first comprehensive genomic characterization for patients with metastatic SCCA provides further rationale for the integration of SCCA into the development of novel targeted therapies across HPV-related cancers. Mol Cancer Res; 15(11); 1542-50. ©2017 AACR.

APOBEC-mediated genomic alterations link immunity and viral infection during human papillomavirus-driven cervical carcinogenesis

Cervical cancer is one of the most frequently diagnosed cancers and is a major cause of death from gynecologic cancers worldwide; the cancer burden from cervical cancer is especially heavy in less developed countries. Most cases of cervical cancer are caused by persistent infection with carcinogenic human papillomavirus (HPV) genotypes 16 and 18. Non-resolving inflammation caused by HPV infection provides a microenvironment that facilitates cancer development. Molecular alterations during the process of HPV-induced carcinogenesis are characterized by DNA methylation within the HPV genome, promoter hypermethylation of tumor suppressor genes in the host genome, as well as genomic instability caused by viral DNA integrating into the host genome. Catalytic polypeptide-like apolipoprotein B mRNA editing enzymes (APOBECs) normally function as part of the innate immune system. APOBEC expression is stimulated upon viral infection and plays an important role in HPV-induced cervical cancer. APOBECs catalyze the deamination of cytosine bases in nucleic acids, which leads to a conversion of target cytosine (C) to uracil (U) and consequently a change in the single-stranded DNA/RNA sequence. APOBEC proteins mediate the complex interactions between HPV and the host genome and link immunity and viral infection during HPV-driven carcinogenesis. Understanding the effects of APOBECs in HPV-induced cervical carcinogenesis will enable the development of better tools for HPV infection control and personalized prevention and treatment strategies.

Diagnosis of Sepsis with Cell-free DNA by Next-Generation Sequencing Technology in ICU Patients

BACKGROUND AND AIMS

Bacteremia is a common serious manifestation of disease in the intensive care unit (ICU), which requires quick and accurate determinations of pathogens to select the appropriate antibiotic treatment. To overcome the shortcomings of traditional bacterial culture (BC), we have adapted next-generation sequencing (NGS) technology to identify pathogens from cell-free plasma DNA.

METHODS

In this study, 78 plasma samples from ICU patients were analyzed by both NGS and BC methods and verified by PCR amplification/Sanger sequencing and ten plasma samples from healthy volunteers were analyzed by NGS as negative controls to define or calibrate the threshold of the NGS methodology.

RESULTS

Overall, 1578 suspected patient samples were found to contain bacteria or fungi by NGS, whereas ten patients were diagnosed by BC. Seven samples were diagnosed with bacterial or fungal infection both by NGS and BC. Among them, two samples were diagnosed with two types of bacteria by NGS, whereas one sample was diagnosed with two types of bacteria by BC, which increased the detectability of bacteria or fungi from 11 with BC to 17 with NGS. Most interestingly, 14 specimens were also diagnosed with viral infection by NGS. The overall diagnostic sensitivity was significantly increased from 12.82% (10/78) by BC alone to 30.77% (24/78) by NGS alone for ICU patients, which provides more useful information for establishing patient treatment plans.

CONCLUSION

NGS technology can be applied to detect bacteria in clinical blood samples as an emerging diagnostic tool rich in information to determine the appropriate treatment of septic patients.

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.

Somatic Variations in Cervical Cancers in Indian Patients

There are very few reports that describe the mutational landscape of cervical cancer, one of the leading cancers in Indian women. The aim of the present study was to investigate the somatic mutations that occur in cervical cancer. Whole exome sequencing of 10 treatment naïve tumour biopsies with matched blood samples, from a cohort of Indian patients with locally advanced disease, was performed. The data revealed missense mutations across 1282 genes, out of 1831 genes harbouring somatic mutations. These missense mutations (nonsynonymous + stop-gained) when compared with pre-existing mutations in the COSMIC database showed that 272 mutations in 250 genes were already reported although from cancers other than cervical cancer. More than 1000 novel somatic variations were obtained in matched tumour samples. Pathways / genes that are frequently mutated in various other cancers were found to be mutated in cervical cancers. A significant enrichment of somatic mutations in the MAPK pathway was observed, some of which could be potentially targetable. This is the first report of whole exome sequencing of well annotated cervical cancer samples from Indian women and helps identify trends in mutation profiles that are found in an Indian cohort of cervical cancer.

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.

Notch pathway activation is essential for maintenance of stem-like cells in early tongue cancer

BACKGROUND

Notch pathway plays a complex role depending on cellular contexts: promotes stem cell maintenance or induces terminal differentiation in potential cancer-initiating cells; acts as an oncogene in lymphocytes and mammary tissue or plays a growth-suppressive role in leukemia, liver, skin, and head and neck cancer. Here, we present a novel clinical and functional significance of NOTCH1 alterations in early stage tongue squamous cell carcinoma (TSCC).

PATIENTS AND METHODS

We analyzed the Notch signaling pathway in 68 early stage TSCC primary tumor samples by whole exome and transcriptome sequencing, real-time PCR based copy number, expression, immuno-histochemical, followed by cell based biochemical and functional assays.

RESULTS

We show, unlike TCGA HNSCC data set, NOTCH1 harbors significantly lower frequency of inactivating mutations (4%); is somatically amplified; and, overexpressed in 31% and 37% of early stage TSCC patients, respectively. HNSCC cell lines over expressing NOTCH1, when plated in the absence of attachment, are enriched in stem cell markers and form spheroids. Furthermore, we show that inhibition of NOTCH activation by gamma secretase inhibitor or shRNA mediated knockdown of NOTCH1 inhibits spheroid forming capacity, transformation, survival and migration of the HNSCC cells suggesting an oncogenic role of NOTCH1 in TSCC. Clinically, Notch pathway activation is higher in tumors of non-smokers compared to smokers (50% Vs 18%, respectively, P=0.026) and is also associated with greater nodal positivity compared to its non-activation (93% Vs 64%, respectively, P=0.029).

CONCLUSION

We anticipate that these results could form the basis for therapeutic targeting of NOTCH1 in tongue cancer.

Divergent viral presentation among human tumors and adjacent normal tissues

We applied a newly developed bioinformatics system called VirusScan to investigate the viral basis of 6,813 human tumors and 559 adjacent normal samples across 23 cancer types and identified 505 virus positive samples with distinctive, organ system- and cancer type-specific distributions. We found that herpes viruses (e.g., subtypes HHV4, HHV5, and HHV6) that are highly prevalent across cancers of the digestive tract showed significantly higher abundances in tumor versus adjacent normal samples, supporting their association with these cancers. We also found three HPV16-positive samples in brain lower grade glioma (LGG). Further, recurrent HBV integration at the KMT2B locus is present in three liver tumors, but absent in their matched adjacent normal samples, indicating that viral integration induced host driver genetic alterations are required on top of viral oncogene expression for initiation and progression of liver hepatocellular carcinoma. Notably, viral integrations were found in many genes, including novel recurrent HPV integrations at PTPN13 in cervical cancer. Finally, we observed a set of HHV4 and HBV variants strongly associated with ethnic groups, likely due to viral sequence evolution under environmental influences. These findings provide important new insights into viral roles of tumor initiation and progression and potential new therapeutic targets.

Characterization of Intra-Type Variants of Oncogenic Human Papillomaviruses by Next-Generation Deep Sequencing of the E6/E7 Region

Different human papillomavirus (HPV) types are characterized by differences in tissue tropism and ability to promote cell proliferation and transformation. In addition, clinical and experimental studies have shown that some genetic variants/lineages of high-risk HPV (HR-HPV) types are characterized by increased oncogenic activity and probability to induce cancer. In this study, we designed and validated a new method based on multiplex PCR-deep sequencing of the E6/E7 region of HR-HPV types to characterize HPV intra-type variants in clinical specimens. Validation experiments demonstrated that this method allowed reliable identification of the different lineages of oncogenic HPV types. Advantages of this method over other published methods were represented by its ability to detect variants of all HR-HPV types in a single reaction, to detect variants of HR-HPV types in clinical specimens with multiple infections, and, being based on sequencing of the full E6/E7 region, to detect amino acid changes in these oncogenes potentially associated with increased transforming activity.

Non-typhoidal Salmonella DNA traces in gallbladder cancer

BACKGROUND

We earlier proposed a genetic model for gallbladder carcinogenesis and its dissemination cascade. However, the association of gallbladder cancer and 'inflammatory stimulus' to drive the initial cascade in the model remained unclear. A recent study suggested infection with Salmonella can lead to changes in the host signalling pathways in gallbladder cancer.

FINDINGS

We examined the whole exomes of 26 primary gall bladder tumour and paired normal samples for presence of 143 HPV (Human papilloma virus) types along with 6 common Salmonella serotypes (S. typhi Ty2, S. typhi CT18, S. typhimurium LT2, S. choleraesuis SCB67, S. paratyphi TCC, and S. paratyphi SPB7) using a computational subtraction pipeline based on the HPVDetector, we recently described. Based on our evaluation of 26 whole exome gallbladder primary tumours and matched normal samples: association of typhoidal Salmonella species were found in 11 of 26 gallbladder cancer samples, and non-typhoidal Salmonella species in 12 of 26 gallbladder cancer, with 6 samples were found co-infected with both.

CONCLUSIONS

We present the first evidence to support the association of non-typhoidal Salmonella species along with typhoidal strains in gallbladder cancer. Salmonella infection in the chronic carrier state fits the role of the 'inflammatory stimulus' in the genetic model for gallbladder carcinogenesis that may play a role in gallbladder cancer analogous to Helicobacter pylori in gastric cancer.