Racial difference of mutational signature in hepatocellular carcinoma

Bai-Feng Zhang, Xin-Yuan Guan

1Department of Clinical Oncology, The University of Hong Kong-Shenzhen Hospital, Shenzhen 518053, GuangDong, China.

2Department of Clinical Oncology, The University of Hong Kong, Hong Kong, China.

3State Key Laboratory for Liver Research, The University of Hong Kong, Hong Kong, China.

Abstract

Keywords: Mutational signatures, aflatox in exposure, tobacco chewing, aristolochic acid exposure

INTRODUCTION

Hepatocellular carcinoma (HCC), the most common type of primary liver cancer, is the fourth most common cause of cancer-related deaths worldwide[1]. Numerous epidemiologic studies demonstrated that the incidence rates and mortality of HCC patients varied by geographical region and socioeconomic status,such as high incidence in East Asia and Mongolia, the intermediate incidence rate in Europe and North America, and low incidence rate in South-Central Asia[2]. The incidence of HCC is closely associated to chronic viral hepatitis B or C infection and toxin exposure such as alcohol, aflatoxin, tobacco, or nonalcoholic steatohepatitis[3-5].

Over the past two decades, large-scale next-generation sequencing (whole genome and exome sequencing)has been performed in enormous studies to explore the genomic landscape of HCC, including recurrent driver mutations, copy number alterations, related dysfunctional biology pathways, and intratumor heterogeneity[4,6,7-12], which provided abundant resources to study the tumorigenesis of HCC. Somatic mutations in cancer genomes accumulated in mutational processes of both exogenous and endogenous factors occurring through the evolution history of cancer cell[13]. Each mutational process is reflective of the tumor intrinsic features including aging and mismatch despair deficiency and previous exposures to exogenous factors such as tobacco and drugs. In total, 94 single-base substitutions signatures with clear etiologies have been identified in previous pan-cancer study of over 30 cancer type[14]. The mutational signature analysis could provide evidence to build the association of exogenous and endogenous factor exposures with specific cancer risks or preferred treatments and further explain the epidemiology difference of HCC patients.

To study whether there is any racial difference of mutational signatures in HCC patients, we collected the genomic mutation profile of 589 HCC patients with different racial origins from previous studies[11,12],including 231 Asian-Korea HCC patients, 156 Asian-TCGA patients, 176 White-TCGA patients, and 16 Black-TCGA patients. Mutational signature extraction analysis was performed to compare their composition difference across races, and their influence on clinical outcomes was further evaluated.

METHODS

Data sources and clinical information

