Distinct immune signature predicts progression of vestibular schwannoma and unveils a possible viral etiology
Journal of Experimental & Clinical Cancer Research volume 41, Article number: 292 (2022)
The management of sub-totally resected sporadic vestibular schwannoma (VS) may include observation, re-resection or irradiation. Identifying the optimal choice can be difficult due to the disease’s variable progression rate.
We aimed to define an immune signature and associated transcriptomic fingerprint characteristic of rapidly-progressing VS to elucidate the underpinnings of rapidly progressing VS and identify a prognostic model for determining rate of progression.
We used multiplex immunofluorescence to characterize the immune microenvironment in 17 patients with sporadic VS treated with subtotal surgical resection alone. Transcriptomic analysis revealed differentially-expressed genes and dysregulated pathways when comparing rapidly-progressing VS to slowly or non-progressing VS.
Rapidly progressing VS was distinctly enriched in CD4+, CD8+, CD20+, and CD68+ immune cells. RNA data indicated the upregulation of anti-viral innate immune response and T-cell senescence. K − Top Scoring Pair analysis identified 6 pairs of immunosenescence-related genes (CD38-KDR, CD22-STAT5A, APCS-CXCR6, MADCAM1-MPL, IL6-NFATC3, and CXCL2-TLR6) that had high sensitivity (100%) and specificity (78%) for identifying rapid VS progression.
Rapid progression of residual vestibular schwannoma following subtotal surgical resection has an underlying immune etiology that may be virally originating; and despite an abundant adaptive immune response, T-cell immunosenescence may be associated with rapid progression of VS. These findings provide a rationale for clinical trials evaluating immunotherapy in patients with rapidly progressing VS.
Vestibular schwannomas (VS) are benign tumors arising from the vestibular division of cranial nerve VIII. The incidence of VS is 1.1 per 100,000 person-years in the United States, showing no gender preponderance [1, 2]. Advances in neuroimaging, particularly magnetic resonance imaging (MRI), have led to increasingly frequent diagnoses of smaller VS, thereby contributing to the increased incidence of VS reported in recent years . Treatment options for VS include watchful waiting with serial MRI, surgical resection or radiotherapy. Given that VS tumors grow at a variable and unpredictable rate, choosing among treatment options remains difficult .
Surgical refinements and surgeon experience have minimized the morbidity and mortality associated with the surgical treatment of VS to 22% and 0.2%, respectively . Although the goal of surgical excision is total tumor removal, adherence to the facial nerve, the brainstem, or local vasculature precludes such an extent of resection to only 29% of patients . When gross total resection is achieved, the recurrence-free survival rates are 96%, 82%, 73%, and 56% at five, 10, 15 and 20 years. With a subtotal resection these recurrence-free survival rates drop to 47%, 17%, and 8% at five, 10, and 15 years . Progression following subtotal resection is three times more likely than with gross total or near total resection and approximates 40% [6,7,8,9]. In such cases of subtotal resection, the increased risk of progression often necessitates further intervention, including subsequent surgical resections or radiosurgery. To date, there are no meaningful clinical or molecular markers that distinguish non-progressing or slowly progressing tumor remnants from those that will progress rapidly.
In the present study, we analyzed transcriptomes of tumors from patients with sporadic VS to identify molecular perturbations that might be used to discriminate slow-growing VS remnants from more aggressive VS remnants. The ability to identify patients whose tumor remnants are likely to progress rapidly would select these patients for multimodal treatment.
Between July 2005 and April 2014 tissue samples were collected during the surgical excision of tumors from 17 patients with sporadic VS. Patients with neurofibromatosis, history of radiation treatment to the head and neck region or lacking follow up data were excluded. All patients had less than total surgical excision, and all samples were collected before radiotherapy or any other adjuvant treatment. Progression was defined as any radiographic evidence of lesion growth following surgery. Progression-free survival was defined as the time from the date of surgery to the date of radiographic evidence of progression. The study was approved by our Institutional Review Board (protocol PA13-0067), and clinical samples and data were collected only after patients provided written informed consent.
Multiplex immunofluorescence staining and image analysis
Our study cohort consisted of 8 patients with early disease progression (< 5 years after surgery), and 9 patients with either no progression or late disease progression (≥ 5 years after surgery). The 5 µm slides of the Formalin-sixed, paraffin-embedded (FFPE) tumor tissue were prepared and were stained for mIF. Multiplex immunofluorescence staining on a single slide was performed with use of the Opal 7-Color Manual IHC Kit (AKOYA Biosciences, #NEL811001KT). The antibodies used included CD4, CD8, CD20, CD68, and CD1A. The slides were scanned (Vectra Polaris, Akoya Biosciences) and the scanned images were analyzed at The University of Texas MD Anderson Cancer Center using AI-based software (VIS Image Analysis, Visiopharm).
Formalin-fixed, paraffin-embedded samples were submitted for analysis using HTG EdgeSeq panels. The area of each sample was measured, and HTG Lysis Buffer was added to obtain a per-well concentration of 6 mm2/35 µL. To improve sample lysis, we added proteinase K to the lysis buffer at a ratio of 1:20, and the samples were incubated at 50°C for 180 min. We added 35 µL of each sample to a single well of a 96-well plate. We also added 25 ng of human universal RNA to 3 wells to serve as a process control.
Inform spectral unmixing
We used the Vectra Polaris 3.0.3 multispectral imaging system (Akoya Biosciences) through the full emission spectrum from 440 to 780 nm, to extract fluorescence intensity information from the images using positive tonsil controls from each run staining to calibrate the spectral image scanner protocol at 20 × magnification (0.5 µm/pixel). Each marker was quantified individually using a spectral signature for each fluorophore obtained by the “spectral unmixing library” using the same algorithm from the InForm 2.4.8 image analysis software (Akoya Biosciences). The percentages representing each marker were calculated by dividing the absolute number of each marker by the absolute number of total nucleated cells (DAPI +) on each core at each time point. Following whole-slide image acquisition, images are analyzed with inForm image analysis software to quantify the cell-level biological features. The inForm software program was developed to integrate multispectral capabilities with image analysis to (1) spectrally unmix and isolate multiple Opal signals and background autofluorescence; (2) detect different tissue architecture (e.g., tumor, stroma, vessels, and necrosis) using a machine learning–based neural network pattern recognition function; (3) segment individual cells starting with nuclei, based on DAPI, and using other markers to detect membranous and cytoplasmic regions of cells; and (4) identify cell types of interest based on marker signal levels and cellular staining pattern using user-trained multinomial logistic regression algorithms. Once the slides were stained, they were scanned on a multispectral digital slide imaging system, the Vectra Polaris. The Vectra Polaris uses multispectral imaging technology to compensate for optical spectral bleed-through among channels and to isolate signal from background autofluorescence. Using inForm, designated library slides are used to isolate the exact spectral signature of each fluorophore to properly unmix each whole-slide composite image, as well as isolate and remove tissue autofluorescence. With multispectral unmixing, residual bleed-through was reduced to < 1% in both cases.
HTG EdgeSeq assay
Samples were run on an HTG EdgeSeq processor using an HTG EdgeSeq Oncology Biomarker Panel that allows for the measurement of expression of 2,549 genes (Supplementary Table S1). The samples were then individually barcoded using a 19-cycle PCR reaction to add adapters and molecular barcodes. Barcoded samples were individually purified using AMPure XP beads and quantitated using a KAPA Library Quantification kit. The library was sequenced on an Illumina MiSeq using a V3 150-cycle kit with 2 index reads. PhiX was spiked into the library at 5%; this spike-in control is standard for Illumina sequencing libraries.
Data were returned from the sequencer in the form of de-multiplexed FASTQ files, with one file per original well of the assay. To collate the data, we used HTG EdgeSeq Reveal software to align the FASTQ files to the probe list. We applied Median Ratio Normalization (MRN) to the aligned raw data prior downstream analysis. MRN data is available in the Supplementary Table S2.
Identification of differentially expressed genes
Gene expression analyses were performed using JMP PRO 15.2.1 software. To identify differentially expressed genes between patients with early (< 5 years) and late (≥ 5 years) progression, we performed unpaired two-sample Wilcoxon test and the p-value correction for multiple testing was done by the Benjamini & Hochberg FDR method. However, considering the small sample size (n = 17), we considered unadjusted p-values for downstream analyses. Gene expression patterns among samples were visualized by hierarchical clustering analysis (HCA) using the Ward minimum variance method for defining distances between clusters, and by principal component analyses (PCA).
Canonical pathway integrative analysis
To determine the potential biological mechanisms associated with the gene expression differences between the groups of schwannoma samples, we initially performed an overrepresentation analysis (ORA) using the KOBAS 3.0 online tool (http://kobas.cbi.pku.edu.cn) . This software evaluates whether a list of genes is statistically enriched by pathways and terms from several databases such as KEGG, Reactome, and Gene Ontology (GO). The gene list used for this analysis was performed by differentially expressed genes (p < 0.01) between early and late/no progression VS cases. ORA results were considered significant at FDR < 0.05 (by Benjamani and Hochberg).
To further explore the biology behind the transcriptional differences between the VS groups, we used the Ingenuity Pathway Analysis (IPA) software (Qiagen) that indicates whether differentially expressed genes are involved in the activation or inhibition of a curated set of canonical pathways and other biological mechanisms. IPA uses a priori knowledge of expected interactions between transcriptional regulators and their target genes stored in the Ingenuity Knowledge Base, a scientific literature–based database. Considering that the IPA analysis is more complex and considers statistical parameters associated with the degree of transcriptional variation between the compared groups, all differentially expressed genes (p < 0.05) between VS groups were used as input. Pathways were considered significantly enriched when the z-score was < 1 or > 1 and p < 0.05. Positive z-scores are indicative of activation, while negative z-scores are suggestive of inhibition of a given pathway.
Additionally, we investigated the existence of protein–protein interaction evidence among genes differentially expressed between VS groups (p < 0.05) using the online tool STRING. Gene interactions were considered only at the highest confidence score (> 0.9) . After filtering out genes without significant interactions, an additional ORA was performed using the STRING database.
Prediction model development
To evaluate potential predictive markers to discriminate between patients with early progression and those with late or no progression, we employed the k − Top Scoring Pair (KTSP) classifier . Briefly, KTSP is based on the Top Scoring Pair approach proposed by Geman et al. , in which gene expression data are converted into binary classifiers based on the expression level difference between 2 genes (e.g., GeneA > GeneB and GeneA < GeneB). The best gene-pair classifiers are those whose expression levels switch more consistently between the 2 groups of interest (e.g., GeneA > GeneB in 90% of group 1 samples, whereas GeneB > GeneA in 100% of group 2 samples). The KTSP approach is based on the same concept but considers a “k” number of gene pairs in the final prediction algorithm. KTSP initially calculates the discriminatory power—called “votes”—for all possible gene pairs in a gene expression set. The prognostic score is obtained by summing the votes among all “k” pairs. The best prognostic marker is that for which the “k” number of pairs provides the highest score . KTSP calculations were performed with the switchBox package in R .
CD8+ T cell isolation, activation, and co-culture with cancer cells (TR6Bc1), and single cell immune analysis
Mouse schwannoma cell line, TR6Bc1, is maintained in our lab with DMEM medium (with additional 2 mM Glutamine, 10% tryptose phosphate broth, and 10% fetal bovine serum) and was plated into a cell culture insert (24 well format) 24 h before co-culture with T cells. CD8+ T cells were then isolated and activated as previously described . Briefly, spleen was harvested from C3H mice and was pressed though cell strainer. Red blood cells (RBC) were lysed with lysis buffer (eBioscience 10X RBC Lysis Buffer (Multi-species), ThermoFisher). The CD8+ T cells were isolated with CD8+ T Cell Isolation Kit (Miltenyi Biotec B.V. & Co. KG) following the company’s protocol and activated for 18 h with CD3 antibody (coated, 2 µL/mL) and CD28 antibody (5 µg/ml) (BioLegend) in RPMI-1640 medium. After activation, the T cells were plated in a 24 well-plate and cultured with or without TR6Bc1 cells, for 24 h. Next, the T cells were harvested and stained with Cell stain 405 and anti-CD8 antibody (Alexa Fluor 647 anti-mouse CD8a (AF647-CD8) (IsoCode Kit, Isoplexis) following the manufacturer’s protocol. After staining, the T cells were re-suspended in RPMI-1640 medium (1 × 106 cells/mL) and 30 µL of the T cells was transferred into the inlet port of the IsoCode chips (IsoCode Kit, Isoplexis). The chips were then loaded into the Isolight instrument. The data for single cell CD8+ immune function were analyzed with IsoSpeak Software (Isoplexis).
Rapidly progressing and slowly/non-progressing VS have distinct immune signatures
Our study cohort consisted of age- and sex-matched patients with VS treated with surgery alone at The University of Texas MD Anderson Cancer Center from 2005–2014 stratified by those with early disease progression (< 5 years after surgery; n = 8) or either no progression or late disease progression (≥ 5 years after surgery; n = 9; Fig. 1A, B, Table 1). From patient tumors, we identified 5 major intratumoral immune cell types: CD4+ T cells, CD8+ T cells, CD20+ B cells, CD68+ macrophages, and CD1A+ dendritic cells. Using inForm spectral unmixing, we classified the phenotypes of 3.3 × 105 immune cells within neural cell adhesion molecule-positive tumor regions. Quantification of the innate immune cells revealed a significant enrichment of CD68+ macrophages and a significant depletion of dendritic cells (CD1a+) in rapidly progressing VS (P < 0.001; Fig. 1C, D). Regarding the adaptive immune cells, there was a significant enrichment in CD4+ and CD8+ T cells in rapidly progressing VS (P < 0.001; Fig. 1C, D). Although the overall number of CD20+ B cells across the tissue samples was low, a spectral isolation algorithm allowed us to confidently identify their phenotypes despite the low expression of CD20 across the tissue samples; this analysis revealed that the distribution of CD20+ B cells differed significantly between rapidly progressing VS and slowly progressing VS (P = 0.008; Fig. 1C, D). These different cellular distributions suggest that rapidly progressing VS and slowly progressing VS have distinct immune signatures in both the innate and adaptive immune compartments.
Innate immune response pathways are enriched in rapidly progressing VS
To interrogate the etiology of rapidly-progressing VS immune cell enrichment, we performed gene expression analyses using HTG EdgeSeq next-generation RNA sequencing. We analyzed the expression of 2,549 genes in 17 samples from the same patient cohort. Rapidly progressing VS samples exhibited 19 downregulated, and 24 upregulated genes, when compared with slowly or non-progressing VS cases (Fig. 2A, B) (Supplementary Table S3). To investigate the functional roles of these differentially expressed genes, we performed overrepresentation analyses (ORA) against several molecular databases (Supplementary Table S4). Considering the results obtained from the KEGG and Reactome databases (Fig. 1C), the differentially expressed genes were involved in immune activity and signal transduction, especially with Ras-ERK and PI3K-mTOR and their downstream signaling pathways. Similar findings were observed considering the enriched terms from the GO database (Fig. 2D). Genes such as NF1 (Neurofibromatosis type 1), RALA (RAS Like Proto-Oncogene A), PIK3C3, and RPS6KA5, upregulated in early VS progressors, were linked to the enrichment of most of the Ras-ERK and PI3K-mTOR cell signaling pathways and their downstream targets. On the other hand, the cytokines CSF2 and CXCL2 and the cytokine receptor IL2RA, were downregulated in early VS progressors, and, among other downregulated genes, were enriched among immune-related pathways and GO terms. The enrichment of pathways associated with immunity corroborates with our microscopic findings. Additionally, it suggests that a gamut of cell signaling pathways are altered among VS cases, which could represent a potential mechanism for our findings.
To improve our perspective on the direction of molecular enrichment changes between the VS groups, we explored our findings using the IPA software that provides an estimated activation or inhibition score (z-score) for significantly enriched pathways. Cellular senescence, cancer, epithelial-mesenchymal transition, and Ras-ERK & PI3K-mTOR-related signaling, were the top enriched pathways displaying a positive z-score (Fig. 3A), suggesting that these pathways are activated among fast-progressing VS. On the other hand, immune, stress response, and tumor suppressor signaling pathways were found significantly inhibited among fast VS progressors. These results corroborate with our other findings, suggesting that early VS progressors have a decreased local immunity associated with immune signaling inhibition. Additionally, this phenotype is accompanied by an increase of signaling pathways associated with malignant neoplasms, especially those pathways directly or indirectly related to Ras-ERK & PI3K-mTOR signaling, which are major modulators of survival and proliferation mechanisms.
Altered adaptive immune response is associated with a possible viral etiology in rapidly progressing VS
The lack of enrichment in adaptive immune response − related pathways in the presence of increased B- and T-cell populations suggests adaptive immune exhaustion or dysfunction in rapidly progressing VS. When taken in the context of well-preserved, upregulated innate immunity − related pathways in the presence of severely altered adaptive immune response, we hypothesized that adaptive immunosenescence underlies rapidly progressing VS. This hypothesis is supported by the identification of neuroinflammatory and cellular senescence pathways enriched in rapidly progressing VS (Fig. 3A, B and Supplementary Tables S5 and S6).
Although immunosenescence is most often associated with aging, it can also result from continuous exposure to external perturbations. In fact, we observed an enrichment of several stress-related pathways upregulated in rapidly progressing VS, such as Senescence, Telomerase Signaling, Hepatic Fibrosis Signaling Pathway, Autophagy, LPS-stimulated MAPK Signaling, Role of BRCA1 in DNA Damage Response, Activation of IRF by Cytosolic Pattern Recognition Receptors, UVB-Induced MAPK Signaling, among others, suggesting that these tumors are under the effect of a stressor. Interestingly, we found that rapidly progressing VS cases exhibited activation of several pathways linked to viral infection, such as the canonical pathway NF-κB Activation by Viruses, and others highlighted in Fig. 3C. These findings may indicate that fast progressing VS are associated with an underlying viral infection, which is associated with the dramatic upregulation of stress-related and proliferation/survival signaling pathways, as well as the significant downregulation of immune signaling.
We used the STRING database to construct a protein–protein interaction network. The protein–protein interaction network consisted of 343 nodes and 559 edges with a significant number of interactions (P < 0.0001). Considering only the genes with known high confidence interactions with each other, we performed an additional ORA using the STRING database. These genes significantly enriched several pathways and GO terms associated with molecular signaling associated with viral infection, especially mechanisms associated with immune innate response against viruses (Fig. 4, Supplementary Table S7).
Immunosenescence gene signature predicts rapid VS progression
Given our finding that rapidly progressing VS has enriched expression of immune response-related genes, we focused our biomarker discovery approach on 438 immune-related genes included in our HTG EdgeSeq Oncology Biomarker panel (designated “Cluster of Differentiation” and “Immuno-Oncology”). k − Top Scoring Pair (KTSP) analysis identified 9 gene pairs (18 genes total) whose combined scores provided the best discriminatory power between rapidly progressing VS and slowly or non-progressing VS. The finding of more CD4+ and CD8+ T cells in rapidly progressing VS in the absence of increased T-cell activation levels suggests immunosenescence.
Immunosenescence is most often described as a state in which either the number or function of a patient’s total T cells is markedly poorer than those of younger or healthier patients. This decrease in available or functional T cells is thought to mediate older individuals’ increased susceptibility to disease . Hence, we further pruned the list of 18 genes based on the impact of these genes’ relative up- or downregulation on immunosenescence. KTSP analysis of the final list of 12 genes revealed several genes that correlate with functional adaptive immune activity—and thus negatively correlate with immunosenescence—that were upregulated in our samples. These genes included CD38, CD22, MADCAM-1, APCS, IL-6, and CXCL2. Conversely, our analysis uncovered several genes associated with T-lymphocyte dysfunction that were downregulated. The downregulation of these genes, which included KDR, CXCR6, STAT5A, MPL, NFATC3, and TLR6, also signifies a healthy adaptive T-cell response. The overexpression of each of these genes typically results in the impaired T-cell activity seen in immunosenescence. This final analysis generated a prognostic gene-set biomarker composed of 6 gene pairs (CD38-KDR, CD22-STAT5A, APCS-CXCR6, MADCAM1-MPL, IL6-NFATC3, and CXCL2-TLR6; Fig. 5A, B), which correctly included all rapidly progressing VS in our cohort and correctly excluded 7 of 9 slowly or non-progressing VS (sensitivity = 1.00, specificity = 0.78; Fig. 5A, B). of note, the two patients that were incorrectly classified as rapidly progressing VS, had relatively large tumors (i.e. 28 and 30 mm); yet neither experienced disease progression during follow up period.
Single cell proteomic analysis shows loss of polyfunctionality in CD8+ T cells cultured with VS cells
To test our hypothesis that VS can alter immune function leading to suppressed immunosurveillance that potentially allows for VS progression, we next sought to characterize the functional impact of VS on effector T cell function. To accomplish this, we analyzed single cell proteomic profiles (Mouse Adaptive Immune panel, IsoPlexis, USA) of activated murine T cells cultured with and without murine vestibular schwannoma TR6Bc1 cells. Principal component analysis showed that CD8+ T cells cultured with vestibular schwannoma cells are functionally distinct from activated T cells cultured alone (Fig. 6A – C). In the presence of VS cells we found a significant downregulation of effector phenotype markers (e.g. Granzyme B, IFNγ and TNFα), chemoattractive molecules (e.g. IP-10, CCL11 and RANTES), inflammatory markers (e.g. IL17A and MCP1) and stimulatory signals (e.g. GM-CSF, IL2 and IL7) (Fig. 6D – E). Taken together, the loss of effector T cells polyfunctionality in the presence of VS supports the notion that schwannoma cells might drive immunosenescence in the tumor microenvironment of patients with rapidly progressing VS.
VS recur and progress at unpredictable rates, complicating disease management especially in the context of less than total surgical resection. Our results shed light on the potential mechanisms involved in VS progression. We found that multiple pathways, particularly immune-associated pathways involved in innate immunity and antigen presentation, are dysregulated in rapidly progressing forms of sporadic VS. Surprisingly, many of the cell replication- and tumorigenesis-associated pathways were not enriched, and some of these pathways (e.g., those related to DNA synthesis and central nervous system solid tumor formation) were significantly downregulated in patients with rapidly progressing VS. Taken together, these results indicate that the drivers of VS progression are external rather than internal genetic perturbations.
The drivers of sporadic VS associated with progression have never been characterized previously, although electron microscopical studies of bovine schwannomas have revealed viral particles in the tumors, suggesting a potential viral etiology [17, 18]. In the present study, we uncovered the involvement of viral response pathways, underscoring the importance of adaptive immune dysfunction in rapidly progressing VS, and perhaps even a role of viruses themselves in VS oncogenesis. Our findings show that in rapidly progressing VS, virus-inducible cellular gene networks, such as interferon gamma signaling, are enriched. The timing of potential viral infection still remains to be elucidated, that is, whether viral infection occurred after initial surgery for VS or whether the presence of such enrichment represents a reactivation of latent viral infection . Yet the interaction between a virus and its host organism, which is innately inflammatory, can modify the tumor immunobiology, and may overwhelm and inhibit specific acquired responses necessary to achieve viral clearance. Similarly, we detected significant dysregulation of the expression of Toll-like receptors and IRF3 and IRF7 in sporadic VS. Upon activation, Toll-like receptors recruit adapter proteins that orchestrate inflammatory responses by the infected cells; immune cells in the tumor microenvironment that have detected a virus may also release anti-viral factors such as interferons.
Our findings showed a recruitment of CD4+, CD8+, and CD20+ cells to the tumor microenvironment in rapidly progressing VS. However, RNA sequencing suggested that these lymphocytes are senescent. Hence, we evaluated the predictive role of these cells’ dysregulated immunosenescence-related genes; these genes, identified and validated in our cohort, comprise the prognostic biomarker for VS progression rate. Our final biomarker is composed of 6 gene pairs, each known to play a distinct and important role in adaptive immune function: CD38-KDR, CD22-STAT5A, APCS-CXCR6, MADCAM1-MPL, IL6-NFATC3, and CXCL2-TLR6. CD38, expressed on activated T-cells is a well-defined surface protein denoting a robust T-cell response . Other identified genes, are also indicative of a healthy T-cell response, although they are expressed on cooperating immune cells . CD22, expressed on activated B-cells, denotes T-cell:B-cell crosstalk . MADCAM-1 is a mucosal adhesion molecule that binds surface proteins on T lymphocytes, thereby guiding them into lymphoid tissue for activation . Although MADCAM-1 has been studied primarily in the context of gastrointestinal lymphoid tissue , queries of the GTEx database reveal robust MADCAM-1 expression in several other organs, including the thyroid and brain. In line with these findings, our data indicates that MADCAM-1 expression was increased in late recurring VS.
Other genes associated with T-cell activity, like APCS, do not directly interact with T lymphocytes. Rather, they interact with important T cell-related proteins, thereby indirectly reflecting T-cell activity. APCS, also known as SAP, is a major acute-phase reactant protein that is dramatically upregulated during IL-6-mediated inflammation [25, 26], indirectly driving T-cells to release further IL-6, IL-17, and CXCL2. Thus, the increased expression of APCS, IL-6, and CXCL2 found in our patient cohort corresponds with improved immune cell function, and subsequent delays in VS tumor recurrence. Additionally, the 6 markers (KDR, CXCR6, STAT5A, MPL, NFATC3, and TLR6) we found to be downregulated in our late-recurring tumor samples are also biologically relevant to functional T-cell activity, which helps explain their expression and function in this context.
Other identified biomarker genes impact T-cell migration and differentiation rather than cytotoxicity and activation. Activated KDR, also known as VEGFR2, acts in an immunosuppressive manner by inhibiting the migration of T lymphocytes [27, 28]. Concordantly, the inhibition of KDR activity with an anti-KDR antibody improves the T-cell response [29, 30]. Thus, the decreased expression of KDR seen in our analysis associates with a stronger immune response and a subsequent delay in tumor recurrence. Another membrane-bound marker of T-cell activity is CXCR6, a chemokine receptor expressed on T lymphocytes that is downregulated upon T-cell activation . Thus, the downregulation of CXCR6 expression is a useful biomarker of functional T-cell response. In our patients, downregulation of CXCR6 associated with better tumor control and late VS recurrence. In concert with this, CXCR6 knockout mice are better able to control infections that are known to depend on a T-cell response. However, these effects appear to be independent of T-cell activity, suggesting the need for further work to elucidate the mechanism underlying these effects .
STAT5A is a multi-functional protein that regulates several immune-related processes, including T-cell differentiation [33, 34]. However, while serving as a potential marker of adaptive immune function, STAT5A activity and overexpression predict both early recurrence and tumor aggression in head and neck carcinoma , as well as prostate cancer patients , matching our findings in VS. Similarly, the protein MPL, also known as CD110, is vital for proper immune function; MPL knockout mice generate only 10% of the megakaryocytes that wild type mice do . However, MPL overexpression is associated with increased aggression and poor prognosis in cancer patients , again mirroring the expression pattern and cancer phenotype seen in our cohort. Given these proteins’ relation to proper immune cell function and the importance of immunosenescence in cancer development, progression, and recurrence, the misregulation of STAT5A and MPL likely mediates increased oncogenesis in VS via their immune-related pathways.
Nuclear factor of activated T cells 3 (NFATC3) is a transcription factor that promotes the expression of several genes that are essential for proper T-cell development and activity. Knockout studies of members of the NFAT gene family demonstrate a functional redundancy between family members . Thus, NFATC3 downregulation may not be particularly damaging to a cell. However, NFATC3 also drives the expression of stem cell-promoting proteins like OCT4 . Thus, NFATC3 overexpression and hyperactivity may be particularly oncogenic. Similarly, TLR6 is a membrane-bound receptor found on the surface of T cells. TLR6 activation, a sign of innate immune activation, results in inflammation at the site of expression . Inflammation is a key immune-related factor in cancer development and progression, and sequence variation in TLR6 predicts inflammation-related cancer development . Moreover, TLR6 overexpression is associated with poor outcomes in patients with certain cancers . These descriptions match the findings of the present study, in which NFATC3 and TLR6 overexpression were associated with early recurrence, whereas decreased NFATC3 and TLR6 expression were associated with late recurrence. Given the links between NFATC3 and TLR6 activity and proper immune system functioning, the overexpression described here likely has both immunosenescent and oncogenic effects. Therefore, similar to the decreased expression of KDR, CXCR6, STAT5A, and MPL, the decreased expression of NFATC3 and TLR6 is an indicator of functional T-cell activity and subsequent late VS recurrence.
In vitro, VS cells induce “deactivation” of CD8+ T cells; these findings might indicate that inhibited immune surveillance results in rapid progression of VS and raise the tantalizing possibility that drugs that target immune cell regulation could be useful in the treatment of VS.
Our study was limited by the low number of patients. Additionally, this is the first analysis of its kind in sporadic VS, and there currently exist no external datasets with available expression data; hence, our work requires further validation, including validation using in vivo models. Still, our well-defined, surgically-treated patient population enabled the detection of signals that allowed us to further explore the possibility of a viral link with progression of residual tumor following the subtotal resection of sporadic VS. Our data confirm the role of innate immune response in, and the potential viral etiology of, VS progression. Our study tested a relatively small number of genes (2,549) and it employed bulk rather than single-cell sequencing. Although the assay we used was not designed for the evaluation of the impact that external perturbations have on the tumor microenvironment, its results clearly demonstrate their role in disease progression. These novel results provide data that support further investigation of the immunobiology of VS. Further studies are needed to identify any potential viral pathogens, prior to considering any antiviral approach. The senescent or exhausted adaptive immune microenvironment in rapidly progressing VS, however, suggests a potential role for immune checkpoint (e.g. anti-PD1) and innate immune cell targeting.
The rapid progression of residual VS following subtotal surgical resection has an underlying immune etiology that may be virally originating. Despite the adaptive immune response, we find evidence that T-cell immunosenescence may be associated with the rapid progression of VS. These findings provide a rationale for clinical trials evaluating immunotherapy in patients with rapidly progressing VS. Here we propose a non-platform dependent (KTSP), transcriptomic signature to allow identification of patients with rapidly progressing VS. Harnessing this signature will allow for rapid transcriptomic sequencing of surgically resected samples to stratify patients based on their risk of rapid progression, thereby identifying candidates for adjuvant multimodal therapy or immunotherapy.
Availability of data and materials
The datasets generated and/or analysed during the current study are not publicly available to maintain individual patient privacy, but are available from the corresponding author on reasonable request.
K-top scoring pair
Kshettry VR, Hsieh JK, Ostrom QT, Kruchko C, Barnholtz-Sloan JS. Incidence of vestibular schwannomas in the United States. J Neurooncol. 2015;124(2):223–8.
Carlson ML, Habermann EB, Wagie AE, Driscoll CL, Van Gompel JJ, Jacob JT, et al. The Changing Landscape of Vestibular Schwannoma Management in the United States - A Shift Toward Conservatism. Otolaryngol Head Neck Surg (United States). 2015;153(3):440–6.
Nakatomi H, Jacob JT, Carlson ML, Tanaka S, Tanaka M, Saito N, et al. Long-term risk of recurrence and regrowth after gross-total and subtotal resection of sporadic vestibular schwannoma. J Neurosurg. 2017;1–7. https://doi.org/10.3171/2016.11.JNS16498.
Sughrue ME, Yang I, Aranda D, Rutkowski MJ, Fang S, Cheung SW, et al. Beyond audiofacial morbidity after vestibular schwannoma surgery: Clinical article. J Neurosurg. 2011;114(2):367–74.
Fabbris C, Gazzini L, Paltrinieri D, Marchioni D. Exclusive surgical treatment for vestibular schwannoma regrowth or recurrence: A meta-analysis of the literature. Clin Neurol Neurosurg. 2020;193(February):105769. https://doi.org/10.1016/j.clineuro.2020.105769.
El-Kashlan HK, Zeitoun H, Arts HA, Hoff JT, Telian SA. Recurrence of acoustic neuroma after incomplete resection. Am J Otol. 2000;21(3):389–92 (http://www.ncbi.nlm.nih.gov/pubmed/10821553 Available from).
Strickland BA, Ravina K, Rennert RC, Jackanich A, Aaron K, Bakhsheshian J, et al. Intentional Subtotal Resection of Vestibular Schwannoma: A Reexamination. J Neurol Surgery Part B Skull Base. 2020;81(2):136–41.
Monfared A, Corrales CE, Theodosopoulos PV, Blevins NH, Oghalai JS, Selesnick SH, et al. Facial nerve outcome and tumor control rate as a function of degree of resection in treatment of large acoustic neuromas: Preliminary report of the Acoustic Neuroma Subtotal Resection Study (ANSRS). Neurosurgery. 2016;79(2):194–200.
Bloch DC, Oghalai JS, Jackler RK, Osofsky M, Pitts LH. The fate of the tumor remnant after less-than-complete acoustic neuroma resection. Otolaryngol Head Neck Surg. 2004;130(1):104–12.
Bu D, Luo H, Huo P, Wang Z, Zhang S, He Z, et al. KOBAS-i: Intelligent prioritization and exploratory visualization of biological functions for gene enrichment analysis. Nucleic Acids Res. 2021;49(W1):W317–25.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Marchionni L, Afsari B, Geman D, Leek JT. A simple and reproducible breast cancer prognostic test. BMC Genomics. 2013;14(1):0–6.
Geman D, D’Avignon C, Naiman DQ, Winslow RL. Classifying gene expression profiles from pairwise mRNA comparisons. Stat Appl Genet Mol Biol. 2004;3(1):1–16.
Afsari B, Fertig EJ, Geman D, Marchionni L. SwitchBox: An R package for k-Top Scoring Pairs classifier development. Bioinformatics. 2015;31(2):273–4.
Lim JF, Berger H, Su IH. Isolation and activation of murine lymphocytes. J Vis Exp. 2016;2016(116):1–8.
Pawelec G. Immunosenescence and cancer. Biogerontology. 2017;18(4):717–21.
Chandra S, Jerva MJ, Clemis JD. Ultrastructural Characteristics of Human Neurilemoma Cell Nuclei. Cancer Res. 1975;35(8):2000–6.
Canfield PJ, Doughty FR. A study of virus-like particles present in bovine nerve sheath tumours. Aust Vet J. 1980;56(6):257–61 (Available from: papers3://publication/uuid/C138351A-2768-45C9-A28F-9C3BCF9764E4).
Franco-Vidal V, Nguyen DQ, Guerin J, Darrouzet V. Delayed facial paralysis after vestibular schwannoma surgery: Role of herpes viruses reactivation - Our experience in eight cases. Otol Neurotol. 2004;25(5):805–10.
Deaglio S, Vaisitti T, Billington R, Bergui L, Omede’ P, Genazzani AA, et al. CD38/CD19: A lipid raft-dependent signaling complex in human B cells. Blood. 2007;109(12):5390–8.
Habtezion A, Nguyen LP, Hadeiba H, Butcher EC. Leukocyte Trafficking to the Small Intestine and Colon. Gastroenterology. 2016;150(2):340–54.
Clark EA, Giltiay NV. CD22: A Regulator of Innate and Adaptive B Cell Responses and Autoimmunity. Front Immunol. 2018;9(September):2235.
Berlin C, Berg EL, Briskin MJ, Andrew DP, Kilshaw PJ, Holzmann B, et al. α4β7 integrin mediates lymphocyte binding to the mucosal vascular addressin MAdCAM-1. Cell. 1993;74(1):185–95.
Buettner M, Dittrich-Breiholz O, Falk CS, Lochner M, Smoczek A, Menzel F, et al. Stromal cells as trend-setters for cells migrating into the lymph node. Mucosal Immunol. 2015;8(3):640–9. https://doi.org/10.1038/mi.2014.97.
Pepys MB, Baltz ML. Acute phase proteins with special reference to C-reactive protein and related proteins (pentaxins) and serum amyloid A protein. Adv Immunol. 1983;34:141–212.
Mortensen RF, Beisel K, Zeleznik NJ, Le PT. Acute-phase reactants of mice II. Strain dependence of serum amyloid P-component (SAP) levels and response to inflammation. J Immunol. 1983;130(2):885–9.
Voron T, Marcheteau E, Pernot S, Colussi O, Tartour E, Taieb J, et al. Control of the immune response by pro-angiogenic factors. Front Oncol. 2014;4(April):1–9.
Voron T, Colussi O, Marcheteau E, Pernot S, Nizard M, Pointet AL, et al. VEGF-A modulates expression of inhibitory checkpoints on CD8+ T cells in tumors. J Exp Med. 2015;212(2):139–48.
Wallin JJ, Bendell JC, Funke R, Sznol M, Korski K, Jones S, et al. Atezolizumab in combination with bevacizumab enhances antigen-specific T-cell migration in metastatic renal cell carcinoma. Nat Commun. 2016;7:1–8.
Tada Y, Togashi Y, Kotani D, Kuwata T, Sato E, Kawazoe A, et al. Targeting VEGFR2 with Ramucirumab strongly impacts effector/activated regulatory T cells and CD8+ T cells in the tumor microenvironment. J Immunother Cancer. 2018;6(1):1–14.
Koprak S, Matheravidathu S, Springer M, Gould S, Dumont FJ. Down-regulation of cell surface CXCR6 expression during T cell activation is predominantly mediated by calcineurin. Cell Immunol. 2003;223(1):1–12.
Ashhurst AS, Flórido M, Lin LCW, Quan D, Armitage E, Stifter SA, et al. CXCR6-deficiency improves the control of pulmonary mycobacterium tuberculosis and influenza infection independent of T-lymphocyte recruitment to the lungs. Front Immunol. 2019;10(MAR):1–16.
Ren Z, Aerts JL, Vandenplas H, Wang JA, Gorbenko O, Chen JP, et al. Phosphorylated STAT5 regulates p53 expression via BRCA1/BARD1-NPM1 and MDM2. Cell Death Dis. 2016;7(12):e2560.
Kagami SI, Nakajima H, Suto A, Hirose K, Suzuki K, Morita S, et al. Stat5a regulates T helper cell differentiation by several distinct mechanisms. Blood. 2001;97(8):2358–65.
Lai SY, Childs EE, Xi S, Coppelli FM, Gooding WE, Wells A, et al. Erythropoietin-mediated activation of JAK-STAT signaling contributes to cellular invasion in head and neck squamous cell carcinoma. Oncogene. 2005;24(27):4442–9.
Li H, Zhang Y, Glass A, Zellweger T, Gehan E, Bubendorf L, et al. Activation of signal transducer and activator of transcription-5 in prostate cancer predicts early recurrence. Clin Cancer Res. 2005;11(16):5863–8.
Storan MJ, Heazlewood SY, Heazlewood CK, Haylock DN, Alexander WS, Neaves RJ, et al. Brief report: Factors released by megakaryocytes thrombin cleave osteopontin to negatively regulate hematopoietic stem cells. Stem Cells. 2015;33(7):2351–7.
Yan Z, Ohuchida K, Zheng B, Okumura T, Takesue S, Nakayama H, et al. CD110 promotes pancreatic cancer progression and its expression is correlated with poor prognosis. J Cancer Res Clin Oncol. 2019;145(5):1147–64. https://doi.org/10.1007/s00432-019-02860-z.
Urso K, Alfranca A, Martínez-Martínez S, Escolano A, Ortega I, Rodríguez A, et al. NFATc3 regulates the transcription of genes involved in T-cell activation and angiogenesis. Blood. 2011;118(3):795–803.
Lee SH, Kieu C, Martin CE, Han J, Chen W, Kim JS, et al. NFATc3 plays an oncogenic role in oral/oropharyngeal squamous cell carcinomas by promoting cancer stemness via expression of OCT4. Oncotarget. 2019;10(23):2306–19.
Wang HF, Li Y, Wang YQ, Li HJ, Dou L. MicroRNA-494-3p alleviates inflammatory response in sepsis by targeting TLR6. Eur Rev Med Pharmacol Sci. 2019;23(7):2971–7.
Sun J, Wiklund F, Zheng SL, Chang B, Bälter K, Li L, et al. Sequence variants in toll-like receptor gene cluster (TLR6-TLR1-TLR10) and prostate cancer risk. J Natl Cancer Inst. 2005;97(7):525–32.
Akhter A, Street L, Ghosh S, Burns BF, Elyamany G, Shabani-Rad MT, et al. Concomitant high expression of toll-like receptor (TLR) and B-cell receptor (BCR) signalling molecules has clinical implications in mantle cell lymphoma. Hematol Oncol. 2017;35(1):79–86.
We would like to thank Joseph Munch for editing this manuscript.
This work was funded in part by the MDACC Head and Neck Cancer Moonshot, (MA), NIH/NCI R37 CA242006-01A1 (MA), the American Heart Association (20PRE35040011) (PJH), BRASS: Baylor Research Advocates for Student Scientists (PJH), the Von Plonski fund for acoustic neuroma research (FDM), and the Department of Neurosurgery, M.D. Anderson Cancer Center (FDM).
Ethics approval and consent to participate
The study was approved by our Institutional Review Board (protocol PA13-0067), and clinical samples and data were collected only after patients provided written informed consent.
Consent for publication
All patients have given consent for this data to be published.
None of the authors have a conflict of interest. MEN reports stock ownership in 3 M, Amgen, Cardinal Health, Johnson & Johnson, Medtronic, and Pfizer.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Amit, M., Xie, T., Gleber-Netto, F.O. et al. Distinct immune signature predicts progression of vestibular schwannoma and unveils a possible viral etiology. J Exp Clin Cancer Res 41, 292 (2022). https://doi.org/10.1186/s13046-022-02473-4