Next Article in Journal
New Insight into Drugs to Alleviate Atopic March via Network Pharmacology-Based Analysis
Next Article in Special Issue
Sonicated Bordetella bronchiseptica Bacterin Can Protect Dendritic Cells from Differential Cytotoxicity Caused by Doxorubicin and Vincristine and Enhance Their Antigen-Presenting Capability
Previous Article in Journal
Antigenotoxic Effect of Ascorbic Acid and Resveratrol in Erythrocytes of Ambystoma mexicanum, Oreochromis niloticus and Human Lymphocytes Exposed to Glyphosate
Previous Article in Special Issue
Effects of Luteolin on Human Breast Cancer Using Gene Expression Array: Inferring Novel Genes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Next-Generation Sequencing and Bioinformatic Methods to Predict New Genes That May Be Regulated by CD47 in Oral Squamous Cell Carcinoma

1
Institute of Medical Science and Technology, National Sun Yat-sen University, Kaohsiung 80424, Taiwan
2
Department of Dentistry, Zuoying Branch of Kaohsiung Armed Forces General Hospital, Kaohsiung 81342, Taiwan
3
College of Ocean Food and Biological Engineering, Jimei University, Xiamen 361021, China
4
Department of Medical Research, E-Da Hospital/E-Da Cancer Hospital, I-Shou University, Kaohsiung 82445, Taiwan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2022, 44(5), 2243-2256; https://doi.org/10.3390/cimb44050152
Submission received: 17 February 2022 / Revised: 13 April 2022 / Accepted: 9 May 2022 / Published: 17 May 2022
(This article belongs to the Special Issue Targeting Tumor Microenvironment for Cancer Therapy)

Abstract

:
Oral squamous cell carcinoma (OSCC) is one of the most common cancers in the world, and the incidence and death rate of OSCC in men is twice that of women. CD47 is a ubiquitous cell surface transmembrane protein, also known as integrin-related protein (IAP). Previous studies have pointed out that CD47 can inhibit the growth of OSCC, but the detailed mechanism is not clear. This study aimed to explore the effect of CD47 gene expression profiles in OSCC. The OSCC cell lines, OECM-1 and OC-2, overexpressed CD47, and the expression profiles of mRNAs were analyzed through next-generation sequencing (NGS) with a bioinformatic approach. A total of 14 differentially expressed genes (DEGs) were listed. In addition, ingenuity pathway analysis (IPA) was used to analyze the molecular function (MF), biological process (BP), and cellular component (CC) network signaling. The human protein atlas (HPA) database was used to analyze gene expression and the survivability of human cancer. The results found that HSPA5, HYOU1, and PDIA4 were involved in the IPA network and when highly expressed, mediated the survivability of cancer. In addition, HSPA5 was positively and significantly correlated with CD47 expression (p < 0.0001) and induced by CD47-overexpression in the OECM-1 and OC-2 OSCC cancer cell lines. These findings provide important insights into possible new diagnostic strategies, including unfolded protein for OSCC-targeting CD47.

1. Introduction

Oral cancer is one of the most common cancers in the world. It is estimated that 657,000 new cases of oral and throat cancer are generated every year, and more than 330,000 people die. The incidence and death rate of oral cancer in men is twice that of women [1]. Normally, oral cancer is a group of head and neck cancers, formed on the mucosal surface of the lips, hard palate, posterior molar triangle, and mandible [2]. The most common oral cancer tissue type is squamous differentiated carcinoma, caused by epithelial mucosa, also known as oral squamous cell carcinoma (OSCC), which accounts for more than 90% of all oral cancers [3]. Oral cancer can be cured in its early stage by surgery, radiotherapy, and chemoradiation. However, most patients with OSCC already have advanced oral cancer when diagnosed. At this time, the treatment effect and prognosis are poor, and the mortality rate is high [4,5,6]. Studies have shown that the occurrence of oral cancer is related to multiple risk factors including tobacco use, drinking alcohol, betel chewing, human papillomavirus (HPV) [7], a weakened immune system, a lack of nutrition in the diet, and a poor lifestyle. In particular, betel chewing, excessive alcohol consumption, and smoking cause an increased incidence and mortality of OSCC [8,9].
CD47 is a ubiquitous cell surface transmembrane protein, also known as integrin-related protein (IAP), which belongs to the immunoglobulin superfamily (IgSF). It can interact with a variety of cell surface receptors including integrins (such as αvβ3, αIIbβ3, and α2β1), thrombospondin (TSP-1), and signal regulatory protein α (SIRPα), for the signal initiation related pathway [10,11]. For example, CD47 is the ligand of SIRPα, which is a transmembrane receptor of the immunoglobulin (Ig) superfamily with extracellular immunoglobulin-like domains and structural antigen receptors. In bladder cancer, CD47 initiates an immune checkpoint through the SIRPα signaling pathway [12]. In addition, CD47 can also be used as a receptor of TSP-1, which is secreted by several non-hematopoietic cells such as platelets, monocytes, and macrophages. CD47 is involved in a variety of physiological functions including leukocyte adhesion and migration, T cell activation, apoptosis, and phagocytosis [13,14,15]. Related studies have shown that the expression of CD47 increases in migrating hematopoietic stem cells (HSCs), which may provide protection from an important feature of macrophage killing [16]. Therefore, the expression level of CD47 can predict the possibility of HSCs being swallowed during circulation. More evidence indicates that CD47 is a common mechanism found in human solid tumor cells [17,18,19]. Through these mechanisms, tumor cells can protect themselves from phagocytosis, which results in proliferation and metastasis [20].
During immunotherapy, we found that chimeric antigen receptor (CAR)-T cells that bind CD47 antigen specifically kill different types of cancer including ovarian and pancreatic cancer cells [21]. Interestingly, a previous study also found that the effect of radiotherapy can be activated by targeted immunotherapy, as silencing of CD47 and HER2 eliminates radioresistant breast cancer cells [22]. In addition, previous studies have found that high expression of CD47 was detected in OSCC cell lines and in vivo results have shown that CD47 was significantly higher in OSCC than in proximal normal tissues. Furthermore, the poor prognosis in patients with OSCC has been connected to the CD47-SIRPα signaling pathway [23].
In the present study, we examined gene expression profiles by NGS and associated these data with gene network signaling through bioinformatic analyses of OSCC cell lines. We found that unfolded proteins HSPA5, HYOU1, and PDIA4 are induced by CD47 within the OSCC IPA network.

