Skip to main content

Characterization of genomic alterations in Chinese colorectal cancer patients with liver metastases



The exploration of genomic alterations in Chinese colorectal liver metastasis (CRLM) is limited, and corresponding genetic biomarkers for patient’s perioperative management are still lacking. This study aims to understand genome diversification and complexity that developed in CRLM.


A custom-designed IDT capture panel including 620 genes was performed in the Chinese CRLM cohort, which included 396 tumor samples from metastatic liver lesions together with 133 available paired primary tumors.


In this Chinese CRLM cohort, the top-ranked recurrent mutated genes were TP53 (324/396, 82%), APC (302/396, 76%), KRAS (166/396, 42%), SMAD4 (54/396, 14%), FLG (52/396, 13%) and FBXW7 (43/396, 11%). A comparison of CRLM samples derived from left- and right-sided primary lesions confirmed that the difference in survival for patients with different primary tumor sites could be driven by variations in the transforming growth factor β (TGF-β), phosphatidylinositol 3-kinase (PI3K) and RAS signaling pathways. Certain genes had a higher variant rate in samples with metachronous CRLM than in samples with simultaneous metastasis. Overall, the metastasis and primary tumor samples displayed highly consistent genomic alterations, but there were some differences between individually paired metastases and primary tumors, which were mainly caused by copy number variations.


We provide a comprehensive depiction of the genomic alterations in Chinese patients with CRLM, providing a fundamental basis for further personalized therapy applications.


Colorectal cancer (CRC) is one of the most common causes of cancer death worldwide [1], and distant metastases are often observed at diagnosis (synchronous metastases), with the liver being the most frequently affected organ [2]. Even if liver metastasis is not detected initially, there is a considerable probability that it will develop later (metachronous metastases). Perioperative chemotherapy combined with hepatectomy is expected to be the best treatment to cure colorectal liver metastasis (CRLM), but prognosis following resection of CRLM is still poor [3, 4]. There is an urgent need for prognostic and predictive biomarkers that can aid the selection of optimal neoadjuvant regimens and postoperative management strategies.

Many studies regarding the characterization of genomic alterations in metastatic colorectal cancer (mCRC) have been conducted in past years to discover therapeutic biomarkers. Targeted sequencing of mCRC samples at Memorial Sloan Kettering (MSK) Cancer Center revealed that varied survival by tumor laterality can be predicted by certain genes [5]. A large-scale integrated omics study in a Chinese CRC cohort demonstrated the ability of the phosphoproteome to distinguish metastasis and to predict drug response [6]. In research specifically focusing on CRLM, high similarities were observed between primary tumors and liver metastases (LMs), whereas ubiquitous private mutations in LMs suggested individual tumor heterogeneity and specific genetic biomarkers that were capable of predicting potential therapeutics for treating LMs [7,8,9]. However, these studies either handled primary tumors, LMs and other metastases identically or were restricted to quite small numbers of samples. To comprehensively characterize the genomic landscape of Chinese CRLM patients and explore potential prognostic biomarkers, we performed targeted sequencing of LM lesions and available paired primary tumors in a large Chinese CRLM cohort.

Materials and methods

Patients and samples

A total of 529 specimens from 396 patients with CRLM were collected from the Hepato-Pancreato-Biliary Surgery Department I, Peking University Cancer Hospital between January 2015 and October 2020, including 396 tissue samples from metastatic liver lesions and 133 paired primary tumor samples. Acquired samples were subjected to next-generation sequencing using a targeted panel assay as were blood samples obtained from all patients. Microsatellite instability was assessed for each patient. Information about the patient’s age, sex, primary tumor location and times to metastasis was collected by reviewing the medical records, and the data are summarized in Table 1. Detailed clinical data is in Additional file 1: Table S1. The study was reviewed and approved by the Ethics Committee of Peking University Cancer Hospital.

Table 1 Patient characteristics

Targeted sequencing

