Journal of Oncology
Volume 2020 (2020), Article ID 3275315, 8 pages
https://doi.org/10.1155/2020/3275315

Research Article

Next-Generation Sequencing Characterizes the Landscape of Somatic Mutations and Pathways in Metastatic Bile Tract Carcinoma

Qi Li,1,2 Qifan Zhang,³ Xiao Cheng,¹ Xie Weng,⁴ Mian Chen,⁵ Xiaoyun Hu,¹ Jing Huang,¹ and Jinzhang Chen1,2

¹State Key Laboratory of Organ Failure Research, Guangdong Provincial Key Laboratory of Viral Hepatitis Research, Department of Infectious Diseases, Nanfang Hospital, Southern Medical University, No. 1838, North Guangzhou Avenue, Guangzhou, Guangdong 510515, China
²Department of Oncology, Nanfang Hospital, Southern Medical University, Guangzhou 510515, China
³Department of Hepatobiliary Surgery, Nanfang Hospital, Southern Medical University, Guangzhou 510515, China
⁴Cancer Center, Integrated Hospital of Traditional Chinese Medicine, Southern Medical University, Guangzhou, Guangdong 510315, China
⁵Transplant Immunology Laboratory, Churchill Hospital, Oxford University Hospitals NHS Trust, Old Road, Headington, Oxford OX3 8HF, UK

Correspondence should be addressed to Jinzhang Chen; jinzhang_chen1@163.com

Received 9 March 2020; Revised 7 July 2020; Accepted 8 July 2020; Published 4 September 2020

Academic Editor: Thomas E. Adrian

Copyright © 2020 Qi Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Purpose. Tumor metastasis remains the leading cause of cancer-related mortality in biliary tract cancer. The etiology and mechanism of bile tract carcinoma metastasis are unclear. Methods. The primary tumor and blood samples of 14 patients with biliary tract cancer were collected, followed by nucleic acid extraction and library construction. Target sequencing with 556 panel genes and WES were performed to detect the hot spot genes variations. Bioinformatics was used to comprehensively analyze the sequencing data of these samples, including the differences of tumor mutation burden and signaling pathways. Results. The results showed that the mutation frequency of TP53 gene was the highest and the mutations of CTNNB1, EPHA7, ARID2, and PIK3CA were only found in metastatic samples. The TMB mean values of metastatic and non-metastatic groups were 12.97 and 10.38 mutations per Mb, respectively. There were significant differences in the enrichment pathways of cellular components between the tumor metastasis and non-metastatic samples. Conclusions. We identified multiple pathway differences, which helps us better understand metastatic biliary tumors and design clinical therapy for personalized medicine.

1. Introduction

Biliary tract carcinoma (BTC) is an invasive adenocarcinoma that originates from the epithelial cells of the biliary tract, including the intrahepatic and extrahepatic bile ducts and the gallbladder [1]. Although BTCs are considered a rare tumor, they account for about 10%–15% of all primary liver cancers [2]. Only 10% of patients have early disease and are considered candidates for surgical excision, and the recurrence rate is high despite therapeutic surgery. Most patients with locally advanced or metastatic BTCs have poor prognosis and overall survival (OS) is less than 12 months [3].

In recent years, advances in sequencing technology have promoted the application of sequencing technology in cancer research, giving scientists an in-depth understanding of these cancers at the molecular level [4, 5]. However, metastases that cause 90% of cancer deaths have been reported in prostate, breast, and colorectal cancers [6–8]. Even with a small number of metastatic cancers in tumors of the biliary tract carcinomas, only a few specific genes have been studied [9]. Some studies have shown that the process of biliary cancer tumor metastasis has complex signal pathway crossover. Neurotensin can enhance EGFR/AKT signal pathway and promote the metastasis of cholangiocarcinoma cells [10]. LncRNA-UCA1 promoted the activation of ERK/MAPK pathway by regulating the expression of miR-122 and its downstream gene mRNA CLIC1 [11]. Therefore, a better understanding of the biological and phenotypic evolution of BTCs and their molecular and genetic mechanisms during the metastatic process is crucial.