2. Materials and Methods

2.1. Cell Culture

Two human OSCC cell lines (OECM-1 and OC-2) were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). The cancer cell lines were cultured in DMEM/F12 (Gibco, Courtaboeuf, France) supplemented with 10% fetal bovine serum (FBS) (Gibco) and 1% penicillin/streptomycin at 5% CO2 and in a 37 °C humidified incubator.

2.2. Plasmid Transfection

For plasmid transfection, the OSCC cell lines were seeded at a density of 1 × 105 cell/well in a 6-well plate. After 24 h, 1 μg of the pcDNA3.1 (empty control), pcDNA3.1-CD47 (Addgene, Watertown, MA, USA, #65473) and HSPA5 plasmid (Addgene, Watertown, MA, USA, #27165) were incubated with TurboFect (Fisher Scientific, Waltham, MA, USA) transfection reagent according to the manufacturer’s instructions. TurboFect transfection reagent and plasmid were mixed and added to the cell.

2.3. RNA Extraction and qPCR

The total RNA was extracted from cell lines by the TRIzol reagent (Invitrogen, Carlsbad, CA, USA), and cDNA was synthesized using a cDNA Synthesis kit (Promega, Madison, WI, USA) according to the manufacturer’s instructions. The cDNA was PCR-amplified with the primers for CD47 (forward: 5′-AGA TCC GGT GGT ATG GAT GAG A-3′; reverse: 5′-GTC ACA ATT AAA CCA AGG CCA GTA G-3′), HSPA5 (forward: 5′-CTG GCA AGA TGA AGC TCT CC-3′; reverse: 5′-AAA ACC CGA CAG AGG GAC AT-3′), HYOU1 (forward: 5′- GAC TTC GGC ATC TGA GTG GT-3′; reverse: 5′-GCT CCC AAG TCC ACC ATT AC-3′), PDIA4 (forward: 5′- GCT CAG CTC CAG GGA GAG -3′; reverse: 5′- GAT GAT CTC CAC CCA CCT GT-3′). β-actin (forward: 5′-ATG ATA TCG CCG CGC TCG TCG TC-3′; reverse: 5′-CGC TCG GCC GTG GTG GTG AA-3′) was used as the reference gene for normalization. qPCR values were analyzed by the ABI-7500 system (Applied Biosystems, Foster City, CA, USA) and the values were calculated using the 2−ΔΔCT method.

2.4. RNA-Seq Quantification

Differential gene expression (DEG) was analyzed by RNA-seq quantification following the protocol of Illumina (Genomics, Taipei, Taiwan). Briefly, Poly-T oligo-attached magnetic beads were used to purify the mRNA and synthesize double-strand cDNA, which was then sequenced on the Illumina NextSeq 500 platform, which obtained approximately 10 million reads for each sample. The gene expression profiles were analyzed and compared between the CD47-overexpressing and control OECM-1 and OC-2 oral cancer cell lines. All genes showing differential expression were further categorized using Gene Ontology (GO) and DEG report analysis.

2.5. Ingenuity Pathway Analysis

The IPA software (Ingenuity Systems, Redwood City, CA, USA) comprehensively explores DEGs to determine their biological function and signaling pathway. The DEGs were mapped to the IPA genetic network, then the score was calculated and sorted. The score was calculated by computing a statistic for each biological function according to a model that assigns random adjustment directions. The IPA software was used to conduct network analysis of CD47 downstream genes and to explore the underlying genetic networks.

2.6. Database for Annotation, Visualization and Integrated Discovery Analysis

DAVID analysis is a powerful tool for classifying genes through their function. It integrates data and calculates similarities using the global annotation map derived from multiple functional annotation databases including the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) biological processes.

2.7. Cell Viability Assay

Cell viability was analyzed using Cell Counting Kit 8 (CCK-8) (Dojindo, Kumamoto, Japan). The 1 × 104 cells were seeded into 96-well plates, 24 h after transfection. Cells were incubated with 10 μL of CCK-8 solution at 37 °C for 1 h. The OD value was measured at 450 nm with a microplate reader (Multiskan GO, Thermo Scientific, Carlsbad, CA, USA).

2.8. Statistical Analysis

The gene expression was compared between cells with CD47-overexpression and normal cells. All statistical analyses were performed using the GraphPad Prism statistical software (GraphPad Software, San Diego, CA, USA).

3. Results

3.1. Gene Expression Profiling and Ovexpression of CD47 in OSCC Cell Lies

Previous studies showed that CD47-overexpression may regulate the growth process of oral cancer [23]. In this study, we mainly focused on Asian patients suffering from OSCC. OECM-1 cells were derived from the surgical resection of a primary tumor from a Taiwanese patient with a history of betel quid chewing. OC-2 cells were derived from a primary tumor of the buccal mucosa of a Chinese patient with a habit of alcohol drinking, betel quid chewing, and cigarette smoking [24,25]. We used Illumina NextSeq 500 sequencing and bioinformatic studies to analyze the genes in the OECM-1 and OC-2 cells and whether they were regulated by CD47 (Figure 1A). In addition, CD47 was overexpressed in OECM-1 and OC-2 cells (Figure 1B). The RNA was extracted from the cells and the gene profile was analyzed by NGS.

3.2. Gene Expression Profile in CD47-Overexpressing Oral Cancer

Figure 2A shows the differentially upregulated (right panel) and downregulated (left panel) gene expression in CD47-overexpressing OECM-1 and OC-2 oral cancer cell lines. Genes with-log10 (p value) and ≥5-fold and ≤0.2-fold changes were chosen for further analyses (Figure 2B). Our results identified that the expression of 16 upregulated and 16 downregulated genes was found to be associated with CD47-overexpression in OECM-1 and OC-2 cells compared to the empty control (Figure 2C) (Table 1).

