Prognostic Models Based on Machine Learning for Clinical Cancer ResearchView this Special Issue
Research Article | Open Access
Ming Bai, Qi Pan, Chen Sun, "Tumor Purity Coexpressed Genes Related to Immune Microenvironment and Clinical Outcomes of Lung Adenocarcinoma", Journal of Oncology, vol. 2021, Article ID 9548648, 13 pages, 2021. https://doi.org/10.1155/2021/9548648
Tumor Purity Coexpressed Genes Related to Immune Microenvironment and Clinical Outcomes of Lung Adenocarcinoma
Purpose. Lung cancer tissue includes tumor tissue, stromal cells, immune cells, and epithelial cells. These nontumor cells dilute the tumor purity in lung cancer tissues. Tumor purity plays an essential role in the immune response to lung cancer. At present, the biological processes related to the purity of lung cancer tumors remains unclear. Methods. We measured tumor purity in 486 lung carcinoma tissues from TCGA-LUAD FPKM by using the “estimate” R package. Lung carcinoma tumor mutation burden was calculated by analyzing TCGA single nucleotide polymorphism data. The immune cell proportion was also evacuated via the CIBERSORT method. Lung carcinoma samples with were considered significant. Based on the tumor purity and lung carcinoma gene matrix, we performed weighted gene coexpression network analysis (WGCNA), and the tumor purity-related module was identified. Then, we analyzed the functions of the factors involved in the module. We screened the coexpressed factors related to clinical outcome and immunophenotype. Finally, expression levels of these factors were measured at tissue and single-cell levels. Results. A lung cancer tumor purity correlated coexpression network was determined. Five coexpressed genes (CD4, CD53, EVI2B, PLEK, and SASH3) were identified as tumor purity coexpressed genes that negatively correlated with tumor purity. Because the factors in the coexpression network often participate in similar biological processes, we found that CD4, CD53, EVI2B, PLEK, and SASH3 were most related to positive regulation of cytokine production and interleukin−2 production through functional enrichment. In a clinical phenotype analysis, we found that these five factors can be used as independent prognostic risk factors. We found that these factors were significantly negatively correlated with tumor purity and positively correlated with the immune score in the immunophenotyping analysis. Using GSEA analysis, we found that the antigen processing and presentation pathway were related to the five tumor coexpressed genes mentioned above. SASH3 and CD53 were used to conduct a prognostic model based on the interaction analysis of the Support Vector Machine and the Least Absolute Shrinkage and Selection Operator. SASH3 was verified to be related to CD8A using a single-cell analysis. Conclusion. Tumor purity-related coexpression factors in the tumor microenvironment have essential clinical, genomic, and biological significance in lung cancer. These coexpression factors (SASH3 and CD53) can be used to classify tumor purity phenotypes and to predict clinical outcomes.
The structure of tumor tissue is complex. In addition to tumor cells, there are also other components such as stromal cells, inflammatory cells, vasculature, and the extracellular matrix .
Tumor microenvironment possesses complexity because of a mixture of growth-promoting and inhibiting growth factors, nutrients, chemokines, and other noncancer types, which interact with each other and associate with tumor growth, disease progression, drug resistance, and especially, infiltrating T lymphocytes and tumor growth . Approximately 1.8 to 8 million people are diagnosed with lung cancer each year, and 1.6 to 6 million die from lung cancer. The 5-year survival rate of lung cancer is about 4–17% . Although studies have revealed the mechanisms involved in cancer malignant characteristics and identified reasonable therapeutic targets , current clinical prediction and treatment outcome of lung cancer are not satisfactory . We are now aware that lung cancer tissues are rich in nontumor cells, of which stromal cells significantly regulate tumor proliferation, invasion, and angiogenesis .
With the development of bioinformatics in recent years and the acquisition of open lung cancer cohorts, it becomes possible to evaluate tumor purity content in the samples according to the estimated infiltration of stromal and immune cells ESTIMATE algorithm . In this study, we constructed a tumor purity coexpression network based on weighted gene coexpression network analysis (WGCNA) . We explored the coexpression factors most related to tumor purity and related biological functions and demonstrated the most relevant biological functions and mechanism of action affecting tumor purity in the lung cancer tumor microenvironment.
2.1. Data Sources
The Cancer Genome Atlas- (TCGA-) LUAD FPKM data containing 486 cancer tissue samples were obtained (http://cancergenome.nih.gov/) . GSE99254  is a single-cell sequencing cohort with 14 samples and 12346 non-small-cell lung cancer cells and was obtained from the GEO (http://www.ncbi.nlm.nih.gov/geo/) database with the GPL16791 platform and GPL20301 platform. Meanwhile, GSE42127  was also downloaded to verify the conclusion.
2.2. Tumor Purity Evaluation
Expression data (ESTIMATE)  were used to evaluate stromal and immune cells in malignant tumor tissues in this study and estimate the proportion of stromal and immune cells in the tumor microenvironment based on the gene matrix. Through ESTIMATE algorithm, tumor purity of each lung cancer sample in TCGA-LUAD was obtained. CIBERSORT algorithm [12, 13] is a method in order to evaluate the cell content in bulk tissue gene expression matrices. Immune cell infiltration levels were calculated based on the LM22 matrix and CIBERSORT algorithm, and samples with were considered significant and taken into this study.
2.3. Tumor Purity-Related Coexpression Factors
Weighted Gene Coexpression Network Analysis (WGCNA) was demonstrated to determine tumor purity coexpressed genes in lung cancer. This method converted tumor purity coexpression correlations into weight values which determined the coexpression factors. As we know, the expression levels of genes were approximately the same as those possessing similar biological functions . In this research, we set the soft threshold as 5, R square = 0.98, and the factors in the minimum module as 30. We uploaded tumor purity scores and immune cell proportions as phenotype files. In this manner, a cluster of tumor purity coexpression genes with similar biological function was determined via WGCNA .
2.4. Protein-Protein Network and Function Analysis
The encoding genes of tumor purity coexpression proteins were identified by the Pearson correlation coefficient >0.4. The coexpression modules of tumor purity were conducted by Cytoscape software. Meanwhile, the tumor purity coexpressed genes were enriched to explore their biological processes in the tumor microenvironment. The Database for Annotation, Visualization, and Integrated Discovery (DAVID, v6.8) is an online database which provides functional annotation analysis [16, 17]. The Kyoto Encyclopedia of Genes and Genomes (KEGG)  (https://www.genome.jp/kegg/) and Gene Ontology (GO)  (http://geneontology.org/) analysis were used to identify the biological function and related regulation pathways in each coexpression module.
2.5. Prognosis Model Based on LASSO and Support Vector Machine Methods
Univariate Cox regression analysis was performed for tumor purity coexpression genes, and the genes with were taken into the feature selection. Subsequently, we screened the characteristic variables and constructed the prognostic model by observing LASSO regression analysis [20–22] and support vector machine (SVM-RFE algorithms).
2.6. Gene Set Enrichment Analysis (GSEA)
GSEA calculates the significance and consistency differences of a predefined dataset between two biological states . The gene matrix in TCGA was divided into high- and low-expression groups, following the median expression level of lung cancer tumor purity-related genes. Through GSEA analysis, we obtained the related pathways of genes which were correlated with tumor purity and prognosis. These pathways are considered related to the immune microenvironment as well as clinical phenotypes.
2.7. Immune Correlation in Other Types of Cancer
The TIMER database (https://cistrome.shinyapps.io/timer/) [24, 25] was applied to show the correlations between SASH3 and immune cell proportion in 33 types of cancers. A Pearson correlation coefficient higher than 0.4 was considered significant.
2.8. Single-Cell Cohort Analysis
We found that the factor with the strongest negative correlation with tumor purity was CD8+ T lymphocytes. Therefore, we aimed to verify this relationship at the single-cell sequencing level. We obtained the GSE99254 single-cell cohort from the GEO database. The Seurat package was then used to filter and standardize the data . Various cell subpopulations were obtained by the TSNE dimensionality reduction clustering method . Finally, the SingleR package was used to annotate the cell types of these subpopulations . We demonstrated the relationship between SASH3 and CD8+ T lymphocyte infiltration by labeling the relationship between SASH3 and CD8+ T lymphocyte infiltration.
2.9. Statistical Analysis
R 3.6.3 (https://www.r-project.org/) was carried out for statistical analysis. Student’s t-tests are applied to show purity differences in various subgroups in the TCGA cohort. Coexpression coefficients of tumor purity protein encoding genes were evaluated based on the Pearson correlation.
3.1. Tumor Purity Coexpression Network
We obtained tumor purity and immune cell content of each person in TCGA-LUAD which is uploaded into Supplementary Table 1. The corresponding immune cell content in each sample is shown in Figure 1(a). The results demonstrated that CD8+ T lymphocytes and CD4+ cells had the highest content. Next, a dimension-reducing cluster was conducted (Figure 1(b)) for the samples of TCGA-LUAD using omics clustering. We obtained 23 coexpression networks (Figure 1(c)) through WGCNA, where each color represented one coexpression network. Furthermore, we examined the correlation between coexpression networks and tumor purity to identify the most relevant ones. The results elucidated that the yellow and green modules had the strongest correlation with tumor purity (Figure 1(d)). WGCNA results have been uploaded in Supplementary Table 2.
3.2. Protein-Protein Network and Function Enrichment
We plotted a scatter plot of the correlation between tumor purity and coexpression modules in the yellow and green modules (Figure 2(a)). The results showed that the correlation between tumor purity and gene coexpression module in the yellow module was the most significant (COR = 0.96; ), whereas the correlation in the green module was lower (COR = 0.73; ). Furthermore, GO enrichment analysis of genes in the yellow module suggested that positive regulation of chemokines and the generation of interleukin-2 were the most significant enriched pathways (Figure 2(b)), while genes in the green coexpression module are associated with the extracellular matrix. Thereby, the protein encoding genes in the yellow coexpression module were selected for subsequent analyses.
3.3. Clinical Phenotype and Immune Phenotype
The survival analysis for factors in the yellow module that can be used as independent prognostic evaluation for overall survival is shown in Figure 3. We then performed clinical and immunophenotypic assessments of these factors (Figure 4). The results suggested that CD4, CD53, EVI2b, PLEK, and SASH3 correlated with tumor purity, immune score, CD8+ T cells, and clinical phenotypes (Figure 4(a)). Low expression level of these genes led to high tumor purity, low immune score, low CD8+ T lymphocyte content, and shorter 5-year survival. The scatter plots of correlations between CD4, CD53, EVI2b, PLEK, SASH3, and tumor environment score are shown in Figures 4(b)–4(e). Results shown in Figure 5(a) indicated negative correlations between the clinical stages and the expression of these five genes (Figure 5(a)).
3.4. GSEA Analysis
GSEA analysis elucidated that chemokine-chemokine receptor interaction and the T-cell receptor signaling pathway were enriched in the high expression group of factors in the prognosis model (Figure 5(b)). These pathways enhanced the immune response and showed antitumor immune response. These factors might reduce tumor purity of lung cancer by elevating lymphocyte proportion.
3.5. LASSO Regression and SVM
We incorporated the protein encoding genes in the yellow coexpression module into the LASSO regression model and identified five significant prognostic survival genes (SASH3, PLEK, EVI2B, CD53, and CD4). Simultaneously, the support vector machine method was used to screen the features of the abovementioned factors, and four feature variables (SASH3, MNDH, CD53, and CD16) were determined. We finally identified SASH3 and CD53 as tumor purity-related prognostic factors (Figure 5(c)). Risk score = −0.004CD53 – 0.014SASH3. We later found a significant survival difference in lung cancer patients in the TCGA-LUAD cohort between the two risk scores (HR = 1.9; ).
3.6. SASH3 Related to CD8+ T Cell and Immunohistochemistry
In the abovementioned study, we found that SASH3 was significantly negatively correlated with tumor purity but positively correlated with CD8+ T lymphocytes. To further verify this positive correlation, we verified this conclusion in 33 TCGA-type cancers. The results showed that SASH3 was positively correlated with the content of CD8+ T lymphocytes in lung cancer, glioma, liver cancer, and other cancers (Figure 6(a)). At the same time, we found that the staining strength of the SASH3 antibody in lung cancer tissues of China Medical University was higher in paracancerous tissues but relatively lower in tumor tissues (Figure 6(b)). Finally, we found that the distribution of SASH3 in the single-cell cohort was like that of the CD8+ T lymphocyte biomarker CD8A (Figure 6(c)). Finally, we added external queue validation to prove the correlation between SASH3 and CD8A in GSE42127 (Supplementary Figure 1).
In this study, we first calculated the tumor purity of lung cancer tissue. Then, we established a coexpression network related to tumor purity of lung cancer, thereby obtaining the two modules with the highest correlation to tumor purity. A PPI network was established for critical genes in the module, five coexpressed genes were identified, and the enriched pathways were calculated. Then, two machine learning methods (LASSO regression and SVM) were used to establish the model. The intersection was taken to screen out that SASH3 and CD53 were tumor purity-related prognostic genes of lung cancer. Clinical phenotype and immune phenotype assessments of the coexpressed genes showed that SASH3 negatively correlated with tumor purity and positively correlated with CD8+ T lymphocytes. This result was verified by single-cell cohort sequencing, pan-cancer analysis, and immunohistochemistry.
T cells were dominant in lung cancer. CD4+ T cells (26%) were the most abundant T cell population, followed by CD8+ T cells (22%) . CD4 encodes the CD4 membrane glycoprotein of T lymphocytes. The CD4 antigen and the T-cell receptor on the T lymphocyte work together to complete the antigen presentation and recognition . Many scientists are interested in the role of CD4 immunity in the efficacy of PD-L1/PD-1 blocking therapy. Kagamu et al. found that immune monitoring of CD4+ T cells in peripheral blood predicted anti-PD-1 treatment responses in lung cancer patients . Preclinical studies in patients and mouse models have demonstrated the importance of CD4 immunity for immunotherapy . Patients who responded to treatment showed a high proportion of CD4+ T cells before treatment. These CD4+ T cells demonstrated proliferation at baseline and responded to PD-1 blockade . These findings support the idea of using vaccination to enhance CD4+ neoantigen-specific T cells in antitumor immunity .
Antitumor immunity is determined by the presence of different immune cells in the tumor microenvironment (TME). Environmental signals transmitted through the plasma membrane determine whether immune cells are activated or suppressed. Tetrantin proteins are a significant component of the plasma membrane because they aggregate immune receptors, enzymes, and signaling molecules into the tetrantin reticulum . CD53 is a four-transmembrane protein, mainly expressed in the myeloid lymphoid system . Yunta and Lazo found that CD53 antigen stimulation may have a protective effect on programmed cell death. CD53 antigen interaction protects against the apoptotic response caused by serum deprivation and contributes to cell survival in the poorly vascularized region of the tumor mass . CD53 is also essential for B-cell function because CD53 promotes BCR-dependent protein kinase C signaling, allowing it to phosphorylate its substrate .
SAM and SH3 domain containing 3 (SASH3) encoded proteins act as signal transduction proteins in lymphocytes . Pleckstrin (PLEK) is a protein found in platelets and white blood cells that acts as a substrate for protein kinase C . The ecotropic viral integration site 2B (EVI2B) gene was in the intron of the neurofibromatosis type 1 (NF1) gene and transcribed in the opposite direction to the NF1 gene . Like the NF1 gene, EVI2B is involved in the differentiation of melanocytes and keratinocytes . There is currently a lack of research on the relationship between EVI2B, SASH3, and PLEK with lung cancer. However, Huang et al. analyzed the genes of colorectal cancer patients with membrane array and direct sequencing and found that EVI2B may be a potential prognostic marker in CRC patients . Other scholars found that there are mutations in EV2B in the mutation spectrum of breast cancer cell line NZBR under conditions of physiological oxygen concentrations .
Although this study integrated relatively multiple bioinformatics analysis and improved immunohistochemical experimental verification, there are some limitations. More large samples are needed to validate our results. In vitro and in vivo experiments should be conducted, and a feasibility study for clinical practice should be contemplated.
In conclusion, two coexpression factors (SASH3 and CD53) help classify tumor purity phenotypes and predict clinical phenotype in lung cancer with the chemokine signaling pathway. The mechanism might provide concepts to modify the curative effect in patients with high tumor purity.
The datasets TCGA-BRCA for this study can be found in the The Cancer Genome Atlas (http://cancergenome.nih.gov/). The datasets GSE99254 in this study can be found in the GEO (http://www.ncbi.nlm.nih.gov/geo/).
Conflicts of Interest
The authors declare that the research was conducted in the absence of any conflicts of interest.
Ming Bai and Chen Sun conceived and designed the experiments, Qi Pan reviewed the drafts of this paper, and Ming Bai, Chen Sun, and Qi Pan approved the final draft.
The authors want to thank the TCGA and GEO for free use. This work was supported by grants from the Doctor Research Initiation Program of Liaoning Province (2020-BS-094).
Supplementary Table 1: the results of tumor purity. Supplementary Table 2: the results of WGCNA analysis. Supplementary Figure 1: the verification of the correlation between CD8A and SASH3 in the GEO database. (Supplementary Materials)
- L. Hui and Y. Chen, “Tumor microenvironment: sanctuary of the devil,” Cancer Letters, vol. 368, no. 1, pp. 7–13, 2015.
- E. Hirata and E. Sahai, “Tumor microenvironment and differential responses to therapy,” Cold Spring Harbor Perspectives in Medicine, vol. 7, no. 7, Article ID 26781, 2017.
- F. R. Hirsch, G. V. Scagliotti, J. L. Mulshine et al., “Lung cancer: current therapies and new targeted treatments,” The Lancet, vol. 389, no. 10066, pp. 299–311, 2017.
- Y. Mao, D. Yang, J. He, and M. J. Krasna, “Epidemiology of lung cancer,” Surgical Oncology Clinics of North America, vol. 25, no. 3, pp. 439–445, 2016.
- F. Nasim, B. F. Sabath, and G. A. Eapen, “Lung cancer,” Medical Clinics of North America, vol. 103, no. 3, pp. 463–473, 2019.
- V. Mittal, T. El Rayes, N. Narula, T. E. McGraw, N. K. Altorki, and M. H. Barcellos-Hoff, “The microenvironment of lung cancer and therapeutic implications,” Lung Cancer and Personalized Medicine: Novel Therapies and Clinical Management, vol. 890, pp. 75–110, 2016.
- K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.
- P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.
- Z. Wang, M. A. Jensen, and J. C. Zenklusen, “A practical guide to the cancer Genome Atlas (TCGA),” Methods in Molecular Biology, vol. 1418, pp. 111–141, 2016.
- X. Guo, Y. Zhang, L. Zheng et al., “Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing,” Nature Medicine, vol. 24, no. 7, pp. 978–985, 2018.
- S. K. Hight, A. Mootz, R. K. Kollipara et al., “An in vivo functional genomics screen of nuclear receptors and their co-regulators identifies FOXA1 as an essential gene in lung tumorigenesis,” Neoplasia, vol. 22, no. 8, pp. 294–310, 2020.
- A. M. Newman, C. B. Steen, C. L. Liu et al., “Determining cell type abundance and expression from bulk tissues with digital cytometry,” Nature Biotechnology, vol. 37, no. 7, pp. 773–782, 2019.
- C. B. Steen, C. L. Liu, A. A. Alizadeh, and A. M. Newman, “Profiling cell type Abundance and expression in bulk tissues with CIBERSORTx,” Methods in Molecular Biology, vol. 2117, pp. 135–157, 2020.
- E. Radulescu, A. E. Jaffe, R. E. Straub et al., “Identification and prioritization of gene sets associated with schizophrenia risk by co-expression network analysis in human brain,” Molecular Psychiatry, vol. 25, no. 4, pp. 791–804, 2020.
- J. Long, S. Huang, Y. Bai et al., “Transcriptional landscape of cholangiocarcinoma revealed by weighted gene coexpression network analysis,” Briefings in Bioinformatics, 2020.
- D. Huang, B. T. Sherman, Q. Tan et al., “The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists,” Genome Biology, vol. 8, no. 9, p. R183, 2007.
- D. W. Huang, B. T. Sherman, Q. Tan et al., “DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists,” Nucleic Acids Research, vol. 35, no. suppl_2, pp. W169–W175, 2007, Web Server issue.
- M. Kanehisa, M. Furumichi, M. Tanabe, Y. Sato, and K. Morishima, “KEGG: new perspectives on genomes, pathways, diseases and drugs,” Nucleic Acids Research, vol. 45, no. D1, pp. D353–D361, 2017.
- M. Ashburner, C. A. Ball, J. A. Blake et al., “Gene Ontology: tool for the unification of biology,” Nature Genetics, vol. 25, no. 1, pp. 25–29, 2000.
- N. Simon and R. Tibshirani, “Standardization and the group LASSO penalty,” Statistica Sinica, vol. 22, no. 3, pp. 983–1001, 2012.
- J. Bien, J. Taylor, and R. Tibshirani, “A LASSO for hierarchical interactions,” Annals of Statistics, vol. 41, no. 3, pp. 1111–1141, 2013.
- S. Bates and R. Tibshirani, “Log‐ratio lasso: scalable, sparse estimation for log‐ratio models,” Biometrics, vol. 75, no. 2, pp. 613–624, 2019.
- S. Hänzelmann, R. Castelo, and J. Guinney, “GSVA: gene set variation analysis for microarray and RNA-seq data,” BMC Bioinformatics, vol. 14, no. 1, p. 7, 2013.
- T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
- T. Li, J. Fu, Z. Zeng et al., “TIMER2.0 for analysis of tumor-infiltrating immune cells,” Nucleic Acids Research, vol. 48, no. W1, pp. W509–W514, 2020.
- A. Gribov, M. Sill, S. Lück et al., “SEURAT: visual analytics for the integrated analysis of microarray data,” BMC Medical Genomics, vol. 3, no. 1, p. 21, 2010.
- G. Dimitriadis, J. P. Neto, and A. R. Kampff, “t-SNE visualization of large-scale neural recordings,” Neural Computation, vol. 30, no. 7, pp. 1750–1774, 2018.
- D. Aran, A. P. Looney, L. Liu et al., “Reference-based analysis of lung single-cell sequencing reveals a transitional profibrotic macrophage,” Nature Immunology, vol. 20, no. 2, pp. 163–172, 2019.
- B. Stankovic, H. A. K. Bjørhovde, R. Skarshaug et al., “Immune cell composition in human non-small cell lung cancer,” Frontiers in Immunology, vol. 9, p. 3101, 2018.
- D. Mucida, M. M. Husain, S. Muroi et al., “Transcriptional reprogramming of mature CD4+ helper T cells generates distinct MHC class II-restricted cytotoxic T lymphocytes,” Nature Immunology, vol. 14, no. 3, pp. 281–289, 2013.
- H. Kagamu, S. Kitano, O. Yamaguchi et al., “CD4+ T-cell immunity in the peripheral blood correlates with response to anti-PD-1 therapy,” Cancer Immunology Research, vol. 8, no. 3, pp. 334–344, 2020.
- M. H. Spitzer, Y. Carmi, N. E. Reticker-Flynn et al., “Systemic immunity is required for effective cancer immunotherapy,” Cell, vol. 168, no. 3, pp. 487–502, 2017, e15.
- M. Zuazo, H. Arasanz, G. Fernández-Hinojal et al., “Functional systemic CD4 immunity is required for clinical responses to PD-L1/PD-1 blockade therapy,” EMBO Molecular Medicine, vol. 11, no. 7, Article ID e10293, 2019.
- J. R. Veatch, B. L. Jesernig, J. Kargl et al., “Endogenous CD4+ T cells recognize neoantigens in lung cancer patients, including recurrent oncogenic KRAS and ERBB2 (Her2) driver mutations,” Cancer Immunology Research, vol. 7, no. 6, pp. 910–922, 2019.
- F. Schaper and A. B. van Spriel, “Antitumor immunity is controlled by tetraspanin proteins,” Frontiers in Immunology, vol. 9, p. 1185, 2018.
- J. Dagher, B. Delahunt, N. Rioux-Leclercq et al., “Clear cell renal cell carcinoma: validation of world health organization/international society of urological pathology grading,” Histopathology, vol. 71, no. 6, pp. 918–925, 2017.
- M. Yunta and P. A. Lazo, “Apoptosis protection and survival signal by the CD53 tetraspanin antigen,” Oncogene, vol. 22, no. 8, pp. 1219–1224, 2003.
- M. Zuidscherwoude, V. E. Dunlock, G. van den Bogaart et al., “Tetraspanin microdomains control localized protein kinase C signaling in B cells,” Science Signaling, vol. 10, no. 478, 2017.
- E. Astoul, A. D. Laurence, N. Totty, S. Beer, D. R. Alexander, and D. A. Cantrell, “Approaches to define antigen receptor-induced serine kinase signal transduction pathways,” Journal of Biological Chemistry, vol. 278, no. 11, pp. 9267–9275, 2003.
- L. Lian, Y. Wang, M. Flick et al., “Loss of pleckstrin defines a novel pathway for PKC-mediated exocytosis,” Blood, vol. 113, no. 15, pp. 3577–3584, 2009.
- R. M. Cawthon, L. B. Andersen, A. M. Buchberg et al., “cDNA sequence and genomic structure of EVI2B, a gene lying within an intron of the neurofibromatosis type 1 gene,” Genomics, vol. 9, no. 3, pp. 446–460, 1991.
- D. Kaufmann, S. Gruener, F. Braun et al., “EVI2B, a gene lying in an intron of the neurofibromatosis type 1 (NF1) gene, is as the NF1 gene involved in differentiation of melanocytes and keratinocytes and is overexpressed in cells derived from NF1 neurofibromas,” DNA and Cell Biology, vol. 18, no. 5, pp. 345–356, 1999.
- M.-Y. Huang, H.-M. Wang, T.-S. Tok et al., “EVI2B, ATP2A2, S100B, TM4SF3, and OLFM4 as potential prognostic markers for postoperative Taiwanese colorectal cancer patients,” DNA and Cell Biology, vol. 31, no. 4, pp. 625–635, 2012.
- E. Y. Leung, M. E. Askarian-Amiri, D. C. Singleton et al., “Derivation of breast cancer cell lines under physiological (5%) oxygen concentrations,” Frontiers in Oncology, vol. 8, p. 425, 2018.
Copyright © 2021 Ming Bai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.