Skip to main content

ORIGINAL RESEARCH article

Front. Cell Dev. Biol., 13 May 2022
Sec. Cancer Cell Biology
This article is part of the Research Topic Infectious Disease Agents and Cancer View all 6 articles

Pan-Cancer Analysis, Reveals COVID-19-Related BSG as a Novel Marker for Treatment and Identification of Multiple Human Cancers

  • 1Department of Cardiothoracic Vascular Surgery, The Affiliated Hospital of Youjiang Medical University for Nationalities, Baise, China
  • 2The First Clinical Medical College, The First Affiliated Hospital of Guangxi Medical University, Nanning, China

Background: Coronavirus disease 2019 (COVID-19) has been a public threat and healthcare concern caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. During the period of the pandemic of COVID-19, cancer patients should be paid more attention as more severe events are found in cancer patients infected with SARS-CoV-2. Basigin (BSG) is an essential factor for the infection and progression of COVID-19 and tumorigenesis of multiple tumors, which may serve as a novel target for the effective treatment against COVID-19 and multiple human cancers.

Methods: A total of 19,020 samples from multiple centers were included in our research for the comprehensive investigation of the differences in BSG expression among human organs, cancer cells, cancer tissues, and normal tissues. Cox regression analysis and Kaplan–Meier curves were utilized to explore the prognosis factor of BSG in cancers. Correlation analyses were used to determine associations of BSG expression with tumor mutational burden, the immune microenvironment, etc. Gene set enrichment analysis was applied to explore the underlying mechanisms of BSG in cancers.

Results: Compared with normal tissues, BSG expression was high in 13 types of cancers (cholangiocarcinoma, etc.) and low in colon adenocarcinoma and rectum adenocarcinoma. BSG expression was related to the prognosis of eight cancers (e.g., invasive breast carcinoma) (p < 0.05). The gene also demonstrated a pronounced effect in identifying 12 cancers (cholangiocarcinoma, etc.) from their control samples (AUC >0.7). The BSG expression was associated with DNA methyltransferases, mismatch repair genes, immune infiltration levels, tumor mutational burden, microsatellite instability, neoantigen, and immune checkpoints, suggesting the potential of BSG as an exciting target for cancer treatment. BSG may play its role in several cancers by affecting several signaling pathways such as drug cytochrome metabolism P450 and JAK-STAT.

Conclusion: BSG may be a novel biomarker for treating and identifying multiple human cancers.

Introduction

Coronavirus disease 2019 (COVID-19) has been a public threat and healthcare concern caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and spread by respiratory droplets through coughing or sneezing among the contacts (Li Q et al., 2020). During the period of the pandemic of COVID-19, cancer patients should be paid increasing attention as more severe events are found in cancer patients infected with SARS-CoV-2 (Madariaga et al., 2020; Moujaess, Kourie, and Ghosn 2020). Exploring strategies to control the epidemic of SARS-CoV-2 in cancer patients will benefit the continue safe treatment for these patients during the pandemic of COVID-19 (Madariaga et al., 2020).

Basigin (BSG) is a transmembrane glycoprotein of the immunoglobulin superfamily, which is also known as cluster of differentiation 147 (CD147), extracellular matrix metalloproteinase inducer (EMMPRIN), and so on. According to recent reports, BSG serves as a target that correlated with SARS-CoV-2 spike protein that can allow SARS-CoV-2 to enter host cells through endocytosis and even initiate the MAPK pathway through the spike protein/CD147/CyPA” signaling axis and induce cytokine storm (Wang K et al., 2020; Fenizia et al., 2021; Geng et al., 2021; Behl et al., 2022). It is also reported that the existence of humanized anti-CD147 antibodies can inhibit SARS-CoV-2 and its variants entering host cells (Geng et al., 2021). Therefore, BSG is a crucial element for the infection and progression of COVID-19 and may serve as a novel target for effective treatment against COVID-19 (Wang K et al., 2020; Behl et al., 2022).

In addition to COVID-19, BSG is also reported to be involved in the tumorigenesis of multiple tumors, such as melanoma (Kim et al., 2021; Lu et al., 2021), liver hepatocellular carcinoma (LIHC) (Wang SJ et al., 2020), and hypopharyngeal squamous cell carcinoma (Suzuki et al., 2021). BSG has been found to induce the degradation of the extracellular matrix through the matrix metalloproteinase proteins and stimulates the vascular endothelial growth factor, thus participating in the invasion and metastasis of tumors (Weidle et al., 2010; Kim et al., 2021). It is also reported that T cells and NK cells carrying a chimeric antigen receptor that recognizes BSG is detrimental to the tumor cell lines in vitro (Tseng et al., 2020). The inhibition of BSG can reduce the growth, migration, and invasion of tumor cells (Hatanaka et al., 2016; Lian et al., 2017; Dana et al., 2020). Therefore, BSG is a potential predictive biomarker for tumors (Hatanaka et al., 2016; Qiao et al., 2020; Tseng et al., 2020), which may serve as a biological modulator of tumor and COVID-19 and improve the management of cancer patients during the period of the COVID-19 pandemic (Varadarajan et al., 2020). However, BSG expression and its relation with clinical relevance and immune infiltration levels in pan-cancer were rarely reported, which requires further investigation.

Based on analyzing thousands of specimens from several high-throughput databases, including Genotype-Tissue Expression (GTEx), Cancer Cell Line Encyclopedia (CCLE), and The Cancer Genome Atlas (TCGA), this study comprehensively measures the BSG expression level and its potential clinical significance in pan-cancer. Moreover, the study also investigated the association of BSG with alternative genes and immune infiltration levels in cancers, contributing to the understanding of BSG as a novel and potential biomarker for the identification and treatment of multiple cancers.

Materials and Methods

Collection of BSG mRNA Expression, Clinical Characteristics, and Prognosis Data

GTEx Consortium (2013) includes plenty of normal tissue specimens of Homo sapiens, and CCLE (Ghandi et al., 2019) contains numerous samples of various cell lines of cancers. Transcriptome data (for assessing the expression of BSG in normal tissues and cancer lines) of GTEx and CCLE were collected from GTEx Portal (8,671 samples) and DepMap Portal (459 specimens), respectively (Supplementary Table S1). TCGA collects transcriptome and clinical data of 32 kinds of cancers, and 9,054 cancers and 727 controls of TCGA were downloaded from the Xena database; among the 32 cancers, twelve types (including 2,075 samples) had less than three corresponding controls and were not utilized for investigating the difference in BSG expression between normal tissues and cancer tissues. Three types of samples of TCGA were eventually included in this study, containing primary tumor tissue, normal solid tissue, and primary blood-derived cancer peripheral blood. All gene expression levels in data obtained from GTEx, CCLE, and TCGA were processed in R (v4.1.0) using log2 (x + 1).

Three kinds of clinical characteristics—American Joint Committee on Cancer (AJCC), age, and gender—were obtained from the Xena database. Four types of survival data—overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI)—of cancer patients were downloaded from the Xena database. In order to avoid the impact of short-term follow-up time and small sample size, the inclusion criteria of prognostic information for survival analysis were set as follows: 1) the prognostic information of patients with follow-up time not less than 30 days and 2) cancers with at least three cancer samples. The included clinical characteristics and prognosis data can be viewed in Supplementary Table S1.

Collection of BSG Protein Level Data

