Comparison of the miRNA profiles in HPV-positive and HPV-negative tonsillar tumors and a model system of human keratinocyte clones

Background Better insights into the molecular changes involved in virus-associated and -independent head and neck cancer may advance our knowledge of HNC carcinogenesis and identify critical disease biomarkers. Here we aimed to characterize the expression profiles in a matched set of well-characterized HPV-dependent and HPV-independent tonsillar tumors and equivalent immortalized keratinocyte clones to define potential and clinically relevant biomarkers of HNC of different etiology. Methods Fresh frozen tonsillar cancer tissues were analyzed together with non-malignant tonsillar tissues and compared with cervical tumors and normal cervical tissues. Furthermore, relative miRNAs abundance levels of primary and immortalized human keratinocyte clones were evaluated. The global quantitation of miRNA gene abundance was performed using a TaqMan Low Density Array system. The confirmation of differentially expressed miRNAs was performed on a set of formalin-fixed paraffin-embedded tumor samples enriched for the tumor cell fraction by macrodissection. Results We defined 46 upregulated and 31 downregulated miRNAs characteristic for the HPV-positive tonsillar tumors and 42 upregulated miRNAs and 42 downregulated miRNAs characteristic for HPV-independent tumors. In comparison with the expression profiles in cervical tumors, we defined miR-141-3p, miR-15b-5p, miR-200a-3p, miR-302c-3p, and miR-9-5p as specific for HPV induced malignancies. MiR-335-5p, miR-579-3p, and miR-126-5p were shared by the expression profiles of HPV-positive tonsillar tumors and of the HPV immortalized keratinocyte clones, whereas miR-328-3p, miR-34c-3p, and miR-885-5p were shared by the miRNA profiles of HPV-negative tonsillar tumors and the HPV-negative keratinocytes. Conclusions We identified the miRNAs characteristic for HPV-induced tumors and tonsillar tumors of different etiology, and the results were compared with those of the model system. Our report presents the basis for further investigations leading to the identification of clinically relevant diagnostic and/or therapeutic biomarkers for tumors of viral and non-viral etiology

The Role of miRNAs in Virus-Mediated Oncogenesis

To date, viruses are reported to be responsible for more than 15% of all tumors worldwide. The oncogenesis could be influenced directly by the activity of viral oncoproteins or by the chronic infection or inflammation. The group of human oncoviruses includes Epstein–Barr virus (EBV), human papillomavirus (HPV), hepatitis B virus (HBV), hepatitis C virus (HCV), human herpesvirus 8 (HHV-8) or polyomaviruses, and transregulating retroviruses such as HIV or HTLV-1. Most of these viruses express short noncoding RNAs called miRNAs to regulate their own gene expression or to influence host gene expression and thus contribute to the carcinogenic processes. In this review, we will focus on oncogenic viruses and summarize the role of both types of miRNAs, viral as well as host’s, in the oncogenesis

Comparison of the miRNA expression profiles in fresh frozen and formalin-fixed paraffin-embedded tonsillar tumors.

MicroRNAs are considered as promising prognostic and diagnostic biomarkers of human cancer since their profiles differ between tumor types. Most of the tumor profiling studies were performed on rarely available fresh frozen (FF) samples. Alternatively, archived formalin-fixed paraffin-embedded (FFPE) tissue samples are also well applicable to larger-scale retrospective miRNA profiling studies. The aim of this study was to perform systematic comparison of the miRNA expression profiles between FF and macrodissected FFPE tonsillar tumors using the TaqMan Low Density Array system, with the data processed by different software programs and two types of normalization methods. We observed a marked correlation between the miRNA expression profiles of paired FF and FFPE samples; however, only 27-38% of the differentially deregulated miRNAs overlapped between the two source systems. The comparison of the results with regard to the distinct modes of data normalization revealed an overlap in 58-67% of differentially expressed miRNAs, with no influence of the choice of software platform. Our study highlights the fact that for an accurate comparison of the miRNA expression profiles from published studies, it is important to use the same type of clinical material and to test and select the best-performing normalization method for data analysis

Correlation of paired samples visualized on scatterplots.

<p>Scatter plots of the normalized RQ values in pair-wise comparison from fresh-frozen (FF) and formalin-fixed paraffin embedded (FFPE) samples.</p

Comparison on fold-change for selected miRNAs between arrays and qPCR.

<p>The treshold of FC was set to 1.33. T-test or Mann Whitney nonparametric test was applied depending on the data distribution. All results were statistically significant (P-value ≤0.05).</p

Comparison of electropherogram outputs.

<p>At the top—ladder, in the middle—fresh frozen sample, at the bottom—formalin-fixed paraffin-embedded sample. Right—lines from virtual gel.</p

Genome-wide miRNA profiling reinforces the importance of miR-9 in human papillomavirus associated oral and oropharyngeal head and neck cancer

Head and neck cancer is the sixth most common malignancy worldwide, predominantly developing from squamous cell epithelia (HNSCC). The main HNSCC risk factors are tobacco, excessive alcohol use, and the presence of human papillomavirus (HPV). HPV positive (+) cancers are etiologically different from other HNSCC and often show better prognosis. The current knowledge regarding HNSCC miRNA profiles is still incomplete especially in the context of HPV+ cancer. Thus, we analyzed 61 freshly collected primary oral (OSCC) and oropharyngeal (OPSCC) SCC samples. HPV DNA and RNA was found in 21% cases. The Illumina whole-genome small-RNA profiling by next-generation sequencing was done on 22 samples and revealed 7 specific miRNAs to HPV+ OSCC, 77 to HPV+ OPSCC, and additional 3 shared with both; 51 miRNAs were specific to HPV− OPSCC, 62 to HPV− OSCC, and 31 shared with both. The results for 9 miRNAs (miR-9, -21, -29a, -100, -106b, -143 and -145) were assessed by reverse transcription-quantitative polymerase chain reaction on the whole study population. The data was additionally confirmed by reanalyzing publicly available miRNA sequencing Cancer Genome Atlas consortium (TCGA) HNSCC data. Cell signaling pathway analysis revealed differences between HPV+ and HPV− HNSCC. Our findings compared with literature data revealed extensive heterogeneity of miRNA deregulation with only several miRNAs consistently affected, and miR-9 being the most likely HPV related miRNA

Differentially expressed miRNAs detected in tonsillar tumors versus non-malignant tonsillar tissue.

<p>The overlap between FF vs FFPE material is highlighted in green. On the left—results using 50th percentile shift normalization for data analysis. On the right—results using global normalization for data analysis. FC—fold change.</p

Differentially expressed miRNA in both types of clinical material by using two different normalization methods.

<p>Differentially expressed miRNA in both types of clinical material by using two different normalization methods.</p