Abstract
Keywords
Introduction
Sinonasal tumours are rare malignancies with an annual worldwide incidence of approximately 1 in 100000. 1 Histologically, sinonasal tumours are divided into sarcomas, olfactory neuroblastomas, mucosal melanomas and carcinomas. Over 55% of sinonasal tumours are carcinomas, which include squamous cell carcinomas (SCC), adenocarcinomas and adenoid cystic carcinoma.2,3 Inverted papillomas (IP) are common benign tumours of the nasal cavity and paranasal sinuses. The malignant transformation rate of IPs is less than 27%, 4 but the histomorphology, immunohistochemistry and biological markers involved in their transformation remain unknown. 5 It is therefore important to explore the molecular mechanisms underlying the malignant transformation of IPs, which could help identify biomarkers for improved diagnosis of sinonasal tumours and better targeted therapy.
Epigenetic regulation, via DNA methylation and microRNA (miRNA) regulation, plays a critical role in a variety of cancers including head and neck SCC.6,7 Hypomethylation and hypermethylation are thought to contribute to carcinogenesis by overexpressing tumour oncogenes and downregulating tumour suppressor genes.
8
A global methylation study revealed that aberrant hypermethylation of gene promoters is involved in the pathogenesis of sinonasal papilloma.
9
A single gene methylation study showed that low
Materials and methods
Study design
This retrospective tissue analysis study compared DNA methylation and gene expression between tissue samples of SNIP and SCC arising in SNIP samples. A comprehensive methylation profiling technique was used to investigate DNA methylation in C-phosphate-G (CpG) islands and gene promoter regions in tissue samples of SNIP and samples of SCC arising in SNIPs. Gene ontology analysis and network analysis were used to analyse gene function and to construct gene networks of aberrantly methylated genes. Expression of these aberrantly methylated genes was further validated using real-time polymerase chain reaction (PCR) or Western blot analysis in tissue samples of SNIP and samples of SCC arising in SNIPs. All of the SCCs arising in SNIP samples showed signs of malignant transformation of the SNIP based on histological analysis.
Study participants
Tissue samples were collected from surgical patients at Beijing Tongren Hospital, Capital Medical University, Beijing, China between January 2013 and December 2016. None of the patients had received chemotherapy, radiotherapy or biological therapy before sample collection. None of the patients had other types of cancer. Cancer was staged according to the 7th edition of the American Joint Committee on Cancer Staging Manual. 11 All of the SCCs were arising in SNIPs and showed malignant progression. All tissues were snap-frozen in liquid nitrogen and stored at –80°C until further use.
Written informed consent was obtained from all patients involved in this study, which was approved by the Ethics Committee of the Beijing Tongren Hospital, Capital Medical University, Beijing, China (no. TRECKY2014-027).
DNA extraction, microarray analysis and analysis of aberrantly methylated genes
Genomic DNA from each cryopreserved tissue was extracted according to the user manual of the QIAamp DNA Mini kit (Qiagen, Valencia, CA, USA). An EZ DNA Methylation™ Kit (ZymoResearch, Orange, CA, USA) was used to bisulfite the purified DNA. The methylation levels of over 485 000 individual cytosines were measured using an Infinium HumanMethylation450K BeadChip kit (Illumina, San Diego, CA, USA). The Infinium HumanMethylation450K BeadChip assay was performed according to the manufacturer’s instructions using an iScan™ SQ System (Illumina). The methylation level for each locus was designated as the β value, which was a quantitative measure of DNA methylation ranging from 0 (no cytosine methylation) to 1 (complete cytosine methylation) using Illumina GenomeStudio V2011.1 Methylation Module 1.9.0 (Illumina). Quality control was conducted by checking and removing the poorly performing probes or samples with the following exclusion criteria: (i) probes with a detection
Gene ontology analysis
The microarray data were subjected to gene ontology (GO) analysis as previously described.
12
The over-represented GO terms were tested and the overlapping probabilities of differentially methylated region datasets were calculated as previously described.13,14 The cellular components, molecular function and biological processes of differentially methylated genes were analysed.
Gene network analysis
Qiagen’s ingenuity pathway analysis was used to identify the gene networks. The Kyoto Encyclopedia of Genes and Genomes database was used to build the networks of genes. 15 Fisher’s exact test was used to select significant pathways.
Real-time PCR of miR-661 mRNA
Total RNA was extracted using TRIzol® reagent (Invitrogen, Carlsbad, CA, USA) and reverse transcription was performed using a One Step PrimeScript® miRNA cDNA Synthesis Kit (Takara Bio, Dalian, China). Real-time quantitative PCR was performed using SYBR® Premix Ex Taq™ II (Takara Bio). 16 The forward primers used for microRNA 661 (miR-661) mRNA amplification were the sequences of mature miR-661 and the reverse primers were provided by the kit.
Western blot analysis of OPA3 and PLEC
Protein extracts were prepared from tissue samples using ice-cold lysis buffer (50 mM Tris-HCl, pH 6.8, 32 mM 2-mercaptoethanol, 2% w/v sodium dodecyl sulphate [SDS], 10% glycerol) supplemented with ethylenediaminetetra-acetic acid-free complete protease inhibitors (Roche, Penzberg, Germany). Following lysis, samples were centrifuged at 12830
Statistical analyses
All statistical analyses were performed using the SPSS® statistical package, version 17.0 (SPSS Inc., Chicago, IL, USA) for Windows®. The data on gene expression are presented as mean ± SE. Statistical significance between groups were evaluated using Student's
Results
This study analysed 27 tissue samples that were collected from surgical patients: 15 SNIP samples and 12 SCCs arising in SNIPs. The 15 SNIPs were collected from nine males and six females (age range, 35–61 years). The 12 SCCs were collected from eight males and four females (age range, 39–65 years).
The methylated DNAs from six SNIP samples and seven SCC samples were profiled using the Infinium HumanMethylation450K BeadChip array. During the process, two SCC samples failed during analysis. Thus, six SNIP and five SCC samples were included in the final analysis. The differentially methylated sites were determined by the DiffScore value. A total of 42 9607 methylation sites were measured and 11 201 methylation sites showed 2-fold differences between SNIP and SCC tissues after log transformation (r2 = 0.9606, r
2
sel = 0.6311). Among the 11 201 differentially methylated sites, six sites were significantly different between the SNIP and SCC samples at
Six differentially methylated sites identified in six samples of sinonasal inverted papilloma (SNIP) and five samples of sinonasal squamous cell carcinoma arising in SNIPs.
A differential score equal to 20, 30 or 40 corresponds to a significant
UCSC, University of California, Santa Cruz;
The expression levels of the aberrantly methylated genes were further validated using either real-time PCR to measure miR-661 mRNA levels or Western blot analysis to measure OPA3 and PLEC protein levels in nine SNIP samples and five SCC arising in SNIP samples. Real-time PCR confirmed that miR-661 mRNA levels were significantly higher in SCCs compared with SNIPs (

Comparison of the expression levels of the aberrantly methylated genes as determined using either real-time polymerase chain reaction (PCR) or Western blot analysis in nine samples of sinonasal inverted papilloma (SNIP) and five samples of sinonasal squamous cell carcinomas (SCC) arising in SNIPs. (a) Real-time PCR of microRNA 661 (miR-661) mRNA levels in six SNIP samples and seven SCC samples. (b) Western blot analysis of outer mitochondrial membrane lipid metabolism regulator (OPA3) and plectin (PLEC) protein levels. β-actin was used as the internal control. (c) Relative OPA3 protein levels in Western blot analysis. (d) Relative PLEC protein levels in Western blot analysis. Data presented as mean ± SE; *
Discussion
Even though the malignant transformation of SNIP into SCC is documented in the literature, the molecular mechanisms responsible for this transformation have not been elucidated. Aberrant methylation is an important mechanism that regulates the expression of tumour suppressive genes or oncogenes and is subsequently involved in the carcinogenesis and development of cancers.
7
The role of methylation in the malignant transformation of SNIP into SCC has not been currently addressed. In this current study, global DNA methylation in CpG islands and the promoter regions of genes was compared between SNIPs and SCCs arising in SNIPs using a comprehensive methylation profiling technique. This current study found six methylation sites that were significantly different between the two groups at
Previous studies have shown that miR-661 functions as a tumour suppressor in malignancy and acts as a tumour suppressor in breast cancer where it inhibits cell proliferation, cell motility and cell invasion.
17
There are some conflicting findings in breast cancer studies.17,18 A recent study showed that low miR-661 expression was correlated with a poor outcome in patients with breast cancers expressing wild-type p53, whereas high miR-661expression promoted invasion of tumour cells harbouring p53 mutations.
19
Thus, miR-661 may either suppress or promote cancer aggressiveness in the same type of tumour depending upon the status of p53 expression.
19
In glioma cells, miR-661 levels were downregulated and inhibited cancer cell proliferation, migration and invasion.
20
In this study, miR-661 could silence human telomerase reverse transcriptase (hTERT) by recognizing and specifically binding to the predicted site of the hTERT mRNA 3' untranslated region.
20
Another study indicated that hsa-circ-0012129 might act as a natural miR-661 sponge and expression of circ-0012129 could be suppressed by miR-661.
21
Several recent studies have shown that miR-661 may participate in the regulation of occurrence and progression in non-small cell lung cancer by directly targeting runt related transcription factor 3
22
or RB transcriptional corepressor 1
23
and interacting with adenocarcinoma predictive long intergenic non-coding RNA.
24
The expression of the
The OPA3 protein is an integral component of the mitochondrial outer membrane and mutations in the
Plectin is a giant multifunctional cytokine protein that helps stabilize and orchestrate the intermediate filament network in cells.
27
Mutations in the human
Recently, gene body methylation (GbM) was found to frequently occur in the transcribed regions of many oncogenic genes and to be actively involved in multiple regulatory processes.32,33 More detailed genome-wide studies have demonstrated that GbM may alter gene expression by silencing alternative promoters, affecting transcription elongation and regulating splicing.34–36 Therefore, GbM could be serve as a novel biomarker or therapeutic target for human cancers. 37 In the current study, all hypermethylation occurred in the gene body and the expression of these three genes were different between SNIPs and SCCs, indicating that GbM might play a role in the malignant transformation of SNIP.
This current study had a number of limitations. First, the small sample size may lead to insufficient statistical power. The study only identified three genes with abnormal methylation at
In conclusion, this current study identified three methylated genes that were differently expressed between SNIPs and SCCs arising in SNIPs. Although the roles of these three genes have not been previously shown to participate in tumorigenesis and progression, this study is the first to highlight their potential involvement in the malignant transformation of SNIP. This current study also demonstrates the ability of genome-wide epigenetic studies to identify potential biomarkers and thus, provide new insights into their involvement in tumorigenesis. Future research will focus on the association of the methylation landscape with clinical outcome of patients with SNIP.
