Mutational Signature and Integrative Genomic Analysis of Human Papillomavirus-Associated Penile Squamous Cell Carcinomas from Latin American Patients
Authors:
- Luisa Matos CantoLCLuisa Matos Canto
- Jenilson Mota da SilvaJSJenilson Mota da Silva
- Patrícia Valèria Castelo-BrancoPCPatrícia Valèria Castelo-Branco
- Ingrid Monteiro da SilvaISIngrid Monteiro da Silva
- Leudivan NogueiraLNLeudivan Nogueira
- Carlos Eduardo Fonseca-AlvesCarlos Eduardo Fonseca-Alves
- André KhayatAKAndré Khayat
- Alexander BirbrairAlexander Birbrair
- Silma Regina PereiraSPSilma Regina Pereira
Open Access
Abstract
High-throughput DNA sequencing has allowed for the identification of genomic alterations and their impact on tumor development, progression, and therapeutic responses. In PSCC, for which the incidence has progressively increased worldwide, there are still limited data on the molecular mechanisms involved in the disease pathogenesis. In this study, we characterized the mutational signature of 30 human papillomavirus (HPV)-associated PSCC cases from Latin Americans, using whole-exome sequencing. Copy number variations (CNVs) were also identified and compared to previous array-generated data. Enrichment analyses were performed to reveal disrupted pathways and to identify alterations mapped to HPV integration sites (HPVis) and miRNA–mRNA hybridization regions. Among the most frequently mutated genes were NOTCH1, TERT, TTN, FAT1, TP53, CDKN2A, RYR2, CASP8, FBXW7, HMCN2, and ITGA8. Of note, 92% of these altered genes were localized at HPVis. We also found mutations in ten novel genes (KMT2C, SMARCA4, PTPRB, AJUBA, CR1, KMT2D, NBEA, FAM135B, GTF2I, and CIC), thus increasing our understanding of the potential HPV-disrupted pathways. Therefore, our study reveals innovative targets with potential therapeutic benefits for HPV-associated PSCCs. The CNV analysis by sequencing (CNV-seq) revealed five cancer-associated genes as the most frequent with gains (NOTCH1, MYC, NUMA1, PLAG1, and RAD21), while 30% of the tumors showed SMARCA4 with loss. Additionally, four cancer-associated genes (CARD11, CSMD3, KDR, and TLX3) carried untranslated regions (UTRs) variants, which may impact gene regulation by affecting the miRNAs hybridization regions. Altogether, these data contribute to the characterization of the mutational spectrum and its impact on cellular signaling pathways in PSCC, thus reinforcing the pivotal role of HPV infection in the molecular pathogenesis of these tumors.Keywords
Funding Information
- Fundação de Amparo à Pesquisa e ao Desenvolvimento Científico e Tecnológico do Maranhão (COOPI-07895/17, IECT-05551/18, 001, 02598/21)
- Coordenação de Aperfeiçoamento Pessoal de Nível Superior (COOPI-07895/17, IECT-05551/18, 001, 02598/21)
- Fundação de Amparo à Pesquisa e ao Desenvolvimento Científico e Tecnológico do Maranhão (FAPEMA) (COOPI-07895/17, IECT-05551/18, 001, 02598/21)
This publication has 82 references indexed in Scilit:
- ClinVar: public archive of relationships among sequence variation and human phenotypeNucleic Acids Research, 2013
- STRING v9.1: protein-protein interaction networks, with increased coverage and integrationNucleic Acids Research, 2012
- Cutadapt removes adapter sequences from high-throughput sequencing readsEMBnet.Journal, 2011
- A framework for variation discovery and genotyping using next-generation DNA sequencing dataNature Genetics, 2011
- Value of p16INK4a in the Pathology of Invasive Penile Squamous Cell CarcinomasAmerican Journal of Surgical Pathology, 2011
- The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing dataGenome Research, 2010
- ANNOVAR: functional annotation of genetic variants from high-throughput sequencing dataNucleic Acids Research, 2010
- Potent effect of target structure on microRNA functionNature Structural & Molecular Biology, 2007
- Fast and effective prediction of microRNA/target duplexesRNA, 2004
- Sfold web server for statistical folding and rational design of nucleic acidsNucleic Acids Research, 2004