The intratumor microbiome varies by geographical location and anatomical site in head and neck squamous cell carcinoma.

Head and Neck Squamous Cell Carcinoma (HNSCC) is a highly heterogeneous cancer that is characterized by distinct phenotypes based on anatomical site and etiological agents. Recently, the intratumor microbiome has been implicated in cancer pathogenesis and progression. Although it is well established that the gut microbiome varies with geographical location and is highly influenced by factors such as diet, environment, and genetics, the intratumor microbiome is not very well characterized. In this review, we aim to characterize the HNSCC intratumor microbiome by geographical location and anatomical site. We conducted a review of primary literature from PubMed and assessed studies based on relevancy and recency. To the best of our knowledge, we are the first to comprehensively examine the tumor microenvironment of HNSCC with respect to these two primary factors on a large scale. Our results suggest that there are unique bacterial and fungal biomarkers for HNSCC for each of the following geographical locations: North America, Asia, Europe, Australia, and Africa. We also identified a panel of microbial biomarkers that are unique to two primary HNSCC anatomic sites, as well as microbial biomarkers associated with various etiological agents of HNSCC. Future study of these microbes may improve HNSCC diagnostic and therapeutic modalities by accounting for differences based on geographic regions and anatomical sites.

Eliminating Algorithmic Racial Bias in Clinical Decision Support Algorithms: Use Cases from the Veterans Health Administration.

The Veterans Health Administration uses equity- and evidence-based principles to examine, correct, and eliminate use of potentially biased clinical equations and predictive models. We discuss the processes, successes, challenges, and next steps in four examples. We detail elimination of the race modifier for estimated kidney function and discuss steps to achieve more equitable pulmonary function testing measurement. We detail the use of equity lenses in two predictive clinical modeling tools: Stratification Tool for Opioid Risk Mitigation (STORM) and Care Assessment Need (CAN) predictive models. We conclude with consideration of ways to advance racial health equity in clinical decision support algorithms.

Development and validation of a standardized blood culture contamination definition and metric dashboard for a large health care system.

OBJECTIVES: Blood culture contamination is a major problem in health care, with significant impacts on both patient safety and cost. Initiatives to reduce blood culture contamination require a reliable, consistent metric to track the success of interventions. The objective of our project was to establish a standardized definition of blood culture contamination suitable for use in a Veterans Health Administration (VHA) national data query, then to validate this definition and query. A secondary objective was to construct a national VHA data dashboard to display the data from this query that could be used in VHA quality improvement projects aimed at reducing blood culture contamination. METHODS: A VHA microbiology expert work group was formed to generate a standardized definition and oversee the validation studies. The standardized definition was used to generate data for calendar year 2021 using a Structured Query Language data query. Twelve VHA hospital microbiology laboratories compared the data from the query against their own locally derived contamination data and recorded those data in a data collection worksheet that all sites used. Data were collated and presented to the work group. RESULTS: More than 50,000 blood culture accessions were in the validation data set, with more than 1,200 contamination events. The overall blood culture contamination rate for the 12 facilities participating was 2.56% with local definitions and data and 2.43% with the standardized definitions and data query. The main differences noted between the 2 data sets were deemed to be issues in local definitions. The query and definition were then converted into a national data dashboard that all VHA facilities can now access. CONCLUSIONS: A standardized definition for blood culture contamination and a national data query were validated for enterprise-wide VHA use. To our knowledge, this represents the first reported standardized, validated, and automated approach for calculating and tracking blood culture contamination. This tool will be key in quality initiatives aimed at reducing contamination events in VHA.

Clinical characteristics and comorbidities associated with SARS-CoV-2 breakthrough infection in the University of California Healthcare Systems.