In order to further study the genetic characteristics of tumor metastasis in the biliary tract cancer system, we conducted high-depth sequencing in 14 patients with biliary tract cancer. Somatic cell mutations, tumor mutation burden, TITV, and molecular signaling pathways of mutated genes in patients with metastasis and non-metastasis were analyzed. Exploring the mechanism of metastasis in biliary tract cancer can provide effective guidance for potential clinical interventions, thus improving the long-term prognosis of patients.

2. Methods

2.1. Patient Specimen Acquisition

Primary tumor tissue and blood from 14 patients with BTCs were collected. According to the digestive system tumor classification of WHO in 2010, tumor samples were identified by hematoxylin and eosin stained slides. The tumor tissue samples included were histologically confirmed to be adenocarcinoma by two molecular pathologists, and the tumor cell content exceeded 70%. This study was conducted in accordance with the Helsinki Declaration and was approved by the first affiliated southern hospital ethics committee [12]. Prior to inclusion in the study, all patients received informed written consent and all groups participating in this study approved this work.

2.2. DNA Extraction and Qualification

DNA was extracted from tissue sections of tumor samples and matched with the DNA of leukocytes as germline mutation. DNeasy Blood and Tissue Kit (69504, Qiagen, Venlo, Netherlands) was used according to the manufacturer’s instructions. The content of DNA was determined by Agilent Bioanalyzer (USA).

2.3. Target Genes Sequencing and Whole-Exome Sequencing