Immunohistochemical staining data for validating the BSG protein levels in cancer and normal tissues were obtained from The Human Protein Atlas (proteinatlas.org) (Uhlen et al., 2015) by using the “HPAanalyze” package (Tran et al., 2019) in R (v4.1.0). The BSG protein levels of cancers with no fewer than two samples in both normal tissues and cancer tissues were chosen in this study. For one cancer with more than one kind of antibody, data on antibodies used to detect more samples of this cancer were included. Eventually, 109 samples were selected for further analysis, and these samples’ identity documents, groups, and tissue types were stored in Supplementary Table S2. Normal tissues and tumor tissues were obtained from various individuals, and no person provided both normal and tumor samples. Additionally, three normal lung specimens were utilized twice (one for lung adenocarcinoma [LUAD] and the other for lung squamous cell carcinoma [LUSC]) for comparing the difference in the BSG protein levels between normal and cancer tissues (Supplementary Table S2). A bladder cancer (BLCA) tissue sample and a PRAD tissue sample were collected from the same patient with these two diseases (Supplementary Table S2). The BSG protein levels were evaluated with the criteria (listed in Supplementary Table S3) for staining intensity and quantity scores. The BSG protein levels were evaluated by the total immunohistochemical score, the product of the staining intensity, and quantity scores. All protein level scores for the 109 samples collected from The Human Protein Atlas are shown in Supplementary Table S2.

Alternations, DNA Methyltransferases, and the Immune Microenvironment of BSG in Cancers

Alternations of BSG in multiple cancers were explored based on the data (processed by MuTect2 software) of the GDC Portal, and this step was performed in Sanger Box (v3.0). DNA methyltransferases (DNMTs) can influence gene expression and genetic performance without changing the DNA sequence. The expression relationships of BSG with three typical DNMTs—DNMT1, DNMT3A, and DNMT3B—were investigated in this study.

Based on the deconvolution method, the TIMER (Li et al., 2017; Li T et al., 2020) algorithm makes it feasible to calculate the immune abundance of six kinds of immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophil, macrophage, and dendritic cells) by whole-gene expression data. The other algorithm, ESTIMATE (Yoshihara et al., 2013), enables researchers to evaluate the immune abundance of immune cells and immune stromal based on gene expression data matrix. With the ESTIMATE algorithm, the abundance of immune stromal, immune cells, and tumor purity are reflected by the stromal (for stromal cells), immune (for immune cells), and ESTIMATE scores (for tumor purity). The TIMER and ESTIMATE data (Supplementary Table S1) for investigating the correlation of BSG expression with the immune microenvironment were collected from TIMER’s official website and SangerBox 3.0, respectively.

The Expression of BSG With Mismatch Repair Genes, Tumor Mutational Burden, Microsatellite Instability, Neoantigen Count, and Immune Checkpoint Genes

Mismatch repair genes (MMRs) participate in repairing errors of DNA replications, and the lack of MMRs can lead to more somatic mutations. The expression levels of five MMRs (MLH1, MSH2, MSH6, PMS2, and EPCAM) and forty-six immune checkpoints (BTLA, etc.) were extracted from TCGA dataset. For cancer patients included in TCGA dataset, data on tumor mutational burden (TMB), microsatellite instability (MSI), and neoantigen count (including single-nucleotide variation and indel neo-antigens) have been published by the previous study (Liu et al., 2018) and were applied in this research.

The Potential Mechanisms of BSG in Pan-Cancer and the Underlying Target Drugs for the Gene

Using the “clusterProfiler” package (Yu et al., 2012), the study investigated potential mechanisms of BSG in 32 cancers based on KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways by gene set enrichment analysis (GSEA). For essential signaling pathways showing statistical significance in multiple cancers, data of their leading-edge genes were extracted to explore these genes’ correlations with BSG expression, which verify the relation of BSG with the essential signaling pathways. Using CellMiner (Reinhold et al., 2015) and based on half-maximal inhibitory concentration (IC50) values, this research also detected the sensitivity of BSG to 57 drugs approved by the Food and Drug Administration of America or validated by clinical tests.

Statistical Analysis

The Kruskal–Wallis test was selected to assess the difference in BSG expression in various normal tissues and cancer cell lines. The Wilcoxon rank-sum test was used to assess the differential expression of BSG (at both mRNA and protein levels) between cancer tissues and their controls. Multiple comparisons by the false discovery rate method were applied to verify the distinct expression levels of BSG mRNA between cancers and controls. The homogeneity between cancer and control groups was evaluated using the chi-square test or Fisher’s exact test.

The Wilcoxon rank-sum test was utilized to explore the relationship between BSG expression and the three clinical features (AJCC stage, age, and gender). Additionally, if the BSG expression was associated with not less than two clinical parameters (e.g., age and gender), the adjustment test would be performed and the correlation between the BSG expression and age would be analyzed hierarchically based on various genders.

Using “survival” and “forestplot” packages in R (v4.1.0), univariate Cox regression analysis and Kaplan–Meier plots were applied to measure the relevance between BSG expression and OS, DSS, DFI, and PFI. The optimal cut point for high- and low-BSG expression levels in each Kaplan–Meier curve was determined by utilizing the maximally selected rank statistics from the “maxstat” and “survminer” R packages. For detecting the accuracy of BSG expression in identifying cancers from controls, the specificity, sensitivity, and the area under the curve (AUC) of the receiver operating characteristics (ROC) were calculated using the “pROC” package (Robin et al., 2011) for R (v4.1.0), and the specificity, sensitivity, and the AUC of a summary ROC were determined in Stata (v15.0).

The Spearman rank correlation analyses were applied to evaluate the correlations of BSG expression with DNMTs, immune environment, TMB, MSI, neoantigen count, immune checkpoints, and leading-edge genes. The Wilcoxon rank-sum test was utilized to determine the difference in abundance of six kinds of immune cells and the difference in IC50 of certain drugs between high- and low-BSG expression groups.

In this study, p-value < 0.05 demonstrated statistical significance, and except for drug-sensitive analysis, all p-values were two-sided. An overflow of the study can be viewed in Figure 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. Overflow of this study. BSG, basigin; “Cox” represents Cox regression analysis; “KM” represents Kaplan–Meier curves; MMRs, mismatch repair genes; TMB, tumor mutational burden; MSI, microsatellite instability.

Results

Differently Expressed BSG in Human Cancers

Based on GTEx data, the BSG expression was various in different organs its high expression was observed in normal tissues of the adrenal gland, heart, uterus, and so on, while significantly low expression of BSG was found in the liver and pancreas (p < 0.05; Figure 2A). An analysis based on CCLE data demonstrated the different distribution of BSG in various cancer cell lines (p < 0.05; Figure 2B). In the observed 20 cancer types, the BSG expression levels of 16 cancers differed from those of the corresponding control samples. In detail, upregulated BSG expression was observed in 14 cancers: BLCA, invasive breast carcinoma (BRCA), cholangiocarcinoma (CHOL), esophageal carcinoma (ESCA), head and neck squamous cell carcinoma (HNSCC), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), LIHC, LUAD, LUSC, prostate adenocarcinoma (PRAD), stomach adenocarcinoma (STAD), thyroid carcinoma (THCA), and uterine corpus endometrial carcinoma (UCEC); downregulated BSG expression was detected in two cancers: colon adenocarcinoma (COAD) and rectum adenocarcinoma (READ) (p < 0.05; Figure 2C). Notably, in COAD, LIHC, and UCEC, the age distribution between cancer and control groups was statistically different. Similarly, in BLCA, there was a difference in gender distribution between the cancer group and the control group (Supplementary Table S4). In order to exclude the potential effects of gender and age covariates in the difference in BSG expression between cancer and control groups, a hierarchical evaluation was performed for these four cancers. Increasing levels of BSG were observed in female BLCA patients but not in male BLCA patients (Figure 2E). Thus, gender distribution may contribute to the difference in BSG expression in BLCA. Increasing levels of BSG were observed in female BLCA patients but not in male BLCA patients (Figure 2E).