BACKGROUND: To evaluate the degree to which clinical comorbidities or combinations of comorbidities are associated with SARS-CoV-2 breakthrough infection. MATERIALS AND METHODS: A breakthrough infection was defined as a positive test at least 14 days after a full vaccination regimen. Logistic regression was used to calculate aORs, which were adjusted for age, sex, and race information. RESULTS: A total of 110,380 patients from the UC CORDS database were included. After adjustment, stage 5 CKD due to hypertension (aOR: 7.33; 95% CI: 4.86-10.69; p<.001; power=1) displayed higher odds of infection than any other comorbidity. Lung transplantation history (aOR: 4.79; 95% CI: 3.25-6.82; p<.001; power= 1), coronary atherosclerosis (aOR: 2.12; 95% CI: 1.77-2.52; p<.001; power=1), and vitamin D deficiency (aOR: 1.87; 95% CI: 1.69-2.06; p<.001; power=1) were significantly correlated to breakthrough infection. Patients with obesity in addition to essential hypertension (aOR: 1.74; 95% CI: 1.51-2.01; p<.001; power=1) and anemia (aOR: 1.80; 95% CI: 1.47-2.19; p<.001; power=1) were at additional risk of breakthrough infection compared to those with essential hypertension and anemia alone. CONCLUSIONS: Further measures should be taken to prevent breakthrough infection for individuals with these conditions, such as acquiring additional doses of the SARS-CoV-2 vaccine to boost immunity.

Nationwide Genomic Surveillance and Response to COVID-19: The VA SeqFORCE and SeqCURE Consortiums.

Background: The US Department of Veterans Affairs (VA) has dedicated significant resources toward countering the COVID-19 pandemic. Sequencing for Research Clinical and Epidemiology (SeqFORCE) and Sequencing Collaborations United for Research and Epidemiology (SeqCURE) were developed as clinical and research consortiums, respectively, focused on the genetic COVID-19 surveillance. Observations: Through genetic sequencing, VA SeqFORCE and SeqCURE collaborations contributed to the COVID-19 pandemic response and scientific understanding. Future directions for each program include the assessment of the unique impact of COVID-19 on the veteran population, as well as the adaptation of these programs to future infectious disease threats. We foresee the use of these established platforms beyond infectious diseases. Conclusions: VA SeqFORCE and SeqCURE were established as clinical and research programs dedicated to sequencing COVID-19 as part of ongoing clinical and surveillance efforts. In the future, we anticipate that having these programs embedded within the largest integrated health care system in the US will enable the study of pathogens and pandemics beyond COVID-19 and at an unprecedented scale. The investment in these programs will form an integral part of our nation's response to emerging infectious diseases, with future applications to precision medicine and beyond.

The renal clear cell carcinoma immune landscape.

A comprehensive evaluation of the clear cell renal cell carcinoma (ccRCC) immune landscape was found using 584 RNA-sequencing datasets from The Cancer Genome Atlas (TCGA), we identified 17 key dysregulated immune-associated genes in ccRCC based on association with clinical variables and important immune pathways. Of the numerous findings from our analyses, we found that several of the 17 key dysregulated genes are heavily involved in interleukin and NF-kB signaling and that somatic copy number alteration (SCNA) hotspots may be causally associated with gene dysregulation. More importantly, we also found that key immune-associated genes and pathways are strongly upregulated in ccRCC. Our study may lend novel insights into the clinical implications of immune dysregulation in ccRCC and suggests potential immunotherapeutic targets for further evaluation.

Analysis of the immune landscape in virus-induced cancers using a novel integrative mechanism discovery approach.

BACKGROUND: The mechanisms of carcinogenesis from viral infections are extraordinarily complex and not well understood. Traditional methods of analyzing RNA-sequencing data may not be sufficient for unraveling complicated interactions between viruses and host cells. Using RNA and DNA-sequencing data from The Cancer Genome Atlas (TCGA), we aim to explore whether virus-induced tumors exhibit similar immune-associated (IA) dysregulations using a new algorithm we developed that focuses on the most important biological mechanisms involved in virus-induced cancers. Differential expression, survival correlation, and clinical variable correlations were used to identify the most clinically relevant IA genes dysregulated in 5 virus-induced cancers (HPV-induced head and neck squamous cell carcinoma, HPV-induced cervical cancer, EBV-induced stomach cancer, HBV-induced liver cancer, and HCV-induced liver cancer) after which a mechanistic approach was adopted to identify pathways implicated in IA gene dysregulation. RESULTS: Our results revealed that IA dysregulations vary with the cancer type and the virus type, but cytokine signaling pathways are dysregulated in all virus-induced cancers. Furthermore, we also found that important similarities exist between all 5 virus-induced cancers in dysregulated clinically relevant oncogenic signatures and IA pathways. Finally, we also discovered potential mechanisms for genomic alterations to induce IA gene dysregulations using our algorithm. CONCLUSIONS: Our study offers a new approach to mechanism identification through integrating functional annotations and large-scale sequencing data, which may be invaluable to the discovery of new immunotherapy targets for virus-induced cancers.