3.3. Biological Process Analysis of Differentially Expressed Genes in CD47-Overexpressing Oral Squamous Cell Carcinoma

Simultaneously, we executed an intergroup comparison analysis of DEGs in each comparison, determined by RNA-seq by expectation-maximization (RSEM) and differential expression analysis at both gene and isoform level using RNA-seq data (EBSeq) (posterior probability of equal expression (PPEE) ≤ 0.05). In total, 14 genes were commonly identified as DEGs in OECM-1 and OC-2 oral cancer cells with CD47-overexpression (Figure 3A). Furthermore, we used a heatmap to display the expression profiles of union sets (Figure 3B,C).
In addition, we used IPA to analyze the DEGs regulated by CD47-overexpression in oral cancer for network analysis including molecular functions (MFs), biological processes (BPs), and cellular components (CCs). Then, we focused on the upregulated genes by CD47-overexpression in the OECM-1 and OC-2 cell lines. We found six genes (SDF2L1, HSPA5, CALMFR, HYOU1, SESN2, and PDIA4) involved in molecular functions, 13 genes (CHAC1, CTH, DDIT3, HERPUD1, HSPA5, SESN2, XBP1, SDF2L1, DERL3, ASNS, HYOU1, CALR, and PDIA4) involved in biological processes, and eight genes (HSPA5, SDF2L1, CALR, SLC3A2, SLC1A5, PDIA4, SLC1A4, and HYOU1) involved in cellular components (Figure 4A). Furthermore, the results indicate that HSPA5, HYOU1, and PDIA4 genes appear in all three categories (Figure 4B).
To further explore the role of CD47, HSPA5, HYOU1, and PDIA4 genes in cancer, we used the human protein atlas (HPA) database [26] to analyze the expression of these target genes in cancer and their impact on the ability of cells to survive. We know that OSCC occupies a high proportion (90%) of head and neck squamous cell carcinomas [27]. Therefore, we selected head and neck squamous cell carcinoma samples from 499 patients for gene expression and survival analysis using the HPA database. The results showed that the fragments per kilobase per million (FPKM) of CD47, HSPA5, HYOU1, and PDIA4 was 22.19 ± 10.51, 224.2 ± 87.32, 25.35 ± 13.45, 84.83 ± 37.3, respectively. The increased expression of HSPA5 (cut off quantile 25%), HYOU1 (cut off quantile 75%), and PDIA4 (cut off quantile 25%) was significantly associated with survival rate in all 499 patients (Figure 4C). In addition, HSPA5 was highly expressed in cancer and positively and significantly correlated with CD47 (R = 0.8432, p < 0.0001, Figure 4D).
To prove that HSPA5 was associated with CD47-overexpression in the OECM-1 and OC-2 cell lines, we used qPCR to quantify gene expression. We found that the expression of HSPA5 was induced by CD47-overexpression in the OECM-1 (Figure 4E) and OC-2 (Figure 4F) cell lines. We also verified the effect of HSPA5 overexpression in OSCC cell lines. HSPA5 plasmid was transfected into the cells and HSPA5 expression was analyzed by qPCR (Figure 4G). In addition, cell viability assay revealed that HSPA5 increased cell growth in the OECM-1 and OC-2 cell lines (Figure 4H,I). Next, we analyzed the CD47 expression in HSPA5-overexpressing OSCC cell lines. The results found that CD47 expression was induced by HSPA5-overexpression in both in cell lines (Figure 4J). These data suggest that HSPA5, HYOU1 and PDIA4 were significantly associated with survival rate and HSPA5 is a suitable prognostic marker with CD47 in cancer survivability.

3.4. The DAVID Analyses of Differentially Expressed Genes following CD47-Overexpression

We also used the DAVID biological process analysis to analyze DEGs with CD47-overexpressing OECM-1 and OC-2 cell lines from RNA sequencing. We concluded that the top 20 biological processes of the DEGs in the OECM-1 and OC-2 cell lines and the top five biological processes of these genes include responses to topologically incorrect proteins (13 genes), endoplasmic reticulum stress (14 genes), and unfolded proteins (12 genes), as well as cellular responses to topologically incorrect proteins (11 genes), endoplasmic reticulum unfolded protein responses (10 genes), and endoplasmic reticulum unfolded protein responses (10 genes) in OECM-1 cells. Meanwhile, in OC-2 cells, these genes include responses to viruses (41 genes), defense responses to other organisms (39 genes), defense responses to the virus (38 genes), the regulation of symbiosis, encompassing mutualism through parasitism (29 genes), and the regulation of viral processes (29 genes) (Figure 5).

4. Discussion