FIGURE 2
www.frontiersin.org

FIGURE 2. BSG expression in normal tissues and pan-cancer at mRNA and protein levels. (A) Distinct BSG expression in various normal tissues of humans; the dashed line represents the average value (equaling to 8.549) of BSG expression in all normal tissues shown in this panel. (B) Different BSG expressions in various cancer cell lines of humans. (C) Different BSG expressions between cancer tissues with their normal tissues. (D–E) Hierarchical analysis of BSG expression differences in some cancers and their control tissues. (F) Elevated BSG levels in several cancers were validated at the protein level. In (A,B), the value in parentheses represents the number of samples. In (C–E), the two values (e.g., “n = 19 vs. 407” for BLCA) in parentheses represent the numbers of controls and cancers, respectively. In (C–F), each p-value in this figure is based on the Wilcoxon rank-sum test, and then p-values in (C–E) are adjusted by multiple comparisons with the false discovery rate method; *p < 0.05, **p < 0.01, ***p < 0.001. BLCA, bladder cancer; ESCA, esophageal carcinoma; HNSCC, head and neck squamous cell carcinoma; KIRC, kidney renal clear cell carcinoma; LIHC, liver hepatocellular carcinoma; NSCLC, non-small-cell lung carcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma; THCA, thyroid carcinoma; UCEC, uterine corpus endometrial carcinoma; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; GBM, glioblastoma multiforme; KICH, kidney chromophobe; KIRP, kidney renal papillary cell carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PCPG, pheochromocytoma and paraganglioma; READ, rectum adenocarcinoma.

Among the 16 cancers listed earlier, BSG protein levels were investigated in ten types as there were at least two cancer samples and two normal samples for the ten. In these ten kinds of cancers, compared to corresponding normal tissues, upregulated BSG protein levels were detected in kidney cancer (for KICH and KIRP), LIHC, THCA, and UCEC, consistent with the results at mRNA level (p < 0.05; Figure 2D). Notably, the THCA and its control groups’ age distribution was statistically different (Supplementary Table S5). However, we failed to perform a hierarchical evaluation for this resulting from a lack of samples, which requires verification in the future. Although there is no statistical significance for other cancers, a trend can be seen that upregulation of BSG protein levels was identified in BLCA, BRCA, LUAD, LUSC, PRAD, and STAD (p > 0.05; Figure 2D). Under the microscope, the staining intensity of the anti-BSG antibodies in the ten cancer tissues was observed to be stronger than that in normal tissues (Figure 3A).

FIGURE 3
www.frontiersin.org

FIGURE 3. Staining intensity of the anti-BSG antibodies in the ten cancers and the correlations between BSG expression and clinical parameters. (A) Under the microscope, the staining intensity of the anti-BSG antibodies in the ten cancer tissues is observed to be stronger than that in normal tissues; images are available from v21.0.proteinatlas.org. In (A), the value of the lower-left corner of the images (e.g., 2047 for BLCA) is the identity document of patients or normal individuals. (B–D) Close correlations between BSG expression with clinical parameters; all p-values in the three panels are based on the Wilcoxon rank-sum test. The two plots with a blue curve of panel (C) are drawn based on the locally weighted regression. The value in parentheses (e.g., “n = 443” for HNSCC) represents the number of samples. BLCA, bladder cancer; BRCA, invasive breast carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PRAD, prostate adenocarcinoma; STAD, stomach adenocarcinoma; THCA, thyroid carcinoma; UCEC, uterine corpus endometrial carcinoma; HNSCC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; PAAD, pancreatic adenocarcinoma; SKCM, skin cutaneous melanoma.

Relevance Between BSG Expression and Clinical Parameters

Different clinical parameters may result in the different prognoses of cancer patients. For example, patients with cervical cancer and with advanced AJCC stages tended to have poor prognoses (Olawaiye et al., 2021). The other example was that elder age, and male gender were risk factors for patients with oral squamous cancers (Oh et al., 2021). Thus, the relationship between BSG expression and the three clinical parameters—AJCC stages, age, and gender—of patients was explored in neoplasms with the clinical parameter data. The immediate relevance of BSG expression with AJCC stages was found in HNSCC, KICH, LIHC, PAAD, SKCM, and THCA (p < 0.05; Figure 3B), while such a phenomenon was not in the other cancers (ACC, etc.) (Supplementary Figure S1). Higher BSG expression was detected in elderly (65 years old or older) BRCA patients (p < 0.05), and based on the LOESS (locally weighted regression) age trajectory, older patients tend to be found increasing the BSG expression levels (Figure 3C). Elevated BSG expression was generally be observed in young (<65 years old) PRAD patients (p < 0.05); however, the LOESS age trajectory suggested some fluctuations in the correlation between age and BSG expression (Figure 3C). There was no correlation between BSG expression and age in other 30 cancers (e.g., ACC; Supplementary Figure S2). In terms of gender, males with BRCA or KIRP had overexpression of BSG, while the contrast phenomenon was found in patients with KICH or LIHC (p < 0.05; Figure 3D); and for other cancers (ACC, etc.), there was no statistical difference in BSG expression between patients with different genders (Supplementary Figure S3).

Additionally, the adjustment tests were carried out for KICH, LIHC, and BRCA as BSG expression was detected related to no less than two clinical parameters in the three cancers. As a result, adjusted for gender, 1) there was no statistical difference in BSG expression between KICH patients with various AJCC stages for both females and males, although that was observed in LIHC (Supplementary Figure S4A); and 2) elevated BSG expression was detected in the old BRCA patients for female but not for male patients (Supplementary Figure S4B). The conclusion that male BRCA patients tend to have upregulated BSG expression was true for not the old but the young (Supplementary Figure S4C). Moreover, after adjusting AJCC stages, differential BSG expression was observed between female and male patients with KICH (for stage II) or LIHC (for stage I) (Supplementary Figure S4D).

Prognosis Value of BSG Expression in Cancers

We investigated the prognosis value of BSG expression in various cancers. Based on the results of univariate Cox analyses and Kaplan–Meier curves, in terms of both OS and DSS, the BSG expression was relevant to shorter survival time in LGG, LUAD, and UCS (hazard ratio [HR] > 1, p < 0.05) and to longer survival time in KIRP (HR = 0.364, p < 0.05) (Figures 4A–D). Moreover, high BSG expression was also a risk factor of DSS for BRCA patients (HR = 1.459, p < 0.05) (Figures 4C,D). In DFI, overexpression of BSG represented unfavorable prognosis for patients with BRCA, LGG, PRAD, and UCS (HR > 1, p < 0.05) (Figures 5A,B). In PFI, increased BSG expression was associated with poor PFI for BRCA, LGG, and UCS (HR > 1, p < 0.05), while it demonstrated a favorable prognosis for PAAD and THYM (Thymoma) patients (HR < 1, p < 0.05) (Figures 5C,D).

FIGURE 4
www.frontiersin.org