Tobacco Smoke and Electronic Cigarette Vapor Alter Enhancer RNA Expression That Can Regulate the Pathogenesis of Lung Squamous Cell Carcinoma.

Tobacco is the primary etiologic agent in worsened lung squamous cell carcinoma (LUSC) outcomes. Meanwhile, it has been shown that etiologic agents alter enhancer RNAs (eRNAs) expression. Therefore, we aimed to identify the effects of tobacco and electronic cigarette (e-cigarette) use on eRNA expression in relation to LUSC outcomes. We extracted eRNA counts from RNA-sequencing data of tumor/adjacent normal tissue and before/after e-cigarette tissue from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO), respectively. Tobacco-mediated LUSC eRNAs were correlated to patient survival, clinical variables, and immune-associated elements. eRNA expression was also correlated to mutation rates through the Repeated Evaluation of Variables Conditional Entropy and Redundance (REVEALER) algorithm and methylated sites through methylationArrayAnalysis. Differential expression analysis was then completed for the e-cigarette data to compare with key tobacco-mediated eRNAs. We identified 684 downregulated eRNAs and 819 upregulated eRNAs associated with tobacco-mediated LUSC, specifically, with the cancer pathological stage. We also observed a decrease in immune cell abundance in tobacco-mediated LUSC. Yet, we found an increased association of eRNA expression with immune cell abundance in tobacco-mediated LUSC. We identified 16 key eRNAs with significant correlations to 8 clinical variables, implicating these eRNAs in LUSC malignancy. Furthermore, we observed that these 16 eRNAs were highly associated with chromosomal alterations and reduced CpG site methylation. Finally, we observed large eRNA expression upregulation with e-cigarette use, which corresponded to the upregulation of the 16 key eRNAs. Our findings provide a novel mechanism by which tobacco and e-cigarette smoke influences eRNA interactions to promote LUSC pathogenesis and provide insight regarding disease progression at a molecular level.

Test Volume Ratio Benchmarking to Identify and Reduce Low-Value Laboratory Utilization.

OBJECTIVES: We analyzed test volume data to identify low-value test utilization. We subsequently tracked the efficacy of interventions to improve test utilization by decreasing low-value testing. METHODS: Test volume data for analytes included in the Choosing Wisely guidelines were analyzed to identify population outliers. Outliers were defined by test volume ratios of either analyte to sodium or paired analytes to correct for variation in patient volumes at each site. Interventions to improve test utilization were targeted to outlier sites. Relative efficacy in reducing low-value testing was tracked at those sites. RESULTS: After appropriate data cleaning, test volume ratios for 17 analytes paired with sodium and 8 pairs of analytes were acquired from 108 national sites. A site with abnormally high Clostridium difficile/sodium ratio was selected for intervention, leading to a 71% decrease in C difficile tests. Two different interventions to decrease creatine kinase MB isoform (CKMB) testing were performed at two unique sites with abnormally high CKMB/troponin ratios. These interventions decreased CKMB by 11% and 98% at the different sites, showing the efficacy of the different kinds of interventions. CONCLUSIONS: Test volume ratio analysis and benchmarking enable identification of low-value test utilization.

Identifying End Users' Preferences about Structuring Pharmacogenetic Test Orders in an Electronic Health Record System.

Pharmacogenetics (PGx) testing can be used for detecting genetic variations that may affect an individual's anticipated metabolism of, or response to, medications. Although several studies have focused on developing tools for delivering results from PGx testing, there is a relative dearth of information about how to design provider-friendly electronic order-entry systems for PGx. The U.S. Department of Veterans Affairs (VA) is preparing to implement a new electronic health records system. In this study, VA PGx test end users were surveyed about their preferences for how electronic test orders for PGx should be structured, including the nomenclature that should be used to search for and identify PGx-test orders, whether to offer single- versus multigene tests, and whether information about test methodology should be included in the order name. Responses were analyzed systematically to identify areas of agreement and disagreement with the survey options, and areas where respondents' opinions diverged. End users endorsed preferences for flexible ways to identify and order PGx tests and multigene panel tests; opinions on whether test methodology should be included in the test name were divergent. The results could be used for both informing the VA's new electronic health records implementation (including how PGx tests are searched for and ordered) and for providing insights for other health systems implementing PGx-testing programs.

Influence of Intratumor Microbiome on Clinical Outcome and Immune Processes in Prostate Cancer.

Although 1 in 9 American men will receive a diagnosis of prostate cancer (PC), most men with this diagnosis will not die from it, as most PCs are indolent. However, there is a subset of patients in which the once-indolent PC becomes metastatic and eventually, fatal. In this study, we analyzed microbial compositions of intratumor bacteria in PC to determine the influence of the microbiome on metastatic growth. Using large-scale RNA-sequencing data and corresponding clinical data, we correlated the abundance of microbes to immune pathways and PC risk factors, identifying specific microbes that either significantly deter or contribute to cancer aggressiveness. Interestingly, most of the microbes we found appeared to play anti-tumor roles in PC. Since these anti-tumor microbes were overrepresented in tumor samples, we believe that microbes thrive in the tumor microenvironment, outcompete cancer cells, and directly mitigate tumor growth by recruiting immune cells. These include Listeria monocytogenes, Methylobacterium radiotolerans JCM 2831, Xanthomonas albilineans GPE PC73, and Bradyrhizobium japonicum, which are negatively correlated with Gleason score, Tumor-Node-Metastasis (TNM) stage, prostate-specific antigen (PSA) level, and Androgen Receptor (AR) expression, respectively. We also identified microbes that contribute to tumor growth and are positively correlated with genomic alterations, dysregulated immune-associated (IA) genes, and prostate cancer stem cells (PCSC) genes.

Tobacco, but Not Nicotine and Flavor-Less Electronic Cigarettes, Induces ACE2 and Immune Dysregulation.

The COVID-19 pandemic caused by the SARS-CoV-2 virus, overlaps with the ongoing epidemics of cigarette smoking and electronic cigarette (e-cig) vaping. However, there is scarce data relating COVID-19 risks and outcome with cigarette or e-cig use. In this study, we mined three independent RNA expression datasets from smokers and vapers to understand the potential relationship between vaping/smoking and the dysregulation of key genes and pathways related to COVID-19. We found that smoking, but not vaping, upregulates ACE2, the cellular receptor that SARS-CoV-2 requires for infection. Both smoking and use of nicotine and flavor-containing e-cigs led to upregulation of pro-inflammatory cytokines and inflammasome-related genes. Specifically, chemokines including CCL20 and CXCL8 are upregulated in smokers, and CCL5 and CCR1 are upregulated in flavor/nicotine-containing e-cig users. We also found genes implicated in inflammasomes, such as CXCL1, CXCL2, NOD2, and ASC, to be upregulated in smokers and these e-cig users. Vaping flavor and nicotine-less e-cigs, however, did not lead to significant cytokine dysregulation and inflammasome activation. Release of inflammasome products, such as IL-1B, and cytokine storms are hallmarks of COVID-19 infection, especially in severe cases. Therefore, our findings demonstrated that smoking or vaping may critically exacerbate COVID-19-related inflammation or increase susceptibility to COVID-19.

The Landscape of Long Non-Coding RNA Dysregulation and Clinical Relevance in Muscle Invasive Bladder Urothelial Carcinoma.

Bladder cancer is one of the most common cancers in the United States, but few advancements in treatment options have occurred in the past few decades. This study aims to identify the most clinically relevant long non-coding RNAs (lncRNAs) to serve as potential biomarkers and treatment targets for muscle invasive bladder cancer (MIBC). Using RNA-sequencing data from 406 patients in The Cancer Genome Atlas (TCGA) database, we identified differentially expressed lncRNAs in MIBC vs. normal tissues. We then associated lncRNA expression with patient survival, clinical variables, oncogenic signatures, cancer- and immune-associated pathways, and genomic alterations. We identified a panel of 20 key lncRNAs that were most implicated in MIBC prognosis after differential expression analysis and prognostic correlations. Almost all lncRNAs we identified are correlated significantly with oncogenic processes. In conclusion, we discovered previously undescribed lncRNAs strongly implicated in the MIBC disease course that may be leveraged for diagnostic and treatment purposes in the future. Functional analysis of these lncRNAs may also reveal distinct mechanisms of bladder cancer carcinogenesis.

Papillary Thyroid Carcinoma Variants are Characterized by Co-dysregulation of Immune and Cancer Associated Genes.

Papillary thyroid carcinoma (PTC) variants exhibit different prognosis, but critical characteristics of PTC variants that contribute to differences in pathogenesis are not well-known. This study aims to characterize dysregulated immune-associated and cancer-associated genes in three PTC subtypes to explore how the interplay between cancer and immune processes causes differential prognosis. RNA-sequencing data from The Cancer Genome Atlas (TCGA) were used to identify dysregulated genes in each variant. The dysregulation profiles of the subtypes were compared using functional pathways clustering and correlations to relevant clinical variables, genomic alterations, and microRNA regulation. We discovered that the dysregulation profiles of classical PTC (CPTC) and the tall cell variant (TCPTC) are similar and are distinct from that of the follicular variant (FVPTC). However, unique cancer or immune-associated genes are associated with clinical variables for each subtype. Cancer-related genes MUC1, FN1, and S100-family members were the most clinically relevant in CPTC, while APLN and IL16, both immune-related, were clinically relevant in FVPTC. RAET-family members, also immune-related, were clinically relevant in TCPTC. Collectively, our data suggest that dysregulation of both cancer and immune associated genes defines the gene expression landscapes of PTC variants, but different cancer or immune related genes may drive the phenotype of each variant.

Identification and characterization of dysregulated P-element induced wimpy testis-interacting RNAs in head and neck squamous cell carcinoma.

It is clear that alcohol consumption is a major risk factor in the pathogenesis of head and neck squamous cell carcinoma (HNSCC); however, the molecular mechanism underlying the pathogenesis of alcohol-associated HNSCC remains poorly understood. The aim of the present study was to identify and characterize P-element-induced wimpy testis (PIWI)-interacting RNAs (piRNAs) and PIWI proteins dysregulated in alcohol-associated HNSCC to elucidate their function in the development of this cancer. Using next generation RNA-sequencing (RNA-seq) data obtained from 40 HNSCC patients, the piRNA and PIWI protein expression of HNSCC samples was compared between alcohol drinkers and non-drinkers. A separate piRNA expression RNA-seq analysis of 18 non-smoker HNSCC patients was also conducted. To verify piRNA expression, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed on the most differentially expressed alcohol-associated piRNAs in ethanol and acetaldehyde-treated normal oral keratinocytes. The correlation between piRNA expression and patient survival was analyzed using Kaplan-Meier estimators and multivariate Cox proportional hazard models. A comparison between alcohol drinking and non-drinking HNSCC patients demonstrated that a panel of 3,223 piRNA transcripts were consistently detected and differentially expressed. RNA-seq analysis and in vitro RT-qPCR verification revealed that 4 of these piRNAs, piR-35373, piR-266308, piR-58510 and piR-38034, were significantly dysregulated between drinking and non-drinking cohorts. Of these four piRNAs, low expression of piR-58510 and piR-35373 significantly correlated with improved patient survival. Furthermore, human PIWI-like protein 4 was consistently upregulated in ethanol and acetaldehyde-treated normal oral keratinocytes. These results demonstrate that alcohol consumption may cause dysregulation of piRNA expression in HNSCC and in vitro verifications identified 4 piRNAs that may be involved in the pathogenesis of alcohol-associated HNSCC.

Metformin Results in Diametrically Opposed Effects by Targeting Non-Stem Cancer Cells but Protecting Cancer Stem Cells in Head and Neck Squamous Cell Carcinoma.

Cancer stem cells (CSCs) have been shown as a distinct population of cancer cells strongly implicated with resistance to conventional chemotherapy. Metformin, the most widely prescribed drug for diabetes, was reported to target cancer stem cells in various cancers. In this study, we sought to determine the effects of metformin on head and neck squamous cell carcinoma (HNSCC). CSCs and non-stem HNSCC cells were treated with metformin and cisplatin alone, and in combination, and cell proliferation levels were measured through MTS assays. Next, potential targets of metformin were explored through computational small molecule binding analysis. In contrast to the reported effects of metformin on CSCs in other cancers, our data suggests that metformin protects HNSCC CSCs against cisplatin in vitro. Treatment with metformin resulted in a dose-dependent induction of the stem cell genes CD44, BMI-1, OCT-4, and NANOG. On the other hand, we observed that metformin successfully decreased the proliferation of non-stem HNSCC cells. Computational drug⁻protein interaction analysis revealed mitochondrial complex III to be a likely target of metformin. Based on our results, we present the novel hypothesis that metformin targets complex III to reduce reactive oxygen species (ROS) levels, leading to the differential effects observed on non-stem cancer cells and CSCs.

Pharmacogenetic testing in the Veterans Health Administration (VHA): policy recommendations from the VHA Clinical Pharmacogenetics Subcommittee.

PURPOSE: The Veterans Health Administration (VHA) Clinical Pharmacogenetics Subcommittee is charged with making recommendations about whether specific pharmacogenetic tests should be used in healthcare at VHA facilities. We describe a process to inform VHA pharmacogenetic testing policy. METHODS: After developing consensus definitions of clinical validity and utility, the Subcommittee identified salient drug-gene pairs with potential clinical application in VHA. Members met monthly to discuss each drug-gene pair, the evidence of clinical utility for the associated pharmacogenetic test, and any VHA-specific testing considerations. The Subcommittee classified each test as strongly recommended, recommended, or not routinely recommended before drug initiation. RESULTS: Of 30 drug-gene pair tests reviewed, the Subcommittee classified 4 (13%) as strongly recommended, including HLA-B*15:02 for carbamazepine-associated Stevens-Johnston syndrome and G6PD for rasburicase-associated hemolytic anemia; 12 (40%) as recommended, including CYP2D6 for codeine toxicity; and 14 (47%) as not routinely recommended, such as CYP2C19 for clopidogrel dosing. CONCLUSION: Only half of drug-gene pairs with high clinical validity received Subcommittee support for policy promoting their widespread use across VHA. The Subcommittee generally found insufficient evidence of clinical utility or available, effective alternative strategies for the remainders. Continual evidence review and rigorous outcomes research will help promote the translation of pharmacogenetic discovery to healthcare.

Functional Genomics Profiling of Bladder Urothelial Carcinoma MicroRNAome as a Potential Biomarker.

Though bladder urothelial carcinoma is the most common form of bladder cancer, advances in its diagnosis and treatment have been modest in the past few decades. To evaluate miRNAs as putative disease markers for bladder urothelial carcinoma, this study develops a process to identify dysregulated miRNAs in cancer patients and potentially stratify patients based on the association of their microRNAome phenotype to genomic alterations. Using RNA sequencing data for 409 patients from the Cancer Genome Atlas, we examined miRNA differential expression between cancer and normal tissues and associated differentially expressed miRNAs with patient survival and clinical variables. We then correlated miRNA expressions with genomic alterations using the Wilcoxon test and REVEALER. We found a panel of six miRNAs dysregulated in bladder cancer and exhibited correlations to patient survival. We also performed differential expression analysis and clinical variable correlations to identify miRNAs associated with tobacco smoking, the most important risk factor for bladder cancer. Two miRNAs, miR-323a and miR-431, were differentially expressed in smoking patients compared to nonsmoking patients and were associated with primary tumor size. Functional studies of these miRNAs and the genomic features we identified for potential stratification may reveal underlying mechanisms of bladder cancer carcinogenesis and further diagnosis and treatment methods for urothelial bladder carcinoma.

Computational methods reveal novel functionalities of PIWI-interacting RNAs in human papillomavirus-induced head and neck squamous cell carcinoma.

Human papillomavirus (HPV) infection is the fastest growing cause of head and neck squamous cell carcinoma (HNSCC) today, but its role in malignant transformation remains unclear. This study aimed to conduct a comprehensive investigation of PIWI-interacting RNA (piRNA) alterations and functionalities in HPV-induced HNSCC. Using 77 RNA-sequencing datasets from TCGA, we examined differential expression of piRNAs between HPV16(+) HNSCC and HPV(-) Normal samples, identifying a panel of 30 HPV-dysregulated piRNAs. We then computationally investigated the potential mechanistic significances of these transcripts in HPV-induced HNSCC, identifying our panel of piRNAs to associate with the protein PIWIL4 as well as the RTL family of retrotransposon-like genes, possibly through direct binding interactions. We also recognized several HPV-dysregulated transcripts for their correlations with well-documented mutations and copy number variations in HNSCC as well as HNSCC clinical variables, demonstrating the potential ability of our piRNAs to play important roles in large-scale modulation of HNSCC in addition to their direct, smaller-scale interactions in this malignancy. The differential expression of key piRNAs, including NONHSAT077364, NONHSAT102574, and NONHSAT128479, was verified in vitro by evaluating endogenous expression in HPV(+) cancer vs. HPV(-) normal cell lines. Overall, our novel study provides a rigorous investigation of piRNA dysregulation in HPV-related HNSCC, and lends critical insight into the idea that these small regulatory transcripts may play crucial and previously unidentified roles in tumor pathogenesis and progression.