This study included 589 hepatocellular carcinoma patients in total [231 patients from Korea cohort[11]and 358 patients from The Cancer Genome Atlas (TCGA) cohort[12]]. Single nucleotide variants in mutation annotation format (MAF) from whole-exome sequencing and clinical information data, including survival time, race, and age, were downloaded from the cBioPortal website (https://www.cbioportal.org/datasets).Patients without clear race definition were removed, and the remaining patients were grouped into Asian-Korea, Asian-TCGA, White-TCGA, and Black-TCGA populations.

Mutation signature analysis

We applied the computational framework SigProfiler to decipher mutational signature profiles of 589 HCC patients and assign the contributions of each signature to each patient, based on the previously described methodology[14-17]. We fed the computational framework with the MAF profile of single nucleotide variants for all patients as input, and then ran the framework with two main steps: The first step, named SigProfilerExtraction, was based on somatic mutations in sequence context and their distributions in each patient, and it used multiple non-negative matrix factorization (NMF) iterations (10,000-1,000,000) to decipher a minimal set of mutational signatures that optimally explains the faction of each mutation context type in each mutational signature and estimate their activity for each sample. The second step, named SigProfilerAttribution, was responsible for accurately estimating the number of somatic mutations related to each extracted mutational signature in each patient. The optimal number of signatures as seven was determined by the trade-off of the mean sample cosine distance and average stability of solutions in the range from 1 to 10 [Figure 1A].

Proposed etiology annotation of mutational signatures

We used 94 known Catalogue of Somatic Mutations in Cancer (COSMIC) mutational signatures(https://cancer.sanger.ac.uk/signatures/sbs/) to perform decomposition analysis of the extracted signatures in this study. The correlation and cosine similarity of original signatures extracted in this study with reconstructed signatures consisting of known COSMIC signatures were calculated, and their proposed etiology was also annotated [Table 1].

综上所述，咪达唑仑静脉麻醉与利多卡因局部麻醉联合的麻醉方式正在被越来越多的医院科室引入临床应用之中。该种新型联合式的麻醉方式的使用，使患者在检查过程中全程呈现全麻醉状态，该种联合型麻醉法除了可以通过使患者的紧张、恐惧情绪获得缓解，从而使检查进行的更为顺利外，其也能够使患者以呛咳、反射性呕吐、气道收缩痉挛为主的各种机体应激所导致的不良现象的发生进行最大程度的减少，并以此使患者获得更好的检查体验。

Unsupervised clustering of patients based on mutation signature components

Based on the contributions of signatures to each patient, we calculated their fraction among HCC patients to create a numerical matrix (patients as columns, signatures as rows) and then applied the NMF “lee”algorithm to perform molecular subtyping of HCC patients. After the manual inspection, we chose the K =7 solution and reported seven subgroups.

Statistical analysis

Two-sided Mann-Whitney tests were performed using Wilcoxon signed-rank test to generate the empiricalPvalues and plotted with the R package “ggpubr”.

Survival analysis

Chi-square test statistics were computed using log-rank tests, and Kaplan-Meier curves were plotted using the R package “survival” and “survminer”.

RESULTS

Identification of mutational signatures in HCC patients

Somatic mutations occur throughout life and display distinctive tumor cell-intrinsic patterns, including DNA replication machinery and defective DNA repair, as well as exogenous processes such as tobacco exposure and chemotherapy treatment. Different mutational processes generate distinct combinations of mutation types, which can be classified into 96 classes constituted by the six base substitutions C ＞ A, C ＞ G,C ＞ T, T ＞ A, T ＞ C, and T ＞ G. In total, 84,585 single nucleotide variants from 589 HCC patients were downloaded from a public database (cBioPortal) and categorized into a 96 mutation classes × 589 patient numeric matrix. We then applied SigProfiler software to extract optimal mutational signatures of HCC patients. Finally, the seven mutational signatures were identified in our study based on the trade-off curve of error and stability with the range from 1 to 10 [Figure 1A], which showed a distinct mutational pattern:SBS96A was characterized predominantly by C ＞ T and C ＞ G mutations; SBS96B was composed predominantly of T ＞ C mutations; SBS96C was characterized primarily by C ＞ A mutations at GCA and GCC trinucleotides (the mutated base is underlined); SBS96D was composed predominantly of C ＞ T at ACG, CCG, and GCG trinucleotides; SBS96E was composed primarily of C ＞ T at ACC, CCC, and TCC trinucleotides; SBS96F included enrichment of T ＞ A mutations at CTG trinucleotides; and SBS96G wascomposed mainly of C ＞ A at CCA and CCG trinucleotides [Figure 1B]. After dividing our study cohorts into four populations of Asian-Korea (n= 231), Asian-TCGA (n= 156), White-TCGA(n= 176), and Black-TCGA (n= 16), we observed that the seven mutational signatures showed diverse contributions to patients of the four different populations in the estimated count or composition fraction level [Figure 1C], indicating the mutational signature heterogeneity of HCC patients.

Table 1. De novo extraction of mutational signatures and their proposed etiology annotation

Figure 1. Deciphering mutational signatures of 589 hepatocellular carcinoma patients: (A) the distribution of error (mean sample cosine distance) and average stability of solutions for a range from 1 to 10; (B) mutational signatures identified by the SigProfiler; and (C) the estimated count and fractions of seven mutational signatures in HCC patients.

To investigate which etiology was proposed for the seven mutational signatures, we performed decomposition analysis with 94 known cosmic signatures [Table 1]. SBS96A was reconstructed with several different known clock-like COSMIC signatures (SBS1, 5.26%; SBS5, 34.14%; SBS40, 60.60%) and could be defined as a clock-like signature. SBS96B and SBS96D consisted mainly of cosmic signatures associated with defective DNA mismatch repair (MMR) and could be characterized as MMR-like signatures. SBS96C was also reconstructed by 43.54% aflatoxin exposure SBS24 and 50.88% tobacco chewing SBS29. SBS96F consisted of 98.88% SBS22, and could be defined as an aristolochic acid-exposure signature. SBS96E consisted of multiple COSMIC signatures related to clock-like, ultraviolet light exposure, and unknown factor. SBS96G consisted of signatures arising from possible sequencing artifacts [Table 1].

Racial difference of mutational signatures in HCC patients

To explore whether there were any racial differences in mutational signatures, we compared the faction of seven mutation signatures among four populations, including Asian-Korea, Asian-TCGA, White-TCGA,and Black-TCGA. In general, the clock-like SBS96A signature was the major signature of HCC patients with the highest number of somatic mutations per megabase, while the SBS96G signatures as the minor signature showed the lowest activity but still existed in 32.4% of HCC patients (191/589) [Figure 2A]. SBS96A clocklike signature, related to the age of patients, and SBS96E signature, with unclear proposed etiology, showed no obvious difference among the four populations, suggesting that these two signatures were common in HCC patients [Figure 2B]. Compared with the Asian-Korea patients, the fractions of both SBS96B and SBS96D, as the two subsets of defective DNA MMR-related signatures, showed no significant changes in Black-TCGA and Asian-TCGA but showed adverse significant changes in White-TCGA (P＜ 0.05),indicating the heterogeneous mechanisms of defective DNA MMR in HCC patients [Figure 2B].Interestingly, the fraction of SBS96F signature showed no obvious difference in the two Asian populations(Asian-Korea and Asian-TCGA) but was significantly lower in White-TCGA patients [Figure 2B]. The higher fraction of aristolochic acid-exposure signature in Asian patients might be caused by the usage of Chinese herbs that produce aristolochic acid[18-22]. Intriguingly, Asian-Korea populations showed a significantly higher fraction of SBS96C associated with tobacco chewing and aflatoxin exposure compared with the three populations from the TCGA cohort [Figure 2B], suggesting high levels of aflatoxin contamination in food and environment. It is also worth noting that the fraction of SBS96F signature arising from possible sequencing artifacts was observed to be significantly higher in the Asian-Korea cohort than in the other three populations from the TCGA cohort [Figure 2B]. These results demonstrate that the composition of mutational signatures was diverse across the HCC populations.

To investigate whether HCC patients could be characterized by the mutational signature components, we applied a non-negative factorized matrix algorithm on subtyping HCC patients based on the fraction change pattern of seven mutational signatures. We identified seven subgroups of HCC patients: S1-SBS96DMMR dominated, S2-SBS96E-unclear dominated, S3-SBS96A-clock-like dominated, S4-SBS96B-MMR dominated, S5-SBS96C-tobacco and aflatoxin dominated, S6-SBS96G-unclear dominated, and S7-SBS96Faristolochic dominated subgroup patients [Figure 2C]. Interestingly, S6-SBS96G-unclear and S5-SBS96Ctobacco, and aflatoxin exposure dominated patients were enriched in the Asian-Korea population[Figure 2C]. These results suggest that the dominant architectures of mutational signature across HCC patients are distinct and could be used to group HCC patients into heterogenous subtypes.

Figure 2. Racial differences in the mutational signatures of HCC. (A) The distribution of somatic mutations per megabase of signatures.(B) Comparison of prevalence of signatures among Asian-Korea, Asian-TCGA, White-TCGA, and Black-TCGA. Wilcoxon signed-rank test and Bonferroni corrected: P ＞ 0.05, 0.01 ＜ P ≤ 0.05, 0.001 ＜ P ≤ 0.01, 0.0001 ＜ P ≤ 0.001, and 0 ＜ P ≤ 0.0001 are shown as ns, ★, ★★,★★★, and ★★★★, respectively. (C) Unsupervised clustering of HCC patients based on the fractions of seven mutational signatures in each patient.

Prognostic analysis of distinct mutational signatures-dominant HCC subgroup patients

Figure 3. Prognostic analysis of seven mutational signatures-dominant HCC subgroup patients: (A) disease-free survival Kaplan-Meier curve; and (B) overall survival Kaplan-Meier curve. Log-rank sum tests.

DISCUSSION

Hepatocellular carcinoma is the primary type of liver cancer and is one of the most common malignant tumors in Asia. The genomic landscape of HCC was well-characterized by numerous previous studies, but their differences across races remain unclear. In this study, multiple significant racial differences of mutational signatures were demonstrated. We identified more prevalence of SBS96F-aristolochic acid exposure signature in HCC patients from Asia, which might be caused by the common usage of Chinese herbs in this area and is consistent with a previous study[18], suggesting that the detection and alternative usage of ingredients generating aristolochic acid will help reduce the incidence of HCC in Asia. We also showed a higher level of SBS96B-MMR at T ＞ C mutations in Asian patients and a higher level of SBS96DMMR at C ＞ T mutations in White patients. MMR was defined as the predictor to sensitize cancer cells to immunotherapy response and showed high efficiency in the immune checkpoint blockade therapy of patients across many cancer types[23-25]. Differential enrichment of two signatures associated with defective DNA MMR mechanisms across races might explain why HCC patients showed a low response rate to immunotherapy and demonstrated the importance of classifying the MMR-like signature subtype before applying immunotherapy on patients with distinct racial origins.

Clustering of HCC patients based on their mutational signature exposures could stratify the patients into different prognostic subgroups. The HCC patients with the dominated prevalence of SBS96C-tobacco chewing and aflatoxin showed worse disease-free survival prognosis and were enriched in Asian-Korea patients, indicating that the level of aflatoxin contamination in food and environment was higher in Korea compared with North America, which might reflect the correlation of the mortality of HCC patients with the comprehensive socioeconomic status of a region. The SBS96G signature was also observed to be significantly higher in Asian-Korea patients, and their dominated subgroup patients showed better prognosis for both disease-free and overall survival probability, which might suggest the potential role of this signature to predict better prognosis of HCC patients. SBS96G signature was still observed in about 32.4% of HCC patients, and its etiology remained unclear. We also could not exclude the possibility of sequencing bias as the SBS96G signature could be reconstructed by two signatures generated by possible sequencing artifacts reported in a previous study[14].

Overall, our study identified multiple mutational signatures that varied substantially among different racial populations and evaluated their potential usage for predicting the clinical outcome of HCC patients.However, further validation in a larger cohort is needed before applying it to the personalized treatment of patients.

DECLARATIONS

Authors’ contributions

Made substantial contributions to conception and design of the study and performed data analysis and interpretation: Zhang BF, Guan XY

Wrote the draft and revised the manuscript: Zhang BF, Guan XY

Availability of data and materials

Mutation and clinical information data could be obtained from the cBioPortal website(https://www.cbioportal.org/datasets, Korea cohort: Ahn et al. Hepatology 2014, TCGA cohort: TCGA, Cell 2018). Supporting data produced by this study were uploaded to a public website(https://github.com/zhangbaifeng/HCC-HR).

Financial support and sponsorship

This work was supported by supported by grants from the Hong Kong Research Grant Council (RGC)including Collaborative Research Funds (C7065-18GF and C7026-18GF), Theme-based Research Scheme(T12-704/16-R), Shenzhen Science and Technology programs (KQTD20180411185028798), Shenzhen Fundamental Research Program (JCYJ20180508153249223). Professor Guan XY is Sophie YM Chan Professor in Cancer Research.

Conflicts of interest

Both authors declared that there are no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2021.