DNA was extracted from formalin-fixed paraffin-embedded (FFPE) tissues and white blood cells using a QIAamp DNA FFPE Tissue Kit (Qiagen, Hilden, Germany) and a Blood Genomic DNA Mini Kit (Cwbiotech, Beijing, China). A custom-designed IDT capture panel (Integrated DNA Technologies, Coralville, IA) was employed to capture the coding regions of 620 genes (Additional file 2: Table S2). Captured DNA fragments were then used for library preparation and quantification guided by KAPA Hyper Prep protocols (Kapa Biosystems, Wilmington, MA), followed by purification with AMPure XP (Beckman Coulter, Brea, CA) and quantification using a Qubit™ dsDNA HS Assay Kit (Thermo Fisher, Waltham, MA). The final library was sequenced with the a NovoSeq 6000 platform (Illumina, San Diego, CA) with a minimum depth of  ×  500.

Variant identification

The sequencing reads were aligned to a human reference genome (hg19) using Burrows-Wheeler Aligner (BWA) after trimming the adapters with Trimmomatic. Duplicated reads were flagged with Picard, and then reads were realigned using the Genome Analysis Tool Kit (GATK). Mutect2 and HaplotypeCaller were used to identify somatic mutations and germline mutations. FACTERA was used to identify genomic fusions and breakpoints. The filter condition is break support  ≥  5. Mutations were annotated with ANNOVAR. Based on the mutation filtering scheme of MSK-IMPACT (, hotspot mutations with allele depths above 8 and variant frequencies above 2% and non-hotspot mutations with allele depths above 10 and variant frequencies above 5% were retained as somatic mutations. Annotation of pathogenic mutations in ClinVar or truncating mutations of genes related to hereditary colorectal cancer was employed to detect germline mutations. Copy number variations (CNVs) were identified using GATK, and genes with copy number ratios below 0.6 or above 2 were subjected to CNV filtering. Tumor mutational burden (TMB) was defined as the number of nonsynonymous somatic mutations per million bases.

Statistical analysis

Statistical analyses were performed with R v3.6.0. The prevalence of mutants was compared between different groups using the chi-square test. Student’s t test was performed to compare TMB and the percentage of shared mutations between paired primary tumors and LMs.


Mutation landscape of CRLM

In 396 metastatic liver samples, 7327 mutations in 620 genes were identified, including 6038 single nucleotide variants and 1289 small insertions/deletions, 20 germline mutations and 15 fusions (Fig. 1A). The top-ranked recurrently mutated genes observed in this cohort were TP53 (324/396, 82%), APC (302/396, 76%), KRAS (166/396, 42%), SMAD4 (54/396, 14%), FLG (52/396, 13%) and FBXW7 (43/396, 11%). Germline mutations were detected in 20 samples (5.1%), and approximately one-third of the germline mutations (6/20) occurred in the MSH6 gene. Druggable gene fusions were found in 2.8% of all samples (11/396), including 3 ROS1, 3 RET, 2 ALK, 2 NTRK1 and 1 BRAF fusion. Somatic CNVs of 620 genes were also assessed, and fifty-two genes were determined to harbor such genomic events (Fig. 1A). The recurrent CNVs involved prominent amplification of the FLT3, AURKA, MYC, BRCA2, AR and ERBB2 genes and dominant deletion of FGFR1, FGFR3 and PTEN, each of which affected between 2 and 6% of all samples.

Fig. 1
figure 1

The genomic landscape of Chinese CRLM. A Recurrent mutations and CNVs in 396 Chinese patients with CRLM. The number and prevalence of altered genes are indicated on the right side of the heatmap. B Mutation prevalence of genes in Chinese and MSK CRLM cohorts. C Mutation prevalence of genes in treatment-naïve Chinese and MSK CRLM cohorts

Mutation data of 274 LM samples were obtained by taking a subset of the mCRC cohort from MSK Cancer Center, and the data were downloaded from cBioPortal ( To investigate the divergence of genomic alterations between different CRLM populations, the mutation rates were compared between the Chinese and MSK cohorts by assessing the somatic mutations of 337 (Additional file 2: Table S2) genes shared by our panel and MSK-IMPACT. Alterations in PIK3CA and BRAF were significantly enriched in the MSK cohort, while the Chinese cohort had a higher prevalence of mutations in genes like FBXW7, FAT1, ARID1A, ARID1B and NCOR1 (Fig. 1B). Considering that the treatment conditions of the two cohorts were different, the frequencies of genic alterations were further compared considering only treatment-naïve patients (40 in the Chinese cohort and 60 in the MSK cohort). Higher mutation rates in the Chinese CRLM were found in four genes, ARID1B, IRS2, AR and CREBBP (Fig. 1C), and only ARID1B was in the overlap with differentially mutated genes of all samples.

Genomic alterations by primary tumor site

Recent studies suggest that shorter survival of mCRC patients with metastatic tumors originating from right-sided primary tumors is potentially caused by underlying features of somatic alterations [5]. We compared the frequencies of gene mutants in patients with different primary tumor sites. In addition to the previously reported genes KRAS and PIK3CA, other genes (such as AMER1, PIK3R1, PRKDC, ERBB3, DIAPH1, BCLAF1, KDM5A and NOTCH1) were also preferentially altered in LMs from right-sided primary tumors, while enrichment of alterations in TP53 and FBXW7 was observed in LMs originating from the left colon (Fig. 2A). Furthermore, comparative analysis of 10 oncogenic pathways derived from TCGA PanCancer Atlas between right- and left-sided primary sites was also conducted. Two pathways were significantly enriched in right-sided primary tumors, the transforming growth factor β (TGF-β) and phosphatidylinositol 3-kinase (PI3K) signaling pathways (Fig. 2B). Although the difference in mutation rates of RTK-RAS between the two primary site groups did not reach a significant level (82.9% vs. 78.2%, p  =  0.48), RAS showed more mutations in the right-sided primary tumors, while the RTK pathway did not (Additional file 4: Figure S1), and these results are consistent with a previous study [5].

Fig. 2
figure 2

Genomic alterations by primary tumor site. A Gene-level mutation prevalence in liver metastases originating from left- and right-sided primary tumors. B Pathway-level mutation prevalence in liver metastases originating from left- and right-sided primary tumors. Significant elements are labeled as red circles

Genomic alterations according to time to metastasis

It has been reported that synchronous CRLM is correlated with poor survival, but the prognostic role of metastasis timing in CRLM remains controversial [10,11,12]. To clarify this contentious issue from the perspective of genomic alterations, we examined whether there is a differential distribution of genomic alterations between patients with different times to metastasis. Our results showed that mutations in FBXW7, FLT3, XIRP2, TSC2, LATS1 and CREBBP were enriched in metachronous LMs, while only mutations in CDK12 were enriched in synchronous LMs (Fig. 3A). Beyond the associations at the gene level, comparative analysis at the pathway level revealed that the Notch and cell cycle pathways were also selectively altered in metachronous CRLM (Fig. 3B). Additionally, a higher TMB in metachronous LMs than in synchronous LMs was also observed (Additional file 5: Figure S2).

Fig. 3
figure 3

Genomic alterations by the timing of metastasis. A Gene-level mutation prevalence in synchronous and metachronous liver metastases. B Pathway-level mutation prevalence in synchronous and metachronous liver metastases. Significant elements are labeled as red circles

Mutations in primary tumors and liver metastases

The concordance of mutations before and after metastasis was assessed in 133 patients with paired primary tumors and LMs (Additional file 3: Table S3). Shared mutations accounted for 56.7% of all mutations (951/1,685), while 19.7% (330/1,685) were private in primary tumors and 23.6% (399/1,685) were identified only in LMs. The median percentage of shared mutations among all detected genes was 50.0%. This number rose to 72.0% when focusing on the top 10 recurrently mutated genes (Fig. 4A). In particular, the proportions of shared mutations in three genes with the highest mutation rate all exceeded 80% (TP53: 83.5%, APC: 81.9%, and KRAS: 93.0%). No significant difference in shared mutations was observed between left- and right-sided primary tumors or samples from patients with different times to metastasis (Fig. 4B).

Fig. 4
figure 4

Concordances of genomic alterations between paired primary tumors and liver metastases (LMs). A Shared and private variants of recurrently mutated genes in 133 Chinese CRLM patients. The percentage of shared alterations for samples and genes is indicated on the top and right side of the heatmap. B Comparison of shared mutations between different primary tumor sites and the timing of metastasis. C Shared and private CNVs in Chinese CRLM patients

It is noteworthy that the occurrence of shared variants was much lower when only considering CNVs (Fig. 4C). The shared CNVs covered no more than 10% of all CNVs (13/148), and most of them were only in LMs (96/148, 64.9%). Specifically, several genes with only private deletions were observed in the metastases, including FGFR3, MAP2K4 and PTEN.


With extensive research on the genomic characteristics of mCRC, the associations between mutations and the prognosis of patients have gradually been disclosed, and the correlations between primary tumors and LMs have been mined. However, in-depth investigation of genomic alterations of CRLM tumors from Chinese patients has not been performed. In this study, we collected a sizable number of eligible samples and identified comprehensive genomic alterations to explore the differences between patients with different characteristics.

The genomic landscape of Chinese patients with CRLM resembled that previously reported for mCRC [5, 6, 13], and TP53, APC and KRAS were the most frequently mutated genes [14, 15]. Alterations of another recurrently altered gene, FLG, were distributed throughout the entire gene, and most mutations were privately enriched in primary tumors or metastases, indicating that they were passenger mutations that emerged with tumor progression. It has also been confirmed that such variants fail to serve as prognostic indicators for CRC [16]. Published studies have shown that highly frequent mutations in members of the ARID gene family are associated with microsatellite instability [17]. Four out of five microsatellite instability-high patients had such mutations in our cohort. Compared to the MSK cohort, the Chinese cohort showed significantly more mutations of the ARID gene family, regardless of whether all samples or only treatment-naïve patient samples were analyzed. Both genetic and environmental factors can contribute to these disparities, and further validation may be warranted.

Recent studies have confirmed the prognostic role of primary tumor location in mCRC [18, 19] and its association with underlying genomic alterations [5]. In our cohort, the differential distribution of variants in the TGF-β, PI3K and RAS pathways based on primary tumor site was consistent with existing research. Apart from KRAS and PIK3CA, other genes preferentially altered in the right-sided primary tumor sites, such as PRKDC, ERBB3, BCLAF1 and NOTCH1, have been reported to be correlated with poor prognosis or migration in CRC [20,21,22,23,24,25,26]. In our study, alterations in FBXW7 were enriched in LMs derived from left-sided CRCs. Since a previous study found that patients with mutations in FBXW7 may be less likely to respond to anti-EGFR therapy than those without mutations in FBXW7 [12], FBXW7 assessment may be useful to include in molecular testing of patients with left‐sided CRC who would receive targeted therapy.

Compared with left-sided primary tumors, right-sided primary tumors showed a higher proportion of synchronous LMs in our cohort (30.4% vs. 15.7%). Synchronous lesions were once expected to predict worse survival than metachronous lesions in hepatic metastases of CRC, but recent studies have expressed controversial opinions [27]. Interestingly, biomarkers that indicate a worse prognosis for mCRC, such as FLT3 amplification, alterations of FBXW7 and high TMB [28, 29], were found to be enriched in metachronous lesions in our cohort. Consequently, further confirmatory work is required to determine the prognostic value of the timing of metastasis.

Just over half of all genomic alterations were detected in both paired primary tumor and metastasis samples in our cohort. However, when only focusing on recurrent driver genes, this proportion increased to a level that was similar to that seen in several small CRLM cohorts (approximately 70%) [7, 8, 30, 31]. This change in proportion may be because systemic spread of colorectal cancer may occur while the tumor is clinically undetectable [31]. A high percentage of private CNVs in LM was observed in our cohort, which might have been caused by the intratumoral evolution of metastatic clones because metastases were more genetically similar in terms of copy number-based phylogenies [32]. A high CNV burden combined with RAS/BRAFV600E/TP53 formed an efficient prognostic marker [33].

Although this work provides a comprehensive characterization of genomic alterations of liver metastatic tumors in a large Chinese CRLM cohort, there are limitations. This study was performed with samples from a single institution; therefore, the risk of selection bias cannot be completely ruled out. Without available clinical follow-up information, the prognostic effect of differential alterations in subgroups cannot be validated and was only supported by recent literature.


In summary, we delineated a comprehensive genomic landscape of Chinese patients with colorectal liver metastases. Differences in alterations between patient subgroups were characterized, and heterogeneity of individual primary tumor and liver metastasis samples was unveiled. These findings could provide a fundamental basis for further personalized therapy applications.

Availability of data and materials

All data used during the study are available from the corresponding author by request.



Colorectal cancer


Colorectal liver metastasis


Copy number variations


Liver metastases


Formalin-fixed paraffin-embedded


Metastatic colorectal cancer


Memorial Sloan Kettering


Phosphatidylinositol 3-kinase


Transforming growth factor β


Tumor mutational burden


  1. Siegel RL, Miller KD, Goding Sauer A, Fedewa SA, Butterly LF, Anderson JC, et al. Colorectal cancer statistics, 2020. CA Cancer J Clin. 2020;70(3):145–64.

    Article  PubMed  Google Scholar 

  2. McMillan DC, McArdle CS. Epidemiology of colorectal liver metastases. Surg Oncol. 2007;16(1):3–5.

    Article  PubMed  Google Scholar 

  3. Kelly ME, Spolverato G, Le GN, Mavros MN, Doyle F, Pawlik TM, et al. Synchronous colorectal liver metastasis: a network meta-analysis review comparing classical, combined, and liver-first surgical strategies. J Surg Oncol. 2015;111(3):341–51.

    Article  CAS  PubMed  Google Scholar 

  4. Dekker E, Tanis PJ, Vleugels JLA, Kasi PM, Wallace MB. Colorectal cancer. Lancet. 2019;394(10207):1467–80.

    Article  PubMed  Google Scholar 

  5. Yaeger R, Chatila WK, Lipsyc MD, Hechtman JF, Cercek A, Sanchez-Vega F, et al. Clinical sequencing defines the genomic landscape of metastatic colorectal cancer. Cancer Cell. 2018;33(1):125-363 e3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Li C, Sun YD, Yu GY, Cui JR, Lou Z, Zhang H, et al. Integrated omics of metastatic colorectal cancer. Cancer Cell. 2020;38(5):734-47 e9.

    Article  CAS  PubMed  Google Scholar 

  7. Mogensen MB, Rossing M, Ostrup O, Larsen PN, Heiberg Engel PJ, Jorgensen LN, et al. Genomic alterations accompanying tumour evolution in colorectal cancer: tracking the differences between primary tumours and synchronous liver metastases by whole-exome sequencing. BMC Cancer. 2018;18(1):752.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Oga T, Yamashita Y, Soda M, Kojima S, Ueno T, Kawazu M, et al. Genomic profiles of colorectal carcinoma with liver metastases and newly identified fusion genes. Cancer Sci. 2019;110(9):2973–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Lang H, Baumgart J, Heinrich S, Tripke V, Passalaqua M, Maderer A, et al. Extended molecular profiling improves stratification and prediction of survival after resection of colorectal liver metastases. Ann Surg. 2019;270(5):799–805.

    Article  PubMed  Google Scholar 

  10. Cossu ML, Ginesu GC, Feo CF, Fancellu A, Pinna A, Vargiu I, et al. Surgical outcomes in patients with hepatic synchronous and metachronous colorectal metastases. Ann Ital Chir. 2017;88:497–504.

    PubMed  Google Scholar 

  11. Engstrand J, Stromberg C, Nilsson H, Freedman J, Jonas E. Synchronous and metachronous liver metastases in patients with colorectal cancer-towards a clinically relevant definition. World J Surg Oncol. 2019;17(1):228.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Colloca GA, Venturino A, Guarneri D. Different variables predict the outcome of patients with synchronous versus metachronous metastases of colorectal cancer. Clin Transl Oncol. 2020;22(8):1399–406.

    Article  CAS  PubMed  Google Scholar 

  13. Mendelaar PAJ, Smid M, van Riet J, Angus L, Labots M, Steeghs N, et al. Whole genome sequencing of metastatic colorectal cancer reveals prior treatment effects and specific metastasis features. Nat Commun. 2021;12(1):574.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Pietrantonio F, Di Nicolantonio F, Schrock AB, Lee J, Tejpar S, Sartore-Bianchi A, et al. ALK, ROS1, and NTRK rearrangements in metastatic colorectal cancer. J Natl Cancer Inst. 2017.

    Article  PubMed  Google Scholar 

  15. Pagani F, Randon G, Guarini V, Raimondi A, Prisciandaro M, Lobefaro R, et al. The landscape of actionable gene fusions in colorectal cancer. Int J Mol Sci. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Skaaby T, Husemoen LL, Thyssen JP, Meldgaard M, Thuesen BH, Pisinger C, et al. Filaggrin loss-of-function mutations and incident cancer: a population-based study. Br J Dermatol. 2014;171(6):1407–14.

    Article  CAS  PubMed  Google Scholar 

  17. Cajuso T, Hanninen UA, Kondelin J, Gylfe AE, Tanskanen T, Katainen R, et al. Exome sequencing reveals frequent inactivating mutations in ARID1A, ARID1B, ARID2 and ARID4A in microsatellite unstable colorectal cancer. Int J Cancer. 2014;135(3):611–23.

    Article  CAS  PubMed  Google Scholar 

  18. Hugen N, Nagtegaal ID. Distinct metastatic patterns in colorectal cancer patients based on primary tumour location. Eur J Cancer. 2017;75:3–4.

    Article  PubMed  Google Scholar 

  19. Holch JW, Ricard I, Stintzing S, Modest DP, Heinemann V. The relevance of primary tumour location in patients with metastatic colorectal cancer: a meta-analysis of first-line clinical trials. Eur J Cancer. 2017;70:87–98.

    Article  PubMed  Google Scholar 

  20. Seligmann JF, Hatch AJ, Richman SD, Elliott F, Jacobs B, Brown S, et al. Association of tumor HER3 messenger RNA expression with panitumumab efficacy in advanced colorectal cancer. JAMA Oncol. 2018;4(4):564–8.

    Article  PubMed  Google Scholar 

  21. Jackstadt R, van Hooff SR, Leach JD, Cortes-Lavaud X, Lohuis JO, Ridgway RA, et al. Epithelial NOTCH signaling rewires the tumor microenvironment of colorectal cancer to drive poor-prognosis subtypes and metastasis. Cancer Cell. 2019;36(3):319–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Brown GT, Cash B, Alnabulsi A, Samuel LM, Murray GI. The expression and prognostic significance of bcl-2-associated transcription factor 1 in rectal cancer following neoadjuvant therapy. Histopathology. 2016;68(4):556–66.

    Article  PubMed  Google Scholar 

  23. Bogoevska V, Wolters-Eisfeld G, Hofmann BT, El Gammal AT, Mercanoglu B, Gebauer F, et al. HRG/HER2/HER3 signaling promotes AhR-mediated Memo-1 expression and migration in colorectal cancer. Oncogene. 2017;36(17):2394–404.

    Article  CAS  PubMed  Google Scholar 

  24. Li C, Liu X, Liu Y, Liu X, Wang R, Liao J, et al. Keratin 80 promotes migration and invasion of colorectal carcinoma by interacting with PRKDC via activating the AKT pathway. Cell Death Dis. 2018;9(10):1009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Arcaroli JJ, Tai WM, McWilliams R, Bagby S, Blatchford PJ, Varella-Garcia M, et al. A NOTCH1 gene copy number gain is a prognostic indicator of worse survival and a predictive biomarker to a Notch1 targeting antibody in colorectal cancer. Int J Cancer. 2016;138(1):195–205.

    Article  CAS  PubMed  Google Scholar 

  26. Loree JM, Bailey AM, Johnson AM, Yu Y, Wu W, Bristow CA, et al. Molecular landscape of ERBB2/ERBB3 mutated COLORECTAL Cancer. J Natl Cancer Inst. 2018;110(12):1409–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Wisneski AD, Jin C, Huang CY, Warren R, Hirose K, Nakakura EK, et al. Synchronous versus metachronous colorectal liver metastasis yields similar survival in modern era. J Surg Res. 2020;256:476–85.

    Article  PubMed  Google Scholar 

  28. Hasegawa H, Taniguchi H, Nakamura Y, Kato T, Fujii S, Ebi H, et al. FMS-like tyrosine kinase 3 (FLT3) amplification in patients with metastatic colorectal cancer. Cancer Sci. 2020.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Samstein RM, Lee CH, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51(2):202–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Tan IB, Malik S, Ramnarayanan K, McPherson JR, Ho DL, Suzuki Y, et al. High-depth sequencing of over 750 genes supports linear progression of primary tumors and metastases in most patients with liver-limited metastatic colorectal cancer. Genome Biol. 2015;16:32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hu Z, Ding J, Ma Z, Sun R, Seoane JA, Scott Shaffer J, et al. Quantitative evidence for early metastatic seeding in colorectal cancer. Nat Genet. 2019;51(7):1113–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Kim TM, Jung SH, An CH, Lee SH, Baek IP, Kim MS, et al. Subclonal genomic architectures of primary and metastatic colorectal cancer based on intratumoral genetic heterogeneity. Clin Cancer Res. 2015;21(19):4461–72.

    Article  CAS  PubMed  Google Scholar 

  33. Berg KCG, Brunsell TH, Sveen A, Alagaratnam S, Bjornslett M, Hektoen M, et al. Genomic and prognostic heterogeneity among RAS/BRAF(V600E) /TP53 co-mutated resectable colorectal liver metastases. Mol Oncol. 2020.

    Article  Google Scholar 

Download references




This research was supported by Grants (No. 81874143) from the National Natural Science Foundation of China, Municipal Natural Science Foundation of Beijing of China (No. 7192035).

Author information

Authors and Affiliations



Conceived and designed the experiments: CBX, WHW, XLY; performed the data collection: WHW, XLY, LJW, WL, KMJ, QB, JL,KW; analyzed the data: SHZ,CHY,WHW,XLY; contributed reagents/materials/analysis tools: SHZ,CHY,WHW,XLY; contributed to the writing of the manuscript: WHW,XLY; revision of the manuscript: BCX,SHZ,CHY; all authors reviewed the final manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Bao-Cai Xing.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Ethics Committee of Beijing Cancer Hospital and Institute. Informed consent was obtained from individuals for the study.

Consent for publication

Informed consent was obtained from all participants for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

 Detailed clinical data of the patients in the current study

Additional file 2: Table S2.

The gene list of our targeted panel

Additional file 3: Table S3.

 The shared mutations between paired primary tumors and LMs

Additional file 4: Figure S1.

 Comparison of the mutation rates of RTK and RAS between primary tumor site

Additional file 5: Figure S2.

 The association between TMB values and the timing of CRLM diagnosis

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, HW., Yan, XL., Wang, LJ. et al. Characterization of genomic alterations in Chinese colorectal cancer patients with liver metastases. J Transl Med 19, 313 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: