Characteristics of Human and Microbiome RNA Profiles in Saliva.

Saliva is a convenient non-invasive source of liquid biopsy to monitor human health and diagnose diseases. In particular, extracellular vesicles (EVs) in saliva can potentially reveal clinically relevant information for systemic health. Recent studies have shown that RNA in saliva EVs could be exploited as biomarkers for disease diagnosis. However, there is no standardized protocol for profiling RNA in saliva EV nor clear guideline on selecting saliva fractions for biomarker analysis. To address these issues, we established a robust protocol for small RNA profiling from fractionated saliva. With this method, we performed comprehensive small RNA sequencing of four saliva fractions, including cell-free saliva (CFS), EV-depleted saliva (EV-D), exosome (EXO), and microvesicle (MV) from ten healthy volunteers. By comparing the expression profiles of total RNA from these fractions, we found that MV was most enriched in microbiome RNA (76.2% of total reads on average), whereas EV-D was notably enriched in human RNA (70.3% of total reads on average). As for human RNA composition, CFS and EV-D were both enriched in snoRNA and tRNA compared with the two EV fractions (EXO and MV, P < 0.05). Interestingly, EXO and MV had highly correlated expression profiles for various noncoding RNAs such as miRNA, tRNA, and yRNA. Our study revealed unique characteristics of circulating RNAs in various saliva fractions, which provides a guideline on preparing saliva samples to study specific RNA biomarkers of interest.

Sinonasal Adenosquamous Carcinoma - Morphology and Genetic Drivers Including Low- and High-Risk Human Papillomavirus mRNA, DEK::AFF2 Fusion, and MAML2 Rearrangement.

BACKGROUND: Sinonasal adenosquamous carcinoma is rare, and there are almost no studies detailing morphology or characterizing their genetic driver events. Further, many authors have termed sinonasal tumors with combined squamous carcinoma and glands as mucoepidermoid carcinoma but none have analyzed for the presence of MAML2 rearrangement. METHODS: Cases from 2014 to 2020 were collected and diagnosed using World Health Organization criteria. They were tested for p16 expression by immunohistochemistry (70% cut-off), DEK::AFF2 fusion by fluorescence in situ hybridization (FISH) and AFF2 immunohistochemistry, MAML2 rearrangement by FISH, and low- and high-risk HPV by RNA ISH and reverse transcription PCR, respectively. Detailed morphology and clinical features were reviewed. RESULTS: There were 7 male (64%) and 4 female (36%) patients with a median age of 69 years, most Caucasian (10 of 11 or 91%). Most had tobacco exposure (8/11, 73%) and most presented with epistaxis, a visible nasal mass, and/or facial pain. Several had a precursor papillomas (3 of 11, 27%). The squamous component had variable keratinization, 5 of 11 (46%) of which would be described as keratinizing, 3 non-keratinizing, and 2 with mixed features. All had gland formation, by definition, and 2 of 11 (18%) had ciliated tumor cells. None of the 11 cases had MAML2 rearrangement and one had DEK::AFF2 fusion with associated positive nuclear AFF2 protein immunostaining. Most were p16 positive (7 of 11, 64%) and all 7 of these were hrHPV positive either by RNA ISH or RT-PCR. Two of the p16-negative tumors were positive for lrHPV by RNA ISH. Treatment included surgery alone (4 of 11, 36%), surgery with adjuvant radiation (5 of 11, 45%), and surgery with radiation and chemotherapy (2 of 11, 18%). Four of 11 patients (36%) suffered disease recurrence, two requiring re-operation and who were disease free at last follow-up, one receiving additional chemotherapy and who was alive with disease. The other elected to undergo palliative therapy and died of disease. CONCLUSION: Sinonasal adenosquamous carcinoma is a somewhat heterogeneous tumor not infrequently arising ex papilloma and having various drivers including high- and low-risk HPV and rarely DEK::AFF2 fusion. The prognosis appears favorable when proper treatment is possible.

Miz1 promotes KRAS-driven lung tumorigenesis by repressing the protocadherin Pcdh10.

Targeting KRAS-mutated non-small-cell lung cancer (NSCLC) remains clinically challenging. Here we show that loss of function of Miz1 inhibits lung tumorigenesis in a mouse model of oncogenic KRAS-driven lung cancer. In vitro, knockout or silencing of Miz1 decreases cell proliferation, clonogenicity, migration, invasion, or anchorage-independent growth in mutant (MT) KRAS murine or human NSCLC cells but has unremarkable impact on non-tumorigenic cells or wild-type (WT) KRAS human NSCLC cells. RNA-sequencing reveals Protocadherin-10 (Pcdh10) as the top upregulated gene by Miz1 knockout in MT KRAS murine lung tumor cells. Chromatin immunoprecipitation shows Miz1 binding on the Pcdh10 promoter in MT KRAS lung tumor cells but not non-tumorigenic cells. Importantly, silencing of Pcdh10 rescues cell proliferation and clonogenicity in Miz1 knockout/knockdown MT KRAS murine or human tumor cells, and rescues allograft tumor growth of Miz1 knockout tumor cells in vivo. Miz1 is upregulated in MT KRAS lung tumor tissues compared with adjacent non-involved tissues in mice. Consistent with this, Miz1 is upregulated while Pcdh10 is downregulated in human lung adenocarcinomas (LUAD) compared with normal tissues, and high Miz1 levels or low Pcdh10 levels are associated with poor survival in lung cancer patients. Furthermore, the Miz1 signature is associated with worse survival in MT but not WT KRAS LUAD, and Pcdh10 is downregulated in MT compared to WT KRAS LUAD. Taken together, our studies implicate the Miz1/Pcdh10 axis in oncogenic KRAS-driven lung tumorigenesis.

Human Papillomavirus-Associated Oral Cavity Squamous Cell Carcinoma: An Entity with Distinct Morphologic and Clinical Features.

BACKGROUND: HPV-associated oral cavity squamous cell carcinoma (SCC) is not well-characterized in the literature, and also has a clinical significance that is poorly understood. METHODS: We gathered a cohort of oral cavity (OC) SCC with nonkeratinizing morphology, either in the invasive or in situ carcinoma (or both), tested for p16 by immunohistochemistry and high risk HPV E6/E7 mRNA by RTPCR (reference standard for transcriptionally-active high risk HPV) and gathered detailed morphologic and clinicopathologic data. RESULTS: Thirteen patients from two institutions were proven to be HPV-associated by combined p16 and high risk HPV mRNA positivity. All 13 patients (100%) were males, all were heavy smokers (average 57 pack/year), and most were active drinkers (9/11 or 81.8%). All 13 (100%) involved the tongue and/or floor of mouth. All had nonkeratinizing features, but maturing squamous differentiation varied widely (0-90%; mean 37.3%). Nonkeratinizing areas had high N:C ratios and larger nests, frequently with pushing borders, and minimal (or no) stromal desmoplasia. The carcinoma in situ, when present, was Bowenoid/nonkeratinizing with cells with high N:C ratios, full thickness loss of maturation, and abundant apoptosis and mitosis. HPV was type 16 in 11 patients (84.6%) and type 33 in two (15.4%). Nine patients had treatment data available. These underwent primary surgical resection with tumors ranging from 1.6 to 5.2 cm. Most had bone invasion (6/9-66.7% were T4a tumors), and most (6/9-66.7%) had extensive SCC in situ with all 6 of these patients having final margins positive for in situ carcinoma. CONCLUSIONS: HPV-associated OCSCC is an uncommon entity that shows certain distinct clinical and pathologic features. Recognition of these features may help pathologic diagnosis and could potentially help guide clinical management.

HPV+ HNSCC-derived exosomal miR-9-5p inhibits TGF-ß signaling-mediated fibroblast phenotypic transformation through NOX4.

Human papillomavirus (HPV) is a significant risk factor for head and neck squamous cell carcinoma (HNSCC). HPV+ HNSCC patients have a higher survival rate, which may be related to its unique tumor microenvironment. Exosomes are emerging as a communication tool between tumor cells and the tumor microenvironment, including cancer-associated fibroblasts (CAFs). In this study, 111 clinical samples tissues and public sequencing data were analyzed. Our study found fewer CAFs infiltrated in HPV+ HNSCC, and poor CAF infiltration level was associated with a good prognosis. HPV+ HNSCC cell-derived exosomes can significantly reduce the phenotypic transformation of fibroblasts. miR-9-5p, as a miRNA enriched in HPV+ HNSCC cell-derived exosomes, can be transferred to fibroblasts. miR-9-5p mimic transfection decreased the expression of NOX4 and the level of intracellular reactive oxygen species (ROS), which inhibited the transforming growth factor beta 1(TGF-ß1)-induced increase of αSMA levels. Therefore, these results indicated that HPV+ HNSCC-derived exosomal miR-9-5p inhibits TGF-ß signaling-mediated fibroblast phenotypic transformation through NOX4, which is related to the excellent prognosis of HPV patients.

Characteristics of B lymphocyte infiltration in HPV+ head and neck squamous cell carcinoma.

Human papillomavirus (HPV) is an important etiological factor of head and neck squamous cell carcinoma (HNSCC). HPV+ HNSCC patients usually have a better prognosis, which probably results from the higher infiltration of B lymphocytes. This study was purposed to detect the infiltration of B lymphocyte subsets and the correlation between B lymphocyte subsets and the prognosis in HPV-related HNSCC. In this study, 124 HPV+ and 513 HPV- HNSCC samples were obtained from the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database for transcriptomic analysis. Infiltration of B lymphocytes subsets was detected with 7 HPV+ HNSCC and 13 HPV- HNSCC tissues through immunohistochemistry and immunofluorescence. One HPV- HNSCC sample was detected with single-cell sequencing for chemokine analysis. In the results, the infiltration of plasma cells (CD19+ CD38+ ) and memory B cells (MS4A1+ CD27+ ) was higher in HPV+ HNSCC samples. High infiltration of plasma cells and memory B cells was related to a better prognosis. High density of B lymphocytes was positively correlated with high CXCL13 production mainly from CD4+ T lymphocytes in HNSCC. These results indicated that a high density of plasma cells and memory B cells could predict excellent prognosis. CD4+ T lymphocytes might affect B lymphocytes and their subsets through the CXCL13/CXCR5 axis in HNSCC.

Comprehensive profiling of extracellular RNA in HPV-induced cancers using an improved pipeline for small RNA-seq analysis.

Extracellular RNAs (exRNAs) have attracted great attention due to their essential role in cell-to-cell communication as well as their potential as non-invasive disease biomarkers. However, at present, there is no consensus on the best method to profile exRNA expression, which leads to significant variability across studies. To address this issue, we established an experimental pipeline for comprehensive profiling of small exRNAs isolated from cell culture. By evaluating six RNA extraction protocols, we developed an improved method for robust recovery of vesicle-bound exRNAs. With this method, we performed small RNA sequencing of exosomes (EXOs), microvesicles (MVs) and source cells from 14 cancer cell lines. Compared to cells, EXOs and MVs were similarly enriched in tRNAs and rRNAs, but depleted in snoRNAs. By miRNA profiling analysis, we identified a subset of miRNAs, most noticeably miR-122-5p, that were significantly over-represented in EXOs and MVs across all 14 cell lines. In addition, we also identified a subset of EXO miRNAs associated with cancer type or human papillomavirus (HPV) status, suggesting their potential roles in HPV-induced cancers. In summary, our work has laid a solid foundation for further standardization on exRNA analysis across various cellular systems.

HPV + HNSCC-derived exosomal miR-9 induces macrophage M1 polarization and increases tumor radiosensitivity.

Human papillomavirus (HPV) is an etiological risk factor for a subset of head and neck squamous cell carcinoma (HNSCC). HPV + HNSCC is more radiosensitive than HPV- HNSCC, however, the mechanism underlying this observation remains unknown. Tumor microenvironment can regulate tumor response to radiation therapy. Secretory exosomes are emerging as crosstalk mediators between tumor cells and the tumor microenvironment. In this study, we attempted to determine the role of HPV + HNSCC exosomes in increased radiation sensitivity. We found that HPV + HNSCC exosomes were able to transform macrophages into the M1 phenotype, which subsequently increased the radiosensitivity of HNSCC. miR-9 was found enriched in HPV + HNSCC exosomes and it could be transported into macrophages, inducing M1 macrophage polarization via downregulation of PPARδ. After incubating with M1 macrophages or macrophages treated with miR-9 mimics, HNSCC had strikingly increased radiosensitivity. The clinical significance of miR-9 in HNSCC was confirmed by using profiling data from The Cancer Genome Atlas. Our data suggest that miR-9-enriched exosomes from HPV + HNSCC can polarize macrophages into M1 phenotype and increase the radiosensitivity of HPV + HNSCC. Hence, miR-9 may be used as a potential treatment for HNSCC.

IL-6 induced M1 type macrophage polarization increases radiosensitivity in HPV positive head and neck cancer.

Radiation is a crucial component of head and neck squamous cell carcinoma (HNSC) treatment. Human papillomavirus-positive (HPV+) HNSC is significantly more radiosensitive than HPV- HNSC, but the mechanism underlying this increased sensitivity is unknown. We investigated the possible involvement of macrophage subpopulations as key mediators of HNSC radiosensitivity linked to HPV status. We collected forty-one clinical HNSC specimens and determined HPV status and radiosensitivity of each sample. We investigated cytokine mediated induction of macrophage polarization by HPV+ and HPV- HNSC cells. Radiosensitive HNSC tissues exhibited greater numbers of infiltrating M1 macrophages than radioresistant tumor tissue samples. Moreover, M1 macrophage numbers were positively correlated with HNSC radiosensitivity. HPV+ and HPV- tumor cells induced macrophage polarization to M1 and M2 type, respectively. HPV+ HNSC cells secreted more IL-6 than HPV- cells. HPV promoted tumor cell secretion of IL-6, thereby increasing radiosensitivity through M1 polarization of macrophages. M1 macrophages represent an important tissue microenvironment factor with implications for HNSC treatment efficacy and may prove valuable as a biomarker of radiation sensitivity.

Association of lnc-IL17RA-11 with increased radiation sensitivity and improved prognosis of HPV-positive HNSCC.

Human papillomavirus (HPV) associated head and neck squamous cell carcinoma (HNSCC) has a far better prognosis than HPV negative HNSCC. Recent studies suggest that long noncoding RNA (lncRNA) moieties may play a role in HPV associated differential HNSCC prognoses. In this study, we examined differential expression of lncRNAs in HPV+ vs HPV- HNSCC using The Cancer Genome Atlas database. LncRNAs were categorized based on expression level and survival analysis. A group of eight lncRNAs was identified in which altered expression was associated with both HPV infection and better prognosis. Subsequently, genes coexpressed with these lncRNAs in HNSCC cells were sorted into corresponding co-expression modules. The lnc-IL17RA-11 coexpression module exhibited the greatest correlation with HPV infection and radiotherapy efficacy. We identified the lnc-IL17RA-11 transcription factor ER-alpha as the most likely HPV infection associated factor promoting increased lnc-IL17RA-11 levels. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis revealed enrichment among lnc-IL17RA-11 co-expressed genes for functions related to DNA replication and cell proliferation. These observations are consistent with a model in which HPV infection upregulates transcription factor ER-alpha, which increases levels of lnc-IL17RA-11 and coexpressed genes that promote HNSCC cell sensitivity to radiotherapy, thereby improving disease prognosis.

DisV-HPV16, versatile and powerful software to detect HPV in RNA sequencing data.

BACKGROUND: The increasing availability of high-throughput sequencing data provides researchers with unprecedented opportunities to investigate viral genetic elements in host genomes that contribute to virus-linked cancers. Almost all of the available computational tools for secondary analysis of sequencing data detect viral infection or genome integration events. However, viral oncogenes expression is likely of importance in carcinoma. We therefore developed a new software, DisV-HPV16, for the evaluation of HPV16 oncogenes expression. RESULTS: HPV16 virus and viral oncogenes expression was detected more rapidly using DisV-HPV16 compared to other software. DisV-HPV16 was proved highly convenient for detecting candidate virus after modification of the reference file. The accuracy of DisV-HPV16 was empirically confirmed in laboratory experiments. DisV-HPV16 exhibited greater reliability than other software. CONCLUSIONS: DisV-HPV16 is a new, dependable software to detect virus and viral oncogenes expression through analysis of RNA sequencing data. Use of DisV-HPV16 can yield deeper, more comprehensive insights into virus infection status and viral and host cell gene expression.

Association of p16 as Prognostic Factors for Oropharyngeal Cancer: Evaluation of p16 in 1470 Patients for a 16 Year Study in Northeast China.

Human papillomavirus (HPV) is an etiological risk factor for oropharyngeal squamous cell carcinomas (OPSCC). Our study investigates the prevalence, prognostic, and clinicopathologic features of HPV-related oropharyngeal cancer in Northeast China and elucidates the involvement of p16 in the tumorigenesis and progression of OPSCC. Specimens from 1470 OPSCC patients collected from 2000 to 2016 were analyzed using the status of HPV by polymerase chain reaction (PCR) and p16 immunohistochemistry. Overexpression of p16 was observed in 81 (5.51%) of the 1470 cases, and HPV positive was present in 78 cases (5.31%) of the 1470 cases. HPV positive and p16 overexpression have a good concordance. However, we found that the etiological fraction of HPV in cancers of the OPSCCs was obviously lower in Northeast China than other cohorts previously reported. Interestingly, nearly 89% of patients with p16 expression were smokers, and nearly 70% of patients with p16 expression had a history of alcohol. Our study also demonstrates that p16 expression is significantly associated with early stage primary OPSCCs and the patients with p16 expression tend to show better survival following surgery and radiotherapy.

