Xubin Dong1,2†
Percy David Papa Akuetteh3†
Jingjing Song4†
Chao Ni5†
Cong Jin1,2
Huihui Li1Wenjie Jiang1Yuhao Si1
Xiaohua Zhang1,2
Qiyu Zhang3*
Guanli Huang2,6*- 1Department of Breast Surgery, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 2Department of Thyroid Surgery, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 3Department of Hepatobiliary Surgery, The First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
- 4Department of Pediatric Allergy and Immunology, The Second Affiliated Hospital and Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, China
- 5Children’s Heart Center, Institute of Cardiovascular Development and Translational Medicine, the Second Affiliated Hospital and Yuying Children’s Hospital of Wenzhou Medical University, Wenzhou, China
- 6Department of Thyroid Surgery, The Quzhou Affiliated Hospital of Wenzhou Medical University, Quzhou People’s Hospital, Quzhou, China
Papillary thyroid cancer (PTC) is the most common malignancy of the endocrine system, with an increase in incidence frequency. Major vault protein (MVP) is the main structural protein of the vault complex that has already been investigated in specific cancers. Yet the underlying biological functions and molecular mechanisms of MVP in PTC still remain considerably uncharacterized. Comprehensive analyses are predicated on several public datasets and local RNA-Seq cohort. Clinically, we found that MVP was upregulated in human PTC than in non-cancerous thyroid tissue and was correlated with vital clinicopathological parameters in PTC patients. MVP expression was associated with BRAFV600E, RAS, TERT, and RET status, and it was correlated with worse progression-free survival in PTC patients. Functionally, enrichment analysis provided new clues for the close relationship between MVP with cancer-related signaling pathways and the immune microenvironment in PTC. In PTC with high MVP expression, we found CD8+ T cells, regulatory T cells, and follicular helper T cells have a higher infiltration level. Intriguingly, MVP expression was positively correlated with multiple distinct phases of the anti-cancer immunity cycle. MVP knockdown significantly suppressed cell viability and colony formation, and promoted apoptosis. In addition, downregulated MVP markedly inhibited the migration and invasion potential of PTC cells. The rescue experiments showed that MVP could reverse the level of cell survival and migration. Mechanistically, MVP exerts its oncogenic function in PTC cells through activating PI3K/AKT/mTOR and MAPK/ERK pathways. These results point out that MVP is a reliable biomarker related to the immune microenvironment and provide a basis for elucidating the oncogenic roles of MVP in PTC progression.
Introduction
The thyroid is an important endocrine organ and a pathogenic target of human autoimmune diseases. Thyroid cancers are histologically subdivided into papillary (PTC), follicular, poorly differentiated, and anaplastic thyroid cancers (Landa et al., 2016). Plenty of studies had shown that the recurrence risk and cancer-related mortality of PTC patients were tightly associated with certain clinicopathological characteristics, such as the age of the first diagnosis, the size of the primary tumor, extrathyroidal invasion, and distant metastasis (Johnson et al., 1988; Asioli et al., 2010; Randolph et al., 2012). The overall prognosis of PTC patients was relatively satisfactory (Ghossein and Livolsi, 2008); however, patients with specific PTC subtypes (hobnail and tall cell variants) had poor prognoses (Yoo et al., 2016). Published studies had found that tumor-infiltrating immune cells (TIICs) could impact the response to immunotherapy, chemoresistance, and clinical outcome (Deschoolmeester et al., 2010). Therefore, analyzing the immunological function of the tumor microenvironment (TME), elucidating the molecular mechanisms of PTC development, and establishing promising immune-related biomarkers for PTC are crucial.
MVP is the main structural protein of the vault complex with a molecular weight of 110 kDa, which exists plentiful in the cytoplasm of eukaryotic cells (Vanzon et al., 2005). It was primarily reported to be involved in chemoresistance and to be significantly associated with the prognosis of various cancers, such as triple-negative breast cancer (TNBC) (Xiao et al., 2019) and liver (Losert et al., 2012; Yu et al., 2020), lung (Shen et al., 2019), and colon cancers (Henríquez-Hernández et al., 2012; Lötsch et al., 2013a; Teng et al., 2017). MVP also plays an essential role in multidrug resistance, autophagy, and signal transduction (van Zon et al., 2003). As MAPK-driven cancer, PTC mainly interacts with the mutually exclusive drivers (BRAFV600E and mutated RAS); PTC could be molecularly subdivided into BRAFV600E-like (BVL) and RAS-like (RL) that represent differential regulation of the MAPK pathway and thyroid differentiation (Agrawal et al., 2014). Additionally, studies had also shown that MVP was correlated with signaling pathways related to drug resistance, such as PI3K/AKT and MAPK pathways (Park, 2012; Liu et al., 2019). However, there has been no related literature reporting the clinical and biological significance of MVP in PTC, and the potential relationship between MVP and TIME is yet to be determined.
In this study, we comprehensively made an inquiry from several aspects into an expression overview of MVP and its relationship with clinicopathological factors of PTC patients. Furthermore, we have conducted an in-depth exploration of TME through the analyses of the latest immune-related algorithms toward RNA-seq profiles. Significantly, we evaluated the potential clinical application value of MVP by analyzing the critical molecular characteristics of PTC. Our research results were broadly consistent in public databases and the local WMU-PTC cohort. In addition, based on functional enrichment analysis, we carried out a series of cellular and molecular assays to explore the connection between MVP and tumor immune-related molecular characteristics and to demonstrate the effects of knockdown and over-expressed MVP on PTC cell lines. The downstream regulatory mechanism of MVP was also preliminarily explored.
Materials and Methods
Datasets and Algorithms
Our study applied various public datasets of thyroid cancer. For the TCGA-THCA cohort, RNA-seq (RSEM normalization, level 3) and clinical profiles of thyroid cancer patients were acquired from the TCGAbiolinks package (Agrawal et al., 2014). This cohort included 501 PTC, 59 adjacent normal tissues (ANT), and 8 metastatic thyroid cancer samples (Supplementary Table S1, S2). The definition and data of progression-free survival (PFS) were acquired from TCGA-Clinical Data Resource (CDR) (Liu et al., 2018). Mutation status was acquired from VarScan2 MAF files of simple nucleotide variation from the GDC portal. Tumor mutation burden (TMB) was determined by the total number of non-silent mutations in specimens. MSI status (MANTIS score) for TCGA-THCA was used from a recent study (Bonneville et al., 20172017). The mRNA expression data and matched clinical information of GSE33630 (Dom et al., 2012), GSE60542 (Tarabichi et al., 2015), and GSE5364 (Yu et al., 2008) datasets were obtained from the GEO. In the present research, the Human Protein Atlas (HPA) immunohistochemistry images were used to explore the distribution and subcellular localization of MVP, comparing protein expression in malignant tissues and ANT (Uhlén et al., 2015; Uhlen et al., 2017). All IHC staining in the HPA project were performed using a standard protocol as described before (Kampf et al., 2012).
The “limma” R package was used to identify the differentially expressed genes (DEGs) between PTC subgroups. t-SNE was performed to explore the distribution of different groups using the “Rtsne” R package. The Spearman correlation analysis was evaluated by cor. test function in R. Genes and genes which had Spearman rho >0.5 or < −0.5, p < 0.001 with MVP were selected. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were completed by clusterprofiler package. The Gene set enrichment analysis (GSEA) was carried out utilizing the hallmark and C6 oncogenic signature gene sets from the Molecular Signature Database (Subramanian et al., 2005). Significantly enriched gene sets were identified when the FDR-adjusted q value <0.25.
ESTIMATE package was used to quantify the immune score, stromal score, microenvironment score, and tumor purity (Yoshihara et al., 2013). For particular TIIC analyses, quanTIseq, CIBERSORT, and MCPcounter algorithms were performed to calculate the relative proportions of different immune components in TME. CIBERSORT is a deconvolution-based method for quantitatively estimating TIME contexture from RNA-seq data (Newman et al., 2019). MCPcounter is a computational approach based on RNA-Seq data of immune-specific marker genes (Becht et al., 2016a). quanTIseq is a deconvolution-based algorism to quantify the proportions of immune components from bulk RNA-sequencing data in TME (Finotello et al., 2019). Tumor ImmunoPhenotype (TIP) was an ssGSEA-based pipeline applied for cancer-immunity cycle profiling (Xu et al., 2018), which evaluates the relative activity of the main 7 steps of the cancer-immunity cycle.
Patients, Thyroid Tissue Specimens, and RNA-Seq
79 pairs of thyroid specimens were derived from the Department of Thyroid Surgery, the First Affiliated Hospital of Wenzhou Medical University. Fresh tissues were snap-frozen in liquid nitrogen at the time of thyroidectomy and then stored at -80°C. Detailed clinicopathological information of WMU-PTC patients is summarized in Supplementary Table S3, S4. Details of the RNA-seq experimental protocol have been described in the previous publication (Dong et al., 2021). Briefly, total RNA was utilized to establish cDNA libraries for high-throughput RNA sequencing. The RNA expression proles were identified in the sequencing libraries obtained from an NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, United States). The clustering of the sample was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumina), and the library was sequenced on an Illumina NovaSeq platform. Due to the overall young age in our local cohort patients, to have more pT3/T4 and Stage IIIIV patients and better demonstrate the potential factors that affect the PTC progression, we evaluate our cohort with the seventh edition of the AJCC/TNM staging system (Lamartina et al., 2018).
Cell Incubation and RNA Interference
HTori-3, BCPAP, TPC-1, and KTC-1 cells were all obtained from Shanghai Cell Biology, Institute of the Chinese Academy of Sciences (Shanghai, China). These cells were all cultivated in the RPMI 1640 medium (Gibco, United States) supplemented with 10% FBS (PAN Biotech, Germany). The genotypes of the PTC cell lines are shown in Supplementary Table S5. Lipofectamine RNAiMAX transfection reagent (Thermo Fisher Scientific, United States) was mixed with siRNA to transfect PTC cells. The sequences of the si-MVP were as follows: siMVP-1 sense: 5′- CCTACATGCTGACCCAGGA -3′, siMVP-2 sense: 5′-ATCATTCGCACTGCTGTC -3′, siMVP-3 sense: 5′- GCAGATGACAGAGGCCATA-3′.
Ectopic Expression
Full-length MVP cDNA was synthesized and inserted into pCDH-GFP +PURO-3xFlag and pCDH-GFP PURO vectors (Genepharma, Shanghai, China). The resulting vector or empty vector was transfected into PTC cells using Lipofectamine 2000 Transfection Reagent (Life Technologies, Carlsbad, CA) according to the manufacturer’s protocol. Infected cells were selected with puromycin (Invivogen) at 1 μg/ml.
qRT-PCR and Primers
Total RNA was isolated from patient tissues and thyroid cancer cell lines by TRIzol reagent (Invitrogen, United States). All RNA specimens were temporarily stored at -80°C. The isolated RNA was analyzed at 260/280 nm, which ranges from 1.81 to 1.97. RNA reverse transcription was accomplished by the ReverTra Ace qPCR RT Kit (Toyobo, Japan). Real-time PCR was carried out and analyzed through an ABI 7500 System (Life Technologies, United States). The relative expression of MVP mRNA was presented using the normalized method of 2−ΔΔCT with the endogenous control GAPDH. The primer sequences were as follows: MVP forward primer, 5′- CCCAACACTGCCCTCCATCTAAAG-3'; MVP reverse primer, 5′- ATCTCCACGACCTCCACTTCCTTC-3'; GAPDH forward primer, 5′- GTCTCCTCTGACTTCAACAGCG-3'; GAPDH reverse primer, 5′- ACCACCCTGTTGCTGTAGCCAA-3'.
Cellular Proliferation, Migration, and Invasion Assay
The CCK-8 (Dojindo, Japan) was employed to assess cell viability. Transfected PTC cells (1–1.5 × 103/well) were seeded into 96-well plates. The optical density (OD) in each well was detected at 450 nm and recorded on a microplate reader (SpectraMax Plus 384, Molecular Devices Corporation, United States). For the colony formation assay, transfected PTC cells (1 × 103/well) were seeded into a 6-well plate and then incubated for 5–8 days till clonogenicity. The plates were then softly washed by PBS and stained with crystal violet. Colony fields were calculated by the ColonyArea software (Guzmán et al., 2014). Transwell plates (Corning, United States) were applied to the migration test. Cells (4 × 104 cells/chamber) were seeded onto the upper chamber, and the growth medium with 10% FBS was supplied to the bottom chamber. After culturing for 24 h in 37°C incubation, non-migrated cells in the upper chamber were cautiously wiped out using a cotton swab. Migrated cells were then fixed with methanol and stained with 0.1% crystal violet. For the invasion assay, the experimental procedure was similar to the migration test; in addition, the chamber is replaced by the Matrigel® invasion chamber (Corning, United States). Migrated or invasive cells were photographed in a 10 × magnification microscope in at least 5 randomly selected fields for each well and counted using ImageJ software. For the wound healing assay, transfected cells were plated in 6-well plates, and the confluent monolayer was scratched after incubation for 24 h.
Flow Cytometry Assay
The PTC cell lines were collected after being transfected for 72 h. Collected cells were washed with PBS in triplicate, and then resuspended with 300l 1binding buffer. Annexin V-fluorescein isothiocyanate (5l) and propidium iodide (5l) (BD Biosciences, United States) were added to 300l of cell suspensions at room temperature for 15 min in the dark. The apoptosis rate was defined as the percentage of Q2+ Q3 and analyzed by FlowJo (Tree Star, United States).
Western Blotting
Transfected PTC cells were lysed in RIPA buffer (Solarbio, China), and phenylmethylsulfonyl chloride was used as a protease inhibitor to stabilize the whole lysate. The extracted proteins were quantified by the bicinchoninic acid assay (Thermo Scientific, United States). Then the proteins were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (BioRad, United States) followed by transferring them to the polyvinylidene difluoride (PVDF) membranes (Millipore, United States). The primary antibodies were as follows: MVP (16478-1-AP, Proteintech), phospho-AKTSer473 (4060S, Cell Signaling Technology), total-AKT (4691S, Cell Signaling Technology), phospho-mTOR (381557, Zen Bioscience), total-mTOR (380411, Zen Bioscience), Phospho-p44/42 (4370T, Cell Signaling Technology), total-p44/42 (4695T, Cell Signaling Technology), phospho-p38 (4511T, Cell Signaling Technology), total-p38 (8690T, Cell Signaling Technology), and -Actin (AP0060, Bioworld Technology). Primary antibodies were used for immunoblotting at 1:1,000 dilution. The membranes were then incubated with a secondary antibody (ab97047 or ab6728, Abcam). Eventually, proteins were detected by the chemiluminescence kit (Thermo Scientific), and images of the protein bands were quantified by ImageJ software (NIH, United States).
Statistical Analysis
Univariate and multivariate Cox analyses were used to search for independent risk factors. The MannWhitney test or Wilcoxon signed-rank test was adopted for comparisons between two groups. The KruskalWallis one-way analysis of variance (ANOVA) was used in multiple groups. The predictive performance of MVP was assessed by the receiver operating characteristic (ROC) curve analysis. The KaplanMeier curve analysis and log-rank test were applied to calculate the prognostic results. The best MVP critical point for KaplanMeier curves was determined by res. cut function. The CCK-8 assay was analyzed by two-way analysis of variance. In all experiments, at least three biological replicates were performed for each group. R 4.0.0 and Graphpad Prism 8.3.0 were recruited in our work.
Results
Major Vault Protein Was Upregulated in Papillary Thyroid Cancer
To identify molecules that are differentially expressed in BRAF-driver PTC, we analyzed TCGA-PTC datasets for gene expression differences between BVL and RL PTCs. The t-SNE analysis revealed that patients in two PTC subgroups were distributed in discrete directions (Supplementary Figure S1A). Notably, MVP mRNA was markedly higher in BVL-PTC than in RL-PTC (Supplementary Figure S1B,C).
For the patients in the TCGA-THCA database, the MVP expression in PTC tissues was significantly higher than that in ANT (Figure 1A, p < 0.0001). This result was consistent with the conclusions drawn from the GSE60542, GSE35570, and GSE33630 cohorts (Figures 1B–D, all p < 0.0001). To verify the MVP expression at the transcriptional level in the local PTC cohort (WMU-PTC), we sequenced 79 matched PTC and ANT, finding that the MVP expression in PTC was remarkably upregulated (Figure 1E, p < 0.0001). All the data mentioned above indicated that MVP was upregulated in PTCs. The representative immunohistochemical images and detailed information also verify that the MVP protein expression in PTC was upregulated (Figure 1F).
FIGURE 1. MVP expression levels in PTC and normal thyroid tissues. The MVP mRNA expression in PTC and ANT obtained from the (A) TCGA-THCA, (B) GSE60542, (C) GSE35570, and (D) GSE33630 datasets. (E) MVP expression in local WMU-PTC samples and matched ANT by RNA-Seq. (F) Representative IHC images and detailed MVP profiles in malignant tissues and ANT were acquired from THPA database.
Major Vault Protein Expression Was Positively Associated With the Progression of Papillary Thyroid Cancer
The following analysis illustrated the relationship between MVP expression at mRNA level in PTC patients and clinicopathological characteristics in the TCGA dataset and our local patients. For the TCGA-THCA cohort, the MVP mRNA expression in PTC was considerably correlated with tumor size (Figure 2A, p = 0.00052), lymph node status (Figure 2B, p = 2.2e-10), extrathyroidal invasion (Figure 2C, p = 9.3e-09), and pathologic stages (Figure 2D, p = 3.7e-07). Moreover, the MVP level of follicular variant subtypes was lower than that of tall cell variant and classic subtypes (Figure 2E, p < 2.2e-16). We further analyzed the sequencing data of our local PTC patients, finding a positive relationship between the MVP level and LNM (Figure 2G, p = 0.0002). Still, there is no statistical difference between MVP expression and T stage (Figure 2F, p = 0.78) as well as disease stage (Figure 2H, p = 0.34).
FIGURE 2. MVP level and relationship with clinical and pathological features of PTC patients. (A-E) For the TCGA cohort, the MVP level in distinct (A) T stage, (B) lymph node status, (C) extrathyroidal extension, (D) pathologic stage, and (E) histological subtypes. (F-H) For WMU-PTC dataset, MVP level in distinct (F) T stage, (G) lymph node status, and (H) pathologic stage.
Association Between Major Vault Protein and Molecular Characteristics
BRAFV600E caused persistent activation of the MAPK signaling pathway, resulting in unlimited cell proliferation and tumor formation (Xing, 2005). Moreover, a large number of studies had also authenticated that genes such as RAS (Liu et al., 2008), TERT (Xing et al., 2014), and RET (Romei et al., 2016) played a vital role in the tumorigenesis and progression of PTC. Therefore, we divided the patients based on driver mutation status and noticed that the MVP level was markedly increased in the BRAFV600E mutated group (Figure 3A, p < 2.2e-16) and TERT mutated group (Figure 3C, p = 0.037). Nevertheless, compared with the RAS mutated group, the RAS wild type has a higher level of MVP (Figure 3B, p < 2.2e-16). Additionally, compared with the RET wild-type subset, the MVP level in the RET fusion group was higher (Figure 3D, p = 0.019). According to 71 gene expression characteristics, PTC can be classified into BVL and RL (Cancer Genome Atlas Research, 2014). The MVP level in the BVL subgroup was considerably higher than that in the RL subgroup (Figure 3E, p < 2.2e-16). Nevertheless, there was no remarkable difference in the MVP mRNA level between BRAFV600E mutated and wild type in our local PTC patients (Supplementary Figure S2A). Interestingly, MVP had a strong negative correlation with VEGFA (Figure 3F, R = -0.65, p < 2.2e-16). On the contrary, MVP had no significant correlation with TMB (Figure 3G, R = 0.05, p = 0.27) and MSI scores (Figure 3H, R = 0.086, p = 0.06).
FIGURE 3. MVP was correlated with molecular characteristics in PTC. (A–E) MVP expression in PTC patients with different genotypic features: (A) BRAFV600E pattern, (B) RAS pattern, (C) TERT pattern, (D) RET pattern, and (E) BVL/RL. (F–H) Spearman’s correlations between MVP expression and (F) VEGFA, (G) TMB, and (H) MSI.
Major Vault Protein Could Function as a Valuable Biomarker for Papillary Thyroid Cancer
Taking the above findings into consideration, we conjectured that MVP might be a whole new biomarker of PTC. AUC values calculated through ROC curve analysis in TCGA, GSE35570, GSE60542, and GSE33630 databases, and our local validated cohort were 0.8940, 0.9750, 0.9667, 0.9501, and 0.8945, respectively (Figure 4A), which indicated that MVP had a diagnostic value. In addition, MVP also had predictive capability for LNM in the TCGA dataset (AUC = 0.6742) and local validation group (AUC = 0.7440, Figure 4B). For tumor size and disease stages, the AUC values of TCGA were 0.6086 and 0.6168, respectively (Figures 4C,D). Our results demonstrated that MVP might be a valuable disease biomarker of PTCs.
FIGURE 4. MVP is a promising diagnostic and prognostic factor in PTC. (A) ROC curve analysis revealed the diagnostic potential of MVP in PTC patients derived from public datasets and local WMU-PTC cohort. (B-D) ROC curve analysis exhibited the prediction capabilities of MVP in (B) LNM status, (C) tumor size, and (D) disease stage. (E, F) Survival curves contrast MVPhigh and MVPlow patients together with the log-rank test. The p-value and hazard ratios (HR) with 95% confidence interval (CI) are calculated. (G) Multivariable Cox regression analyses of PFS in PTC patients. (H, I) Kaplan–Meier curves contrast MVPhigh and MVPlow patients with different (H) BRAFV600E phenotype and (I) lymph node status. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Then, we further speculated that MVP might be a promising prognostic predictor in PTC patients. Survival analysis revealed that high-level MVP was dramatically correlated with worse PFS in PTC patients (Figures 4E,F, HR: 1.8; 95% CI: 1.02–3.17; p = 0.033). By univariate Cox regression analysis, we determined that pathologic stages, T stage, extra invasion, and lymph node status were remarkably correlated with PFS (Supplementary Table S6). Yet the multivariate Cox regression analysis of PFS showed no independent prognostic factors (Figure 4G). Moreover, further survival analysis exhibited that high-level MVP was markedly correlated with shorter PFS in PTCs with BRAFV600E mutation or lymph node metastasis (Figure 4H, HR: 2.54; 95% CI: 1.22–5.25; p = 0.02; Figure 4I, HR: 2.41; 95% CI: 1.09–5.29; p = 0.014), but not in the BRAFV600E wild-type group or PTCs with negative lymph node metastasis. In our local WMU-PTC cohorts, survival analysis indicated that high MVP was associated with worse PFS (Supplementary Figure S2B, p = 0.13). The results indicated that MVP might affect the outcome of aggressive PTC patients.
Predicted Functions and Pathways of Major Vault Protein
In particular biological circumstances, genes usually act as drivers to manipulate potential downstream pathways and trigger specific molecular functions (Jia et al., 2016). On this basis, we performed the co-expression analysis using the data of the TCGA cohort to clarify the genes that have intimately interacted with MVP. Subsequently, we picked out certain genes with significant co-expression correlation with MVP for GO and KEGG analyses (Figures 5A,B; Supplementary Table S7, S8). GO ontology and KEGG enrichment analyses showed that MVP was principally related to immune-related functions and the activation, adhesion, and proliferation of T cell. Additionally, GSEA analysis uncovered that the high-level MVP was mainly enriched in distinct oncogenic pathways and immune-related terms like IL-2 and IL-6 signaling, INF-α and INF-γ response, allograft rejection, and P53 pathway (Figures 5C,D; Supplementary Table S9). These results revealed that MVP might play an essential role in the progression and tumor-immunity process of PTC.
FIGURE 5. Summary of MVP co-expressed genes and functional enrichment analyses. Gene ontology (A) and KEGG pathways (B) were considerably related to MVP expression. Genes strongly co-expressed with MVP are exhibited in order of Spearman’s correlation coefficient. (C) GSEA presented crucial hallmarker terms enriched in high (up) MVP expression subgroup. (D) GSEA exhibited key oncogenic terms enriched in high (up) and low (down) MVP expression subgroup. The cutoff value between high (up) and low (down) level MVP is the median value in the cohort. Entire lists are shown in Supplementary Table S5–S7.
Major Vault Protein Interrelated With Tumor-Infiltrating Immune Cells and Cancer-Immunity Cycle
It has been expounded that TIICs widely existed in the TME of PTC and affected pathological processes, including tumorigenicity. Stromal and immune scores were used to describe the proportion of stromal and immune cells in the TME, respectively (Mao et al., 2013). PTCs in the TCGA database were subdivided into MVPhigh and MVPlow sets according to the MVP median level. We quantitatively evaluated the immune and stromal infiltration by the ESTIMATE algorithm. We found MVP was correlated with the immune score (R = 0.65), TME score (R = 0.57), and tumor purity (R = -0.57) (Figure 6A). Similarly, utilizing another algorithm (xCell), we also found an apparent correlation between MVP expression and immune score (R = 0.65), as well as the TME score (R = 0.39) (Figure 6B). Interestingly, MVP expression was positively correlated with the ESTIMATE stromal score (R = 0.37) but negatively correlated with the xCell stromal score (R = -0.42). To sum up, the above analysis indicated that MVP had a positive relationship with the PTC immune microenvironment. With the purpose of further investigating the potential association between MVP expression and TIME, we established two different algorithms, CIBERSORT and quanTIseq, to calculate the abundance of TIICs. The clinicopathological features and infiltration level of TIICs based on the two algorithms are shown in Figure 6C.
FIGURE 6. MVP interacted with TIME phenotype. (A) Spearman’s correlations between MVP level and immune, stromal, microenvironment scores, and tumor purity by ESTIMATE. (B) Spearman’s correlations between MVP level and immune, stromal, and microenvironment scores by xCell. (C) Profiles of TIICs in PTC immune microenvironment assessed through CIBERSORT and quanTIseq algorithms.
Next, we attempted to determine whether the TIME changed in PTC samples with different expression levels of MVP (Figures 7A–C). By quantifying the immune infiltrating cellular component using three different algorithms (quanTIseq, MCPcounter, and CIBERSORT), we found that PTCs with high MVP expression were dramatically correlated with the high infiltrating level of anti-cancer CD8+ T cells and pro-cancer regulatory T cells. Simultaneously, more anti-cancer follicular helper T cells existed in MVPhigh PTC.
FIGURE 7. Profiles of immune cell infiltration in MVPhigh and MVPlow PTC from TCGA dataset. Compare the infiltrating immunocyte scores calculated by (A) quanTIseq, (B) MCP-counter, and (C) CIBERSORT in MVPhigh and MVPlow groups. (D) Compare each step of activity of cancer-immunity cycle between MVPhigh and MVPlow groups. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
The cancer-immunity cycle could trigger successive molecular behaviors, which could initiate the anti-cancer immunologic reaction to executing cancer cells efficiently (Chen and Mellman, 2013). As exhibited in Figure 7D, killed tumor cells first release antigens, and the antigens are subsequently captured by antigen-presenting cells, which could present specific antigens to T cells and induce the priming and activation of effector T-cell responses (steps 1 to step 3). MVPhigh PTCs were more immunocompetent than MVPlow in step 3. We noticed that MVPhigh PTCs had a significantly higher immune score in trafficking and infiltration of T cells to tumors (step 4 to step5). The MVPhigh group was significantly more active in the recruitment of CD4+ T cells, CD8+ T cells, Th1, dendritic cell, Th22, macrophage, monocyte, neutrophil, NK cells, eosinophil, basophil, B cells, Th2, and MDSC cells. Interestingly, we noticed that the process of recognition of tumor cells (step 6) was enhanced in MVPlow PTC. Taken together, high MVP expression was connected with anti-cancer immune activity and was instrumental in the regulation of the cancer-immunity cycle.
Deregulation of Major Vault Protein Promotes Papillary Thyroid Cancer Cell Proliferation, Migration, Invasion, and Apoptosis In Vitro
As a validation of the above bioinformatic analyses, a series of cell biological experiments were designed and performed in three PTC cell lines in vitro. We observed that the MVP expression of PTC cell lines (BCPAP, KTC-1, and TPC-1) was significantly increased by contrast with normal HTori-3 cells (Figure 8A). After that, we attempted to elucidate the effects caused by deregulating MVP.
FIGURE 8. MVP was over-expressed in PTC cell lines, promoting the capability of proliferation, colony, migration, and invasion but attenuated apoptosis in vitro. (A) Relative transcript level of MVP was substantially higher in PTC cell lines (BCPAP, TPC-1, and KTC-1) than that in normal HTori-3 cell line. (B) Transcript levels of MVP in si-NC, siMVP-1, siMVP-2, and siMVP-3 groups were determined by qRT-PCR. (C) Protein levels of MVP in si-NC and siMVP-1 PTC cells were determined by Western blotting. (D) Cell proliferative capacity was detected through CCK-8 assays. (E) Clonogenic ability was determined through Colony formation assays. (F) Cell apoptosis was detected through flow cytometry. Early and late apoptotic cells were respectively characterized as PI–/Annexin V+ and PI+/Annexin V+. (G) Barplots showed repeated experimental results in Colony formation assays and flow cytometry. (H) Migration and invasion abilities were detected through Transwell assays. (I) Migration capability was determined by wound healing assays. All assays were independently repeated at least three times. Data are presented as the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
On the one hand, we already had demonstrated that MVP was correlated with tumor size and disease stage in PTC (Figures 2A,D,F,H). Aiming to figure out the role of MVP during the development of PTC, MVP knockdown experiments were performed, and three PTC cell lines were processed by siRNA targeting MVP (si-MVP). Through qRT-PCR evaluation at the mRNA level, we determined that the MVP transcriptional level was reduced in the PTC cell lines knocked down by si-MVP#1 (Figure 8B). Similarly, we determined that MVP protein activity in the PTC cell lines knocked down by si-MVP#1 was significantly decreased by Western blotting (Figure 8C). We observed that knockdown of MVP attenuated cell proliferation and defects in colony formation (Figures 8D,E,G). Using the Annexin V-PI assay to detect early and late apoptotic cells in the PTC cell population, we found an increased apoptosis rate in PTC cells after silencing MVP (Figures 8F,G).
On the other hand, we had demonstrated that the MVP level was correlated with LNM status and extrathyroidal extension in PTC (Figures 2B,C,G). Therefore, to explore the effect of deregulated MVP toward metastasis, a transwell assay was run to evaluate the migratory and invasive activity of PTC cells, which demonstrates that the migration and invasion abilities of PTC cells were hindered under MVP depletion (Figure 8H). Furthermore, it was proven by wound healing assays that silencing MVP restrained the migratory ability of PTC cells (Figure 8I). Our data revealed that MVP participated in regulating proliferation, migration, and invasion abilities of PTC cells in vitro.
Major Vault Protein Promotes the Activity of MAPK/ERK and PI3K/AKT/mTOR Cascades
Earlier research on mechanisms has shown that MAPK/ERK and PI3K/AKT/mTOR cascades functioned as regulatory effects in the development and pathogenesis of thyroid cancer. To illustrate MVP molecular functions on the above signal pathways, we further evaluated the activities of pathways through quantitative analysis of the images obtained by Western blotting. On the one hand, knocking down MVP in PTC cells inhibited activities of the MAPK/ERK pathway, characterized by reduced phosphorylation of ERK1/2 (p44/42) and p38 MAPK (Figure 9A). On the other hand, MVP knocking down in PTC cells reduced phosphorylation of AKT at Ser473 and mTOR (Figure 9B). These results suggested that MVP functions as an oncogene through promoting the activities of both MAPK/ERK and PI3K/AKT/mTOR signaling pathways.
FIGURE 9. MVP promotes the activities of the MAPK/ERK and PI3K/AKT/mTOR cascades. The antibodies against MVP, phospho-p44/42 (p-p44/42), total p44/42 (t-p44/42), phospho-p38 (p-p38), total p38 (t-p38), phospho-AKTSer473 (p-AKTSer473), total AKT (t-AKT), phospho-mTOR (p-mTOR), and total mTOR (t-mTOR) were selected to evaluate the impact of MVP knockdown on the activities of (A) MAPK/ERK and (B) PI3K/AKT/mTOR cascades. GAPDH was used as the loading control. Representative images of Western blotting are showed in the upper panel. Quantitative analysis of Western blotting is displayed in the lower panel. All assays were independently repeated at least three times. Data are presented as the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Rescue Experiment Indicated Major Vault Protein Maintains the Survival of Papillary Thyroid Cancer Cells
To further confirm the oncogenic role of MVP, the expression of MVP was rescued in BCPAP cells which had MVP silenced. First, qPCR and Western blotting assays showed both transcriptional (Figure 10A) and protein levels of MVP (Figures 10B,C) were rescued in MVP silenced BCPAP cells. Second, forced overexpression of MVP could significantly unmask the inhibition of cell survival (CCK-8 and colony formation assays) (Figures 10D,E). Forced overexpression of MVP in these initially silenced BCPAP cells significantly reversed the level of migration compared with the control group (Figure 10F). Taken together, the rescue experiments showed that MVP could maintain PTC cell survival and induce cells migration. Rescue experiments prove that si-MVP#1 is specific and there are no off-target effects.
FIGURE 10. Rescue experiment indicated MVP maintains the survival and migration of PTC cells. (A–C) MVP expression was rescued in MVP silenced BCPAP cells. (A) qPCR and (B, C) Western blotting assays was performed to detect relative expression level of MVP. (D, E) CCK-8 and colony formation assays were performed to test cell proliferation of these BCPAP cells. (F) Transwell assays were performed to test cell migration activity of these BCPAP cells with restored expression of MVP. All assays were independently repeated at least three times. Data are presented as the mean ± SD. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Discussion
The prognosis of most PTC patients is favorable with surgery and postoperative thyroxine replacement therapy. However, for dedifferentiated, recrudescent, and metastatic lesions, their median survival time is less than 6 months as radiotherapy and chemotherapy are the few treatment modalities available. A study had clarified a strong connection between TIME and immunotherapy in PTC (Na and Choi, 2018). Therefore, to ameliorate the prognosis of highly invasive thyroid cancer and develop more optimized treatment strategies, exploring the underlying signaling mechanisms becomes extremely urgent.
Vault is the largest cellular ribonucleoprotein complex with a hollow barrel structure, which comprises 3 protein components, namely, major vault protein (MVP), telomerase-associated protein 1 (TEP1), and vault poly ADP ribose polymerase (VPARP) (Park, 2012). Recent research works illustrate the vault is instrumental in a series of cellular processes, covering nuclear pore assembly, subcellular transport, signal transduction, and interferon response (Steiner et al., 2006; Berger et al., 2009; Ryu and Park, 2009; Vollmar et al., 2009). Recently, the literature has shown that MVP played an important part in cancer progression. MVP upregulation promotes oncogenesis and development of multiple tumor types (Xiao et al., 2019). Contradictively, a recent study showed that MVP inhibited lung cancer cell proliferation by suppressing the STAT3 pathway (Bai et al., 2019). In addition, MVP promotes glioblastoma survival and migration (Lötsch et al., 2013a), and inhibits apoptosis of human senescent diploid fibroblasts (Ryu et al., 2008) and human colorectal cancer cells (Ikeda et al., 2008). According to the information we have, this is the first work to dissect the clinical application value and potential molecular role of MVP in PTC from an entirely new perspective.
In this study, based on the analysis of the available data of 658 PTC patients derived from four public databases, we identified that MVP was significantly upregulated in PTC as compared to ANT. Subsequently, we took 79 paired PTC specimens and their corresponding ANT as a validation group and got the same conclusion by analyzing their RNA-Seq data. PTCs with a high pathologic grade or poorly differentiated still involve a certain amount of risk for disease recurrence and metastasis, and even endanger their lives (Farahati et al., 2004). Analysis of the TCGA database revealed that the increased MVP expression was consistent with more malignant clinicopathological features, involving larger tumor, more extensive LNM, and higher tumor stage, and extrathyroidal invasion implied that MVP might be bound up with the progression of PTC. In the local cohort, we also found that MVP is related to LNM.
According to the latest clinical diagnosis and treatment guidelines (ATA 2015) (Haugen et al., 2015), tissue/serum markers are not recommended for distinguishing benign and malignant thyroid nodules. The clinical differentiation of benign and malignant thyroid nodules mainly relies on pathological biopsy and tissue/cytological biopsy. A machine learning model has established a similar model on the TCGA database, and the discrimination between normal tissue and thyroid tissue exceeds 0.95 (Park et al., 2020). Our research group also published a model with good discrimination performance in identifying normal thyroid and thyroid cancer tissues (Wang et al., 2016). ROC curve analysis indicated that MVP had an accurate distinguishing ability for tumors in multi-center datasets (all AUC >0.89), but its actual ability to be used in the clinic still needs further research such as prospective clinical trials. In addition, MVP also had a specific predictive value for T-stage, LNM status, and disease stages of PTC patients in the TCGA cohort. In short, the above evidence revealed that MVP had the potential to serve as a promising biomarker of PTC.
BRAFV600E is a principal human oncogene originally reported in melanoma and incrementally discovered in other cancers, especially in PTC. Apart from clinicopathological characteristics, several genic features (such as BRAF, TERT, and RET) can be applied for risk stratification and prognostic scoring systems of thyroid cancer patients (Xing et al., 2005; Lee et al., 2007; Kim et al., 2012). PTC predominantly involves contradictory driving factors, and its comparative monotonicity of the genome brought about piles of studies which in detail analyzed the fundamental difference in genomics, epigenomics, and proteomics between BVL-PTCs and RL-PTCs. A study exposited that the existence of BRAFV600E mutation was the sole clinicopathological predictor of persistent disease after following up for 5 years (Elisei et al., 2012). BRAFV600E mutation is the most frequent mutation of BRAF and signifies high diagnostic accurate PTC, which can independently activate the downstream MER-ERK pathway and then lead to high MAPK signaling. A single mutation in RAS tends to indicate low risks, including benign nodules, NIFTP, and low-danger cancer. However, combined mutations in RAS usually mean increased malignancy. We found that MVP expression was elevated in BRAFV600E mutated, RAS wild-type, and BVL genotypes. We did not find the difference in the MVP level between BRAFV600E mutated and wild-type groups in our local PTC patients, which may be due to the low sample number in the wild-type group. TERT promoter mutations readily occurred in invasive thyroid cancer and were a strong predictor of poor clinical outcomes in thyroid tumor (Liu and Xing, 2016). RET fusions generally presented the BVL phenotype and upregulated the MAPK signaling (Mitsutake et al., 2006). We observed that the MVP level was substantially elevated in TERT mutated and RET fusion groups. Vascular endothelial growth factor (VEGF) promotes angiogenesis and endothelial cell proliferation, and is correlated with tumor growth and aggressive behavior of thyroid cancer (Salajegheh et al., 2011; Li et al., 2014). Correlation analysis of available RNA-Seq profiles suggested a strong negative correlation between MVP expression and VEGFA. Recently, Ronald J Hause et al. found MSI existed in varieties of cancers through whole-exome sequencing of 5,930 genomes from 18 types of cancers, particularly with regard to colorectal, endometrial, gastric, and thyroid cancers (Hause et al., 2016; Genutis et al., 2019). Earlier literature illustrated that cancers with high-frequency MSI often showed up as poor differentiation, expansive growth, histological heterogeneity, and high-level TIICs (Zheng et al., 2020). Interestingly, although SNPs in the BRAF and RAS were strongly associated with MVP expression, no significant correlation was found between MVP level and TMB and MSI scores in PTCs. Recently, research based on proteomic and bioinformatic analyses identified MVP as a prognostic biomarker for fatal prostate cancer (Ramberg et al., 2021). There were molecular characteristics that proved to be associated with MVP expression, and a high MVP level was associated with poor prognosis. However, we found MVP was not an independent prognostic factor of PTC. It might attribute to the favorable outcome of PTC patients, which causes the number of progression events to be small. These clinical analyses suggested that MVP may involve in the occurrence and progression of PTC.
Based on the transcriptional-level analysis, we then investigated the potential molecular function of MVP in PTC. Co-expressed genes commonly function coordinatively in biological molecular signals regulated by various factors and play a positive role in adaptive evolution (Niehrs and Pollet, 1999). GO ontology and KEGG enrichment analyses showed that MVP was correlated with the enriched immune-related genes. Hallmark terms related to immunologic processes or tumorigenesis pathways (e.g., IL-2 and IL-6 signaling, INF-α and INF-γ response, and P53 pathway) were enriched in GSEA analyses and associated with the high MVP level. A comprehensive dissection of TIME landscapes could contribute to identifying novel immunotherapeutic targets in PTC patients. Na et al. summarized immune cell abundances as an immune score to describe the profiles of the TIME in PTC (Na and Choi, 2018). Therefore, relying on two algorithms (ESTIMATE and xCell), we made use of transcriptional profiles of PTC tissues to quantify the immune and stromal proportions so that we can figure out the relationship between MVP and TME. Our research observed that the MVP expression level was positively related to immune scores, while stromal scores exhibited diametrically opposite results in two algorithms. In the light of the results of functional enrichment analysis, we found MVP might participate in regulating immune cell infiltration in PTC immune microenvironment.
Deconvolution algorithms are novel for decomposing heterogeneous cellular admixtures in the TME (Gentles et al., 2015; Orhan et al., 2020). In addition, MCP-counter quantifies the absolute abundance of 8 immune and 2 stromal cell clusters in heterogeneous tissues using transcriptome information (Becht et al., 2016b). Because of no criterion for assessing immune infiltration from RNA-seq profiles, we employed the deconvolution-based CIBERSORT and quanTIseq algorithms, and the marker gene-based MCP-counter algorithm to calculate the level of TIIC components quantitatively. Now there is evidence that high infiltration of CD8+ T cells was related to increased disease-free survival in thyroid cancer (Schreiber et al., 2011; Cunha et al., 2012). However, another retrospective study showed an opposite result that infiltration of CD8+ T cells was positively related to recurrence risk in differentiated thyroid cancer (Cunha et al., 2015). Recently, a study elucidated that immunotherapy was more effective against cancers with high CD8+ T-cell infiltration (Farhood et al., 2019). Our research found that MVPhigh status was associated with the increased number of anti-cancer CD8+ T cells, pro-cancer regulatory T cells, and anti-cancer follicular helper T cells in PTC.
The cancer-immunity cycle is a multi-step system in which the immune system distinguishes and eliminates tumor cells, and each step has corresponding positive or negative regulatory factors (Chen and Mellman, 2013). Immunotherapy is to activate the entire immune cycle by targeting these regulators at different phases and ultimately achieve the purpose of treatment. For example, the CTLA-4 inhibitor mainly acts on priming and activation (step 3), while PD-L1 and PD-1 inhibitors target step 7 (Chen and Mellman, 2017). Notably, our results showed that MVPhigh PTC exhibited higher overall anti-cancer immune scores in the cancer-immunity cycle. Furthermore, a variety of immune infiltrated cell components were increasingly identified in PTC with a high level of MVP. We also noticed that MVPhigh PTC was more vigorous in recruiting CD8+ T cells, CD4+ T cells, Th1, dendritic cells, Th22, macrophage, monocyte, neutrophil, NK cells, eosinophil, basophil, B cells, Th2, and MDSC cells, which was consistent with the alterations of TIICs in TIME. In summary, our findings elucidate that MVP may affect and regulate the cancer-immunity cycle.
Some studies reported that TIICs in TME might impact the proliferation, migration, and drug resistance of tumor cells (Gocheva et al., 2010; Joyce and Fearon, 2015). So far, there is no relevant research to determine the biological functions of MVP in PTC. Considering that transfection efficiency and genotypes varied among different cancer cells, we conducted a series of in vitro experiments after the knockdown of the MVP in three different PTC cell lines and expectedly elaborated that downregulated MVP inhibited the proliferation, migration, and invasion abilities of PTC cells. In the gain-of-function experiments, overexpression of MVP increases the proliferation and migration abilities of PTC cells. The effects of MVP knockdown on cell proliferation and migration were rescued in part by overexpression of MVP. In a nutshell, these lines of evidence demonstrated that MVP played a cancer-promoting role in PTCs.
Recently, the mechanisms of genetic events and cancer-related signaling pathways in thyroid cancer have attracted increasing attention. It has been clarified that MAPK/ERK and PI3K/AKT/mTOR cascades robustly participate in the tumorigenesis and development of thyroid cancer (Nikiforov and Nikiforova, 2011; Xing, 2013). Our findings showed that BRAF mutation was significantly correlated with MVP expression, whereas RAS mutation was connected with MVPlow status. Several studies had demonstrated that MAPK/ERK and PI3K/AKT/mTOR pathways played critical roles in the tumorigenesis and progression of PTC. BRAF mutations primarily affected downstream signaling of the MAPK/ERK pathway, while RAS mutations mainly influenced the MAPK and PI3K/AKT/mTOR cascades (Nikiforov and Nikiforova, 2011; Xing, 2013). Therefore, to explore the cancer-promoting effect of MVP more deeply, we further discussed the impact of MVP in the above pathways. After the knockdown of MVP, we observed that the levels of p44/42 and p38 were obviously lower, indicating that MVP exerted its biological functions through the MAPK/ERK signaling pathway. A previous study on TNBC reached a similar conclusion as well (Xiao et al., 2019). Similarly, we observed that protein p-AKTser473 and p-mTOR decreased after the knockdown of MVP, revealing that MVP may exert its effects by activating the PI3K/AKT/mTOR cascade. Our results were consistent with the relevant reports of MVP in other cancers (Lötsch et al., 2013b; Xiao et al., 2019). The further rescue experiments also validated that MVP could maintain PTC cells’ survival and induce PTC cells’ migration.
Several limitations in the current research still required to be considered. First, larger sample size information needs to be collected for verification of the association between MVP and clinicopathological tumor parameters, and also the prognostic value of MVP. After that, this study demonstrated that MVP interacted with anti-cancer immune response, neoplasm formation, and metastasis of PTC. However, potential MVP-related molecular functions in PTC progression require further exploration. In vivo experiments are necessary to further complement the regulatory mechanisms of MVP. Moreover, the gain-of-function experiments are also needed to explore if overexpressed MVP increases the proliferation and migration abilities in the HTori-3 cell line, which can help us understand the role of the MVP in PTC progression. In the end, the crosstalk between mutation drivers and altered tumor immune landscape in MVP-induced PTC lesions needs to be determined.
Conclusion
Taken together, the biological effects and molecule mechanisms of MVP in PTC seem to be complicated. This work is the first comprehensive research to determine the expression level, relationship with clinical parameters, pathological characteristics, molecular hallmarks, TIICs, and cancer-immunity cycle functions of MVP in PTC. Our work first revealed that MVP is a reliable immune microenvironment-related biomarker, which also functions as a potential oncogene in PTC through activating the PI3K/AKT/mTOR and MAPK/ERK cascades. More in-depth studies of the MVP will contribute to our understanding of the mechanism of PTC progression.
Data Availability Statement
The data that support the findings of this study have been deposited into CNGB Sequence Archive (CNSA, https://db.cngb.org/cnsa/) of China National GeneBank DataBase (CNGBdb) with accession number CNP0002696. The original contributions generated for this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
All research protocols have been approved and implemented through the ethical standards of the Institutional Review Board of the First Affiliated Hospital of Wenzhou Medical University (Approval No. 2012-57). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
XD contributed to study design, bioinformatic analysis, and manuscript draft. XD, PA, JS, and CN contributed to molecular biology experiments. CJ contributed to the manuscript draft. HL, WJ, and YS collected and analyzed clinical samples for the study. XZ, QZ, and GH discussed the results and participated in the critical review of the manuscript. All authors read and approved the final manuscript.
Funding
This study was supported by the funding of the National Natural Science Foundation of China (No. 81802328), Major Science and Technology Projects of Zhejiang Province (2015C03052), Wenzhou Science and Technology Plan Project (Y20180846, Y20190388, and Y20210258), Young Talents Program of the First Affiliated Hospital of Wenzhou Medical University (No. qnyc094), and Scientific Research Incubator Project of the First Affiliated Hospital of Wenzhou Medical University (FHY2019060).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We appreciate all team members of TCGA-THCA, GSE33630, GSE60542, and GSE35570 for providing available public data.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.688370/full#supplementary-material
Supplementary Figure 1 | Identification of MVP as an important BRAF-driven oncogene in PTC. (A) t-SNE analysis of the TCGA cohort. (B) Volcano plot of differentially expressed genes between BVL and RL PTCs. (C) Heatmap of top 20 differentially expressed genes between BVL and RL PTCs.
Supplementary Figure 2 | Association of MVP level with BRAFV600E and prognosis in WMU-PTC patients. (A) MVP level in PTC patients stratified by BRAFV600E status. (B) Survival curves comparing patients with distinct MVP expression for progression-free survival.
Abbreviations
MAPK, mitogen-activated protein kinase; TME, tumor microenvironment; TIME, tumor immune microenvironment; TIICs, tumor-infiltrating immune cells; ESTIMATE, estimation of stromal and immune cells in malignant tumors using expression data; CIBERSORT, cell-type identification by estimating relative subsets of known RNA transcripts; MCP-counter, microenvironment cell population-counter; ssGSEA, single-sample gene set enrichment analysis; MSI, microsatellite instability; SNP, single-nucleotide polymorphism; TIP, tumor immunophenotype; PBS, phosphate-buffered saline; PI, propidium iodide; PCR, polymerase chain reaction; NC, negative control; CCK-8, Cell Counting Kit-8; LNM, lymph node metastasis; AUC, area under the curve; RSEM, RNA-seq by expectation-maximization; VEGFA, vascular endothelial growth factor A; MDSC, myeloid-derived suppressor cells; NIFTP, non-invasive follicular thyroid neoplasm with papillary-like nuclear features.
References
Agrawal, N., Akbani, R., Aksoy, B. A., Ally, A., Arachchi, H., Asa, S. L., et al. (2014). Integrated Genomic Characterization of Papillary Thyroid Carcinoma. Cell 159 (3), 676–690. doi:10.1016/j.cell.2014.09.050
Asioli, S., Erickson, L. A., Sebo, T. J., Zhang, J., Jin, L., Thompson, G. B., et al. (2010). Papillary Thyroid Carcinoma with Prominent Hobnail Features: a New Aggressive Variant of Moderately Differentiated Papillary Carcinoma. A Clinicopathologic, Immunohistochemical, and Molecular Study of Eight Cases. Am. J. Surg. Pathol. 34 (1), 44–52. doi:10.1097/PAS.0b013e3181c46677
Bai, H., Wang, C., Qi, Y., Xu, J., Li, N., Chen, L., et al. (2019). Major Vault Protein Suppresses Lung Cancer Cell Proliferation by Inhibiting STAT3 Signaling Pathway. BMC Cancer 19 (1), 454. doi:10.1186/s12885-019-5665-6
Becht, E., Giraldo, N. A., Lacroix, L., Buttard, B., Elarouci, N., Petitprez, F., et al. (2016). Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol. 17 (1), 218–220. doi:10.1186/s13059-016-1070-5
Becht, E., Giraldo, N. A., Lacroix, L., Buttard, B., Elarouci, N., Petitprez, F., et al. (2016). Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol. 17 (1), 218. doi:10.1186/s13059-016-1070-5
Berger, W., Steiner, E., Grusch, M., Elbling, L., and Micksche, M. (2009). Vaults and the Major Vault Protein: Novel Roles in Signal Pathway Regulation and Immunity. Cell. Mol. Life Sci. 66 (1), 43–61. doi:10.1007/s00018-008-8364-z
Bonneville, R., Krook, M. A., Kautto, E. A., Miya, J., Wing, M. R., Chen, H-Z., et al. (20172017). Landscape of Microsatellite Instability across 39 Cancer Types. JCO Precis Oncol. doi:10.1200/PO.17.00073
Cancer Genome Atlas Research, N. (2014). Integrated Genomic Characterization of Papillary Thyroid Carcinoma. Cell 159 (3), 676–690. doi:10.1016/j.cell.2014.09.050
Chen, D. S., and Mellman, I. (2017). Elements of Cancer Immunity and the Cancer-Immune Set point. Nature 541 (7637), 321–330. doi:10.1038/nature21349
Chen, D. S., and Mellman, I. (2013). Oncology Meets Immunology: the Cancer-Immunity Cycle. Immunity 39 (1), 1–10. doi:10.1016/j.immuni.2013.07.012
Cunha, L. L., Marcello, M. A., Nonogaki, S., Morari, E. C., Soares, F. A., Vassallo, J., et al. (2015). CD8+ Tumour-Infiltrating Lymphocytes and COX2 Expression May Predict Relapse in Differentiated Thyroid Cancer. Clin. Endocrinol. 83 (2), 246–253. doi:10.1111/cen.12586
Cunha, L. L., Morari, E. C., Guihen, A. C. T., Razolli, D., Gerhard, R., Nonogaki, S., et al. (2012). Infiltration of a Mixture of Immune Cells May Be Related to Good Prognosis in Patients with Differentiated Thyroid Carcinoma. Clin. Endocrinol. 77 (6), 918–925. doi:10.1111/j.1365-2265.2012.04482.x
Deschoolmeester, V., Baay, M., Van Marck, E., Weyler, J., Vermeulen, P., Lardon, F., et al. (2010). Tumor Infiltrating Lymphocytes: an Intriguing Player in the Survival of Colorectal Cancer Patients. BMC Immunol. 11 (1), 19. doi:10.1186/1471-2172-11-19
Dom, G., Tarabichi, M., Unger, K., Thomas, G., Oczko-Wojciechowska, M., Bogdanova, T., et al. (2012). A Gene Expression Signature Distinguishes normal Tissues of Sporadic and Radiation-Induced Papillary Thyroid Carcinomas. Br. J. Cancer 107 (6), 994–1000. doi:10.1038/bjc.2012.302
Dong, X., Yang, Q., Gu, J., Lv, S., Song, D., Chen, D., et al. (2021). Identification and Validation of L Antigen Family Member 3 as an Immune-Related Biomarker Associated with the Progression of Papillary Thyroid Cancer. Int. Immunopharmacol 90, 107267. doi:10.1016/j.intimp.2020.107267
Elisei, R., Viola, D., Torregrossa, L., Giannini, R., Romei, C., Ugolini, C., et al. (2012). TheBRAFV600E Mutation Is an Independent, Poor Prognostic Factor for the Outcome of Patients with Low-Risk Intrathyroid Papillary Thyroid Carcinoma: Single-Institution Results from a Large Cohort Study. J. Clin. Endocrinol. Metab. 97 (12), 4390–4398. doi:10.1210/jc.2012-1775
Farahati, J., Geling, M., Mäder, U., Mörtl, M., Luster, M., Müller, J. G., et al. (2004). Changing Trends of Incidence and Prognosis of Thyroid Carcinoma in Lower Franconia, Germany, from 1981-1995. Thyroid 14 (2), 141–147. doi:10.1089/105072504322880382
Farhood, B., Najafi, M., and Mortezaee, K. (2019). CD8 + Cytotoxic T Lymphocytes in Cancer Immunotherapy: A Review. J. Cel Physiol 234 (6), 8509–8521. doi:10.1002/jcp.27782
Finotello, F., Mayer, C., Plattner, C., Laschober, G., Rieder, D., Hackl, H., et al. (2019). Molecular and Pharmacological Modulators of the Tumor Immune Contexture Revealed by Deconvolution of RNA-Seq Data. Genome Med. 11 (1), 34–20. doi:10.1186/s13073-019-0638-6
Gentles, A. J., Newman, A. M., Liu, C. L., Bratman, S. V., Feng, W., Kim, D., et al. (2015). The Prognostic Landscape of Genes and Infiltrating Immune Cells across Human Cancers. Nat. Med. 21 (8), 938–945. doi:10.1038/nm.3909
Genutis, L. K., Tomsic, J., Bundschuh, R. A., Brock, P. L., Williams, M. D., Roychowdhury, S., et al. (2019). Microsatellite Instability Occurs in a Subset of Follicular Thyroid Cancers. Thyroid 29 (4), 523–529. doi:10.1089/thy.2018.0655
Ghossein, R., and Livolsi, V. A. (2008). Papillary Thyroid Carcinoma Tall Cell Variant. Thyroid 18 (11), 1179–1181. doi:10.1089/thy.2008.0164
Gocheva, V., Wang, H.-W., Gadea, B. B., Shree, T., Hunter, K. E., Garfall, A. L., et al. (2010). IL-4 Induces Cathepsin Protease Activity in Tumor-Associated Macrophages to Promote Cancer Growth and Invasion. Genes Dev. 24 (3), 241–255. doi:10.1101/gad.1874010
Guzmán, C., Bagga, M., Kaur, A., Westermarck, J., and Abankwa, D. (2014). ColonyArea: an ImageJ Plugin to Automatically Quantify colony Formation in Clonogenic Assays. PLoS One 9 (3), e92444. doi:10.1371/journal.pone.0092444
Haugen, B. R., Alexander, E. K., Bible, K. C., Doherty, G. M., Mandel, S. J., Nikiforov, Y. E., et al. (2015). 2015 American Thyroid Association Management Guidelines for Adult Patients with Thyroid Nodules and Differentiated Thyroid Cancer: The American Thyroid Association Guidelines Task Force on Thyroid Nodules and Differentiated Thyroid Cancer. Thyroid 26 (1), 1–133. doi:10.1089/thy.2015.0020
Hause, R. J., Pritchard, C. C., Shendure, J., and Salipante, S. J. (2016). Classification and Characterization of Microsatellite Instability across 18 Cancer Types. Nat. Med. 22 (11), 1342–1350. doi:10.1038/nm.4191
Henríquez-Hernández, L. A., Moreno, M., Rey, A., Lloret, M., and Lara, P. C. (2012). MVP Expression in the Prediction of Clinical Outcome of Locally Advanced Oral Squamous Cell Carcinoma Patients Treated with Radiotherapy. Radiat. Oncol. 7, 147. doi:10.1186/1748-717x-7-147
Ikeda, R., Iwashita, K.-i., Sumizawa, T., Beppu, S.-i., Tabata, S., Tajitsu, Y., et al. (2008). Hyperosmotic Stress Up-Regulates the Expression of Major Vault Protein in SW620 Human colon Cancer Cells. Exp. Cel Res. 314 (16), 3017–3026. doi:10.1016/j.yexcr.2008.07.001
Jia, Z., Liu, Y., Guan, N., Bo, X., Luo, Z., and Barnes, M. R. (2016). Cogena, a Novel Tool for Co-expressed Gene-Set Enrichment Analysis, Applied to Drug Repositioning and Drug Mode of Action Discovery. BMC genomics 17 (1), 414. doi:10.1186/s12864-016-2737-8
Johnson, T. L., Lloyd, R. V., Thompson, N. W., Beierwaltes, W. H., and Sisson, J. C. (1988). Prognostic Implications of the Tall Cell Variant of Papillary Thyroid Carcinoma. Am. J. Surg. Pathol. 12 (1), 22–27. doi:10.1097/00000478-198801000-00003
Joyce, J. A., and Fearon, D. T. (2015). T Cell Exclusion, Immune Privilege, and the Tumor Microenvironment. Science 348 (6230), 74–80. doi:10.1126/science.aaa6204
Kampf, C., Olsson, I., Ryberg, U., Sjöstedt, E., and Pontén, F. (2012). Production of Tissue Microarrays, Immunohistochemistry Staining and Digitalization within the Human Protein Atlas. J. Vis. Exp. (63). doi:10.3791/3620
Kim, T. H., Park, Y. J., Lim, J. A., Ahn, H. Y., Lee, E. K., Lee, Y. J., et al. (2012). The Association of the BRAFV600E Mutation with Prognostic Factors and Poor Clinical Outcome in Papillary Thyroid Cancer. Cancer 118 (7), 1764–1773. doi:10.1002/cncr.26500
Lamartina, L., Grani, G., Arvat, E., Nervo, A., Zatelli, M. C., Rossi, R., et al. (2018). 8th Edition of the AJCC/TNM Staging System of Thyroid Cancer: what to Expect (ITCO#2). Endocrine-Related Cancer 25 (3), L7–L11. doi:10.1530/erc-17-0453
Landa, I., Ibrahimpasic, T., Boucai, L., Sinha, R., Knauf, J. A., Shah, R. H., et al. (2016). Genomic and Transcriptomic Hallmarks of Poorly Differentiated and Anaplastic Thyroid Cancers. J. Clin. Invest. 126 (3), 1052–1066. doi:10.1172/jci85271
Lee, J.-H., Lee, E.-S., and Kim, Y.-S. (2007). Clinicopathologic Significance ofBRAF V600E Mutation in Papillary Carcinomas of the Thyroid. Cancer 110 (1), 38–46. doi:10.1002/cncr.22754
Li, J., Teng, L., and Jiang, H. (2014). Relationship between Preoperative Serum TSH Levels and Expression of VEGF in Papillary Thyroid Carcinoma. Asia-pac J. Clin. Oncol. 10 (2), 149–152. doi:10.1111/ajco.12075
Liu, J., Lichtenberg, T., Hoadley, K. A., Poisson, L. M., Lazar, A. J., Cherniack, A. D., et al. (2018). An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell 173 (2), 400–e11. doi:10.1016/j.cell.2018.02.052
Liu, R., and Xing, M. (2016). TERT Promoter Mutations in Thyroid Cancer. Endocr. Relat. Cancer 23 (3), R143–R155. doi:10.1530/erc-15-0533
Liu, Z., Hou, P., Ji, M., Guan, H., Studeman, K., Jensen, K., et al. (2008). Highly Prevalent Genetic Alterations in Receptor Tyrosine Kinases and Phosphatidylinositol 3-kinase/akt and Mitogen-Activated Protein Kinase Pathways in Anaplastic and Follicular Thyroid Cancers. J. Clin. Endocrinol. Metab. 93 (8), 3106–3116. doi:10.1210/jc.2008-0273
Liu, Z., Zhang, W., Phillips, J. B., Arora, R., McClellan, S., Li, J., et al. (2019). Immunoregulatory Protein B7-H3 Regulates Cancer Stem Cell Enrichment and Drug Resistance through MVP-Mediated MEK Activation. Oncogene 38 (1), 88–102. doi:10.1038/s41388-018-0407-9
Losert, A., Lötsch, D., Lackner, A., Koppensteiner, H., Peter-Vörösmarty, B., Steiner, E., et al. (2012). The Major Vault Protein Mediates Resistance to Epidermal Growth Factor Receptor Inhibition in Human Hepatoma Cells. Cancer Lett. 319 (2), 164–172. doi:10.1016/j.canlet.2012.01.002
Lötsch, D., Steiner, E., Holzmann, K., Spiegl-Kreinecker, S., Pirker, C., Hlavaty, J., et al. (2013). Major Vault Protein Supports Glioblastoma Survival and Migration by Upregulating the EGFR/PI3K Signalling axis. Oncotarget 4 (11), 1904–1918. doi:10.18632/oncotarget.1264
Mao, Y., Keller, E. T., Garfield, D. H., Shen, K., and Wang, J. (2013). Stromal Cells in Tumor Microenvironment and Breast Cancer. Cancer Metastasis Rev. 32 (1-2), 303–315. doi:10.1007/s10555-012-9415-3
Mitsutake, N., Miyagishi, M., Mitsutake, S., Akeno, N., Mesa, C., Knauf, J. A., et al. (2006). BRAF Mediates RET/PTC-induced Mitogen-Activated Protein Kinase Activation in Thyroid Cells: Functional Support for Requirement of the RET/PTC-RAS-BRAF Pathway in Papillary Thyroid Carcinogenesis. Endocrinology 147 (2), 1014–1019. doi:10.1210/en.2005-0280
Na, K. J., and Choi, H. (2018). Immune Landscape of Papillary Thyroid Cancer and Immunotherapeutic Implications. Endocr. Relat. Cancer 25 (5), 523–531. doi:10.1530/erc-17-0532
Newman, A. M., Steen, C. B., Liu, C. L., Gentles, A. J., Chaudhuri, A. A., Scherer, F., et al. (2019). Determining Cell Type Abundance and Expression from Bulk Tissues with Digital Cytometry. Nat. Biotechnol. 37 (7), 773–782. doi:10.1038/s41587-019-0114-2
Niehrs, C., and Pollet, N. (1999). Synexpression Groups in Eukaryotes. Nature 402 (6761), 483–487. doi:10.1038/990025
Nikiforov, Y. E., and Nikiforova, M. N. (2011). Molecular Genetics and Diagnosis of Thyroid Cancer. Nat. Rev. Endocrinol. 7 (10), 569–580. doi:10.1038/nrendo.2011.142
Orhan, A., Vogelsang, R. P., Andersen, M. B., Madsen, M. T., Hölmich, E. R., Raskov, H., et al. (2020). The Prognostic Value of Tumour-Infiltrating Lymphocytes in Pancreatic Cancer: a Systematic Review and Meta-Analysis. Eur. J. Cancer 132, 71–84. doi:10.1016/j.ejca.2020.03.013
Park, K. S., Kim, S. H., Oh, J. H., and Kim, S. Y. (2020). Highly Accurate Diagnosis of Papillary Thyroid Carcinomas Based on Personalized Pathways Coupled with Machine Learning. Brief. Bioinform. doi:10.1093/bib/bbaa336
Park, K. (2012). The Role of Major Vault Protein (MVP) in Drug Resistance. J. Controlled Release 163 (2), 266. doi:10.1016/j.jconrel.2012.10.005
Ramberg, H., Richardsen, E., de Souza, G. A., Rakaee, M., Stensland, M. E., Braadland, P. R., et al. (2021). Proteomic Analyses Identify Major Vault Protein as a Prognostic Biomarker for Fatal Prostate Cancer. Carcinogenesis 42 (5), 685–693. doi:10.1093/carcin/bgab015
Randolph, G. W., Duh, Q.-Y., Heller, K. S., LiVolsi, V. A., Mandel, S. J., Steward, D. L., et al. (2012). The Prognostic Significance of Nodal Metastases from Papillary Thyroid Carcinoma Can Be Stratified Based on the Size and Number of Metastatic Lymph Nodes, as Well as the Presence of Extranodal Extension. Thyroid 22 (11), 1144–1152. doi:10.1089/thy.2012.0043
Romei, C., Ciampi, R., and Elisei, R. (2016). A Comprehensive Overview of the Role of the RET Proto-Oncogene in Thyroid Carcinoma. Nat. Rev. Endocrinol. 12 (4), 192–202. doi:10.1038/nrendo.2016.11
Ryu, S. J., An, H. J., Oh, Y. S., Choi, H. R., Ha, M. K., and Park, S. C. (2008). On the Role of Major Vault Protein in the Resistance of Senescent Human Diploid Fibroblasts to Apoptosis. Cell Death Differ 15 (11), 1673–1680. doi:10.1038/cdd.2008.96
Ryu, S. J., and Park, S. C. (2009). Targeting Major Vault Protein in Senescence-Associated Apoptosis Resistance. Expert Opin. Ther. Targets 13 (4), 479–484. doi:10.1517/14728220902832705
Salajegheh, A., Smith, R. A., Kasem, K., Gopalan, V., Nassiri, M. R., William, R., et al. (2011). Single Nucleotide Polymorphisms and mRNA Expression of VEGF-A in Papillary Thyroid Carcinoma: Potential Markers for Aggressive Phenotypes. Eur. J. Surg. Oncol. (Ejso) 37 (1), 93–99. doi:10.1016/j.ejso.2010.10.010
Schreiber, R. D., Old, L. J., and Smyth, M. J. (2011). Cancer Immunoediting: Integrating Immunity's Roles in Cancer Suppression and Promotion. Science 331 (6024), 1565–1570. doi:10.1126/science.1203486
Shen, W., Qiu, Y., Li, J., Wu, C., Liu, Z., Zhang, X., et al. (2019). IL‐25 Promotes Cisplatin Resistance of Lung Cancer Cells by Activating NF‐κB Signaling Pathway to Increase of Major Vault Protein. Cancer Med. 8 (7), 3491–3501. doi:10.1002/cam4.2213
Steiner, E., Holzmann, K., Pirker, C., Elbling, L., Micksche, M., Sutterlüty, H., et al. (2006). The Major Vault Protein Is Responsive to and Interferes with Interferon-γ-Mediated STAT1 Signals. J. Cel Sci 119 (Pt 3), 459–469. doi:10.1242/jcs.02773
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene Set Enrichment Analysis: a Knowledge-Based Approach for Interpreting Genome-wide Expression Profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Tarabichi, M., Saiselet, M., Trésallet, C., Hoang, C., Larsimont, D., Andry, G., et al. (2015). Revisiting the Transcriptional Analysis of Primary Tumours and Associated Nodal Metastases with Enhanced Biological and Statistical Controls: Application to Thyroid Cancer. Br. J. Cancer 112 (10), 1665–1674. doi:10.1038/bjc.2014.665
Teng, Y., Ren, Y., Hu, X., Mu, J., Samykutty, A., Zhuang, X., et al. (2017). MVP-mediated Exosomal Sorting of miR-193a Promotes colon Cancer Progression. Nat. Commun. 8, 14448. doi:10.1038/ncomms14448
Uhlen, M., Zhang, C., Lee, S., Sjöstedt, E., Fagerberg, L., Bidkhori, G., et al. (2017). A Pathology Atlas of the Human Cancer Transcriptome. Science 357 (6352), 357. doi:10.1126/science.aan2507
Uhlén, M., Fagerberg, L., Hallström, B. M., Lindskog, C., Oksvold, P., Mardinoglu, A., et al. (2015). Tissue-based Map of the Human Proteome. Science 347 (6220), 1260419. doi:10.1126/science.1260419
van Zon, A., Mossink, M. H., Scheper, R. J., Sonneveld, P., and Wiemer, E. A. C. (2003). The Vault Complex. Cell Mol. Life Sci. (Cmls) 60 (9), 1828–1837. doi:10.1007/s00018-003-3030-y
Vanzon, A., Mossink, M., Houtsmuller, A., Schoester, M., Scheffer, G., Scheper, R., et al. (2005). Vault Mobility Depends in Part on Microtubules and Vaults Can Be Recruited to the Nuclear Envelope. Exp. Cel Res. 312 (3), 245–255. doi:10.1016/j.yexcr.2005.10.016
Vollmar, F., Hacker, C., Zahedi, R.-P., Sickmann, A., Ewald, A., Scheer, U., et al. (2009). Assembly of Nuclear Pore Complexes Mediated by Major Vault Protein. J. Cel Sci 122 (Pt 6), 780–786. doi:10.1242/jcs.039529
Wang, Q. X., Chen, E. D., Cai, Y. F., Li, Q., Jin, Y. X., Jin, W. X., et al. (2016). A Panel of Four Genes Accurately Differentiates Benign from Malignant Thyroid Nodules. J. Exp. Clin. Cancer Res. 35 (1), 1–8. doi:10.1186/s13046-016-0447-3
Xiao, Y.-S., Zeng, D., Liang, Y.-K., Wu, Y., Li, M.-F., Qi, Y.-Z., et al. (2019). Major Vault Protein Is a Direct Target of Notch1 Signaling and Contributes to Chemoresistance in Triple-Negative Breast Cancer Cells. Cancer Lett. 440-441, 156–167. doi:10.1016/j.canlet.2018.09.031
Xing, M. (2005). BRAF Mutation in Thyroid Cancer. Endocr. Relat. Cancer 12 (2), 245–262. doi:10.1677/erc.1.0978
Xing, M., Liu, R., Liu, X., Murugan, A. K., Zhu, G., Zeiger, M. A., et al. (2014). BRAF V600E and TERT Promoter Mutations Cooperatively Identify the Most Aggressive Papillary Thyroid Cancer with Highest Recurrence. Jco 32 (25), 2718–2726. doi:10.1200/jco.2014.55.5094
Xing, M. (2013). Molecular Pathogenesis and Mechanisms of Thyroid Cancer. Nat. Rev. Cancer 13 (3), 184–199. doi:10.1038/nrc3431
Xing, M., Westra, W. H., Tufano, R. P., Cohen, Y., Rosenbaum, E., Rhoden, K. J., et al. (2005). BRAF Mutation Predicts a Poorer Clinical Prognosis for Papillary Thyroid Cancer. J. Clin. Endocrinol. Metab. 90 (12), 6373–6379. doi:10.1210/jc.2005-0987
Xu, L., Deng, C., Pang, B., Zhang, X., Liu, W., Liao, G., et al. (2018). TIP: A Web Server for Resolving Tumor Immunophenotype Profiling. Cancer Res. 78 (23), 6575–6580. doi:10.1158/0008-5472.Can-18-0689
Yoo, S.-K., Lee, S., Kim, S.-j., Jee, H.-G., Kim, B.-A., Cho, H., et al. (2016). Comprehensive Analysis of the Transcriptional and Mutational Landscape of Follicular and Papillary Thyroid CancersCentral PMCID: PMCPMC4975456 along with Any Other Relevant Declarations Relating to Employment, Consultancy, Patents, Products in Development or Marketed Products Etc, Do Not Alter the Authors' Adherence to PLOS GENETICS Policies on Sharing Data and Materials. Plos Genet. 12 (8), e1006239. doi:10.1371/journal.pgen.1006239
Yoshihara, K., Shahmoradgoli, M., Martínez, E., Vegesna, R., Kim, H., Torres-Garcia, W., et al. (2013). Inferring Tumour Purity and Stromal and Immune Cell Admixture from Expression Data. Nat. Commun. 4 (1), 1–11. doi:10.1038/ncomms3612
Yu, H., Li, M., He, R., Fang, P., Wang, Q., Yi, Y., et al. (2020). Major Vault Protein Promotes Hepatocellular Carcinoma through Targeting Interferon Regulatory Factor 2 and Decreasing P53 Activity. Hepatology 72 (2), 518–534. doi:10.1002/hep.31045
Yu, K., Ganesan, K., Tan, L. K., Laban, M., Wu, J., Zhao, X. D., et al. (2008). A Precisely Regulated Gene Expression Cassette Potently Modulates Metastasis and Survival in Multiple Solid Cancers. Plos Genet. 4 (7), e1000129. doi:10.1371/journal.pgen.1000129
Keywords: major vault protein, papillary thyroid cancer, biomarker, tumor microenvironment, BRAFV600E, RAS
Citation: Dong X, Akuetteh PDP, Song J, Ni C, Jin C, Li H, Jiang W, Si Y, Zhang X, Zhang Q and Huang G (2022) Major Vault Protein (MVP) Associated With BRAFV600E Mutation Is an Immune Microenvironment-Related Biomarker Promoting the Progression of Papillary Thyroid Cancer via MAPK/ERK and PI3K/AKT Pathways. Front. Cell Dev. Biol. 9:688370. doi: 10.3389/fcell.2021.688370
Received: 30 March 2021; Accepted: 23 December 2021;
Published: 31 March 2022.
Edited by:
Hernandes F Carvalho, State University of Campinas, BrazilReviewed by:
Smrithi Rajendiran, University of California, United StatesShiyang Liu, Huazhong University of Science and Technology, China
Copyright © 2022 Dong, Akuetteh, Song, Ni, Jin, Li, Jiang, Si, Zhang, Zhang and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiyu Zhang, qiyu_zhang@wmu.edu.cn; Guanli Huang, huangguanli@wzhospital.cn
†These authors have contributed equally to this work and share first authorship