The targeted capture pulldown and exon-wide libraries from genomic DNA were generated through the xGen® Exome Research Panel (Integrated DNA Technologies, Inc., Illinois, USA) and the TruePrep DNA Library Prep Kit V2 for Illumina (#TD501, Vazyme, Nanjing, China). The captured libraries sequencing was performed as paired-end reads on the Illumina NoveSeq platform.

2.4. Variant Annotation Analysis and Visual Mapping

Split-reads were re-aligned to hg19/GRch37 using blast to exclude false positive results. We used Picard tool (http://broadinstitute.github.io/picard/command-line-overview. html#CollectHsMetrics) on the computational efficiency of WES data capture. The tool is used to filter low-quality bases, repetitive bases, bases that deviate from the target, and ends that are overlapped to read at both ends due to short inserts. This tool calculated the strictest depth distribution. If the mean bait coverage of sequencing data was more than 100, we defined the data quality control as qualified. Further local rearrangements were performed with SpeedSeq to improve the alignment of individual reads [13]. Somatic mutations identification and indels were annotated through Mutect [14] and Somatic Indel Detector software [15]. ANNOVAR [16] and Oncotator [17] annotated the variant data in use and converted it to a MAF file through the maf tool [18]. The cancer driver genes were analyzed using Intogen [19], including Oncodrive FM and Oncodrive CLUST. The landscape of top driver mutation spectrum predicted by Intogen for tumors was visualized via R Script, including mutation rate and mutation subclass/subtypes (ONCODRIVEFM value ≤ 0.1).

2.5. Pathway Enrichment Analysis

Wayne map was used to show that there were differences in the frequency of somatic gene mutation between metastatic and non-metastatic group, because some mutation genes exist only in the metastatic group or in the non-metastatic group, while some genes were mutated in both. GO or Gene Ontology canonical pathways with R packages: clusterProfiler were used to analyze the genes that contain single nucleotide variants (SNV) or SV [20]. The calculation of value was based on hypergeometric distribution, and Benjamin and Hochberg methods were used for FDR correction [21].

2.6. Statistical Analyses

All the correlate clinical and biological variables were employed using the SPSS Statistics 22.0 package and ggpubr package [22] in R [23] by methods of a non-parametric test (such as Welch’s t test and Wilcoxon) when necessary.

3. Result

3.1. Patient Characteristics

To better understand the molecular mechanism of BTC metastasis, fresh-frozen tumor tissue and corresponding blood specimens from 14 BTC patients were selected for bioinformatics analysis based on stringent criteria. S08 and S09 samples were sequenced by whole exome sequencing (WES), the other 12 samples sequenced by 556 panel genes. A total of three patients had a history of alcohol consumption, and half had HBV. According to the clinical stages, there were two patients in stage I, one patient in stage IIA, three patients in stage IIIB, two patients in stage IVA, and six patients in stage IVB. In addition, nine samples including S01, S03, S04, S05, S06, S08, S09, S10, and S14 were assigned to the metastatic group, while the remaining 5 samples of S02, S07, S11, S12, and S13 were in the non-metastatic group. The detailed clinical characteristics of the patients are shown in Table 1 and Supplementary Table 1.

Table 1: Patient characteristics.

Characteristic
No. of cases
Proportion (%)

Total number
n = 14

Age, years (mean)
53.28 (30–70)

Sex

    Male
7
50.0

    Female
7
50.0

Drinking history

    Drinker
3
16.7

    Nondrinker
11
83.3

Metastatic
9
64.3

Non-metastatic
5
35.8

Tumor staging

    I
2
14.3

    IIA
1
7.1

    IIIB
3
21.4

    IVA
2
14.3

    IVB
6
42.9

State

    Alive
6
42.9

    Dead
8
57.1

3.2. Identification of Somatic Mutations in 14 Patients with BTCs

We performed NGS sequencing on DNA from 14 primary tumor tissues along with matched blood and annotated several somatic mutations using Mutect and Somatic Indel Detector: WES depth with a mean depth of 200.2X and 556 panel sequencing with a mean depth of 2559.95X Overall. We detected 119 nonsynonymous single nucleotide variations (SNVs) in all samples, including 31 nonsynonymous SNVs detected in WES and 88 nonsynonymous SNVs detected in 556 panel sequencing, which were all missense mutations (Supplementary Table 2). As shown in Figure 1, case S08 had the most SNVs, followed by case S07. We listed the top 29 genes based on the frequency of somatic mutations. Among them, TP53 (64%) was the highest mutation gene. Missense mutation was the most common type of mutation, along with nonsense, multi hit, and so on (Figure 1). Driver genes, such as CTNNB1 (22%), EP300 (22%), KMT2C (22%), and IDH1 (22%), were mutated only in the metastatic group, while XPO1 (40%) mutation was found only in non-metastatic group.

Figure 1: Somatic mutation in biliary tumors. The top figure shows the variation number of each sample. The middle picture shows the mutated genes and mutation types of the samples. The bottom chart shows the clinical information such as age, gender, clinical diagnosis, drinking history, survival status, etc. The red bars on the right represent the mutation frequency of somatic cell genes in patients with cholangiocarcinoma metastasis, and the blue bars represent the mutation frequency of somatic genes in patients without metastasis.

We also calculated tumor mutation burden (TMB) using only somatic nonsynonymous mutations. On the whole, the TMB mean values of metastatic and non-metastatic groups were 12.97 and 10.38 mutations per Mb, respectively, and there was no significant difference between these two groups,  = 0.6885 (Welch’s t test) (Figures 2(a) and 2(b)). In addition, the TMB value of this study (8.635829) was significantly higher than that of TCGA-CHOL database (1.433684), with statistical significance (Wilcoxon,  = 0.000032) (Figure 2(c)).

(a)

(b)

(c)

Figure 2: Comparative analysis of TMB difference. (a) Tumor mutation burden in metastatic and non-metastatic samples; (b) comparison of TMB between metastatic and non-metastatic samples; (c) comparison of TMB between the data in this study and the TCGA-CHOL.

3.3. Characteristics of Signaling Pathways in BTCs

In order to further characterize the functions of mutational genes and pathways involved in BTC, we used PANTHER classification system [24], an ontology-based pathway database coupled with data analysis tools. According to the metastasis of the patient, we used a Venn diagram to divide the mutant genes detected by 556 gene panel sequencing into three clusters, including metastatic genes cluster (44 genes), non-metastatic genes cluster (50 genes), and intersection genes cluster (9 genes) (Figure 3(a)).

(a)

(b)

(c)

(d)

Figure 3: Analysis of pathways involving nonsynonymous genes in BTCs. (a) Venn diagram of gene distribution in three clusters. (b) The cellular component pathway of mutant genes in the metastatic group. (c) The cellular component pathway of mutant genes in the non-metastatic group. (d) The cellular component pathway of mutant genes in the intersection group.

In the metastatic cluster, the mutant somatic genes were mainly enriched in the transferase complex, SWI/SNF complex, ATPase complex, and nBAF complex pathway (Figure 3(b)). However, the altered somatic genes in the non-metastatic cluster were mainly enriched in extrinsic component of membrane, chromosome-telomeric region, cytoplasmic ribonucleoprotein granule, DNA repair complex, and lateral plasma membrane in the pathway of cellular components (Figure 3(c)). In addition, the cellular components of mutant genes in the intersection of metastatic and non-metastatic clusters were mainly in histone removal factor TFIID complex, components of synaptic membrane, intrinsic components of presynaptic membrane, postsynaptic membrane, and activity of postsynaptic membrane enzyme modulator (Figure 3(d)).

4. Discussion

BTC is a rare malignant tumor [25]. At the initial diagnosis stage, about 60% or 70% of cholangiocarcinomas are pathologically advanced [26], and the overall survival time is between 8 and 715 days, with an average of 302 days, which is consistent with the previous statement that OS is less than 12 months [27]. Metastatic clones can occur in the early and late stages of primary tumors [8, 28]. Better understanding of the genetic characteristics of metastatic diseases may reveal the differences between the treatment weaknesses of primary and metastatic tumors and provide insights into the biology of metastasis. The metastasis of biliary malignancies is often accompanied by mutations or changes in expression levels of multiple pathway genes, including abnormal activation of proto-oncogenes and inactivation of tumor suppressor genes. We designed 556 genes panel of tumor mutation hot spot gene and sequenced it in high depth. We also found high frequency mutated genes in tumors of the biliary system that have been reported, such as TP53, ARID1B, CTNNB1, EPHA7, IDH1, and so on [29–31], while CTNNB1, EPHA7, ARID2, and PIK3CA were mutated in patients with metastasis (Figure 1). It has been reported that, in intrahepatic biliary cell carcinoma, mutations in CTNNB1 are mostly associated with alcohol intake risk factors, while mutations in TP53 are caused by HBV risk factors [32]. Frequent mutations of TP53 and ARID1A associated with chromatin remodeling and chromosomal organization may be involved in the carcinogenesis and development of intrahepatic cholangiocarcinoma [33]. We also analyzed 51 samples from the TCGA-CHOL database and found that the mutation frequency of TP53 was only (4/51) 8%, far lower than the targeted sequencing (9/14) 62% in this study (Supplementary Figure 1). This indicates that the target gene sequencing method is more accurate in detecting target genes.

Previous studies have reported that, in univariate analysis, TP53 and ARID1A were predictors of poor prognosis in cholangiocarcinoma [29]. We calculated the relationship between TP53 mutation and overall survival time, indicating that no significant difference was found ( = 0.325), which may be largely due to the low sample size (Supplementary Figure 2).

Currently, there are dozens of biomarkers related to checkpoint inhibitors, among which TMB, PD-L1, and MSI/dMMR have been verified by phase III clinical trials and are widely used in clinical practice [34–36]. TMB is a biomarker for predicting PD-1/PD-L1 immune response [37, 38]. Even though it has been reported that TMB-H alone is not suitable for predicting the immunotherapy effect of solid tumor type [39], we found that there was a significant difference in TMB between this study and TCGA, but the TMB did not exceed 20 mutations per Mb (Figure 2). For different cancer types, the setting of high TMB threshold may need more clinical studies and a large number of patient information statistics. A Japanese paper suggests that TMBs are overestimated in targeted sequencing compared to WES [40]. Therefore, we compared the TMB value of BTC samples in this study with TCGA-CHOL and found that the depth of targeted sequencing was significantly higher than that of WES.

Signaling pathways are characterized by mutations in a single gene or changes in expression, which usually involve simultaneous changes in multiple pathways, such as angiogenesis and notch signaling pathway [41]. During the process of tumor metastasis, cellular components such as the transferase complex, SWI/SNF complex, ATPase complex, and nBAF complex pathway are all adjusted accordingly. SWI/SNF, including BRM or BRG-1, ATPase subunits, controls many aspects of normal cellular function [42]. Ribonucleoprotein promotes tumor metastasis by induction of genes involved in extracellular matrix, cell movement, and angiogenesis [43]. Functional recovery is promising in cancer therapy because epigenetic inhibition regulates the expression of SWI/SNF components at least in some cases. More research is needed to unblock the role of SWI/SNF in cancer and determine how it affects tumor metastasis. This is an exciting but poorly understood molecule that may play a role in causing cancer.

5. Conclusion

Our study identified somatic mutations, TMB in metastatic and non-metastatic groups. The GO analysis showed that the metastatic and non-metastatic groups were completely different in terms of cellular composition pathways. Although the sequencing depth of 556 panel gene was very high, the detected genes were known or have been reported, and new site mutations cannot be found. In addition, our 556 panel was designed for the locus with high mutation hotspots in the Asian population, which was suitable for the detection of mutated genes in the whole cancer species, not specifically for the cholangiocarcinoma, so the detection range of mutated genes is still limited. In addition, the small sample size and the lack of matched BTC transfer samples were also limitations of this study. In conclusion, the current findings could help identify specific pathways and hot spots that altered during metastasis and provide a direction for targeted therapy of these tumors.

Data Availability

All the related software and scripts are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that no conflicts of interest exist.

Authors’ Contributions

Qi Li and Qifan Zhang contributed equally to this work

Acknowledgments

The authors acknowledge the support of Shanghai Tongshu Biotechnology Co., Ltd. This research was supported by the Natural Science Foundation of Guangdong Province (Grant no. 2017A030313645).

Supplementary Materials

Supplementary figure 1: mutational landscape of somatic alterations in TCGA-CHOL database. Supplementary figure 2: Kaplan–Meier survival analysis using the log-rank test (mantel-cox test) for comparing the differences between mutant TP53 and wild-type TP53 genes. Supplementary Table 1: patient clinical information. Supplementary Table 2: NGS results of our cohort. Support files legends. Supplementary table a: variant allele frequency (VAF) of 556 genes in target sequencing data and WES data. Supplementary table b: 556 somatic mutant genes were divided into TCGA-CHOL recognized genes and unrecognized genes in this study. Supplementary table c: mutation screening of 44 transferred genes in the TCGA-CHOL dataset. Supplementary table d: list of genes involved in 556 genes panel sequencing in this study. Supplementary table e: the list of genes related to the pathways in the three clusters. Supplementary figure a: Go-CC pathway of transferred genes in TCGA database. (Supplementary Materials)

References

P. Letelier, I. Riquelme, A. Hernández, N. Guzmán, J. Farías, and J. Roa, “Circulating microRNAs as biomarkers in biliary tract cancers,” International Journal of Molecular Sciences, vol. 17, no. 5, p. 791, 2016.

F. Eckel, T. Brunner, and S. Jelic, “Biliary cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up,” Annals of Oncology, vol. 22, no. 5, pp. v65–v69, 2011.

A. F. Hezel, V. Deshpande, and A. X. Zhu, “Genetics of biliary tract cancers and emerging targeted therapies,” Journal of Clinical Oncology, vol. 28, no. 21, pp. 3531–3540, 2010.

J. N. Weinstein, J. B M. Gordon, E. A. Collisson et al., “The cancer genome atlas pan-cancer analysis project,” Nature Genetics, vol. 45, no. 10, pp. 1113–1120, 2013.

X. Ma, Y. Liu, Y. Liu et al., “Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours,” Nature, vol. 555, no. 7696, pp. 371–376, 2018.

G. Gundem, K. Barbara, P. Van Loo et al., “The evolutionary history of lethal metastatic prostate cancer,” Nature, vol. 520, no. 7547, pp. 353–357, 2015.

K. Naxerova, J. G. Reiter, E. Brachtel et al., “Origins of lymphatic and distant metastases in human colorectal cancer,” Science, vol. 357, no. 6346, pp. 55–60, 2017.

L. R. Yates, S. Knappskog, D. Wedge et al., “Genomic evolution of breast cancer metastasis and relapse,” Cancer Cell, vol. 32, no. 2, pp. 169–184, 2017.

J. Ghosh, S. Ganguly, D. Dabkara et al., “Metachronous muscle metastasis in a case of metastatic gallbladder cancer with TP35 gene mutation: a rare case report,” South Asian Journal of Cancer, vol. 8, no. 4, p. 240, 2019.

Z. J. Su, X. Y. Liu, J. H. Zhang, S. Y. Ke, and H. J. Fei, “Neurotensin promotes cholangiocarcinoma metastasis via the EGFR/AKT pathway,” Gene, vol. 687, pp. 143–150, 2019.

L. Kong, Q. Wu, and L. Zhao, “Upregulated lncRNA-UCA1 contributes to metastasis of bile duct carcinoma through regulation of miR-122/CLIC1 and activation of the ERK/MAPK signaling pathway,” Cell Cycle, vol. 18, pp. 1–17, 2019.

World Medical Association, “Declaration of Helsinki: ethical principles for medical research involving human subjects,” Jahrbuch für Wissenschaft und Ethik, vol. 14, no. 1, pp. 233–238, 2009.

C. Chiang, R. M. Layer, G. G. Faust et al., “Speedseq: ultra-fast personal genome analysis and interpretation,” Nature Methods, vol. 12, no. 10, pp. 966–968, 2015.

K. Cibulskis, M. S. Lawrence, S. L. Carter et al., “Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples,” Nature Biotechnology, vol. 31, no. 3, pp. 213–219, 2013.

M. A. Depristo, E. Banks, R. Poplin et al., “A framework for variation discovery and genotyping using next-generation DNA sequencing data,” Nature Genetics, vol. 43, no. 5, pp. 491–498, 2011.

K. Wang, M. Li, and H. Hakonarson, “Annovar: functional annotation of genetic variants from high-throughput sequencing data,” Nucleic Acids Research, vol. 38, no. 16, p. e164, 2010.

A. H. Ramos, L. Lichtenstein, M. Gupta et al., “Oncotator: cancer variant annotation tool,” Human Mutation, vol. 36, no. 4, pp. E2423–E2429, 2015.

A. Mayakonda, D.-C. Lin, Y. Assenov, C. Plass, and H. P. Koeffler, “Maftools: efficient and comprehensive analysis of somatic variants in cancer,” Genome Research, vol. 28, no. 11, pp. 1747–1756, 2018.

A. Gonzalez-Perez, C. Perez-Llamas, J. Deu-Pons et al., “Intogen-mutations identifies cancer drivers across tumor types,” Nature Methods, vol. 10, no. 11, pp. 1081-1082, 2013.

G. Yu, L.-G. Wang, Y. Han, and Q.-Y. He, “Clusterprofiler: an r package for comparing biological themes among gene clusters,” Omics: A Journal of Integrative Biology, vol. 16, no. 5, pp. 284–287, 2012.

A. Li and R. F. Barber, “Multiple testing with the structure-adaptive benjamini-hochberg algorithm,” Journal of the Royal Statistical Society: Series B (Statistical Methodology), vol. 81, no. 1, pp. 45–74, 2019.

A. Kassambara, “ggpubr: “ggplot2” based publication ready plots,” 2018, http://www.sthda.com/english/rpkgs/ggpubr.

R Core Team, “R: a language and environment for statistical computing,” 2018, https://www.R-project.org/.

M. Huaiyu, A. Muruganujan, J. T. Casagrande, and P. D. Thomas, “Large-scale gene function analysis with the panther classification system,” Nature Protocols, vol. 8, no. 8, pp. 1551–1566, 2013.

R. Sumera, S. A. Khan, and C. L. Hallemeier, “Cholangiocarcinoma—evolving concepts and therapeutic strategies,” Nature Reviews. Clinical Oncology, vol. 15, no. 2, pp. 95–111, 2018.

M. M. Kirstein and A. Vogel, “Epidemiology and risk factors of cholangiocarcinoma,” Visceral Medicine, vol. 32, no. 6, pp. 395–400, 2016.

J. Weigt and P. Malfertheiner, “Cisplatin plus gemcitabine versus gemcitabine for biliary tract cancer,” Expert Review of Gastroenterology & Hepatology, vol. 4, no. 4, pp. 395–397, 2010.

S. Yachida, S. Jones, I. Bozic et al., “Distant metastasis occurs late during the genetic evolution of pancreatic cancer,” Nature, vol. 467, no. 7319, pp. 1114–1117, 2010.

M. Simbolo, C. Vicentini, A. Ruzzenente et al., “Genetic alterations analysis in prognostic stratified groups identified TP53 and ARID1A as poor clinical performance markers in intrahepatic cholangiocarcinoma,” Scientific Reports, vol. 8, no. 1, 2018.

S. Conrad, “Mesenchymal stem cell-derived extracellular vesicles as mediators of anti-inflammatory effects,” in Stem Cell Transplantation for Autoimmune Diseases and Inflammation, pp. 89–123, Springer, Berlin, Germany, 2019.

C. Guo, Y. Shen, and Q. Zhang, “Prediction of risk factors for postoperative pancreatic fistula using a nomogram based on the updated definition,” HPB, vol. 20, pp. S279–S280, 2018.

K. Schulze, S. Imbeaud, E. Letouzé et al., “Exome sequencing of hepatocellular carcinomas identifies new mutational signatures and potential therapeutic targets,” Nature Genetics, vol. 47, no. 5, pp. 505–511, 2015.

Z. H. Liu, “Whole-exome mutational and transcriptional landscapes of combined hepatocellular cholangiocarcinoma and intrahepatic cholangiocarcinoma reveal molecular diversity,” Biochimica et Biophysica Acta, vol. 1864, no. 6, 2018.

Q. Yang, Z. Xu, L. Zheng, L. Zhang, Q. You, and J. Sun, “Multimodal detection of PD-L1: reasonable biomarkers for immune checkpoint inhibitor,” American Journal of Cancer Research, vol. 8, no. 9, pp. 1689–1696, 2018.

L. Feng, S. Hong, J. Gao, and J. Li, “Whole-exome sequencing characterized the landscape of somatic mutations and pathways in colorectal cancer liver metastasis,” Journal of Oncology, vol. 2019, Article ID 2684075, 8 pages, 2019.

A. Henriksen, “Checkpoint inhibitors in pancreatic cancer,” Cancer Treatment Reviews, vol. 78, 2019.

B.-D. Qin, X.-D. Jiao, and Y.-S. Zang, “Tumor mutation burden to tumor burden ratio and prediction of clinical benefit of anti-PD-1/PD-L1 immunotherapy,” Medical Hypotheses, vol. 116, pp. 111–113, 2018.

S. Lu, J. E. Stein, D. L. Rimm et al., “Comparison of biomarker modalities for predicting response to PD-1/PD-L1 checkpoint blockade,” JAMA Oncology, vol. 5, no. 8, pp. 1195–1204, 2019.

Z. Xu, J. Dai, D. Wang et al., “Assessment of tumor mutation burden calculation from gene panel sequencing data,” OncoTargets and Therapy, vol. 12, pp. 3401–3409, 2019.

K. Hatakeyama, T. Nagashima, K. Urakami et al., “Tumor mutational burden analysis of 2,000 Japanese cancer genomes using whole exome and targeted gene panel sequencing,” Biomedical Research, vol. 39, no. 3, pp. 159–167, 2018.

R. Lim, “Deubiquitinase USP10 regulates notch signaling in the endothelium,” Science, vol. 364, no. 6436, pp. 188–193, 2019.

G. M. Halliday, “SWI/SNF: a chromatin-remodelling complex with a role in carcinogenesis,” International Journal of Biochemistry & Cell Biology, vol. 41, no. 4, p. 728, 2009.

R. Gao, Y. Yu, A. Inou et al., “Heterogeneous nuclear ribonucleoprotein k (hnRNP-k) promotes tumor metastasis by induction of genes involved in extracellular matrix, cell movement, and angiogenesis,” Journal of Biological Chemistry, vol. 288, no. 21, pp. 15046–15056, 2013.