FIGURE 4. Relation of BSG expression with overall survival ( (A–B)) and disease-specific survival of cancer patients ( (C–D)). (A,C) Forest plots produced by univariate Cox regression analysis. (B,D) Kaplan–Meier curves with p-values based on log-rank tests. The optimal cut point for high and low expression levels of BSG in each Kaplan–Meier curve was identified using the maximally selected rank statistics from the “maxstat” and “survminer” R packages. In A and C, the value in parentheses (e.g., “n = 77” for ACC) represents the number of samples. ACC, adrenocortical carcinoma; BLCA, bladder cancer; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSCC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma.

FIGURE 5
www.frontiersin.org

FIGURE 5. Relation of BSG expression with disease-free survival (A–B) and progression-free survival of cancer patients (C–D). (A,C) Forest plots produced by univariate Cox regression analysis. (B,D) Kaplan–Meier curves with p-values based on log-rank tests. The optimal cut point for high and low expression levels of BSG in each Kaplan–Meier curve was identified using the maximally selected rank statistics from the “maxstat” and “survminer” R packages. In (A,C), the value in parentheses (e.g., “n = 44” for ACC) represents the number of samples. ACC, adrenocortical carcinoma; BLCA, bladder cancer; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSCC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma.

Identification Value of BSG Expression in Cancers

It is of great significance for the clinical management of cancers to be timely for determining the cancer status of patients. In 12 of 20 cancers in this study, BSG demonstrated pronounced effects in identifying cancer tissues from their control counterparts (AUC >0.7; Figure 6A). Especially in CHOL, KICH, and LIHC, the AUC values were more than 0.95 (Figure 6A), suggesting that BSG expression made it feasible to distinguish these cancers and implying the significant potential of BSG expression to predict these cancers. Moreover, in an overview of the 20 cancers, BSG expression also had the potential to distinguish the cancer patients from normal individuals (AUC = 0.89 [0.86–0.91]; Figure 6B). Thus, BSG may serve as a marker for predicting the disease status of certain cancers.

FIGURE 6
www.frontiersin.org

FIGURE 6. Ability of BSG expression to distinguish these cancer tissues from their normal tissues and a landscape of BSG mutations in pan-cancer. (A) In receiver operating characteristic curves, BSG demonstrated pronounced effects in identifying cancer tissues from their control counterparts based on the area under the curve (AUC). The two values (e.g., “n = 19 vs. 407” for BLCA) in parentheses represent the numbers of controls and cancers, respectively. (B) In the summary receiver operating characteristic curve, the BSG expression can accurately identify the 20 types of cancers samples from their controls. (C) Missense mutation is the most common alternation for BSG in pan-cancer; the panel plot was obtained from sangerBox 3.0; Del, deletion; Ins, insertion. BLCA, bladder cancer; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSCC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; STAD, stomach adenocarcinoma; THCA, thyroid carcinoma; UCEC, uterine corpus endometrial carcinoma; LGG, brain lower grade glioma.

The Landscape of BSG Mutations and the Expression Correlation of BSG With DNMTs, MMRs, and Immune Microenvironment

In the landscape of BSG mutation, missense mutation was the most common alternation, and UCEC possessed the most significant number (3.4%) of mutation (Figure 6C). The expression of three DNMTs (DNMT1, DNMT3A, and DNMT3B) was positively correlated with BSG, especially in GBM, LIHC, and LUSC (Figure 7A). The four proteins encoded by MLH1, MSH2, MSH6, and PMS2 have prominent roles in the MMR procession. The four genes and another MMR gene, EPCAM, are the characteristics of Lynch syndrome (a cancer-prone syndrome) (Baretti and Le 2018; Latham et al., 2019). Thus, the five MMRs (MLH1, MSH2, MSH6, PMS2, and EPCAM) were investigated in cancer-related studies (Deshpande et al., 2020; Liu et al., 2020), which were also selected in our research. As a result, BSG expression showed that it was significantly related to the five MMRs in multiple cancers, especially in UVM and GBM (Figure 7B), implying that BSG may affect the expression of MMRs in certain cancers, which needs experimental validation.

FIGURE 7
www.frontiersin.org

FIGURE 7. Correlations of BSG expression with DNA methylation, mismatch repair genes expression, and infiltration levels of all the six immune cells. (A) BSG expression was positively correlated with the three DNA methyltransferases. (B) BSG expression is significantly related to the five mismatch repair genes in multiple cancers. (C) BSG expression is negatively associated with nearly all six types of immune cells in the three cancers. The Spearman correlation coefficient follows the letter “ρ.” In (A–C), the value in parentheses (e.g., “n = 77” for ACC) represents the number of samples. ACC, adrenocortical carcinoma; BLCA, bladder cancer; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma.

In the aspect of immune infiltration levels, BSG expression showed a negative correlation with nearly all six types of immune cells (B cells, CD4+ T cells, CD8+ T cells, neutrophil, macrophage, and dendritic cells) in THCA, CHOL, and SKCM (p < 0.05, Figure 7C) rather than other cancers (Supplementary Figure S5). For nine kinds of cancers (BLCA, BRCA, KICH, LIHC, LUAD, LUSC, PRAD, THCA, and UCEC), both low infiltration levels of at least one of six types of immune cells and upregulated BSG expression were detected in cancer tissues (Supplementary Table S6), suggesting that the less infiltration of specific immune cells may be due to high expression of BSG. Notably, decreasing levels of all six types of cells and highly expressed BSG can be observed in LUSC (Supplementary Table S6), highlighting that BSG was likely to affect the six immune cell levels in LUSC patients. As for the ESTIMATE scores, BSG expression appeared negatively related to all of the immune, stromal, and estimate scores in THCA, CHOL, GBM, PAAD, and KICH (p < 0.05, Figure 8A). However, that was not conspicuous for other cancers (Supplementary Figure S6). In short, the biological functions of BSG in specific cancers may attribute to its negative correlation with the immune microenvironment.

FIGURE 8
www.frontiersin.org

FIGURE 8. Relevance of BSG expression with immune microenvironment scores, tumor mutational burden, microsatellite instability, and neoantigen. (A) BSG expression was positively correlated with all five cancers’ immune, stromal, and estimate scores. (B–C) BSG expression demonstrates distinct (positive or negative) correlation with tumor mutational burden and microsatellite in various cancers. (D) BSG expression is negatively associated with the count of neoantigens in specific cancers. The Spearman correlation coefficient follows the letter “ρ.” In (A,C), the value in parentheses (e.g., “n = 503” for THCA in A, and “77” for AA in (B)) represents the number of samples. ACC, adrenocortical carcinoma; BLCA, bladder cancer; BRCA, invasive breast carcinoma; CHOL, cholangiocarcinoma; COAD, colon adenocarcinoma; DLBC, lymphoid neoplasm diffuse large B-cell lymphoma; ESCA, esophageal carcinoma; GBM, glioblastoma multiforme; HNSCC, head and neck squamous cell carcinoma; KICH, kidney chromophobe; KIRC, kidney renal clear cell carcinoma; KIRP, kidney renal papillary cell carcinoma; LAML, acute myeloid leukemia; LGG, brain lower grade glioma; LIHC, liver hepatocellular carcinoma; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; MESO, mesothelioma; OV, ovarian serous cystadenocarcinoma; PAAD, pancreatic adenocarcinoma; PCPG, pheochromocytoma and paraganglioma; PRAD, prostate adenocarcinoma; READ, rectum adenocarcinoma; SARC, sarcoma; SKCM, skin cutaneous melanoma; STAD, stomach adenocarcinoma; TGCT, testicular germ cell tumors; THCA, thyroid carcinoma; THYM, thymoma; UCEC, uterine corpus endometrial carcinoma; UCS, uterine carcinosarcoma; UVM, uveal melanoma.