Prevalence and Prognostic Significance of HPV in Laryngeal Squamous Cell Carcinoma in Northeast China.

BACKGROUND/AIMS: Human papillomavirus (HPV) is an etiological risk factor for a subset of head and neck squamous cell carcinomas. HPV has been proven to be a powerful prognostic biomarker for oropharyngeal cancer, but its role in the larynx has not been explored in depth. Here, we sought to evaluate the prevalence and genotype distribution of HPV in patients with laryngeal squamous cell carcinoma (LSCC) in northeast China. METHODS: HPV DNA in specimens from 211 patients diagnosed with LSCC was analyzed by the polymerase chain reaction and in situ hybridization, and p16 overexpression was evaluated by immunohistochemistry. p16 expression was scored positive if strong and diffuse nuclear and cytoplasmic staining was present in > 75% of tumor cells. RESULTS: In this study, infection with HPV and p16 expression were not absolutely consistent. Among all patients, 132 (62.6%) were positive for HPV DNA (HPV+), while 23 (10.9%) were inconsistent for HPV and p16. Multivariate analysis indicated that HPV, but not p16, is an independent prognostic factor for overall survival in LSCC. Overall survival was significantly improved in HPV+ LSCC patients compared with the HPV-negative group (hazard ratio, 0.395; 95% confidence interval, 0.185-0.843; p = 0.016). Among the 132 HPV+ patients, 28 (21.2%) were HPV-16 single infection. CONCLUSION: This study indicates that HPV DNA is a more reliable surrogate marker than p16 for the prediction of survival in patients with LSCC.

Immunological network analysis in HPV associated head and neck squamous cancer and implications for disease prognosis.

Human papillomavirus-positive (HPV+) head and neck squamous cell cancer (HNSCC) exhibits a better prognosis than HPV-negative (HPV-) HNSCC. This difference may in part be due to enhanced immune activation in the HPV+ HNSCC tumor microenvironment. To characterize differences in immune activation between HPV+ and HPV- HNSCC tumors, we identified and annotated differentially expressed genes based upon mRNA expression data from The Cancer Genome Atlas (TCGA). Immune network between immune cells and cytokines was constructed by using single sample Gene Set Enrichment Analysis and conditional mutual information. Multivariate Cox regression analysis was used to determine the prognostic value of immune microenvironment characterization. A total of 1673 differentially expressed genes were functionally annotated. We found that genes upregulated in HPV+ HNSCC are enriched in immune-associated processes. And the up-regulated gene sets were validated by Gene Set Enrichment Analysis. The microenvironment of HPV+ HNSCC exhibited greater numbers of infiltrating B and T cells and fewer neutrophils than HPV- HNSCC. These findings were validated by two independent datasets in the Gene Expression Omnibus (GEO) database. Further analyses of T cell subtypes revealed that cytotoxic T cell subtypes predominated in HPV+ HNSCC. In addition, the ratio of M1/M2 macrophages was much higher in HPV+ HNSCC. The infiltration of these immune cells was correlated with differentially expressed cytokine-associated genes. Enhanced infiltration of B cells and CD8+ T cells were identified as independent protective factors, while high neutrophil infiltration was a risk enhancing factor for HPV+ HNSCC patients. A schematic model of immunological network was established for HPV+ HNSCC to summarize our findings.

miR-34a and its novel target, NLRC5, are associated with HPV16 persistence.

Persistent infection with human papillomavirus (HPV), particularly type 16, is causally associated with cervical cancer and its precursors. The role of miRNAs in HPV16 persistence currently remains unclear. Preliminary analysis of miRNA profile demonstrated that HPV16 infection caused a striking downregulation of miR-34a. Through bioinformatics analysis and dual-luciferase assay with site-directed mutagenesis strategy, NLRC5, a negative regulator of NF-κB signaling, was identified to be a novel interactor of miR-34a. Transfection of miR-34a mimic strikingly downregulated NLRC5 in the HPV16-positive cervical cells, which might result in the nuclear accumulation of NF-κB p65. However, transfection of miR-34a inhibitor exhibited an opposite effect. The antagonistic expressions of NLRC5 and miR-34a were also observed in keratinocytes harboring HPV16 genome as well as in human cervical samples with persistent infection of HPV16. Our data uncover a previously unknown connection among HPV16 persistence, miR-34a and its interactor NLRC5.