- 1Sports Medicine Center, Southwest Hospital, Third Military Medical University, Chongqing, China
- 2Department of Stomatology, The 970th Hospital of the Joint Logistics Support Force, Yantai, China
- 3Department of Stomatology, Southwest Hospital, Third Military Medical University, Chongqing, China
Oral squamous cell carcinoma (OSCC), a kind of malignant cancer, is associated with increasing morbidity and mortality. Patients with different genetic ancestries may respond differently to clinical treatment. The limited understanding of the influence of genetic ancestry and genetic characteristics on OSCC impedes the development of precision medicine. To provide a reference for clinical treatment, this study comprehensively analyzed multigenomic differences in OSCC patients with different genetic ancestries and their impact on prognosis. An analysis of data from OSCC patients with different genetic ancestries in The Cancer Genome Atlas (TCGA) showed that the overall survival (OS) of African (AFR) patients was lower than that of primarily European (EUR) patients, and differences were also observed in the tumor–stroma ratio (TSR) and tumor-infiltrating lymphocytes (TILs), which are associated with prognosis. FAT1 is a key mutant gene in OSCC, and it has inconsistent effects on clinical evolution for patients with diverse genetic characteristics. PIKfyve and CAPN9 showed a significant difference in mutation frequency between EUR and AFR; PIKfyve was related to Ki-67 expression, suggesting that it could promote tumor proliferation, and CAPN9 was related to the expression of Bcl-2, promoting tumor cell apoptosis. A variant methylation locus, cg20469139, was correlated with the levels of PD-L1 and Caspase-7 and modulated tumor cell apoptosis. A novel ceRNA model was constructed based on genetic ancestries, and it could accurately evaluate patient prognosis. More importantly, although T cell dysfunction scores could determine the potential of tumor immune escape, the efficacy was obviously affected by patients’ genetic ancestries. To provide patients with more precise, personalized therapy and to further improve their quality of life and 5-year survival rate, the influence of genetic ancestry should be fully considered when selecting treatments.
Introduction
Oral squamous cell carcinoma (OSCC) is one of the most common malignant tumors in the oral cavity or head and neck region and is listed as the eighth most common cancer type in the world (Carnielli et al., 2018; Wang et al., 2021), accounting for approximately 90% of oral tumors and 3% of systemic malignant tumors (Zhao et al., 2020; Lu et al., 2021). According to the latest Global Cancer Incidence, the GLOBOCAN 2020 estimates of cancer incidence and mortality produced by the International Agency for Research on Cancer (IARC), newly developed OSCC originating from the alveolar ridge, buccal mucosa, bottom of the oral cavity, upper jaw, tongue, and other parts of the mouth accounted for 377,713 cases worldwide, with a death toll of 177,757 (Sung et al., 2021), which is much higher than the 300,000 new cases and 145,000 deaths reported globally in 2012 (Ferlay et al., 2015). The morbidity and mortality of OSCC are increasing year by year, and the 5-year survival rate of OSCC continues to be lower than 60% (Ho et al., 2019; Bosetti et al., 2020; Ferreira et al., 2021). Even for patients who received standard therapy with proper combinations of surgical and nonsurgical treatments, the recurrence rate of OSCC was still relatively high in the range of 18–76% (da Silva et al., 2012; Wang et al., 2020). Therefore, identifying new treatment methods and more comprehensive prognostic indicators is particularly important.
People of different genetic ancestries often show different reproductive genetic characteristics (Price et al., 2006; Huynh-Le et al., 2021; Zhang et al., 2021). Due to exposure to different living environments or different pathogenic or nonpathogenic factors, their tumor morbidities, outcomes, prognoses, and pathogenic molecular characteristics differ (Yang et al., 2011; Vineis and Fecht, 2018; Ou et al., 2020). Previous studies have mainly focused on molecular differences between normal oral tissues and tumor tissues and their effects (Chen et al., 2020; Yin et al., 2020). No report is available on differences in molecular-genomic characteristics of OSCC patients of different genetic ancestries. In this study, The Cancer Genome Atlas (TCGA) was used as the data source. We intend to use a comprehensive cross-platform and multigenomic analysis, including hematoxylin–eosin (HE) staining, somatic mutations, methylation, RNA expression, immune infiltration, and immune response data, to analyze the differences in the prognoses of OSCC patients of various genetic ancestries, especially primarily European (EUR) and African (AFR). We sought to exploit these results to improve our understanding of the molecular and cellular effects of ancestry across clinical evolution and the relationship between ancestry and clinical treatment and prognosis to provide new options for the effective precision treatment of OSCC patients (Figure 1).
Materials and Methods
Data Collection
The hematoxylin–eosin (HE) staining results, single-nucleotide variant data, methylation data, RNA sequencing data, and clinical data of OSCC patients were downloaded from the TCGA database (https://cancergenome.nih.gov/). The genetic ancestry–related data of OSCC patients were obtained from Carrot-Zhang et al. (2020). In addition, the reverse-phase protein array data of patients were downloaded from The Cancer Proteome Atlas (TCPA) database (https://www.tcpaportal.org/tcpa/) (Li et al., 2013). T cell dysfunction scores of OSCC patients were downloaded from the Tumor Immune Dysfunction and Exclusion (TIDE) database (http://tide.dfci.harvard.edu) (Fu et al., 2020).
Pathological and Survival Analyses of Patients With Different Genetic Ancestries
Based on HE staining results, we evaluated the pathological differentiation degree (well, moderately, poorly) of OSCC patients and the tumor–stroma ratio (TSR) according to the method reported in the literature, which defines a TSR ≥50% as high and a TSR <50% as low (van Pelt et al., 2018; Hagenaars et al., 2021). Tumor-infiltrating lymphocytes (TILs) were also evaluated. TILs ≥10% was considered high, and TILs <10% was considered low (Salgado et al., 2015; Veatch et al., 2021). Differences in the TSR and TILs between OSCC patients of different genetic ancestries were analyzed. In a Kaplan-Meier survival analysis on primarily European (EUR) and African (AFR) OSCC patients, to eliminate the effects of old age and short follow-up times, OSCC patients older than 80 years and those with a follow-up time less than 90 days were excluded.
Analysis of Somatic Mutations and Methylation
We analyzed differences in single-nucleotide variant frequencies between OSCC patients with different genetic ancestries in the TCGA database. Kaplan-Meier survival analysis was performed to compare the effects of mutations between EUR and AFR OSCC patients. Then, we used the cBioPortal database (http://www.cbioportal.org/) (Gao et al., 2013) to perform and visualize a pancancer statistical analysis of FAT1 mutations.
The difference in locus methylation between EUR and AFR OSCC patients was analyzed, where | log2FC | > 1 and a false discovery rate (FDR) < 0.05 were the thresholds for differential methylation. We performed a principal component analysis (PCA) on the differentially methylated sites to determine whether they could effectively differentiate between EUR and AFR OSCC patients. At the same time, survival-related methylation sites were identified, and a correlation analysis between these sites and protein levels was conducted to explore the functions of survival-related methylation sites.
Construction of a ceRNA Risk Regression Model
We, respectively constructed mRNA, lncRNA, and miRNA risk regression models of OSCC patients. Kaplan-Meier survival analysis, multiclinical information heat maps, and receiver operator characteristic (ROC) curves were used to evaluate the effects of the three models. We integrated the three models, yielding the ceRNA risk regression model, and compared it with the three component models.
An independent prognostic analysis was conducted to assess the effects of genetic ancestry and the risk scores from the ceRNA models on the prognosis of OSCC patients. The ceRNA risk regression model was also applied to the pancancer analysis to determine whether it had specificity for OSCC. Finally, we constructed the ceRNA network of mRNA-miRNA-lncRNA to further investigate the relationship between RNAs. We also used the STRING database (http://string-db.org/cgi/input.pl) (Szklarczyk et al., 2019) to construct the protein interaction network for genetic ancestry.
Analysis of Immune Cell Infiltration
We used the CIBERSORTx algorithm (Newman et al., 2019) to calculate the infiltration of 22 kinds of immune cells in OSCC patients and analyzed differences in the immune cell distribution between patients of different genetic ancestries. We compared the effect of using T cell dysfunction scores to evaluate prognoses between OSCC patient groups of different genetic ancestries.
Statistical Analysis
Bilateral tests were performed for all statistical tests. A p value less than 0.05 was considered statistically significant. R software version 4.0.0 (https://www.r-project.org/) was used for analyses. Some R packages were used in this study, including “limma,” “RColorBrewe,” “ggplot2,” “survival,” “maftools,” “corrplot,” “pheatmap,” “survivalROC,” “reshape2,” and “forestplot.” The ceRNA network was visualized using Cytoscape_3.7.2 (https://cytoscape.org) (Shannon, 2003).
Results
OSCC Patients and Their Genetic Ancestry Distributions
We analyzed the genetic ancestries of OSCC patients in the TCGA database. Among 434 patients, primarily European ancestry (EUR; n = 370) accounted for the highest proportion at 85.25%, African ancestry (AFR; n = 32) accounted for 7.37%, Native/Latin American ancestry (AMR; n = 6) accounted for 1.38%, East Asian ancestry (EAS; n = 6) accounted for 1.38%, South Asian ancestry (SAS; n = 3) accounted for 0.69%, and admixed descent ancestry (ADMIX; n = 17) accounted for 3.92% of the sample (Figure 2A). As more than two-thirds of TCGA donors were from the United States (Carrot-Zhang et al., 2020), the numbers of OSCC patients with AMR, EAS, SAS, and ADMIX genetic ancestries in the database were relatively small and may not accurately represent the corresponding overall OSCC population. Therefore, we only discuss differences between the EUR and AFR.
FIGURE 2. Characteristics of OSCC patients of EUR or AFR genetic ancestry. (A) The percentage of OSCC patients with different genetic ancestries in the TCGA database. (B) The OS of patients of AFR ancestry was significantly lower than that of patients of EUR ancestry. (C) The complete remission rates of patients of EUR and AFR ancestries following initial treatment were significantly different. (D) The results of HE staining were used to divide EUR and AFR ancestry patients into TSR ≥50% and TSR <50% groups and TILs ≥10% and TILs <10% groups (scale bar: 50 μM, magnification: 40X; the solid arrowhead indicates tumor cells, and the dotted arrowhead indicates lymphocytes). (E) HE staining results were used to divide EUR and AFR ancestry patients into well differentiated, moderately differentiated, and poorly differentiated tumor groups (scale bar: 50 μM, magnification: 40X; the solid arrowhead indicates tumor cells).
First, the Kaplan-Meier survival analysis indicated that the overall survival (OS) of AFR patients was significantly lower than that of EUR patients (p = 0.026) (Figure 2B). When we compared the initial treatment responses of EUR and AFR OSCC patients, the percentage of complete remission/complete response (CR) in EUR was significantly higher than that in AFR. The CR rates of EUR and AFR OSCC patients were 86.2 and 65.2%, respectively (p = 0.014) (Figure 2C). Analysis of the HE staining results showed that the TSR and TILs in EUR and AFR OSCC patients were significantly different (Figure 2D; Table 1) (p < 0.05), but the pathological differentiation was not significantly different between the two groups of patients (Figure 2E; Table 1).
TABLE 1. Differential analysis of the differentiation degree, the TSR, and TILs in OSCC patients of EUR and AFR ancestries (53 cases without corresponding data among patients of EUR ancestry and 1 case without corresponding data among patients of AFR ancestry).
Somatic Mutation Differences Between OSCC Patients With Different Genetic Ancestries
We analyzed somatic mutations in OSCC patients in the TCGA database. The mutation frequency of FAT1 in OSCC patients was 21%, rendering it the third-most often mutated gene in this patient set (Figure 3A). We performed a pancancer analysis of FAT1 mutations. As shown in Figure 3B, FAT1 had the highest mutation frequency in head and neck squamous cell carcinoma, suggesting that FAT1 mutation might play an important role in the occurrence and development of OSCC (Nishikawa et al., 2011; Hayes et al., 2016; Hsu et al., 2019; Kuo et al., 2019). We analyzed the effect of FAT1 mutation on the survival of patients of EUR and AFR. The results showed that among patients of AFR, the OS of FAT1-mutated patients was significantly better than that of wild-type patients (Figure 3C), but this difference was not significant in EUR (Figure 3D). We next compared genes with the top 10 mutation frequencies in AFR and EUR (Figures 3E–G). UNC13C only appeared in the top 10 genes of AFR and showed a significant difference in mutation frequency between the EUR and AFR (p < 0.01) (Figure 3H). Considering the impact of clinical stage, we also obtained the same results in N0 and N1-3 patients, respectively (Figure 3I). Among the other differentially mutated genes, the mutation frequency of PIKfyve in AFR was significantly higher than that in EUR (AFR = 9.38%, EUR = 1.27%, p < 0.05). PIKfyve mutations mainly lead to an increased expression of the same gene (Figure 3J). Hou et al. (2018) showed that PIKfyve can promote tumor cell proliferation. Therefore, we analyzed the level of Ki-67, a cell proliferation–related marker, and found that the Ki-67 level was significantly higher in patients with PIKfyve mutations (Figure 3K), suggesting that PIKfyve mutation promoted tumor proliferation. The mutation frequency of CAPN9 in AFR was also significantly higher than that in EUR (AFR = 6.25%, EUR = 0.75%, p < 0.05), and CAPN9 mutations were associated with higher expression of the same gene (Figure 3L). Peng et al. (2016) found that overexpression of CAPN9 promoted tumor cell apoptosis. We analyzed the level of Bcl-2, an apoptosis-related marker, and found that Bcl-2 was significantly elevated in patients with CAPN9 mutations (Figure 3M), suggesting that CAPN9 mutations promote tumor cell apoptosis by increasing this gene’s expression, thus playing a role in cancer inhibition.
FIGURE 3. Somatic mutation differences between EUR and AFR OSCC patients. (A) Waterfall plots of the genes with the top 10 single-nucleotide variant frequencies among all OSCC patients in the TCGA database. (B) The pancancer mutation frequency of FAT1, which shows that the mutation frequency was highest in head and neck squamous cell carcinoma. (C) The survival of the FAT1-mutated patients was significantly better than that of the wild-type patients in the AFR ancestry. (D) No significant difference in survival was found between FAT1-mutated and wild-type patients in the EUR ancestry. (E) The genes with the top 10 mutation frequencies in the AFR ancestry. (F) The genes with the top 10 mutation frequencies in the EUR ancestry. (G) Differences in the genes with the highest mutation frequencies between EUR and AFR ancestries. (H) The mutation frequency of UNC13C was significantly different between the EUR and AFR ancestries. (I) The mutation frequency of UNC13C in N0 and in N1-3 stages was significantly different between the EUR and AFR ancestries. (J) PIKfyve mutation caused increased expression of this gene. (K) OSCC patients with PIKfyve mutation had higher levels of Ki-67 expression. (L) CAPN9 mutation caused increased expression of this gene. (M) OSCC patients with CAPN9 mutation had higher Bcl-2 expression.
DNA Methylation Differences Between OSCC Patients With Different Genetic Ancestries
When the methylation information of OSCC patients in the TCGA database was analyzed, the results showed that 75 methylated sites significantly differed between EUR and AFR (Figure 4A). PCA showed that these sites could effectively differentiate patients of EUR and AFR (Figure 4B). Univariate Cox analysis showed that five methylation sites were associated with the prognosis of OSCC patients (Figure 4C). Correlation analysis between the methylation data of these sites and the protein expression data in the TCPA database showed that the degree of methylation at site cg20469139 was correlated with the levels of two apoptosis-related proteins, PDL-1 (Figure 4D) and caspase-7 (Figure 4E), suggesting that this site may affect the prognosis of patients by regulating tumor cell apoptosis.
FIGURE 4. Methylation differences between EUR and AFR OSCC patients. (A) Heat map of differentially methylated sites between the EUR and AFR ancestries. (B) Differentially methylated sites can effectively differentiate patients of EUR and AFR ancestries. (C) Five differentially methylated sites were associated with the prognosis of OSCC patients. (D) The degree of methylation of the cg20469139 site was significantly correlated with the PDL-1 level. (E) The degree of methylation at the cg20469139 site was significantly correlated with the caspase-7 level. (F) Correlations between the methylation levels of the five survival-related loci.
We further performed a correlation analysis of the five methylation sites, and the results showed that the correlation between these methylation sites was relatively low, indicating that the information carried by these sites had almost no redundancy or crossover; therefore, these sites might affect tumor progression in different manners (Figure 4F). In addition, we found that cg10218876 was located in the TNNT1 gene. TNNT1 promotes the proliferation of breast cancer cells by promoting the G1/S transition (Shi et al., 2018). cg20469139 is located in the RNF135 gene. The RNF135 protein has been reported to inhibit the proliferation of tongue squamous carcinoma cell line SCC25 and participate in PTEN signaling transduction (Jin et al., 2016). cg15864601 is located in the C17orf97 gene. Coexpression of CK20 (also known as C17orf97) and Ki-67 may play an important role in the progression of bladder cancer and can be used as a prognostic indicator (Ye et al., 2010). cg05937630 is located in the AKR7A3 gene. AKR7A3 inhibits the tumorigenicity and chemoresistance of hepatic carcinoma by attenuating ERK, c-Jun, and NF-κB signaling pathways (Chow et al., 2017). cg18912855 is located in the PACS2 gene. PACS2 can promote ErbB signal transduction by regulating the recovery of the metalloproteinase ADAM17, thus affecting the biological behavior of tumors (Dombernowsky et al., 2015).
ceRNA Expression Differences Between OSCC Patients With Different Genetic Ancestries
We performed differential expression analyses of mRNA, lncRNA, and miRNA between OSCC patients of different genetic ancestries and constructed respective risk regression models based on differentially expressed mRNA, differentially expressed lncRNA, and differentially expressed miRNA. From the risk heat map, we found that as the risk value increased, the proportion of deaths significantly increased, and the survival time also decreased (Figure 5A). The results of the Kaplan-Meier survival analysis showed that all three risk regression models accurately assessed the prognoses of OSCC patients of all genetic ancestries, including EUR and AFR. The OS and progression-free survival (PFS) of the low-risk group were superior to those of the high-risk group (Figure 5B).
FIGURE 5. Differences in ceRNA expression between EUR and AFR OSCC patients. (A) The multiclinical trait heat map of the mRNA, lncRNA, and miRNA risk regression models (from left to right). (B) The mRNA, lncRNA, miRNA, and ceRNA risk regression models were used to predict the prognoses of OSCC patients of EUR, AFR, and all genetic ancestries. (C) ROC curves for the prognosis at 1, 3, and 5 years from the four risk regression models.
To further evaluate the prognosis of OSCC patients, we integrated the three models to establish a ceRNA risk regression model. From the ROC curve (Figure 5C), we found that the predictive abilities of the ceRNA risk regression model for the prognoses of OSCC patients at one, three, and 5 years were all better than those of the mRNA, lncRNA, and miRNA risk regression models alone, indicating a satisfied predictive effect. Based on the risk score calculated by the ceRNA model, the patients were divided into 10 groups from low to high. As shown in Figure 6A, with the risk score increased, the proportion of patients of AFR increased, and their prognoses were poor, which is consistent with the results of the Kaplan-Meier survival analysis showing that the OS of patients of AFR was significantly lower than that of patients of EUR.
FIGURE 6. Evaluation effectiveness and specificity of the ceRNA risk regression model. (A) The percentages of EUR and AFR OSCC patients in different risk stratification groups. (B) In the ceRNA risk regression model, the relapse rate of the OSCC patients in the high-risk group was significantly higher than that in the low-risk group. (C) Univariate Cox analysis of the ceRNA risk regression model and related clinical information such as genetic ancestry. (D) Multivariate Cox analysis of the ceRNA risk regression model and related clinical information such as genetic ancestry. (E) No significant difference in the prognoses of all cancer patients was found when using the ceRNA risk regression model. (F) The ceRNA network of differentially expressed RNAs between EUR and AFR ancestries. (G) The interaction network of differentially expressed proteins between the EUR and AFR ancestries.
Recurrence is an important factor affecting the survival of OSCC patients. The recurrence plot (Figure 6B) showed that the relapse rate was significantly higher in the high-risk patients than in the low-risk patients (p < 0.01), indicating that the ceRNA risk regression model can not only predict the survival of patients but also accurately predict the risk of tumor recurrence. To further evaluate the effects of different genetic ancestries on prognosis, we performed an independent prognostic analysis. Univariate Cox analysis suggested that genetic ancestry and ceRNA risk scores were correlated with prognosis (Figure 6C). However, multivariate Cox analysis showed that only risk scores could be used as an independent prognostic factor for OSCC patients (Figure 6D).
We have used this ceRNA risk regression model in pancancer prognostic evaluations and found that it could not accurately assess the prognosis of patients with other cancers (Figure 6E), indicating that this model has strong specificity for OSCC patients.
Lastly, we constructed a ceRNA and protein interaction network to analyze interactions between differentially expressed RNAs in OSCC patients with different genetic ancestries. As shown in Figures 6F,G, these mRNAs, lncRNAs, and miRNAs had extensive and close associations.
Immunity Differences Between OSCC Patients With Different Genetic Ancestries
We used gene set enrichment analysis (GSEA) to perform an enrichment analysis on the differentially expressed genes between EUR and AFR, and found that they were significantly enriched in the GO_REGULATION_OF_TYPE_2_IMMUNE_RESPONSE pathway (Figure 7A). Metascape was used to perform an enrichment analysis on the genes used to build the ceRNA risk regression model (Figure 7B), which also showed that these genes were associated with the immune response. Therefore, we speculated that OSCC patients with different genetic ancestries have differences in immune infiltration and immune response, which affected their prognoses.
FIGURE 7. Differences in immunity between EUR and AFR OSCC patients. (A) GSEA of differentially expressed genes between EUR and AFR ancestries in the GO_REGULATION_OF_TYPE_2_IMMUNE_RESPONSE pathway. (B) Enrichment analysis of genes used to build the ceRNA risk regression model using Metascape showed that these genes were associated with the immune response. (C) Histogram of the degree of infiltration of different immune cell types in EUR and AFR ancestries. (D) Violin plot of the differences in the degree of infiltration of different immune cell types between EUR and AFR ancestries. (E) Scatter plot of the differences in the degree of infiltration of different immune cell types between EUR and AFR ancestries. (F) Relationships between T cell dysfunction scores and preliminary treatment outcomes in EUR and AFR OSCC patients.
We used the CIBERSORTx algorithm to calculate the infiltration levels of 22 immune cells in EUR and AFR (Figure 7C). The violin plot (Figure 7D) and the scatter plot (Figure 7E) showed that the degree of plasma cell infiltration in EUR was significantly higher than that in AFR (p < 0.01). Gentles et al. (2015) showed that the plasma cell content is a significant predictor of the survival of patients with solid tumors, which is consistent with our finding that the prognosis of EUR OSCC patients was better than that of AFR. In EUR, the T cell dysfunction scores of patients with complete remission after initial treatment was significantly higher than that of patients with incomplete remission (p < 0.05); however, the difference between the two patient groups was not significant in the population with AFR (Figure 7F).
We further explored the relationship between the infiltration of different immune cells and survival. Naive B cells, plasma cells, and mast cells were associated with the OS of EUR but not AFR. CD8+ T cells, M1 macrophages, and M2 macrophages were associated with the PFS of AFR but not EUR, further confirming the differences in immune infiltration between the patients of the two genetic ancestries (Figures 8A,B).
FIGURE 8. Relationship between immune cells and survival with different genetic ancestries. (A) The relationship between immune cell infiltration and survival in EUR and AFR OSCC patients. (B) P-value heat map of the relationship between immune cell infiltration and survival in EUR and AFR OSCC patients.
Discussion
OSCC refers to a group of heterogeneous tumors originating from the mucosal lining (Meng et al., 2021). Most of these tumors are associated with lifestyle factors, such as smoking, excessive drinking, and betel nut chewing (Adil et al., 2021; Weiße et al., 2021). Cetuximab, which targets epidermal growth factor receptors, was approved for OSCC treatment in 2006 and remains the only available targeted molecular therapy for OSCC to date (Chai et al., 2020; Johnson et al., 2020). Therefore, exploring more effective targets and accurately assessing the prognosis are key to improving the 5-year survival rate of OSCC patients.
With the rapid development of precision medicine, a wide range of personal data, including clinical, lifestyle, genetic, and biological marker information, must be fully considered when selecting treatments (Konig et al., 2017; Elmore et al., 2021). Different genetic ancestries play important roles in the occurrence and progression of tumors. For example, Jiagge et al. (2018) found that the mortality rate of breast cancer among African-American women was significantly higher than that among white American women. Gentles et al. (2015) showed that some genomic differences may lead to higher prostate cancer morbidity and mortality in African-American men. In this study, we comprehensively analyzed the molecular characteristics of OSCC in patients between EUR and AFR. Our results showed that OSCC patients with different genetic ancestries had significant differences (p < 0.05) not only in initial treatment responses and overall survival but also in the TSR and TILs, two variables that affect prognosis.
At the molecular level, somatic mutations in genes such as FAT1 are more common in OSCC than in other tumors (Chai et al., 2020; Sequeira et al., 2020). We found that FAT1 mutation had different effects on prognosis in OSCC patients with different genetic ancestries. In addition, the mutation frequencies of UNC13C, PIKfyve, and CAPN9 were significantly different between EUR and AFR (p < 0.05) and might affect prognosis by regulating tumor proliferation and apoptosis. Five methylation sites differed between OSCC patients with different genetic ancestries and affected their prognoses (p < 0.05). Our ceRNA risk regression model revealed differences in survival and prognosis between AFR and EUR OSCC patients (p < 0.05). The degree of plasma cell infiltration in EUR was significantly higher than that in AFR (p < 0.05). The prognosis of EUR was better than that of AFR, which is consistent with the results of Gentles et al. (2015). Jiang et al. (2018) found that T cell dysfunction scores could predict the responses of patients to immunotherapy. We found differences in complete remission and incomplete remission between OSCC patients with different genetic ancestries, suggesting that the predictive ability of T cell dysfunction scores might be different in patients with different genetic ancestries and that the effects of different genetic ancestries should be considered when using T cell dysfunction scores.
The present study has some limitations. TCGA cohort consists of convenience samples, so it may not represent a general cancer population (Wang et al., 2018). Only almost 15% OSCC patients in our study were at least partially non-EUR, a robust distinction of tumor-specific ancestral associations will require comparative analyses of more tissue samples from patients of different genetic ancestries, especially non-EUR OSCC patient tissues. However, we found fairly strong evidence suggesting that genetic ancestry–based differences in somatic mutations, methylation, ceRNA expression, immune infiltration, and immune responses are associated with cancer prognosis. In particular, the differences in OSCC variables between the EUR and AFR populations suggest that genetic ancestry should be considered when contemplating disease-causing factors and related treatments. Our ceRNA risk regression model can be used to evaluate the prognosis of OSCC patients with different genetic ancestries and is specific for OSCC. These results provide new ideas for the treatment and prevention of OSCC in patients with different genetic ancestries.
This study explored the effect of genetic ancestry on the clinical evolution of OSCC patients and the molecular mechanism through a comprehensive analysis of multigenomic differences between OSCC patients with different genetic ancestries. The genetic ancestry of a patient must be considered when selecting a treatment strategy. Such considerations will be conducive to achieving precision treatment and further improving therapeutic effects and 5-year survival rate.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: https://cancergenome.nih.gov/.
Author Contributions
Conceptualization, JG and WZ; methodology, JG, and WZ; software, JG; validation, JG and XL; formal analysis, JG and XL; investigation, TL and KT; resources, YZ; data curation, YZ; writing—original draft preparation, JG and JG; writing—review and editing, JG and WZ; visualization, JG; supervision, WZ; project administration, TL; funding acquisition, KT. All authors have read and agreed to the published version of the manuscript. All of the changes are highlighted in red.
Funding
This study was supported by the Personalization Training Program for the Training Object of the Outstanding Talents of Army Medical University (2019) and Chief medical expert of Chongqing (2019).
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.
Abbreviations
OSCC, oral squamous cell carcinoma; TCGA, the cancer genome atlas; HE, hematoxylin–eosin; EUR, primarily European; AFR, African; IARC, international agency for research on cancer; TCPA, the cancer proteome atlas; TIDE, tumor immune dysfunction and exclusion; TSR, tumor–stroma ratio; TILs, tumor-infiltrating lymphocytes; PCA, principal component analysis; ROC, receiver operator characteristic; OS, overall survival; CR, complete remission/response; PFS, progression-free survival; LIHC, liver hepatocellular carcinoma; BRCA, breast invasive carcinoma; COAD, colon adenocarcinoma; OV, ovarian serous cystadenocarcinoma; KIRC, kidney renal clear cell carcinoma; PAAD, pancreatic adenocarcinoma; BLCA, bladder urothelial carcinoma; STAD, stomach adenocarcinoma.
References
Adil, N., Ali, H., Siddiqui, A. J., Ali, A., Ahmed, A., El-Seedi, H. R., et al. (2021). Evaluation of Cytotoxicity of Areca Nut and its Commercial Products on normal Human Gingival Fibroblast and Oral Squamous Cell Carcinoma Cell Lines. J. Hazard. Mater. 403, 123872. doi:10.1016/j.jhazmat.2020.123872
Bosetti, C., Carioli, G., Santucci, C., Bertuccio, P., Gallus, S., Garavello, W., et al. (2020). Global Trends in Oral and Pharyngeal Cancer Incidence and Mortality. Int. J. Cancer 147 (4), 1040–1049. doi:10.1002/ijc.32871
Carnielli, C. M., Macedo, C. C. S., De Rossi, T., Granato, D. C., Rivera, C., Domingues, R. R., et al. (2018). Combining Discovery and Targeted Proteomics Reveals a Prognostic Signature in Oral Cancer. Nat. Commun. 9 (1), 3598. doi:10.1038/s41467-018-05696-2
Carrot-Zhang, J., Chambwe, N., Damrauer, J. S., Knijnenburg, T. A., Robertson, A. G., Yau, C., et al. (2020). Comprehensive Analysis of Genetic Ancestry and its Molecular Correlates in Cancer. Cancer Cell 37 (5), 639–654. doi:10.1016/j.ccell.2020.04.012
Chai, A. W. Y., Lim, K. P., and Cheong, S. C. (2020). Translational Genomics and Recent Advances in Oral Squamous Cell Carcinoma. Semin. Cancer Biol. 61, 71–83. doi:10.1016/j.semcancer.2019.09.011
Chen, Z. F., Wang, Y., Sun, L. L., Ding, S. Y., and Jinag, H. (2020). LncRNA SNHG20 Enhances the Progression of Oral Squamous Cell Carcinoma by Regulating the miR-29a/DIXDC1/Wnt Regulatory axis. Eur. Rev. Med. Pharmacol. Sci. 24 (10), 5436–5445. doi:10.26355/eurrev_202005_21328
Chow, R. K. K., Tsz-Kwan Sin, S., Liu, M., Li, Y., Man Chan, T. H., Song, Y., et al. (2017). AKR7A3 Suppresses Tumorigenicity and Chemoresistance in Hepatocellular Carcinoma through Attenuation of ERK, C-Jun and NF-Κb Signaling Pathways. Oncotarget 8 (48), 83469–83479. doi:10.18632/oncotarget.12726
da Silva, S. D., Hier, M., Mlynarek, A., Kowalski, L. P., and Alaoui-Jamali, M. A. (2012). Recurrent Oral Cancer: Current and Emerging Therapeutic Approaches. Front. Pharmacol. 3, 149. doi:10.3389/fphar.2012.00149
Dombernowsky, S. L., Samsøe-Petersen, J., Petersen, C. H., Instrell, R., Hedegaard, A.-M. B., Thomas, L., et al. (2015). The Sorting Protein PACS-2 Promotes ErbB Signalling by Regulating Recycling of the Metalloproteinase ADAM17. Nat. Commun. 6, 7518. doi:10.1038/ncomms8518
Elmore, L. W., Greer, S. F., Daniels, E. C., Saxe, C. C., Melner, M. H., Krawiec, G. M., et al. (2021). Blueprint for Cancer Research: Critical Gaps and Opportunities. CA A. Cancer J. Clin. 71 (2), 107–139. doi:10.3322/caac.21652
Ferlay, J., Soerjomataram, I., Dikshit, R., Eser, S., Mathers, C., Rebelo, M., et al. (2015). Cancer Incidence and Mortality Worldwide: Sources, Methods and Major Patterns in GLOBOCAN 2012. Int. J. Cancer 136 (5), E359–E386. doi:10.1002/ijc.29210
Ferreira, A., Carvalho, S., Granville-Garcia, A., Sarmento, D., Agripino, G., Abreu, M., et al. (2021). Survival and Prognostic Factors in Patients with Oral Squamous Cell Carcinoma. Med. Oral 26 (3), e387–e392. doi:10.4317/medoral.24242
Fu, J., Li, K., Zhang, W., Wan, C., Zhang, J., Jiang, P., et al. (2020). Large-scale Public Data Reuse to Model Immunotherapy Response and Resistance. Genome Med. 12 (1), 21. doi:10.1186/s13073-020-0721-z
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. Sci. Signal. 6 (269), pl1. doi:10.1126/scisignal.2004088
Gentles, A. J., Newman, A. M., Liu, C. L., Bratman, S. V., Feng, W., Kim, D., et al. (2015). The Prognostic Landscape of Genes and Infiltrating Immune Cells across Human Cancers. Nat. Med. 21 (8), 938–945. doi:10.1038/nm.3909
Hagenaars, S. C., Groot, S., Cohen, D., Dekker, T. J. A., Charehbili, A., Meershoek‐Klein Kranenbarg, E., et al. (2021). Tumor‐stroma Ratio Is Associated with Miller‐Payne Score and Pathological Response to Neoadjuvant Chemotherapy in HER2 ‐negative Early Breast Cancer. Int. J. Cancer 149 (5), 1181–1188. doi:10.1002/ijc.33700
Hayes, T. F., Benaich, N., Goldie, S. J., Sipilä, K., Ames-Draycott, A., Cai, W., et al. (2016). Integrative Genomic and Functional Analysis of Human Oral Squamous Cell Carcinoma Cell Lines Reveals Synergistic Effects of FAT1 and CASP8 Inactivation. Cancer Lett. 383 (1), 106–114. doi:10.1016/j.canlet.2016.09.014
Ho, S. Y., Wu, W. S., Lin, L. C., Wu, Y. H., Chiu, H. W., Yeh, Y. L., et al. (2019). Cordycepin Enhances Radiosensitivity in Oral Squamous Carcinoma Cells by Inducing Autophagy and Apoptosis through Cell Cycle Arrest. Int. J. Mol. Sci. 20 (21), 5366. doi:10.3390/ijms20215366
Hou, J. Z., Xi, Z. Q., Niu, J., Li, W., Wang, X., Liang, C., et al. (2018). Inhibition of PIKfyve Using YM201636 Suppresses the Growth of Liver Cancer via the Induction of Autophagy. Oncol. Rep. 41 (3), 1971–1979. doi:10.3892/or.2018.6928
Hsu, T.-N., Huang, C.-M., Huang, C.-S., Huang, M.-S., Yeh, C.-T., Chao, T.-Y., et al. (2019). Targeting FAT1 Inhibits Carcinogenesis, Induces Oxidative Stress and Enhances Cisplatin Sensitivity through Deregulation of LRP5/WNT2/GSS Signaling Axis in Oral Squamous Cell Carcinoma. Cancers 11 (12), 1883. doi:10.3390/cancers11121883
Huynh-Le, M. P., Fan, C. C., Karunamuni, R., Thompson, W. K., Martinez, M. E., Eeles, R. A., et al. (2021). Polygenic hazard Score Is Associated with Prostate Cancer in Multi-Ethnic Populations. Nat. Commun. 12 (1), 1236. doi:10.1038/s41467-021-21287-0
Jiagge, E., Chitale, D., and Newman, L. A. (2018). Triple-Negative Breast Cancer, Stem Cells, and African Ancestry. Am. J. Pathol. 188 (2), 271–279. doi:10.1016/j.ajpath.2017.06.020
Jiang, P., Gu, S., Pan, D., Fu, J., Sahu, A., Hu, X., et al. (2018). Signatures of T Cell Dysfunction and Exclusion Predict Cancer Immunotherapy Response. Nat. Med. 24 (10), 1550–1558. doi:10.1038/s41591-018-0136-1
Jin, J., Zhao, L., and Li, Z. (2016). The E3 Ubiquitin Ligase RNF135 Regulates the Tumorigenesis Activity of Tongue Cancer SCC25 Cells. Cancer Med. 5 (11), 3140–3146. doi:10.1002/cam4.832
Johnson, D. E., Burtness, B., Leemans, C. R., Lui, V. W. Y., Bauman, J. E., and Grandis, J. R. (2020). Head and Neck Squamous Cell Carcinoma. Nat. Rev. Dis. Primers 6 (1), 92. doi:10.1038/s41572-020-00224-3
König, I. R., Fuchs, O., Hansen, G., von Mutius, E., and Kopp, M. V. (2017). What Is Precision Medicine? Eur. Respir. J. 50 (4), 1700391. doi:10.1183/13993003.00391-2017
Kuo, T. M., Nithiyanantham, S., Lee, C. P., Hsu, H. T., Luo, S. Y., Lin, Y. Z., et al. (2019). Arecoline N‐oxide Regulates Oral Squamous Cell Carcinoma Development through NOTCH1 and FAT1 Expressions. J. Cel Physiol 234 (8), 13984–13993. doi:10.1002/jcp.28084
Li, J., Lu, Y., Akbani, R., Ju, Z., Roebuck, P. L., Liu, W., et al. (2013). TCPA: a Resource for Cancer Functional Proteomics Data. Nat. Methods 10 (11), 1046–1047. doi:10.1038/nmeth.2650
Lu, Y., Zheng, Z., Yuan, Y., Pathak, J. L., Yang, X., Wang, L., et al. (2021). The Emerging Role of Exosomes in Oral Squamous Cell Carcinoma. Front. Cel Dev. Biol. 9, 628103. doi:10.3389/fcell.2021.628103
Meng, X., Lou, Q. Y., Yang, W. Y., Wang, Y. R., Chen, R., Wang, L., et al. (2021). The Role of Non‐coding RNAs in Drug Resistance of Oral Squamous Cell Carcinoma and Therapeutic Potential. Cancer Commun. 41 (10), 981–1006. doi:10.1002/cac2.12194
Miyazaki, T., Miyazaki, T., Nakashiro, K., Yamagata, H., Isokane, M., Goda, H., et al. (2011). Human FAT1 Cadherin Controls Cell Migration and Invasion of Oral Squamous Cell Carcinoma through the Localization of β-catenin. Oncol. Rep. 26 (3), 587–592. doi:10.3892/or.2011.1324
Newman, A. M., Steen, C. B., Liu, C. L., Gentles, A. J., Chaudhuri, A. A., Scherer, F., et al. (2019). Determining Cell Type Abundance and Expression from Bulk Tissues with Digital Cytometry. Nat. Biotechnol. 37 (7), 773–782. doi:10.1038/s41587-019-0114-2
Ou, J. Y., Kirchhoff, A. C., and Hanson, H. A. (2020). Air Pollution across the Cancer Continuum: Extending Our Understanding of the Relationship between Environmental Exposures and Cancer. Cancer Epidemiol. Biomarkers Prev. 29 (10), 1876–1879. doi:10.1158/1055-9965.Epi-19-1588
Peng, P., Wu, W., Zhao, J., Song, S., Wang, X., Jia, D., et al. (2016). Decreased Expression of Calpain-9 Predicts Unfavorable Prognosis in Patients with Gastric Cancer. Sci. Rep. 6, 29604. doi:10.1038/srep29604
Price, A. L., Patterson, N. J., Plenge, R. M., Weinblatt, M. E., Shadick, N. A., and Reich, D. (2006). Principal Components Analysis Corrects for Stratification in Genome-wide Association Studies. Nat. Genet. 38 (8), 904–909. doi:10.1038/ng1847
Salgado, R., Denkert, C., Demaria, S., Sirtaine, N., Klauschen, F., Pruneri, G., et al. (2015). The Evaluation of Tumor-Infiltrating Lymphocytes (TILs) in Breast Cancer: Recommendations by an International TILs Working Group 2014. Ann. Oncol. 26 (2), 259–271. doi:10.1093/annonc/mdu450
Sequeira, I., Rashid, M., Tomás, I. M., Williams, M. J., Graham, T. A., Adams, D. J., et al. (2020). Genomic Landscape and Clonal Architecture of Mouse Oral Squamous Cell Carcinomas Dictate Tumour Ecology. Nat. Commun. 11 (1), 5671. doi:10.1038/s41467-020-19401-9
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Shi, Y., Zhao, Y., Zhang, Y., AiErken, N., Shao, N., Ye, R., et al. (2018). TNNT1 Facilitates Proliferation of Breast Cancer Cells by Promoting G1/S Phase Transition. Life Sci. 208, 161–166. doi:10.1016/j.lfs.2018.07.034
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71 (3), 209–249. doi:10.3322/caac.21660
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47 (D1), D607–d613. doi:10.1093/nar/gky1131
van Pelt, G. W., Kjær-Frifeldt, S., van Krieken, J. H. J. M., Al Dieri, R., Morreau, H., Tollenaar, R. A. E. M., et al. (2018). Scoring the Tumor-Stroma Ratio in colon Cancer: Procedure and Recommendations. Virchows Arch. 473 (4), 405–412. doi:10.1007/s00428-018-2408-z
Veatch, J. R., Simon, S., and Riddell, S. R. (2021). Tumor-infiltrating Lymphocytes Make Inroads in Non-small-cell Lung Cancer. Nat. Med. 27 (8), 1339–1341. doi:10.1038/s41591-021-01445-z
Vineis, P., and Fecht, D. (2018). Environment, Cancer and Inequalities-The Urgent Need for Prevention. Eur. J. Cancer 103, 317–326. doi:10.1016/j.ejca.2018.04.018
Wang, J., Wang, Y., Kong, F., Han, R., Song, W., Chen, D., et al. (2020). Identification of a Six‐gene Prognostic Signature for Oral Squamous Cell Carcinoma. J. Cel Physiol 235 (3), 3056–3068. doi:10.1002/jcp.29210
Wang, S.-H., Chen, Y.-L., Hsiao, J.-R., Tsai, F.-Y., Jiang, S. S., Lee, A. Y.-L., et al. (2021). Insulin-like Growth Factor Binding Protein 3 Promotes Radiosensitivity of Oral Squamous Cell Carcinoma Cells via Positive Feedback on NF-κB/IL-6/ROS Signaling. J. Exp. Clin. Cancer Res. 40 (1), 95. doi:10.1186/s13046-021-01898-7
Wang, X., Steensma, J. T., Bailey, M. H., Feng, Q., Padda, H., and Johnson, K. J. (2018). Characteristics of the Cancer Genome Atlas Cases Relative to U.S. General Population Cancer Cases. Br. J. Cancer 119 (7), 885–892. doi:10.1038/s41416-018-0140-8
Weiße, J., Rosemann, J., Müller, L., Kappler, M., Eckert, A. W., Glaß, M., et al. (2021). Identification of Lymphocyte Cell-specific Protein-Tyrosine Kinase (LCK) as a Driver for Invasion and Migration of Oral Cancer by Tumor Heterogeneity Exploitation. Mol. Cancer 20 (1), 88. doi:10.1186/s12943-021-01384-w
Yang, J. J., Cheng, C., Devidas, M., Cao, X., Fan, Y., Campana, D., et al. (2011). Ancestry and Pharmacogenomics of Relapse in Acute Lymphoblastic Leukemia. Nat. Genet. 43 (3), 237–241. doi:10.1038/ng.763
Ye, Y.-k., Bi, X.-c., He, H.-c., Han, Z.-d., Dai, Q.-s., Liang, Y.-x., et al. (2010). CK20 and Ki-67 as Significant Prognostic Factors in Human Bladder Carcinoma. Clin. Exp. Med. 10 (3), 153–158. doi:10.1007/s10238-009-0088-3
Yin, Y., Tan, Y., Yao, Y., Lu, N., and Zhang, F. (2020). SNHG12/miR‐326/E2F1 Feedback Loop Facilitates the Progression of Oral Squamous Cell Carcinoma. Oral Dis. 26, 1631–1639. doi:10.1111/odi.13458
Zhang, J., Chu, M., Tian, Z., Xie, K., Cui, Y., Liu, L., et al. (2021). Clinical Profile of Fatal Familial Insomnia: Phenotypic Variation in 129 Polymorphisms and Geographical Regions. J. Neurol. Neurosurg. Psychiatry, 1–7. doi:10.1136/jnnp-2021-327247
Keywords: OSCC, genetic ancestry, genomic characteristics, TCGA, prognosis
Citation: Guo J, Liu X, Zeng Y, Liang T, Tang K, Guo J and Zheng W (2021) Comprehensive Analysis of the Effects of Genetic Ancestry and Genetic Characteristics on the Clinical Evolution of Oral Squamous Cell Carcinoma. Front. Cell Dev. Biol. 9:678464. doi: 10.3389/fcell.2021.678464
Received: 09 March 2021; Accepted: 22 November 2021;
Published: 07 December 2021.
Edited by:
Alberto Danese, University of Ferrara, ItalyReviewed by:
Griselda Escobedo, University of Guadalajara, MexicoSrinivas Patnaik, KIIT University, India
Copyright © 2021 Guo, Liu, Zeng, Liang, Tang, Guo and Zheng. 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: Junfeng Guo, Z3VvamZrcUAxNjMuY29t; Weiwei Zheng, b3ZlcmZvZXJAdmlwLnNpbmEuY29t