Novel age-dependent targets in vestibular schwannomas
© Toren et al.; licensee BioMed Central Ltd. 2014
Received: 7 May 2014
Accepted: 15 June 2014
Published: 30 June 2014
Schwannomas are the most common neurofibromatosis type 2 (NF2)-associated tumors with significant phenotypic heterogeneity in patients. The most severe subtype has an early and rapid progression and the mild type has a later onset and a less aggressive course. The aim of this study was to elucidate the underlying molecular differences between these groups. We compared the gene expression pattern between patients with early to late age of onset.
A gene signature of 21 genes was constructed to differentiate between early-onset and late-onset patients. We confirmed these results by real-time PCR for SNF1LK2, NGFRAP1L1 (BEX 5), GMNN, and EPHA2.
Genes identified here may be additional aberrations in merlin-depleted cells that govern the disease onset. A significant number of these genes have been suggested as having a role in carcinogenesis and are used as biomarkers for prognosis in several other cancers. The role of these genes in NF2 carcinogenesis and their potential as biomarkers or drug target are worthwhile exploring.
KeywordsNeurofibromatosis Schwannomas Targeted therapy
Neurofibromatosis type 2 (NF2) is an autosomal dominant disorder that predisposes affected individuals to several types of cerebral and spinal tumors . Epidemiological studies place the incidence of NF2 at 1 in 25,000 live births . The current treatment of NF2 is surgery and/or radiation therapy. However, there is no current safe therapeutic agent available at this time. Bilateral vestibular schwannomas (BVS) are the most common of NF2-associated tumors . The consequences of a vestibular schwannomas (VS) are numerous, including dizziness, imbalance, tinnitus, hearing loss progressing to deafness, facial nerve paralysis, brainstem compression, and, if left untreated, death.
Mutations of the NF2 gene are thought to be the primary cause of this disease. Various types of mutations have been identified, among them single-base substitutions, insertions, and deletions . The NF2 gene product, merlin, functions as a signal transduction pathway regulating cell-to-cell and cell matrix interactions. Merlin has been demonstrated to exert its antiproliferative effect by impact on several intracellular signaling pathways: inhibition of phosphatidylinositol 3-kinase (PI3-kinase), through binding to PI3-kinase enhancer, long form (PIKE-L) , modulation of the RAS/RAC oncogenic pathways , inhibition of the mitogen-activated protein kinase (MAPK) pathways  and repression of PAK-induced cyclin D1 expression . Merlin regulates the MAPK signaling pathway by blocking the Cdc42-MLK3 interaction and MLK3 activity inhibition . Furthermore, inactivation of merlin induces oncogenic gene expression, hyperproliferation, and tumorigenicity by releasing the inhibition of E3 ubiquitin ligase CRL4 (DCAF1) in the nucleus .
Since the disease has a variable presentation, with the severe subtype having an early and rapid progression and the milder type having a later onset and a less aggressive course, the diagnosis and management of these patients present a unique challenge . Studies have indicated that a truncating NF2 mutation (nonsense and frameshift) is linked with the more severe form of NF2 . The more severe form is particularly virulent, with unrelenting growth of schwannomas and meningiomas from childhood, resulting in blindness, deafness, paralysis, and death by the age of 40. The mild form of NF2 is less debilitating. The schwannomas may remain the same size for years, few meningiomas will develop. The patient may not develop symptoms until later in life and will have fewer disabilities .
There is considerable variability in growth rates of VS in patients with NF2, but they tend to be higher in patients who are younger at onset of signs or symptoms . The most important predictors of disease severity are age at diagnosis and age at onset of symptoms. The average age of onset is 18 to 24 years. Almost all affected individuals develop bilateral vestibular schwannomas by age 30 years. Vestibular schwannomas are much less common in children than adults . There are age-specific risks of NF2 to an offspring of a patient. Patient presenting with BVS <20 years 29.3% and patient presenting with asymmetric disease after 40 years 5.5% . Patients of ages 35 years and below show higher p53 phosphorylation compared to the tumors of older patients .
The development of microarray tools assists effective, current, and comprehensive understanding of complex diseases including cancer . Differences in gene expression between schwannomas or between schwannoma cells from NF2 patients and normal human primary schwann cells previously examined do not confer age impact [17–19]. Here, we report whole genome gene expression variation in NF tumors differing in age of onset. Improved understanding of these molecular differences may allow elucidation of age-dependent onset etiology, novel molecular markers for prognosis, and new therapeutic targets for neurofibromatosis.
Results and discussion
List of age-dependent differentially expressed genes
Role in cancer
Overexpressed in pancreatic carcinoma
Neuron growth, regeneration
Zinc finger protein
Associated with prostate cancer
Protein sorting and membrane trafficking
Overexpressed in multiple cancers
Nuclear protein transport export
Correlated with poorer clinical outcome
Regulation of Wnt pathway
Expression increases with ascending glioma proliferation activity
Overexpressed in many types of cancer
Correlated with ovarian cancer poor survival
Nerve growth factor receptor-associated
May be a tumor suppressor in VS
Suppress canonical Wnt signaling
ELL translocations in malignancies.
Expressed abnormally in many cancer types
Glucose transport inhibitor
Suppresses cell proliferation
Lost in breast and prostate cancer
Neuronal cell death
Proposed applicable uses of differentially expressed genes
To be defined
Ten of the 21 genes we identified fall into four broad functional categories (Table 1): protein transport (genes PSCD3, XPO1, and SYTL3), protein phosphorylation and dephosphorylation (XPO1, EPHA2, and SNF1LK2), cell growth (NAV3, GMNN, and NGFRAP1L1), and cell death (ZNF331 and RNF19A). The latter two are particularly relevant to cancerous growth and the first two may underpin it as well. Additional cellular functions and cancer phenotypes are also listed in Table 1. Furthermore, some, at least 14 (EPA, XPO, SNF1LK, NAV3, KCNK17, ZNF331, GNMM, ZNF718, WHSC1, ANLN, PSCD3, RGN, ELL2, and PTPN1) of the 21 genes, have been reported to be involved in cancers (Table 1).
We further grouped the 21 genes we identified into 4 subgroups according to their potential role in NF2: in carcinogenesis, as biomarker, as drug target, or potential for clinical outcome prediction. This was according to other studies findings for each gene above. Another group was constructed from genes that their role in cancer has not been indicated yet (Table 2).
In summary, we report here a gene signature comprising 21 genes that characterize the clinically important age at onset in NF2 patients. Interestingly, many of these genes have known or at least inferred functions which may lead to a more complete molecular understanding of the age-based pathogenic differences in NF2 in the future. Such data may, in turn, lead to improved pre-symptomatic diagnoses and personalized age-focused treatment.
There is great phenotypic difference between early and late onset in NF2. Our data indicate that these groups have distinct gene expression signatures. Most of the differentiating genes have already been shown to have important roles in other cancers. Thus, they may also have a role in VS NF2 patients. Moreover, these genes may be involved in pathogenesis, used as potential age-oriented biomarkers and treatment targets. Our gene lists will provide fertile ground for testing novel hypotheses on disease mechanism, pathogenesis, and therapeutic targets.
Thirty-seven samples from human tumor (schwannomas) tissues were obtained at the time of surgery in accordance with applicable human ethics regulations (USC IRB protocol # 97–157) and snap frozen in liquid nitrogen. Tumor specimens were divided into two groups (young and old) according to the patient's age of onset. The young-onset group included 21 specimens from patients whose age of onset was less than 16 years old (ranging from 7 to 16 years old). The old-onset group included 16 specimens from patients whose age of onset was 30 years and older (ranging from 30 to 67 years old). The average age of the young group was 12 years old. The average age of the old group was 43 years. RNA was extracted and cRNA were hybridized to Human Genome U133 Plus 2.0 Array (Affymetrix) for analysis of over 47,000 transcripts. The detailed protocol for the sample preparation and microarray processing is available from Affymetrix at http://www.affymetrix.com (Santa Clara, CA, USA). For real-time PCR, total RNA was extracted and RT-qPCR was performed in the DNA Engine Opticon 2 Continuous Fluorescence Detection System (Bio-Rad, Hercules, CA, USA).
The Bioconductor limma package  was used to preprocess the expression data from Human Genome U133 Plus 2.0 Affymetrix Array. To remove batch effects, the complete dataset were normalized by RMA. Before doing statistical analysis, we filter the probe sets by present/absent calls using the Wilcoxon signed rank-based algorithm. Out of 54,675 probes, 42,339 were kept and the remaining ones were removed as ‘Absent’ to reduce false positive rate. We identified differentially expressed gene list by fitting linear models to the normalized expression values. The empirical Bayes shrinkage was applied to t statistics using the limma package in Bioconductor. Because patient subtypes were still observed after RMA normalization, balanced samples from the same subtype were selected to identify significantly differentially expressed genes and to reduce the potential of small-sample artifacts. In general, selected genes have greater than 1.5 fold change and less than 0.01 p value.
The analysis of RT-PCR is as follows: Sample RNA was processed in triplicate with serial dilutions of Human Universal Reference RNA (Stratagene, La Jolla, CA, USA) on the same 96-well plate. Standard curves were constructed from the Universal RNA wells with arbitrary units of 1 unit equivalent to 1 pg of Universal RNA. Results were plotted as units per nanogram of RNA. Graphically, the middle line of the bar represents mean, bar size corresponds to the difference between lower quartile and upper quartile, and error bars represent the range of the raw data. t test p values were calculated for each differentiating gene.
Ethical conduct of research
The authors state that they have obtained the University of Southern California institutional review board approval and have followed the principles outlined in the Declaration of Helsinki for all human experimental investigations. In addition, informed consent has been obtained from the participants involved.
- Gutmann DH, Giordano MJ, Fishback AS, Guha A: Loss of merlin expression in sporadic meningiomas, ependymomas and schwannomas. Neurology. 1997, 49 (1): 267-270. 10.1212/WNL.49.1.267.View ArticlePubMedGoogle Scholar
- Asthagiri AR, Parry DM, Butman JA, Kim HJ, Tsilou ET, Zhuang Z, Lonser RR: Neurofibromatosis type 2. Lancet. 2009, 373 (9679): 1974-1986. 10.1016/S0140-6736(09)60259-2.View ArticlePubMedGoogle Scholar
- Evans DG, Trueman L, Wallace A, Collins S, Strachan T: Genotype/phenotype correlations in type 2 neurofibromatosis (NF2): evidence for more severe disease associated with truncating mutations. J Med Genet. 1998, 35 (6): 450-455. 10.1136/jmg.35.6.450.PubMed CentralView ArticlePubMedGoogle Scholar
- Sughrue ME, Yeung AH, Rutkowski MJ, Cheung SW, Parsa AT: Molecular biology of familial and sporadic vestibular schwannomas: implications for novel therapeutics. J Neurosurg. 2009, 114 (2): 359-366.View ArticlePubMedGoogle Scholar
- Rong R, Tang X, Gutmann DH, Ye K: Neurofibromatosis 2 (NF2) tumor suppressor merlin inhibits phosphatidylinositol 3-kinase through binding to PIKE-L. Proc Natl Acad Sci USA. 2004, 101 (52): 18200-18205. 10.1073/pnas.0405971102.PubMed CentralView ArticlePubMedGoogle Scholar
- Yohay KH: The genetic and molecular pathogenesis of NF1 and NF2. Semin Pediatr Neurol. 2006, 13 (1): 21-26. 10.1016/j.spen.2006.01.007.View ArticlePubMedGoogle Scholar
- Lim JY, Kim H, Jeun SS, Kang SG, Lee KJ: Merlin inhibits growth hormone-regulated Raf-ERKs pathways by binding to Grb2 protein. Biochem Biophys Res Commun. 2006, 340 (4): 1151-1157. 10.1016/j.bbrc.2005.12.122.View ArticlePubMedGoogle Scholar
- Xiao GH, Gallagher R, Shetler J, Skele K, Altomare DA, Pestell RG, Jhanwar S, Testa JR: The NF2 tumor suppressor gene product, merlin, inhibits cell proliferation and cell cycle progression by repressing cyclin D1 expression. Mol Cell Biol. 2005, 25 (6): 2384-2394. 10.1128/MCB.25.6.2384-2394.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhan Y, Chadee DN: Inhibition of Cdc42-mediated activation of mixed lineage kinase 3 by the tumor suppressor protein merlin. Small Gtpases. 2010, 1 (3): 183-186. 10.4161/sgtp.1.3.14935.PubMed CentralView ArticlePubMedGoogle Scholar
- Li W, Giancotti FG: Merlin's tumor suppression linked to inhibition of the E3 ubiquitin ligase CRL4 (DCAF1). Cell Cycle. 2010, 9 (22): 4433-4436. 10.4161/cc.9.22.13838.PubMed CentralView ArticlePubMedGoogle Scholar
- Uppal S, Coatesworth AP: Neurofibromatosis type 2. Int J Clin Pract. 2003, 57 (8): 698-703.PubMedGoogle Scholar
- Mautner VF, Baser ME, Thakkar SD, Feigen UM, Friedman JM, Kluwe L: Vestibular schwannoma growth in patients with neurofibromatosis type 2: a longitudinal study. J Neurosurg. 2002, 96 (2): 223-228. 10.3171/jns.2002.96.2.0223.View ArticlePubMedGoogle Scholar
- Evans DG, Huson SM, Donnai D, Neary W, Blair V, Newton V, Harris R: A clinical study of type 2 neurofibromatosis. Q J Med. 1992, 84 (304): 603-618.PubMedGoogle Scholar
- Evans DG, Wallace A: An update on age related mosaic and offspring risk in neurofibromatosis 2 (NF2). J Med Genet. 2009, 46 (11): 792-10.1136/jmg.2009.070342.View ArticlePubMedGoogle Scholar
- Dayalan AH, Jothi M, Keshava R, Thomas R, Gope ML, Doddaballapur SK, Somanna S, Praharaj SS, Ashwathnarayanarao CB, Gope R: Age dependent phosphorylation and deregulation of p53 in human vestibular schwannomas. Mol Carcinog. 2006, 45 (1): 38-46. 10.1002/mc.20150.View ArticlePubMedGoogle Scholar
- Shai RM: Microarray tools for deciphering complex diseases. Front Biosci. 2006, 11: 1414-1424. 10.2741/1892.View ArticlePubMedGoogle Scholar
- Hanemann CO, Bartelt-Kirbach B, Diebold R, Kampchen K, Langmesser S, Utermark T: Differential gene expression between human schwannoma and control Schwann cells. Neuropathol Appl Neurobiol. 2006, 32 (6): 605-614. 10.1111/j.1365-2990.2006.00769.x.View ArticlePubMedGoogle Scholar
- Aarhus M, Bruland O, Saetran HA, Mork SJ, Lund-Johansen M, Knappskog PM: Global gene expression profiling and tissue microarray reveal novel candidate genes and down-regulation of the tumor suppressor gene CAV1 in sporadic vestibular schwannomas. Neurosurgery. 2010, 67 (4): 998-1019. 10.1227/NEU.0b013e3181ec7b71. discussion 1019View ArticlePubMedGoogle Scholar
- Martinez-Glez V, Franco-Hernandez C, Alvarez L, De Campos JM, Isla A, Vaquero J, Lassaletta L, Casartelli C, Rey JA: Meningiomas and schwannomas: molecular subgroup classification found by expression arrays. Int J Oncol. 2009, 34 (2): 493-504.PubMedGoogle Scholar
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3 (1): ISSN (Online) 1544-6115, doi:10.2202/1544-6115.1027. Article3Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.