Based on previous studies, we know that CD47 is abundantly expressed in many tumors, and it plays an important role in macrophages, affecting the development of cancer [28,29]. In the current study, we identified that the differentially expressed genes in OSCC cells are enriched in the biological functions of CD47, including regulating the response of unfolded proteins by NGS and systematic bioinformatic analysis. The results showed that the unfolded proteins HSPA5, HYOU1, and PDIA4 were associated with survivability in cancer, and HSPA5 was positively and significantly correlated with CD47, and induced by CD47 in OSCC cells lines.
Heat shock protein family A (Hsp70) member 5 (HSPA5) is a member of the molecular chaperones, which regulates the response of unfolded proteins [30]. HSPA5 promotes cancer progression, drug resistance, and metastasis and is a poor prognosis marker [31,32]. Hypoxia and nutrient deprivation in the environment activate the UPR and upregulate HSPA5 expression, leading to the rapid growth of solid tumors [33]. In addition, knockout mouse models have demonstrated that HSPA5 is required for cell proliferation and is critical for embryonic cell growth and survival [34]. In vitro and in vivo studies found that silencing of HSPA5 expression in cancer cells inhibits tumor cell metastasis [35]. Therefore, we believe that HSPA5 plays an important role in cancer progression.
Hypoxia-upregulated protein 1 (HYOU1) belongs to the heat shock protein 70 family, which is overexpressed in many tumors. Highly expressed HYOU1 induces tumor invasion, results in a poor prognosis, and inhibits apoptosis [36,37,38]. A previous study found that heat shock protein modulation primary macrophages are polarized into classic (M1) macrophages and alternative (M2) macrophages. The Hsp70 family might be strictly involved in the effector stages of macrophage activity [39]. Previous research found that anti-CD47 treatment enhanced the phenotype of macrophages converted to M1 subtype, and produced a higher phagocytosis rate through M1 macrophages [40]. In addition, CD47 secreted significant levels of IL-10 via M2 macrophages differentiation to promote cancer cell migration [41]. According to the above results, we believe that CD47 plays an important role in the regulation of HSP. CD47 has the same ability to promote cancer progression [42], drug resistance [43], poor prognosis marker [44] and macrophage polarization [40] as HSP. In addition, CD47 as a transcription factor [45] can mediate gene expression. Our data also found that overexpression of CD47 induced the expression of HSPA5, and overexpression of HSPA5 induced the expression of CD47 in the OECM-1 and OC-2 cells. A previous study also demonstrated that inhibition of glucose-regulated protein 78 (GRP78), a member of the HSP70 family, downregulated CD47 expression in tumor cells [46]. Therefore, we believe that CD47 and HSPA5 have a feedback loop in OSCC cell lines, and suggest that CD47 mediates cancer progression through the HSPA5 signaling pathway in oral squamous cell carcinoma. In the future, we will further explore the mechanisms and regulatory capabilities of CD47 and HSPA5.
PDIA4, as a member of the protein disulfide isomerases (PDI) family, is a key factor for protein folding in the endoplasmic reticulum. A previous study found that PDIA4 participates in tumor progression by influencing cell apoptosis and DNA repair mechanisms [47]. Meanwhile, a recent study pointed out that PDIA4 may be regarded as a biomarker for the treatment of ovarian cancer [48]. Interestingly, PDIA4 has similar properties to HSP proteins, and both have the function of affecting protein folding. Protein folding is important for CD47 evasion of macrophage phagocytosis [49]. However, the detailed mechanism is unclear at present.

5. Conclusions

These results provide important insights into possible new diagnostic genes for OSCC-targeting CD47. Although the mechanism by which CD47 regulates unfolded proteins is still unclear, the detailed signaling pathway will be discussed in more depth in the future.

Author Contributions