Correlation of BSG Expression With TMB, MSI, Neoantigen Count, and Immune Checkpoint Genes

TMB was a measurement of mutation levels in tumors. It was positively related to BSG expression in MESO and THYM and negatively correlated in BRCA and PAAD by the Spearman rank correlation analysis (p < 0.05; Figure 8B). MSI was a pattern of hypermutation caused by several mechanisms such as MMRs, and it was positively associated with BSG expression in BLCA, GBM, LGG, LUAD, LUSC, PRAD, SARC, STAD, and THCA (p < 0.05) and negatively related with DLBC and LIHC (p < 0.05) (Figure 8C). Higher TMB and MSI levels may contribute to elevated neoantigens (particularly tumor antigens) and thus trigger immune responses of the body (Jardim et al., 2021) (Schumacher, Scheper, and Kvistborg 2019). In our research, neoantigens demonstrated mild to moderate correlation with the count of neoantigens in DLBC and UCEC (p < 0.05, Figure 8D).

Previous studies have shown that tumor cells can suppress body immune response by motivating the expression of the immune checkpoints (Toor et al., 2020). Based on the correlation analysis, BSG expression was negatively (ρ < -0.3, p < 0.05) correlated with at least ten checkpoints in seven cancers (CHOL, KICH, KIRC, PAAD, THCA, THYM, and UVM) (Supplementary Figure S7).

GSEA and Drug Susceptibility of BSG in Cancers

The study investigated potential mechanisms of BSG in 32 cancers based on KEGG signaling pathways by GSEA. BSG was likely an essential factor in the occurrence and development of ESCA, KICH, KIRC, KIRP, PAAD, and THCA by complex mechanisms as it was identified to involve no fewer than five signaling pathways (Figure 9A). Totally, 39 KEGG signaling pathways were found for BSG in 32 cancers in this study (Supplementary Table S7). In at least five cancers, BSG was shown closely related to three pathways—olfactory transduction, metabolism of xenobiotics by cytochrome P450, and hematopoietic cell lineage (Supplementary Table S7) indicating the gene may play its roles in multiple human cancers. Moreover, some pathways such as drug metabolism cytochrome P450 and JAK-STAT signaling pathways were also found in the GSEA of BSG (Supplementary Table S7). For the top seven essential signaling pathways (olfactory transduction, etc.) presenting statistical significance in at least four cancers, correlations of their leading-edge genes with BSG expression were explored. For instance, for the olfactory transduction pathway demonstrating statistical significance in 15 cancers (KIRC, etc.), 21 leading-edge genes were found in the pathway. The relationship between them with BSG expression was explored in the 15 cancers. As a result, the significant (positive or negative) relevance of BSG expression with the 21 leading-edge genes was identified in certain types of the 15 cancers based on the Spearman rank correlation analysis (p < 0.05; Supplementary Figure S8A). A similar analysis was carried out for the other six essential pathways—hematopoietic cell lineage, metabolism of xenobiotics by cytochrome p450, drug metabolism cytochrome p450, JAK-STAT signaling pathway, Leishmania infection, and pentose and glucuronate interconversions; as a result, BSG expression primarily showed negative correlations with leading-edge genes of the six signaling pathways (p < 0.05; Supplementary Figures S8B–G), indicating the roles of BSG played in specific cancers may be due to the disorder of these pathways.

FIGURE 9
www.frontiersin.org

FIGURE 9. Potential mechanisms of BSG in pan-cancer and the drug sensitivity for the gene. (A) Gene set enrichment analysis suggests that BSG is identified to involve no fewer than five signaling pathways. (B) Drug sensitivity of BSG; the high-BSG expression group represented low-IC50 (half maximal inhibitory concentration) for the 12 drugs shown in the panel, implying cancer cells with a high-BSG expression may be sensitive to these drugs. *p < 0.05, **p < 0.01; p values are based on the one-sided Wilcoxon rank-sum test.

Furthermore, by using CellMiner, this research detected the sensitivity of BSG to 57 drugs approved by the Food and Drug Administration of America or validated by clinical tests. Notably, high-BSG expression was sensitive to 12 (e.g., KPT-8602 and brigatinib) of the 57 drugs (Figure 9B), suggesting the potential of these drugs in treating cancer patients with elevated BSG expression.

Discussion

Nowadays, COVID-19 has been a severe global public health problem since December 2019 (Li Q et al., 2020). Because of the susceptibility of cancer patients to the infection of SARS-CoV-2 and the occurrence of severe complications in the cancer patients, the prevalence of COVID-19 poses significant challenges to the management of cancer patients (Liang et al., 2020; Moujaess, Kourie, and Ghosn 2020). BSG is a vital entrance of SARS-CoV-2 into human cells and is upregulated in both COVID-19 patients and tumor patients (Monteiro et al., 2014; Xin et al., 2016; Bergsneider et al., 2021). Thus, there may be a linkage between COVID-19 and cancers (Varadarajan et al., 2020). However, although the other essential receptor for SARS-CoV-2 invading normal cells—angiotensin-converting enzyme 2—has received extensive attention in tumor-related research, effort on BSG is relatively limited, which needs to be further analyzed.

To promote the understanding of BSG in multiple human cancers, a systematic analysis of COVID-19-related BSG in pan-cancer was performed based on 19,020 samples in this study. Upregulated and downregulated BSG expression was observed in 14 cancers (BRCA, etc.) and two cancers (COAD and READ), respectively. The increasing BSG protein levels were also observed in kidney cancer (for KICH and KIRP), LIHC, THCA, and UCEC. BSG served as a risk role for the prognosis of patients with BRCA, LGG, LUAD, PRAD, and UCS in at least one of OS, DSS, DFI, and PFI, while it also represented a protective factor in the prognosis of individuals with KIRP, PAAD, and THYM in OS, DSS, or PFI. BSG expression showed a pronounced effect in identifying multiple cancers (e.g., CHOL) from their control samples. Moreover, the associations of BSG expression with TMB, MSI, neoantigen, and immune checkpoints suggested the potential of BSG as an exciting target for cancer treatment.

Generally, overexpression of BSG is common in many cancers. Our study revealed differences in the BSG expression levels in various human organs and cancer cells. In terms of tissues, previous reports have identified differently expressed BSG in a few cancers, such as melanoma (Kim et al., 2021; Lu et al., 2021) and LIHC (Wang SJ et al., 2020). In our study, compared with normal tissues, BSG mRNA was highly expressed in cancers containing BLCA, BRCA, CHOL, ESCA, HNSCC, KICH, KIRP, LIHC, LUAD, and LUSC, PRAD, STAD, THCA, and UCEC, while it was lowly expressed in COAD and READ. In most of these cancers, increasing (exceptionally, decreasing for COAD) BSG expression had been detected before (Han et al., 2010; Xue, Lu, and Sun 2011; Zhong et al., 2012; Liu et al., 2013; Xu et al., 2013; Huang et al., 2019; Yu et al., 2019; He et al., 2021). These reports mainly focused on protein levels of BSG, which were consistent with the mRNA results in our study. Indeed, we also validated elevated BSG protein levels in kidney cancer (for KICH and KIRP), LIHC, THCA, and UCEC. Collectively, differential (mainly increasing) BSG expression was observed in various cancers.

BSG is not only differently expressed in most tumors but represents different prognosis results and identifies cancer status in some cancers. Previously, Suzuki et al. (2021) revealed CD147 (alias BSG) as an indicator of the poor prognosis of patients with hypopharyngeal cancer. Miyazaki et al. (2020) reported that the upregulated EMMPRIN (another name for BSG) expression was also associated with a poor external auditory canal carcinoma prognosis. Our study discussed several types of survival information for further exploring the roles BSG played in the prognosis of multiple human cancers. In OS, DSS, DFI, and/or PFI, overexpressed BSG expression was associated with an unfavorable prognosis in patients with BRCA, LGG, LUAD, PRAD, and UCS. At the same time, it showed a favorable prognosis for patients with KIRP, PAAD, and THYM. Moreover, BSG made it feasible to differentiate 12 cancer (particularly CHOL, KICH, and LIHC) tissues from their normal tissues, implying its potential in screening patients with these cancers from normal individuals. This finding, to our knowledge, had not been revealed before, indicating the novelty of our study. Thus, the prognosis value and identification effect of BSG expression in cancers were conspicuous.

BSG may be an important potential target for treating various cancers. Constantly, the function of immunocytes facilitates the body’s antivirus and anti-tumor biological process. Thus, decreasing levels of immune cells will cause deterioration in COVID-19 patients and cancer patients. Indeed, previous research suggests lymphocyte reduction is an indicator of COVID-19 patients and various cancer patients (Li M et al., 2020; Tan et al., 2020). Furthermore, immunosuppression and weakened immune system were significant causes of severe disease course and high mortality in cancer patients with COVID-19 (Xia and Dubrovska 2020). Interestingly, based on this study and previous reports, BSG played multiple biological roles: 1) one of the significant known receptors of SARS-CoV-2 (Wang K et al., 2020; Fenizia et al., 2021; Geng et al., 2021; Behl et al., 2022); 2) a factor negatively related to the infiltration levels of innate and adaptive immune cells; 3) a factor conversely correlated with immune matrix levels, and 4) a prognostic indicator for multiple cancers. These results supported each other and suggested that BSG may be involved in the negative regulation of immune response. In addition, this study not only found that BSG expression was related to the DNMT expression, MMRs expression, TMB level of cancer (KICH, etc.) patients, and MSI level of cancer (BLCA, etc.) patients but also found that the expression level of BSG was related to neoantigens (for DLBC and UCEC) and the expression levels of multiple immune checkpoints. These results may promote the understanding of the potential role of BSG in tumor immunology.

This study also preliminarily explored the potential molecular mechanism of BSG in various human cancers and discussed the sensitivity of some drugs to BSG. In the GSEA, BSG was closely related to olfactory transduction, xenobiotics metabolism by cytochrome P450, hematopoietic cell lineage, drug cytochrome metabolism P450, JAK-STAT signaling pathways, Leishmania infection, and pentose and glucuronate interconversions. Indeed, previous single tumor studies had demonstrated these findings, such as Knutti et al. (2019), revealing that BSG affects the malignancy of BRCA cells by affecting Wnt and JAK/STAT signaling pathways. Another example was that BSG was also considered to be a vital molecule causing multiple cancer resistance (Landras and Mourah 2020); combined with the GSEA results of this study, this may be related to BSG participating in drug metabolism cytochrome P450 and other pathways. However, the mechanism of BSG in tumors was complex (e.g., 39 KEGG signaling pathways were found associated with BSG), which still needs further experimental studies. Based on drug sensitivity, our research also revealed that some drugs might be applied to cancer patients with high expression of BSG, but more effort is also needed to confirm these results further.

Some limitations to this study should be of concern. Above all, we failed to collect samples from patients diagnosed with both COVID-19 and cancer, and thus, we could not directly identify whether these patients with differently expressed BSG expression had poor prognosis in our study. Specimens utilized for exploring BSG protein levels are relatively limited. The potential mechanisms of BSG in multiple cancers require experimental validation.

This study revealed that COVID-19-related BSG was differently expressed in various human cancers. The gene is a prognostic marker for various cancers and has the potential to screen several cancers. In short, BSG may be a novel target for cancer treatment and identification.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

TH: data collection, data analysis, and manuscript writing. W-YH: conception, design, data analysis, manuscript revision, and supervision.

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 are deeply grateful to reviewers for their detailed work and thoughtful suggestions, which have contributed greatly to the improvement of this article. We also acknowledge Guangxi Key Pathology Laboratory for its technical support. The results shown in the article are in part based upon data generated by the GTEx, CCLE, TCGA, HPA, and Sanger Box 3.0.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.876180/full#supplementary-material

References

Baretti, M., and Le, D. T. (2018). DNA Mismatch Repair in Cancer. Pharmacol. Ther. 189, 45–62. doi:10.1016/j.pharmthera.2018.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Behl, T., Kaur, I., Aleya, L., Sehgal, A., Singh, S., Sharma, N., et al. (2022). CD147-spike Protein Interaction in COVID-19: Get the ball Rolling with a Novel Receptor and Therapeutic Target. Sci. Total Environ. 808, 152072. doi:10.1016/j.scitotenv.2021.152072

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergsneider, B., Bailey, E., Ahmed, Y., Gogineni, N., Huntley, D., and Montano, X. (2021). Analysis of SARS-CoV-2 Infection Associated Cell Entry Proteins ACE2, CD147, PPIA, and PPIB in Datasets from Non SARS-CoV-2 Infected Neuroblastoma Patients, as Potential Prognostic and Infection Biomarkers in Neuroblastoma. Biochem. Biophys. Rep. 27, 101081. doi:10.1016/j.bbrep.2021.101081

PubMed Abstract | CrossRef Full Text | Google Scholar

Dana, P., Saisomboon, S., Kariya, R., Okada, S., Obchoei, S., Sawanyawisuth, K., et al. (2020). CD147 Augmented Monocarboxylate Transporter-1/4 Expression through Modulation of the Akt-FoxO3-NF-Κb Pathway Promotes Cholangiocarcinoma Migration and Invasion. Cell Oncol. 43 (2), 211–222. doi:10.1007/s13402-019-00479-3

CrossRef Full Text | Google Scholar

Deshpande, M., Romanski, P. A., Rosenwaks, Z., and Gerhardt, J. (2020). Gynecological Cancers Caused by Deficient Mismatch Repair and Microsatellite Instability. Cancers 12 (11), 3319. doi:10.3390/cancers12113319

PubMed Abstract | CrossRef Full Text | Google Scholar

Fenizia, C., Galbiati, S., Vanetti, C., Vago, R., Clerici, M., Tacchetti, C., et al. (2021). SARS-CoV-2 Entry: At the Crossroads of CD147 and ACE2. Cells 10 (6), 1434. doi:10.3390/cells10061434

PubMed Abstract | CrossRef Full Text | Google Scholar

Geng, J., Chen, L., Yuan, Y., Wang, K., Wang, Y., Qin, C., et al. (2021). CD147 Antibody Specifically and Effectively Inhibits Infection and Cytokine Storm of SARS-CoV-2 and its Variants delta, Alpha, Beta, and Gamma. Sig Transduct Target. Ther. 6 (1), 347. doi:10.1038/s41392-021-00760-8

CrossRef Full Text | Google Scholar

Ghandi, M., Huang, F. W., Jané-Valbuena, J., Kryukov, G. V., Lo, C. C., McDonald, E. R., et al. (2019). Next-generation Characterization of the Cancer Cell Line Encyclopedia. Nature 569 (7757), 503–508. doi:10.1038/s41586-019-1186-3

