{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,9,30]],"date-time":"2026-09-30T22:08:25Z","timestamp":1790806105037,"version":"4.1.0"},"reference-count":43,"publisher":"Springer Science and Business Media LLC","issue":"1","license":[{"start":{"date-parts":[[2021,12,20]],"date-time":"2021-12-20T00:00:00Z","timestamp":1639958400000},"content-version":"tdm","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"},{"start":{"date-parts":[[2021,12,20]],"date-time":"2021-12-20T00:00:00Z","timestamp":1639958400000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0"}],"funder":[{"DOI":"10.13039\/100000002","name":"U.S. Department of Health & Human Services | National Institutes of Health","doi-asserted-by":"publisher","award":["Intramural"],"award-info":[{"award-number":["Intramural"]}],"id":[{"id":"10.13039\/100000002","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["link.springer.com"],"crossmark-restriction":false},"short-container-title":["Commun Biol"],"abstract":"<jats:title>Abstract<\/jats:title><jats:p>Recurrent respiratory papillomatosis (RRP) is a debilitating neoplastic disorder of the upper aerodigestive tract caused by chronic infection with low-risk human papillomavirus types 6 or 11. Patients with severe RRP can require hundreds of lifetime surgeries to control their disease and pulmonary papillomatosis can be fatal. Here we report the comprehensive genomic and transcriptomic characterization of respiratory papillomas. We discovered and characterized distinct subtypes with transcriptional resemblance to either a basal or differentiated cell state that associate with disease aggressiveness and differ in key molecular, immune and APOBEC mutagenesis profiles. Through integrated comparison with high-risk HPV-associated head and neck squamous cell carcinoma, our analysis revealed divergent molecular and immune papilloma subtypes that form independent of underlying genomic alterations. Cumulatively our results support the development of dysregulated cellular proliferation and suppressed anti-viral immunity through distinct programs of squamous cell differentiation and associated expression of low-risk HPV genes. These analyses provide insight into the pathogenesis of respiratory papillomas and provide a foundation for the development of therapeutic strategies.<\/jats:p>","DOI":"10.1038\/s42003-021-02942-0","type":"journal-article","created":{"date-parts":[[2021,12,20]],"date-time":"2021-12-20T11:13:41Z","timestamp":1639998821000},"update-policy":"https:\/\/doi.org\/10.1007\/springer_crossmark_policy","source":"Crossref","is-referenced-by-count":19,"title":["Comprehensive multiomic characterization of human papillomavirus-driven recurrent respiratory papillomatosis reveals distinct molecular subtypes"],"prefix":"10.1038","volume":"4","author":[{"given":"Cem","family":"Sievers","sequence":"first","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Yvette","family":"Robbins","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Ke","family":"Bai","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Xinping","family":"Yang","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Paul E.","family":"Clavijo","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Jay","family":"Friedman","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Andrew","family":"Sinkoe","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Scott M.","family":"Norberg","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Christian","family":"Hinrichs","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Carter","family":"Van Waes","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0001-6586-5804","authenticated-orcid":false,"given":"Clint T.","family":"Allen","sequence":"additional","affiliation":[],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"297","published-online":{"date-parts":[[2021,12,20]]},"reference":[{"key":"2942_CR1","doi-asserted-by":"crossref","first-page":"340","DOI":"10.1002\/lio2.545","volume":"6","author":"JJ Benedict","year":"2021","unstructured":"Benedict, J. J. & Derkay, C. S. Recurrent respiratory papillomatosis: a 2020 perspective. Laryngoscope Investig. Otolaryngol. 6, 340\u2013345 (2021).","journal-title":"Laryngoscope Investig. Otolaryngol."},{"key":"2942_CR2","doi-asserted-by":"crossref","first-page":"2225","DOI":"10.1002\/lary.26662","volume":"127","author":"SR Best","year":"2017","unstructured":"Best, S. R., Mohr, M. & Zur, K. B. Systemic bevacizumab for recurrent respiratory papillomatosis: a national survey. Laryngoscope 127, 2225\u20132229 (2017).","journal-title":"Laryngoscope"},{"key":"2942_CR3","doi-asserted-by":"crossref","first-page":"119","DOI":"10.1186\/s40425-019-0603-3","volume":"7","author":"CT Allen","year":"2019","unstructured":"Allen, C. T. et al. Safety and clinical activity of PD-L1 blockade in patients with aggressive recurrent respiratory papillomatosis. J. Immunother. Cancer 7, 119 (2019).","journal-title":"J. Immunother. Cancer"},{"key":"2942_CR4","doi-asserted-by":"publisher","unstructured":"Meites, E. et al. Significant declines in juvenile onset recurrent respiratory papillomatosis following HPV vaccine introduction in the United States. Clin. Infect. Dis. https:\/\/doi.org\/10.1093\/cid\/ciab171 (2021).","DOI":"10.1093\/cid\/ciab171"},{"key":"2942_CR5","doi-asserted-by":"crossref","first-page":"398","DOI":"10.1093\/infdis\/jiy107","volume":"218","author":"H Artemchuk","year":"2018","unstructured":"Artemchuk, H. et al. Seroprevalences of antibodies to 11 human papillomavirus (HPV) types mark cumulative HPV exposure. J. Infect. Dis. 218, 398\u2013405 (2018).","journal-title":"J. Infect. Dis."},{"key":"2942_CR6","doi-asserted-by":"crossref","first-page":"e1000318","DOI":"10.1371\/journal.ppat.1000318","volume":"5","author":"D Pyeon","year":"2009","unstructured":"Pyeon, D., Pearce, S. M., Lank, S. M., Ahlquist, P. & Lambert, P. F. Establishment of human papillomavirus infection requires cell cycle progression. PLoS Pathog. 5, e1000318 (2009).","journal-title":"PLoS Pathog."},{"key":"2942_CR7","doi-asserted-by":"crossref","first-page":"153","DOI":"10.1016\/j.virol.2011.03.019","volume":"414","author":"GA Maglennon","year":"2011","unstructured":"Maglennon, G. A., McIntosh, P. & Doorbar, J. Persistence of viral DNA in the epithelial basal layer suggests a model for papillomavirus latency following immune regression. Virology 414, 153\u2013163 (2011).","journal-title":"Virology"},{"key":"2942_CR8","volume":"11","author":"B Zhu","year":"2020","unstructured":"Zhu, B. et al. Mutations in the HPV16 genome induced by APOBEC3 are associated with viral clearance. Nat. Commun. 11, 886 (2020).","journal-title":"Nat. Commun."},{"key":"2942_CR9","doi-asserted-by":"crossref","first-page":"970","DOI":"10.1038\/ng.2702","volume":"45","author":"SA Roberts","year":"2013","unstructured":"Roberts, S. A. et al. An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers. Nat. Genet. 45, 970\u2013976 (2013).","journal-title":"Nat. Genet."},{"key":"2942_CR10","doi-asserted-by":"crossref","first-page":"54","DOI":"10.1016\/j.immuni.2020.06.014","volume":"53","author":"E Bartok","year":"2020","unstructured":"Bartok, E. & Hartmann, G. Immune sensing mechanisms that discriminate self from altered self and foreign nucleic acids. Immunity 53, 54\u201377 (2020).","journal-title":"Immunity"},{"key":"2942_CR11","doi-asserted-by":"crossref","first-page":"576","DOI":"10.1038\/nature14129","volume":"517","author":"Cancer Genome Atlas, N.","year":"2015","unstructured":"Cancer Genome Atlas, N. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 517, 576\u2013582 (2015).","journal-title":"Nature"},{"key":"2942_CR12","doi-asserted-by":"crossref","first-page":"378","DOI":"10.1038\/nature21386","volume":"543","author":"Cancer Genome Atlas Research, N.","year":"2017","unstructured":"Cancer Genome Atlas Research, N. et al. Integrated genomic and molecular characterization of cervical cancer. Nature 543, 378\u2013384 (2017).","journal-title":"Nature"},{"key":"2942_CR13","doi-asserted-by":"crossref","first-page":"3475","DOI":"10.1038\/s41388-018-0657-6","volume":"38","author":"VL Cannataro","year":"2019","unstructured":"Cannataro, V. L. et al. APOBEC-induced mutations and their cancer effect size in head and neck squamous cell carcinoma. Oncogene 38, 3475\u20133487 (2019).","journal-title":"Oncogene"},{"key":"2942_CR14","doi-asserted-by":"crossref","DOI":"10.1038\/s41598-020-72927-2","volume":"10","author":"J Fernandez-Mateos","year":"2020","unstructured":"Fernandez-Mateos, J. et al. Oncogenic driver mutations predict outcome in a cohort of head and neck squamous cell carcinoma (HNSCC) patients within a clinical trial. Sci. Rep. 10, 16634 (2020).","journal-title":"Sci. Rep."},{"key":"2942_CR15","doi-asserted-by":"crossref","first-page":"448","DOI":"10.1093\/jnci\/djx243","volume":"110","author":"Q Mo","year":"2018","unstructured":"Mo, Q. et al. Prognostic power of a tumor differentiation gene signature for bladder urothelial carcinomas. J. Natl Cancer Inst. 110, 448\u2013459 (2018).","journal-title":"J. Natl Cancer Inst."},{"key":"2942_CR16","doi-asserted-by":"crossref","first-page":"119","DOI":"10.1016\/j.virusres.2016.12.017","volume":"231","author":"N Egawa","year":"2017","unstructured":"Egawa, N. & Doorbar, J. The low-risk papillomaviruses. Virus Res. 231, 119\u2013127 (2017).","journal-title":"Virus Res"},{"key":"2942_CR17","doi-asserted-by":"crossref","first-page":"6758","DOI":"10.1038\/sj.onc.1209943","volume":"25","author":"MS Hayden","year":"2006","unstructured":"Hayden, M. S., West, A. P. & Ghosh, S. NF-kappaB and the immune response. Oncogene 25, 6758\u20136780 (2006).","journal-title":"Oncogene"},{"key":"2942_CR18","doi-asserted-by":"crossref","first-page":"1236","DOI":"10.1097\/MLG.0b013e31816a7135","volume":"118","author":"CS Derkay","year":"2008","unstructured":"Derkay, C. S. & Wiatrak, B. Recurrent respiratory papillomatosis: a review. Laryngoscope 118, 1236\u20131247 (2008).","journal-title":"Laryngoscope"},{"key":"2942_CR19","doi-asserted-by":"crossref","first-page":"684","DOI":"10.1001\/archotol.133.7.684","volume":"133","author":"JK Stephen","year":"2007","unstructured":"Stephen, J. K. et al. An epigenetically derived monoclonal origin for recurrent respiratory papillomatosis. Arch. Otolaryngol. Head. Neck Surg. 133, 684\u2013692 (2007).","journal-title":"Arch. Otolaryngol. Head. Neck Surg."},{"key":"2942_CR20","doi-asserted-by":"crossref","first-page":"19055","DOI":"10.1073\/pnas.1906184116","volume":"116","author":"SB Drutman","year":"2019","unstructured":"Drutman, S. B. et al. Homozygous NLRP1 gain-of-function mutation in siblings with a syndromic form of recurrent respiratory papillomatosis. Proc. Natl Acad. Sci. USA 116, 19055\u201319063 (2019).","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2942_CR21","doi-asserted-by":"crossref","first-page":"55","DOI":"10.1016\/j.tim.2014.09.007","volume":"23","author":"IY Chen","year":"2015","unstructured":"Chen, I. Y. & Ichinohe, T. Response of host inflammasomes to viral infection. Trends Microbiol 23, 55\u201363 (2015).","journal-title":"Trends Microbiol"},{"key":"2942_CR22","doi-asserted-by":"crossref","first-page":"E320","DOI":"10.1002\/lary.24602","volume":"124","author":"R Rodman","year":"2014","unstructured":"Rodman, R., Mutasa, S., Dupuis, C., Spratt, H. & Underbrink, M. Genetic dysregulation in recurrent respiratory papillomatosis. Laryngoscope 124, E320\u2013E325 (2014).","journal-title":"Laryngoscope"},{"key":"2942_CR23","doi-asserted-by":"crossref","first-page":"289","DOI":"10.1177\/000348940511400407","volume":"114","author":"R Rahbar","year":"2005","unstructured":"Rahbar, R. et al. Role of vascular endothelial growth factor\u2014a in recurrent respiratory papillomatosis. Ann. Otol. Rhinol. Laryngol. 114, 289\u2013295 (2005).","journal-title":"Ann. Otol. Rhinol. Laryngol."},{"key":"2942_CR24","doi-asserted-by":"crossref","first-page":"608","DOI":"10.2119\/2008-00060.DeVoti","volume":"14","author":"JA DeVoti","year":"2008","unstructured":"DeVoti, J. A. et al. Immune dysregulation and tumor-associated gene changes in recurrent respiratory papillomatosis: a paired microarray analysis. Mol. Med. 14, 608\u2013617 (2008).","journal-title":"Mol. Med."},{"key":"2942_CR25","doi-asserted-by":"publisher","first-page":"2078","DOI":"10.1093\/bioinformatics\/btp352","volume":"25","author":"H Li","year":"2009","unstructured":"Li, H. et al. The sequence alignment\/map format and SAMtools. Bioinformatics 25, 2078\u20132079 (2009).","journal-title":"Bioinformatics"},{"key":"2942_CR26","doi-asserted-by":"crossref","first-page":"1297","DOI":"10.1101\/gr.107524.110","volume":"20","author":"A McKenna","year":"2010","unstructured":"McKenna, A. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297\u20131303 (2010).","journal-title":"Genome Res"},{"key":"2942_CR27","volume":"10","author":"M Wang","year":"2020","unstructured":"Wang, M. et al. SomaticCombiner: improving the performance of somatic variant calling based on evaluation tests and a consensus approach. Sci. Rep. 10, 12898 (2020).","journal-title":"Sci. Rep."},{"key":"2942_CR28","doi-asserted-by":"crossref","first-page":"11189","DOI":"10.1093\/nar\/gks918","volume":"40","author":"A Wilm","year":"2012","unstructured":"Wilm, A. et al. LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res. 40, 11189\u201311201 (2012).","journal-title":"Nucleic Acids Res."},{"key":"2942_CR29","volume":"17","author":"Y Fan","year":"2016","unstructured":"Fan, Y. et al. MuSE: accounting for tumor heterogeneity using a sample-specific error model improves sensitivity and specificity in mutation calling from sequencing data. Genome Biol. 17, 178 (2016).","journal-title":"Genome Biol."},{"key":"2942_CR30","doi-asserted-by":"crossref","first-page":"591","DOI":"10.1038\/s41592-018-0051-x","volume":"15","author":"S Kim","year":"2018","unstructured":"Kim, S. et al. Strelka2: fast and accurate calling of germline and somatic variants. Nat. Methods 15, 591\u2013594 (2018).","journal-title":"Nat. Methods"},{"key":"2942_CR31","doi-asserted-by":"crossref","first-page":"1220","DOI":"10.1093\/bioinformatics\/btv710","volume":"32","author":"X Chen","year":"2016","unstructured":"Chen, X. et al. Manta: rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics 32, 1220\u20131222 (2016).","journal-title":"Bioinformatics"},{"key":"2942_CR32","doi-asserted-by":"crossref","DOI":"10.1186\/s13059-016-0974-4","volume":"17","author":"W McLaren","year":"2016","unstructured":"McLaren, W. et al. The Ensembl variant effect predictor. Genome Biol. 17, 122 (2016).","journal-title":"Genome Biol."},{"key":"2942_CR33","doi-asserted-by":"crossref","first-page":"1747","DOI":"10.1101\/gr.239244.118","volume":"28","author":"A Mayakonda","year":"2018","unstructured":"Mayakonda, A., Lin, D. C., Assenov, Y., Plass, C. & Koeffler, H. P. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 28, 1747\u20131756 (2018).","journal-title":"Genome Res."},{"key":"2942_CR34","doi-asserted-by":"crossref","first-page":"2847","DOI":"10.1093\/bioinformatics\/btw313","volume":"32","author":"Z Gu","year":"2016","unstructured":"Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847\u20132849 (2016).","journal-title":"Bioinformatics"},{"key":"2942_CR35","doi-asserted-by":"crossref","first-page":"e1007603","DOI":"10.1371\/journal.pcbi.1007603","volume":"16","author":"C Flensburg","year":"2020","unstructured":"Flensburg, C., Sargeant, T., Oshlack, A. & Majewski, I. J. SuperFreq: Integrated mutation detection and clonal tracking in cancer. PLoS Comput. Biol. 16, e1007603 (2020).","journal-title":"PLoS Comput. Biol."},{"key":"2942_CR36","volume":"12","author":"B Li","year":"2011","unstructured":"Li, B. & Dewey, C. N. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323 (2011).","journal-title":"BMC Bioinform."},{"key":"2942_CR37","doi-asserted-by":"crossref","first-page":"D499","DOI":"10.1093\/nar\/gkw879","volume":"45","author":"K Van Doorslaer","year":"2017","unstructured":"Van Doorslaer, K. et al. The Papillomavirus Episteme: a major update to the papillomavirus sequence database. Nucleic Acids Res. 45, D499\u2013D506 (2017).","journal-title":"Nucleic Acids Res."},{"key":"2942_CR38","doi-asserted-by":"crossref","first-page":"448","DOI":"10.1101\/gr.257246.119","volume":"31","author":"S Uhrig","year":"2021","unstructured":"Uhrig, S. et al. Accurate and efficient detection of gene fusions from RNA sequencing data. Genome Res. 31, 448\u2013460 (2021).","journal-title":"Genome Res."},{"key":"2942_CR39","doi-asserted-by":"crossref","first-page":"15","DOI":"10.1093\/bioinformatics\/bts635","volume":"29","author":"A Dobin","year":"2013","unstructured":"Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15\u201321 (2013).","journal-title":"Bioinformatics"},{"key":"2942_CR40","volume":"15","author":"MI Love","year":"2014","unstructured":"Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).","journal-title":"Genome Biol."},{"key":"2942_CR41","doi-asserted-by":"crossref","first-page":"2084","DOI":"10.1093\/bioinformatics\/bty895","volume":"35","author":"A Zhu","year":"2019","unstructured":"Zhu, A., Ibrahim, J. G. & Love, M. I. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 35, 2084\u20132092 (2019).","journal-title":"Bioinformatics"},{"key":"2942_CR42","doi-asserted-by":"crossref","first-page":"D481","DOI":"10.1093\/nar\/gkv1351","volume":"44","author":"A Fabregat","year":"2016","unstructured":"Fabregat, A. et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 44, D481\u2013D487 (2016).","journal-title":"Nucleic Acids Res"},{"key":"2942_CR43","doi-asserted-by":"crossref","first-page":"284","DOI":"10.1089\/omi.2011.0118","volume":"16","author":"G Yu","year":"2012","unstructured":"Yu, G., Wang, L. G., Han, Y. & He, Q. Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284\u2013287 (2012).","journal-title":"OMICS"}],"container-title":["Communications Biology"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.nature.com\/articles\/s42003-021-02942-0.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/www.nature.com\/articles\/s42003-021-02942-0","content-type":"text\/html","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/www.nature.com\/articles\/s42003-021-02942-0.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2022,12,3]],"date-time":"2022-12-03T20:02:10Z","timestamp":1670097730000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.nature.com\/articles\/s42003-021-02942-0"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,12,20]]},"references-count":43,"journal-issue":{"issue":"1","published-online":{"date-parts":[[2021,12]]}},"alternative-id":["2942"],"URL":"https:\/\/doi.org\/10.1038\/s42003-021-02942-0","relation":{},"ISSN":["2399-3642"],"issn-type":[{"value":"2399-3642","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,12,20]]},"assertion":[{"value":"29 June 2021","order":1,"name":"received","label":"Received","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"2 December 2021","order":2,"name":"accepted","label":"Accepted","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"20 December 2021","order":3,"name":"first_online","label":"First Online","group":{"name":"ArticleHistory","label":"Article History"}},{"value":"The authors declare no competing interests.","order":1,"name":"Ethics","group":{"name":"EthicsHeading","label":"Competing interests"}}],"article-number":"1416"}}