Conceptualization, T.-H.H. and C.-C.T.; methodology, S.-Y.H. and C.-H.T.; validation, C.-W.W.; writing—original draft preparation, T.-H.H. and C.-W.W.; writing—review and editing, T.-H.H.; supervision, T.-H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of Taiwan (grant number MOST 110-2314-B-650-004), and by Kaohsiung Armed Forces General Hospital ZuoYing RESEARCH PROJECT (KAFGH-ZY-A-109023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hung, L.C.; Kung, P.T.; Lung, C.H.; Tsai, M.H.; Liu, S.A.; Chiu, L.T.; Huang, K.H.; Tsai, W.C. Assessment of the Risk of Oral Cancer Incidence in A High-Risk Population and Establishment of A Predictive Model for Oral Cancer Incidence Using A Population-Based Cohort in Taiwan. Int. J. Environ. Res. Public Health 2020, 17, 665. [Google Scholar] [CrossRef] [Green Version]
  2. Dhanuthai, K.; Rojanawatsirivej, S.; Thosaporn, W.; Kintarak, S.; Subarnbhesaj, A.; Darling, M.; Kryshtalskyj, E.; Chiang, C.P.; Shin, H.I.; Choi, S.Y.; et al. Oral cancer: A multicenter study. Med. Oral. Patol. Oral Cir. Bucal. 2018, 23, e23–e29. [Google Scholar] [CrossRef]
  3. Manikandan, M.; Deva Magendhra Rao, A.K.; Arunkumar, G.; Manickavasagam, M.; Rajkumar, K.S.; Rajaraman, R.; Munirajan, A.K. Oral squamous cell carcinoma: MicroRNA expression profiling and integrative analyses for elucidation of tumourigenesis mechanism. Mol. Cancer 2016, 15, 28. [Google Scholar] [CrossRef] [Green Version]
  4. Pignon, J.P.; le Maitre, A.; Maillard, E.; Bourhis, J.; Group, M.-N.C. Meta-analysis of chemotherapy in head and neck cancer (MACH-NC): An update on 93 randomised trials and 17,346 patients. Radiother. Oncol. 2009, 92, 4–14. [Google Scholar] [CrossRef] [PubMed]
  5. Dzebo, S.; Mahmutovic, J.; Erkocevic, H.; Foco, F. Frequency of Depression and its Correlation with Quality of Life of Patients with Oral Cavity Cancer. Mater. Sociomed. 2017, 29, 97–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Manasa, V.G.; Kannan, S. Impact of microRNA dynamics on cancer hallmarks: An oral cancer scenario. Tumour. Biol. 2017, 39, 1010428317695920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yete, S.; D’Souza, W.; Saranath, D. High-Risk Human Papillomavirus in Oral Cancer: Clinical Implications. Oncology 2018, 94, 133–141. [Google Scholar] [CrossRef] [PubMed]
  8. Santos, H.B.; dos Santos, T.K.; Paz, A.R.; Cavalcanti, Y.W.; Nonaka, C.F.; Godoy, G.P.; Alves, P.M. Clinical findings and risk factors to oral squamous cell carcinoma in young patients: A 12-year retrospective analysis. Med. Oral. Patol. Oral Cir. Bucal. 2016, 21, e151–e156. [Google Scholar] [CrossRef] [PubMed]
  9. Polesel, J.; Franceschi, S.; Talamini, R.; Negri, E.; Barzan, L.; Montella, M.; Libra, M.; Vaccher, E.; Franchin, G.; La Vecchia, C.; et al. Tobacco smoking, alcohol drinking, and the risk of different histological types of nasopharyngeal cancer in a low-risk population. Oral Oncol. 2011, 47, 541–545. [Google Scholar] [CrossRef]
  10. Brown, E.J.; Frazier, W.A. Integrin-associated protein (CD47) and its ligands. Trends Cell. Biol. 2001, 11, 130–135. [Google Scholar] [CrossRef]
  11. Van Beek, E.M.; Cochrane, F.; Barclay, A.N.; van den Berg, T.K. Signal regulatory proteins in the immune system. J. Immunol. 2005, 175, 7781–7787. [Google Scholar] [CrossRef] [Green Version]
  12. Kiss, B.; van den Berg, N.S.; Ertsey, R.; McKenna, K.; Mach, K.E.; Zhang, C.A.; Volkmer, J.P.; Weissman, I.L.; Rosenthal, E.L.; Liao, J.C. CD47-Targeted Near-Infrared Photoimmunotherapy for Human Bladder Cancer. Clin. Cancer Res. 2019, 25, 3561–3571. [Google Scholar] [CrossRef] [Green Version]
  13. Nakaishi, A.; Hirose, M.; Yoshimura, M.; Oneyama, C.; Saito, K.; Kuki, N.; Matsuda, M.; Honma, N.; Ohnishi, H.; Matozaki, T.; et al. Structural insight into the specific interaction between murine SHPS-1/SIRP alpha and its ligand CD47. J. Mol. Biol. 2008, 375, 650–660. [Google Scholar] [CrossRef]
  14. Oldenborg, P.A.; Sundqvist, I.M.; Sehlin, J. Different effects of glucose on extracellular and intracellular respiratory burst response in normal human neutrophils activated with the soluble agonist fMet-Leu-Phe. Diabet. Med. 2000, 17, 532–537. [Google Scholar] [CrossRef]
  15. Bornstein, P. Diversity of function is inherent in matricellular proteins: An appraisal of thrombospondin 1. J. Cell. Biol. 1995, 130, 503–506. [Google Scholar] [CrossRef] [Green Version]
  16. Jaiswal, S.; Jamieson, C.H.; Pang, W.W.; Park, C.Y.; Chao, M.P.; Majeti, R.; Traver, D.; van Rooijen, N.; Weissman, I.L. CD47 is upregulated on circulating hematopoietic stem cells and leukemia cells to avoid phagocytosis. Cell 2009, 138, 271–285. [Google Scholar] [CrossRef] [Green Version]
  17. Imam, R.; Chang, Q.; Black, M.; Yu, C.; Cao, W. CD47 expression and CD163(+) macrophages correlated with prognosis of pancreatic neuroendocrine tumor. BMC Cancer 2021, 21, 320. [Google Scholar] [CrossRef]
  18. Yuan, J.; Shi, X.; Chen, C.; He, H.; Liu, L.; Wu, J.; Yan, H. High expression of CD47 in triple negative breast cancer is associated with epithelial-mesenchymal transition and poor prognosis. Oncol. Lett. 2019, 18, 3249–3255. [Google Scholar] [CrossRef] [Green Version]
  19. Barrera, L.; Montes-Servin, E.; Hernandez-Martinez, J.M.; Garcia-Vicente, M.L.A.; Montes-Servin, E.; Herrera-Martinez, M.; Crispin, J.C.; Borbolla-Escoboza, J.R.; Arrieta, O. CD47 overexpression is associated with decreased neutrophil apoptosis/phagocytosis and poor prognosis in non-small-cell lung cancer patients. Br. J. Cancer 2017, 117, 385–397. [Google Scholar] [CrossRef] [Green Version]
  20. Willingham, S.B.; Volkmer, J.P.; Gentles, A.J.; Sahoo, D.; Dalerba, P.; Mitra, S.S.; Wang, J.; Contreras-Trujillo, H.; Martin, R.; Cohen, J.D.; et al. The CD47-signal regulatory protein alpha (SIRPa) interaction is a therapeutic target for human solid tumors. Proc. Natl. Acad. Sci. USA 2012, 109, 6662–6667. [Google Scholar] [CrossRef] [Green Version]
  21. Golubovskaya, V.; Berahovich, R.; Zhou, H.; Xu, S.; Harto, H.; Li, L.; Chao, C.C.; Mao, M.M.; Wu, L. CD47-CAR-T Cells Effectively Kill Target Cancer Cells and Block Pancreatic Tumor Growth. Cancers 2017, 9, 139. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Candas-Green, D.; Xie, B.; Huang, J.; Fan, M.; Wang, A.; Menaa, C.; Zhang, Y.; Zhang, L.; Jing, D.; Azghadi, S.; et al. Dual blockade of CD47 and HER2 eliminates radioresistant breast cancer cells. Nat. Commun. 2020, 11, 4591. [Google Scholar] [CrossRef] [PubMed]
  23. Pai, S.; Bamodu, O.A.; Lin, Y.K.; Lin, C.S.; Chu, P.Y.; Chien, M.H.; Wang, L.S.; Hsiao, M.; Yeh, C.T.; Tsai, J.T. CD47-SIRPalpha Signaling Induces Epithelial-Mesenchymal Transition and Cancer Stemness and Links to a Poor Prognosis in Patients with Oral Squamous Cell Carcinoma. Cells 2019, 8, 1658. [Google Scholar] [CrossRef] [Green Version]
  24. Wong, D.Y.; Chang, K.W.; Chen, C.F.; Chang, R.C. Characterization of two new cell lines derived from oral cavity human squamous cell carcinomas--OC1 and OC2. J. Oral. Maxillofac. Surg. 1990, 48, 385–390. [Google Scholar] [CrossRef]
  25. Yang, C.Y.; Meng, C.L. Regulation of PG synthase by EGF and PDGF in human oral, breast, stomach, and fibrosarcoma cancer cell lines. J. Dent. Res. 1994, 73, 1407–1415. [Google Scholar] [CrossRef] [PubMed]
  26. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjostedt, E.; Asplund, A.; et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  27. Singh, P.; Rai, A.; Verma, A.K.; Alsahli, M.A.; Rahmani, A.H.; Almatroodi, S.A.; Alrumaihi, F.; Dev, K.; Sinha, A.; Sankhwar, S.; et al. Survival-Based Biomarker Module Identification Associated with Oral Squamous Cell Carcinoma (OSCC). Biology 2021, 10, 760. [Google Scholar] [CrossRef]
  28. Sakakura, K.; Takahashi, H.; Kaira, K.; Toyoda, M.; Murata, T.; Ohnishi, H.; Oyama, T.; Chikamatsu, K. Relationship between tumor-associated macrophage subsets and CD47 expression in squamous cell carcinoma of the head and neck in the tumor microenvironment. Lab. Invest. 2016, 96, 994–1003. [Google Scholar] [CrossRef]
  29. Vaeteewoottacharn, K.; Kariya, R.; Pothipan, P.; Fujikawa, S.; Pairojkul, C.; Waraasawapati, S.; Kuwahara, K.; Wongkham, C.; Wongkham, S.; Okada, S. Attenuation of CD47-SIRPalpha Signal in Cholangiocarcinoma Potentiates Tumor-Associated Macrophage-Mediated Phagocytosis and Suppresses Intrahepatic Metastasis. Transl. Oncol. 2019, 12, 217–225. [Google Scholar] [CrossRef]
  30. Lee, A.S. Glucose-regulated proteins in cancer: Molecular mechanisms and therapeutic potential. Nat. Rev. Cancer 2014, 14, 263–276. [Google Scholar] [CrossRef] [Green Version]
  31. Chang, Y.W.; Tseng, C.F.; Wang, M.Y.; Chang, W.C.; Lee, C.C.; Chen, L.T.; Hung, M.C.; Su, J.L. Deacetylation of HSPA5 by HDAC6 leads to GP78-mediated HSPA5 ubiquitination at K447 and suppresses metastasis of breast cancer. Oncogene 2016, 35, 1517–1528. [Google Scholar] [CrossRef]
  32. Luo, X.; Yao, J.; Nie, P.; Yang, Z.; Feng, H.; Chen, P.; Shi, X.; Zou, Z. FOXM1 promotes invasion and migration of colorectal cancer cells partially dependent on HSPA5 transactivation. Oncotarget 2016, 7, 26480–26495. [Google Scholar] [CrossRef]
  33. Lee, A.S. GRP78 induction in cancer: Therapeutic and prognostic implications. Cancer Res. 2007, 67, 3496–3499. [Google Scholar] [CrossRef] [Green Version]
  34. Luo, S.; Mao, C.; Lee, B.; Lee, A.S. GRP78/BiP is required for cell proliferation and protecting the inner cell mass from apoptosis during early mouse embryonic development. Mol. Cell. Biol. 2006, 26, 5688–5697. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, J.; Jiang, Y.; Jia, Z.; Li, Q.; Gong, W.; Wang, L.; Wei, D.; Yao, J.; Fang, S.; Xie, K. Association of elevated GRP78 expression with increased lymph node metastasis and poor prognosis in patients with gastric cancer. Clin. Exp. Metastasis 2006, 23, 401–410. [Google Scholar] [CrossRef]
  36. Stojadinovic, A.; Hooke, J.A.; Shriver, C.D.; Nissan, A.; Kovatich, A.J.; Kao, T.C.; Ponniah, S.; Peoples, G.E.; Moroni, M. HYOU1/Orp150 expression in breast cancer. Med. Sci. Monit. 2007, 13, BR231–BR239. [Google Scholar]
  37. Slaby, O.; Sobkova, K.; Svoboda, M.; Garajova, I.; Fabian, P.; Hrstka, R.; Nenutil, R.; Sachlova, M.; Kocakova, I.; Michalek, J.; et al. Significant overexpression of Hsp110 gene during colorectal cancer progression. Oncol. Rep. 2009, 21, 1235–1241. [Google Scholar] [CrossRef] [Green Version]
  38. Kusaczuk, M.; Cechowska-Pasko, M. Molecular chaperone ORP150 in ER stress-related diseases. Curr. Pharm. Des. 2013, 19, 2807–2818. [Google Scholar] [CrossRef]
  39. Fagone, P.; Di Rosa, M.; Palumbo, M.; De Gregorio, C.; Nicoletti, F.; Malaguarnera, L. Modulation of heat shock proteins during macrophage differentiation. Inflamm. Res. 2012, 61, 1131–1139. [Google Scholar] [CrossRef]
  40. Zhang, M.; Hutter, G.; Kahn, S.A.; Azad, T.D.; Gholamin, S.; Xu, C.Y.; Liu, J.; Achrol, A.S.; Richard, C.; Sommerkamp, P.; et al. Anti-CD47 Treatment Stimulates Phagocytosis of Glioblastoma by M1 and M2 Polarized Macrophages and Promotes M1 Polarized Macrophages In Vivo. PLoS ONE 2016, 11, e0153550. [Google Scholar] [CrossRef]
  41. Zhang, Y.; Sime, W.; Juhas, M.; Sjolander, A. Crosstalk between colon cancer cells and macrophages via inflammatory mediators and CD47 promotes tumour cell migration. Eur. J. Cancer 2013, 49, 3320–3334. [Google Scholar] [CrossRef]
  42. Liu, R.; Wei, H.; Gao, P.; Yu, H.; Wang, K.; Fu, Z.; Ju, B.; Zhao, M.; Dong, S.; Li, Z.; et al. CD47 promotes ovarian cancer progression by inhibiting macrophage phagocytosis. Oncotarget 2017, 8, 39021–39032. [Google Scholar] [CrossRef] [Green Version]
  43. Rastgoo, N.; Wu, J.; Liu, A.; Pourabdollah, M.; Atenafu, E.G.; Reece, D.; Chen, W.; Chang, H. Targeting CD47/TNFAIP8 by miR-155 overcomes drug resistance and inhibits tumor growth through induction of phagocytosis and apoptosis in multiple myeloma. Haematologica 2020, 105, 2813–2823. [Google Scholar] [CrossRef] [Green Version]
  44. Xu, Y.; Li, J.; Tong, B.; Chen, M.; Liu, X.; Zhong, W.; Zhao, J.; Wang, M. Positive tumour CD47 expression is an independent prognostic factor for recurrence in resected non-small cell lung cancer. ESMO Open 2020, 5, e000823. [Google Scholar] [CrossRef] [PubMed]
  45. Soto-Pantoja, D.R.; Kaur, S.; Roberts, D.D. CD47 signaling pathways controlling cellular differentiation and responses to stress. Crit. Rev. Biochem. Mol. Biol. 2015, 50, 212–230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Cook, K.L.; Soto-Pantoja, D.R. “UPRegulation” of CD47 by the endoplasmic reticulum stress pathway controls anti-tumor immune responses. Biomark Res. 2017, 5, 26. [Google Scholar] [CrossRef] [PubMed]
  47. Wang, Z.; Zhang, H.; Cheng, Q. PDIA4: The basic characteristics, functions and its potential connection with cancer. Biomed. Pharmacother. 2020, 122, 109688. [Google Scholar] [CrossRef]
  48. Chanjiao, Y.; Chunyan, C.; Xiaoxin, Q.; Youjian, H. MicroRNA-378a-3p contributes to ovarian cancer progression through downregulating PDIA4. Immun. Inflamm. Dis. 2021, 9, 108–119. [Google Scholar] [CrossRef]
  49. Gagnon, E.; Duclos, S.; Rondeau, C.; Chevet, E.; Cameron, P.H.; Steele-Mortimer, O.; Paiement, J.; Bergeron, J.J.; Desjardins, M. Endoplasmic reticulum-mediated phagocytosis is a mechanism of entry into macrophages. Cell 2002, 110, 119–131. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Investigation of CD47 expression in OECM-1 and OC-2 cells. (A) The flow chart shows the method of studying various gene expression changes with overexpressed CD47 in OECM-1 and OC-2 cells, through a gene expression array. The clusterProfiler tool was used to analyze bioinformatic approaches including gene ontology (GO) and DEGs. (B) Overexpression of CD47 was performed by transfecting cells with CD47 plasmid (1 μg), and CD47 expression was detected with qPCR and compared to the empty control. The data are presented as the mean ± standard deviation (SD) of three experiments. The data are presented as the mean ± SD of three experiments. * p < 0.05 vs. empty control; two-tailed Student’s t-test.
Figure 1. Investigation of CD47 expression in OECM-1 and OC-2 cells. (A) The flow chart shows the method of studying various gene expression changes with overexpressed CD47 in OECM-1 and OC-2 cells, through a gene expression array. The clusterProfiler tool was used to analyze bioinformatic approaches including gene ontology (GO) and DEGs. (B) Overexpression of CD47 was performed by transfecting cells with CD47 plasmid (1 μg), and CD47 expression was detected with qPCR and compared to the empty control. The data are presented as the mean ± standard deviation (SD) of three experiments. The data are presented as the mean ± SD of three experiments. * p < 0.05 vs. empty control; two-tailed Student’s t-test.
Cimb 44 00152 g001
Figure 2. Fold changes in mRNA expression in CD47-overexpressing oral cancer cells. (A) A volcano plot showing the differentially upregulated (right panel) and downregulated (left panel) genes in CD47-overexpressing OECM-1 and OC-2 cell lines. (B) Heatmap visualization of the RNA sequencing analysis results showing differentially expressed genes with-log10 (p value) and ≥5-fold and ≤ 0.2-fold changes. (C) Venn diagram of genes shows that 617 genes were upregulated and 653 genes were downregulated due to CD47-overexpression in oral cancer.
Figure 2. Fold changes in mRNA expression in CD47-overexpressing oral cancer cells. (A) A volcano plot showing the differentially upregulated (right panel) and downregulated (left panel) genes in CD47-overexpressing OECM-1 and OC-2 cell lines. (B) Heatmap visualization of the RNA sequencing analysis results showing differentially expressed genes with-log10 (p value) and ≥5-fold and ≤ 0.2-fold changes. (C) Venn diagram of genes shows that 617 genes were upregulated and 653 genes were downregulated due to CD47-overexpression in oral cancer.
Cimb 44 00152 g002
Figure 3. The function of differentially expressed genes in CD47-overexpressing oral cancer cells. (A) The number of DEGs in each comparison was determined by RSEM and EBSeq (PPEE ≤ 0.05). Fourteen genes were found to be involved in molecular functions, biological processes, and cellular components. (B) Heatmap demonstrating the expression profiles of union sets. The color bar on the right side indicates the log2 fold changes. (C) The list presents the 14 associated genes, including the DEGs in OECM-1 and OC-2 cells and their functions in biological processes. NMD: nonsense-mediated mRNA decay; Readthrough: continuous transcription; NA: not applicable.
Figure 3. The function of differentially expressed genes in CD47-overexpressing oral cancer cells. (A) The number of DEGs in each comparison was determined by RSEM and EBSeq (PPEE ≤ 0.05). Fourteen genes were found to be involved in molecular functions, biological processes, and cellular components. (B) Heatmap demonstrating the expression profiles of union sets. The color bar on the right side indicates the log2 fold changes. (C) The list presents the 14 associated genes, including the DEGs in OECM-1 and OC-2 cells and their functions in biological processes. NMD: nonsense-mediated mRNA decay; Readthrough: continuous transcription; NA: not applicable.
Cimb 44 00152 g003
Figure 4. The biofunction pathway analysis of the differentially expressed genes regulated by CD47-overexpression in oral cancer cells. (A) The network analysis report of upregulated genes by CD47-overexpression in OECM-1 and OC-2 cells. Molecular focus gene was involved in top diseases and function. (B) Three identified genes, HSPA5, HYOU1, and PDIA4, were found to be involved in molecular functions, biological processes, and cellular components in OSCC. (C,D) CD47, HSPA5, HYOU1, and PDIA4 gene expression and survival was analyzed using the HPA database. Cut off FPKM of 25%, 50%, and 75% were used to find the survival rate, respectively. R; correlation coefficient; two-tailed Student’s t-test. (E,F) The expression of HSPA5 was examined in CD47-overexpressing OECM-1 and OC-2 cells using qPCR. (G,J) HSPA5 and CD47 gene expression was analyzed by qPCR and (H,I) cell growth was analyzed by CCK-8 assay. The data are presented as the mean ± SD of three experiments. * p < 0.05 vs. empty control; two-tailed Student’s t-test.
Figure 4. The biofunction pathway analysis of the differentially expressed genes regulated by CD47-overexpression in oral cancer cells. (A) The network analysis report of upregulated genes by CD47-overexpression in OECM-1 and OC-2 cells. Molecular focus gene was involved in top diseases and function. (B) Three identified genes, HSPA5, HYOU1, and PDIA4, were found to be involved in molecular functions, biological processes, and cellular components in OSCC. (C,D) CD47, HSPA5, HYOU1, and PDIA4 gene expression and survival was analyzed using the HPA database. Cut off FPKM of 25%, 50%, and 75% were used to find the survival rate, respectively. R; correlation coefficient; two-tailed Student’s t-test. (E,F) The expression of HSPA5 was examined in CD47-overexpressing OECM-1 and OC-2 cells using qPCR. (G,J) HSPA5 and CD47 gene expression was analyzed by qPCR and (H,I) cell growth was analyzed by CCK-8 assay. The data are presented as the mean ± SD of three experiments. * p < 0.05 vs. empty control; two-tailed Student’s t-test.
Cimb 44 00152 g004
Figure 5. RNA sequencing showing the biofunction of DEGs in CD47-overexpression in the OECM-1 and OC-2 cell lines according to DAVID biological process analysis. The numbers represent the number of genes involved in the biological process.
Figure 5. RNA sequencing showing the biofunction of DEGs in CD47-overexpression in the OECM-1 and OC-2 cell lines according to DAVID biological process analysis. The numbers represent the number of genes involved in the biological process.
Cimb 44 00152 g005
Table 1. Genes with significant changes induced by the overexpression of CD47 in OECM-1 and OC-2 cell lines.
Table 1. Genes with significant changes induced by the overexpression of CD47 in OECM-1 and OC-2 cell lines.
Official Gene SymbolGene NameOECM-1-CD47/OECM-1-Empty Control
FOLD-Change
OC2-CD47/OC2-Empty Control
Fold-Change
GENE Expression
BAALC-AS1BAALC antisense RNA 112.156.26Up
CCDC184Coiled–coil domain containing 1845.468.88Up
CRB1Crumbs 1, cell polarity complex component5.466.26Up
HPCAHippocalcin7.698.88Up
LINC01099Long intergenic non-protein coding RNA 10997.736.26Up
LINC01512Long intergenic non-protein coding RNA 15125.466.26Up
LOC102723703Uncharacterized LOC1027237037.558.88Up
NMBRNeuromedin B receptor5.466.26Up
PCP2Purkinje cell protein 25.056.26Up
RASGRP1RAS guanyl releasing protein 15.468.88Up
SMC2-AS1SMC2 antisense RNA 1 (head to head)9.926.26Up
SMN2Survival of motor neuron 2, centromeric190.2522.01Up
TRIM6-TRIM34TRIM6-TRIM34 readthrough117.7848.70Up
ZC3H12BZinc finger CCCH-type containing 12B5.486.26Up
ZNF32-AS2ZNF32 antisense RNA 27.696.26Up
ZNF664-FAM101ANA10.038.99Up
C16orf71Chromosome 16 open reading frame 710.190.12Down
CELF5CUGBP Elav-like family member 50.140.16Down
FAM53B-AS1FAM53B antisense RNA 10.190.12Down
GAL3ST2Galactose-3-O-sulfotransferase 20.140.16Down
KATNBL1P6Katanin regulatory subunit B1-like 1 pseudogene 60.110.11Down
LOC100506136NA0.140.05Down
LOC101926897NA0.200.16Down
LOC101927178Uncharacterized LOC1019271780.150.12Down
LOC101927666Uncharacterized LOC1019276660.200.09Down
MOGMyelin oligodendrocyte glycoprotein0.110.16Down
NRN1LNeuritin 1 like0.190.12Down
PROSER2-AS1PROSER2 antisense RNA 10.140.12Down
RALGAPA1PRal GTPase activating protein catalytic alpha subunit 1 pseudogene 10.010.02Down
RGPD6RANBP2-like and GRIP domain containing 60.120.20Down
RIMS4Regulating synaptic membrane exocytosis 40.190.16Down
THRILTNF- and HNRNPL-related immunoregulatory long non-coding RNA0.190.16Down
Head-to-head: a genomic locus where two adjacent genes are differentially transcribed from opposite DNA strands. Readthrough: continuous transcription. NA: not applicable.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tseng, C.-C.; Tsou, C.-H.; Huang, S.-Y.; Wu, C.-W.; Hsieh, T.-H. Using Next-Generation Sequencing and Bioinformatic Methods to Predict New Genes That May Be Regulated by CD47 in Oral Squamous Cell Carcinoma. Curr. Issues Mol. Biol. 2022, 44, 2243-2256. https://doi.org/10.3390/cimb44050152

AMA Style

Tseng C-C, Tsou C-H, Huang S-Y, Wu C-W, Hsieh T-H. Using Next-Generation Sequencing and Bioinformatic Methods to Predict New Genes That May Be Regulated by CD47 in Oral Squamous Cell Carcinoma. Current Issues in Molecular Biology. 2022; 44(5):2243-2256. https://doi.org/10.3390/cimb44050152

Chicago/Turabian Style

Tseng, Chung-Chih, Chen-Han Tsou, Shi-Ying Huang, Chia-Wei Wu, and Tsung-Hua Hsieh. 2022. "Using Next-Generation Sequencing and Bioinformatic Methods to Predict New Genes That May Be Regulated by CD47 in Oral Squamous Cell Carcinoma" Current Issues in Molecular Biology 44, no. 5: 2243-2256. https://doi.org/10.3390/cimb44050152

APA Style

Tseng, C. -C., Tsou, C. -H., Huang, S. -Y., Wu, C. -W., & Hsieh, T. -H. (2022). Using Next-Generation Sequencing and Bioinformatic Methods to Predict New Genes That May Be Regulated by CD47 in Oral Squamous Cell Carcinoma. Current Issues in Molecular Biology, 44(5), 2243-2256. https://doi.org/10.3390/cimb44050152

Article Metrics

Back to TopTop