PubMed Abstract | CrossRef Full Text | Google Scholar

GTEx Consortium (2013). The Genotype-Tissue Expression (GTEx) Project. Nat. Genet. 45 (6), 580–585. doi:10.1038/ng.2653

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, Z.-d., He, H.-c., Bi, X.-c., Qin, W.-j., Dai, Q.-s., Zou, J., et al. (2010). Expression and Clinical Significance of CD147 in Genitourinary Carcinomas. J. Surg. Res. 160 (2), 260–267. doi:10.1016/j.jss.2008.11.838

PubMed Abstract | CrossRef Full Text | Google Scholar

Hatanaka, M., Higashi, Y., Kawai, K., Su, J., Zeng, W., Chen, X., et al. (2016). CD147-targeted siRNA in A375 Malignant Melanoma Cells Induces the Phosphorylation of EGFR and Downregulates cdc25C and MEK Phosphorylation. Oncol. Lett. 11 (4), 2424–2428. doi:10.3892/ol.2016.4267

PubMed Abstract | CrossRef Full Text | Google Scholar

He, C., Hua, X., Sun, S., Li, S., Wang, J., and Huang, X. (2021). Integrated Bioinformatic Analysis of SARS-CoV-2 Infection Related Genes ACE2, BSG and TMPRSS2 in Aerodigestive Cancers. Jir 14, 791–802. doi:10.2147/JIR.S300127

CrossRef Full Text | Google Scholar

Huang, W. T., Yang, X., He, R. Q., Ma, J., Hu, X. H., Mo, W. J., et al. (2019). Overexpressed BSG Related to the Progression of Lung Adenocarcinoma with High-Throughput Data-Mining, Immunohistochemistry, In Vitro Validation and In Silico Investigation. Am. J. Transl Res. 11 (8), 4835–4850.

PubMed Abstract | Google Scholar

Jardim, D. L., Goodman, A., de Melo Gagliato, D., and Kurzrock, R. (2021). The Challenges of Tumor Mutational Burden as an Immunotherapy Biomarker. Cancer Cell 39 (2), 154–173. doi:10.1016/j.ccell.2020.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H.-S., Kim, H. J., Lee, M. R., and Han, I. (2021). EMMPRIN Expression Is Associated with Metastatic Progression in Osteosarcoma. BMC Cancer 21 (1), 1059. doi:10.1186/s12885-021-08774-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Knutti, N., Huber, O., and Friedrich, K. (2019). CD147 (EMMPRIN) Controls Malignant Properties of Breast Cancer Cells by Interdependent Signaling of Wnt and JAK/STAT Pathways. Mol. Cel Biochem 451 (1-2), 197–209. doi:10.1007/s11010-018-3406-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Landras, A., and Mourah, S. (2020). CD147 : rôle et ciblage thérapeutique dans les cancers. Med. Sci. (Paris)hors Serie N Degrees 361, 47–49. doi:10.1051/medsci/2020196

PubMed Abstract | CrossRef Full Text | Google Scholar

Latham, A., Srinivasan, P., Kemel, Y., Shia, J., Bandlamudi, C., Mandelker, D., et al. (2019). Microsatellite Instability Is Associated with the Presence of Lynch Syndrome Pan-Cancer. J. Clin. Oncol. 37 (4), 286–295. doi:10.1200/JCO.18.00283

PubMed Abstract | CrossRef Full Text | Google Scholar

Li M, M., Guo, W., Dong, Y., Wang, X., Dai, D., Liu, X., et al. (2020). Elevated Exhaustion Levels of NK and CD8+ T Cells as Indicators for Progression and Prognosis of COVID-19 Disease. Front. Immunol. 11, 580237. doi:10.3389/fimmu.2020.580237

PubMed Abstract | CrossRef Full Text | Google Scholar

Li Q, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., et al. (2020). Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N. Engl. J. Med. 382 (13), 1199–1207. doi:10.1056/NEJMoa2001316

PubMed Abstract | CrossRef Full Text | Google Scholar

Li T, T., Fu, J., Zeng, Z., Cohen, D., Li, J., Chen, Q., et al. (2020). TIMER2.0 for Analysis of Tumor-Infiltrating Immune Cells. Nucleic Acids Res. 48 (W1), W509–W514. doi:10.1093/nar/gkaa407

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.Can-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

Lian, C., Guo, Y., Zhang, J., Chen, X., and Peng, C. (2017). Targeting CD147 Is a Novel Strategy for Antitumor Therapy. Curr. Pharm. Des. 23 (29), 4410–4421. doi:10.2174/1381612823666170710144759

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, W., Guan, W., Chen, R., Wang, W., Li, J., Xu, K., et al. (2020). Cancer Patients in SARS-CoV-2 Infection: a Nationwide Analysis in China. Lancet Oncol. 21 (3), 335–337. doi:10.1016/s1470-2045(20)30096-6

PubMed Abstract | CrossRef Full Text | Google Scholar

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. e11. doi:10.1016/j.cell.2018.02.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Zhang, S., Dai, W., Xie, C., and Li, J.-C. (2020). A Comprehensive Prognostic and Immune Analysis of SLC41A3 in Pan-Cancer. Front. Oncol. 10, 586414. doi:10.3389/fonc.2020.586414

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Xin, T., Jiang, Q.-y., Huang, D.-y., Shen, W.-x., Li, L., et al. (2013). CD147, MMP9 Expression and Clinical Significance of Basal-like Breast Cancer. Med. Oncol. 30 (1), 366. doi:10.1007/s12032-012-0366-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, L., Zhang, J., Gan, P., Wu, L., Zhang, X., Peng, C., et al. (2021). Novel Functions of CD147 in the Mitochondria Exacerbates Melanoma Metastasis. Int. J. Biol. Sci. 17 (1), 285–297. doi:10.7150/ijbs.52043

PubMed Abstract | CrossRef Full Text | Google Scholar

Madariaga, A., McMullen, M., Sheikh, S., Kumar, R., Liu, F.-F., Zimmermann, C., et al. (2020). COVID-19 Testing in Patients with Cancer: Does One Size Fit All? Clin. Cancer Res. 26 (18), 4737–4742. doi:10.1158/1078-0432.Ccr-20-2224

PubMed Abstract | CrossRef Full Text | Google Scholar

Miyazaki, M., Aoki, M., Okado, Y., Koga, K., Hamasaki, M., Nakagawa, T., et al. (2020). Highly Expressed Tumoral Emmprin and Stromal CD73 Predict a Poor Prognosis for External Auditory Canal Carcinoma. Cancer Sci. 111 (8), 3045–3056. doi:10.1111/cas.14508

PubMed Abstract | CrossRef Full Text | Google Scholar

Monteiro, L. S., Delgado, M. L., Ricardo, S., Garcez, F., Amaral, B. d., Pacheco, J. J., et al. (2014). EMMPRIN Expression in Oral Squamous Cell Carcinomas: Correlation with Tumor Proliferation and Patient Survival. Biomed. Res. Int. 2014, 1–9. doi:10.1155/2014/905680

CrossRef Full Text | Google Scholar

Moujaess, E., Kourie, H. R., and Ghosn, M. (2020). Cancer Patients and Research during COVID-19 Pandemic: A Systematic Review of Current Evidence. Crit. Rev. Oncology/Hematology 150, 102972. doi:10.1016/j.critrevonc.2020.102972

PubMed Abstract | CrossRef Full Text | Google Scholar

Oh, L. J., Asher, R., Veness, M., Smee, R., Goldstein, D., Gopalakrishna Iyer, N., et al. (2021). Effect of Age and Gender in Non-smokers with Oral Squamous Cell Carcinoma: Multi-Institutional Study. Oral Oncol. 116, 105210. doi:10.1016/j.oraloncology.2021.105210

PubMed Abstract | CrossRef Full Text | Google Scholar

Olawaiye, A. B., Baker, T. P., Washington, M. K., and Mutch, D. G. (2021). The New (Version 9) American Joint Committee on Cancer Tumor, Node, Metastasis Staging for Cervical Cancer. CA A. Cancer J. Clin. 71 (4), 287–298. doi:10.3322/caac.21663

CrossRef Full Text | Google Scholar

Qiao, X., Gu, Y., Yu, J., Wang, J., Liu, X., Gu, M., et al. (2020). The Combination of CD147 and MMP-9 Serum Levels Is Identified as Novel Chemotherapy Response Markers of Advanced Non-small-cell Lung Cancer. Dis. Markers 2020, 1–10. doi:10.1155/2020/8085053

PubMed Abstract | CrossRef Full Text | Google Scholar

Reinhold, W. C., Sunshine, M., Varma, S., Doroshow, J. H., and Pommier, Y. (2015). Using CellMiner 1.6 for Systems Pharmacology and Genomic Analysis of the NCI-60. Clin. Cancer Res. 21 (17), 3841–3852. doi:10.1158/1078-0432.CCR-15-0335

PubMed Abstract | CrossRef Full Text | Google Scholar

Robin, X., Turck, N., Hainard, A., Tiberti, N., Lisacek, F., Sanchez, J.-C., et al. (2011). pROC: an Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinformatics 12, 77. doi:10.1186/1471-2105-12-77

PubMed Abstract | CrossRef Full Text | Google Scholar

Schumacher, T. N., Scheper, W., and Kvistborg, P. (2019). Cancer Neoantigens. Annu. Rev. Immunol. 37, 173–200. doi:10.1146/annurev-immunol-042617-053402

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, S., Toyoma, S., Kawasaki, Y., Nanjo, H., and Yamada, T. (2021). CD147 Promotes Invasion and MMP-9 Expression through MEK Signaling and Predicts Poor Prognosis in Hypopharyngeal Squamous Cell Carcinoma. Adv. Clin. Exp. Med. 30 (1), 41–48. doi:10.17219/acem/128228

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, L., Wang, Q., Zhang, D., Ding, J., Huang, Q., Tang, Y.-Q., et al. (2020). Lymphopenia Predicts Disease Severity of COVID-19: a Descriptive and Predictive Study. Sig Transduct Target. Ther. 5 (1), 33. doi:10.1038/s41392-020-0148-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Toor, S. M., Sasidharan Nair, V., Decock, J., and Elkord, E. (2020). Immune Checkpoints in the Tumor Microenvironment. Semin. Cancer Biol. 65, 1–12. doi:10.1016/j.semcancer.2019.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, A. N., Dussaq, A. M., Kennell, T., Willey, C. D., and Hjelmeland, A. B. (2019). HPAanalyze: an R Package that Facilitates the Retrieval and Analysis of the Human Protein Atlas Data. BMC Bioinformatics 20 (1), 463. doi:10.1186/s12859-019-3059-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Tseng, H.-c., Xiong, W., Badeti, S., Yang, Y., Ma, M., Liu, T., et al. (2020). Efficacy of Anti-cd147 Chimeric Antigen Receptors Targeting Hepatocellular Carcinoma. Nat. Commun. 11 (1), 4810. doi:10.1038/s41467-020-18444-2

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Varadarajan, S., Balaji, T. M., Sarode, S. C., Sarode, G. S., Sharma, N. K., Gondivkar, S., et al. (2020). EMMPRIN/BASIGIN as a Biological Modulator of Oral Cancer and COVID-19 Interaction: Novel Propositions. Med. Hypotheses 143, 110089. doi:10.1016/j.mehy.2020.110089

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang K, K., Chen, W., Zhang, Z., Deng, Y., Lian, J.-Q., Du, P., et al. (2020). CD147-spike Protein Is a Novel Route for SARS-CoV-2 Infection to Host Cells. Sig Transduct Target. Ther. 5 (1), 283. doi:10.1038/s41392-020-00426-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang Sj, S.-J., Chao, D., Wei, W., Nan, G., Li, J.-Y., Liu, F.-L., et al. (2020). CD147 Promotes Collective Invasion through Cathepsin B in Hepatocellular Carcinoma. J. Exp. Clin. Cancer Res. 39 (1), 145. doi:10.1186/s13046-020-01647-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Weidle, U. H., Scheuer, W., Eggle, D., Klostermann, S., and Stockinger, H. (2010). Cancer-related Issues of CD147. Cancer Genomics Proteomics 7 (3), 157–169.

PubMed Abstract | Google Scholar

Xia, P., and Dubrovska, A. (2020). Tumor Markers as an Entry for SARS‐CoV‐2 Infection? FEBS J. 287 (17), 3677–3680. doi:10.1111/febs.15499

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, X., Zeng, X., Gu, H., Li, M., Tan, H., Jin, Z., et al. (2016). CD147/EMMPRIN Overexpression and Prognosis in Cancer: A Systematic Review and Meta-Analysis. Sci. Rep. 6, 32804. doi:10.1038/srep32804

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X.-y., Lin, N., Li, Y.-m., Zhi, C., and Shen, H. (2013). Expression of HAb18G/CD147 and its Localization Correlate with the Progression and Poor Prognosis of Non-small Cell Lung Cancer. Pathol. - Res. Pract. 209 (6), 345–352. doi:10.1016/j.prp.2013.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, Y.-J., Lu, Q., and Sun, Z.-X. (2011). CD147 Overexpression Is a Prognostic Factor and a Potential Therapeutic Target in Bladder Cancer. Med. Oncol. 28 (4), 1363–1372. doi:10.1007/s12032-010-9582-4

PubMed Abstract | CrossRef Full Text | Google Scholar

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, 2612. doi:10.1038/ncomms3612

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, B., Zhang, Y., Wu, K., Wang, L., Jiang, Y., Chen, W., et al. (2019). CD147 Promotes Progression of Head and Neck Squamous Cell Carcinoma via NF-Kappa B Signaling. J. Cel Mol Med 23 (2), 954–966. doi:10.1111/jcmm.13996

CrossRef Full Text | Google Scholar

Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, W.-d., Liang, Y.-x., Lin, S. X., Li, L., He, H.-c., Bi, X.-c., et al. (2012). Expression of CD147 Is Associated with Prostate Cancer Progression. Int. J. Cancer 130 (2), 300–308. doi:10.1002/ijc.25982

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: cancer, prognosis, prediction, immunology, target therapy

Citation: Huang T and He W-Y (2022) Pan-Cancer Analysis, Reveals COVID-19-Related BSG as a Novel Marker for Treatment and Identification of Multiple Human Cancers. Front. Cell Dev. Biol. 10:876180. doi: 10.3389/fcell.2022.876180

Received: 15 February 2022; Accepted: 25 March 2022;
Published: 13 May 2022.

Edited by:

José Díaz-Chávez, Instituto Nacional de Cancerología (INCAN), Mexico

Reviewed by:

Paula Araya, University of Colorado Denver, United States
Fuhai Li, Washington University in St. Louis, United States

Copyright © 2022 Huang and He. 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: Wei-Ying He, 934333178@qq.com

Disclaimer: 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.