- Research article
- Open Access
Highly multiplexed quantifications of 299 somatic mutations in colorectal cancer patients by automated MALDI-TOF mass spectrometry
BMC Medical Genomics volume 13, Article number: 143 (2020)
Detection of somatic mutations in tumor tissues helps to understand tumor biology and guide treatment selection. Methods such as quantitative PCR can analyze a few mutations with high efficiency, while next generation sequencing (NGS) based methods can analyze hundreds to thousands of mutations. However, there is a lack of cost-effective method for quantitatively analyzing tens to a few hundred mutations of potential biological and clinical significance.
Through a comprehensive database and literature review we selected 299 mutations associated with colorectal cancer. We then designed a highly multiplexed assay panel (8-wells covering 299 mutations in 109 genes) based on an automated MADLI-TOF mass spectrometry (MS) platform. The multiplex panel was tested with a total of 319 freshly frozen tissues and 92 FFPE samples from 229 colorectal cancer patients, with 13 samples also analyzed by a targeted NGS method covering 532 genes.
Multiplex somatic mutation panel based on MALDI-TOF MS detected and quantified at least one somatic mutation in 142 patients, with KRAS, TP53 and APC being the most frequently mutated genes. Extensive validation by both capillary sequencing and targeted NGS demonstrated high accuracy of the multiplex MS assay. Out of 35 mutations tested with plasmid constructs, sensitivities of 5 and 10% mutant allele frequency were achieved for 19 and 16 mutations, respectively.
Automated MALDI-TOF MS offers an efficient and cost-effective platform for highly multiplexed quantitation of 299 somatic mutations, which may be useful in studying the biological and clinical significance of somatic mutations with large numbers of cancer tissues.
The cancer genomes carry many different genomic alterations including somatic mutations. Most notably, driver mutations are causally implicated in oncogenesis. They have conferred a growth advantage on the cancer cells . Increasing number of mutations are being characterized for better treatment selection and prognosis . In non-small cell lung cancer, for example, patients with epidermal growth factor receptor (EGFR) mutation are likely to benefit from tyrosine kinase inhibitors (TKIs) . On the opposite, colorectal cancer patients with KRAS or NRAS mutations should not be treated with either EGFR-inhibitor Cetuximab  or Panitumumab .
Several methods are available for the detection of tumor somatic mutations. Next generation sequencing (NGS) is widely adopted as a discovery tool for somatic mutation analysis in cancer . Using NGS, The Cancer Genome Atlas (TCGA) presents genomic alternations identified from over 11,000 tumor samples from 33 cancer types, providing valuable insights into new targets for drug development, treatment selection, or even combination therapies for personalized medicine . While NGS methods can analyze hundreds to tens of thousands mutations , methods such as qPCR or digital PCR are ideal for analyzing a few clinically important mutations with high efficiency and cost-effectiveness [9, 10]. However, as biologically and clinically important mutations for many cancer types continue to be identified, there is a growing unmet need to cost-effectively and quantitatively analyze tens to a few hundred of mutations, in highly heterogenous cancer tissues where mutant allele frequency can be highly variable.
The MassARRAY platform, based on automated MALDI-TOF mass spectrometry, is suitable to meet this growing need due to its high multiplex capability, flexibility for both sample size and mutation number, quantification capability, and automation in sample processing and data analysis . The MassARRAY can detect up to 40 mutations in a single well on a 96 or 384-well plate.
In this report, we developed and validated an 8-well panel covering 299 most common somatic mutations in colorectal cancer (CRC). Such multiplex level may be sufficient for detecting virtually all clinically relevant mutations in a cancer type, with cost-effectiveness and turnaround time desirable in real-life clinical settings.
Patients diagnosed with colorectal cancer were recruited with informed consent between July 2015 and June 2019 from the First Affiliated Hospital and the Second Affiliated Hospital of Wenzhou Medical University, China. Exclusion criteria include more than two pathological types, metastasis, treated with neoadjuvant chemotherapy or immunotherapy before surgery, or any cancer within the past 5 years. Tumor stages were determined according to the 8th edition of the American Joint Committee on Cancer (AJCC). The study was approved by the ethics committee of Wenzhou Medical University and its affiliated hospital.
Sample collection and DNA extraction
Resections or biopsies of primary solid tumors and adjacent normal tissues located 2-cm away from the tumor tissue, were taken immediately after surgery, snap frozen in liquid nitrogen, and stored at − 80 °C. DNA was extracted from tissue using the QIAamp DNA mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s instruction, and stored at − 20 °C for further use. The DNA concentration was quantified by NanoDrop One (Thermo Fisher Scientific, Foster City, CA, USA).
FFPE tumor samples were analyzed by hematoxylin-and-eosin (H&E) staining to determine tumor purity. For each FFPE tissue, 20 slices of 10-μm-thick sections were used for DNA extraction with the QIAamp DNA FFPE Tissue kit (Qiagen). Concentration of DNA was determined by a Qubit 3.0 Fluorometer (Thermo Fisher Scientific).
CRC mutation panel design
The mutations were selected for either of the following three reasons: 1) present in at least two of the following three sources for CRC samples: TCGA (https://cancergenome.nih.gov/); International Cancer Genome Consortium (ICGC) database (https://icgc.org/); and publication by Muzny DM et al. ; 2) potential resistance to CRC targeted therapy based on My Cancer Genome (http://www.mycancergenome.org), or 3) reported in a commercial cancer panel .
The CRC mutation panel was designed with Assay Designer software (MassARRAY Typer, Version 4.0, Agena Biosciences, San Diego, CA, USA), with primer sequences not overlapping known single nucleotide polymorphisms whenever possible.
CRC mutation detection and quantification by MALDI-TOF MS
CRC Mutation detection was optimized and performed by MassARRAY Analyzer 4 with CPM (Agena Biosciences). All primers were purchased from Integrated DNA Technologies (Integrated DNA Technologies, Coralville, IA, USA). All reagents were purchased from Agena Bioscience unless otherwise specified. Briefly, tissue genomic DNA was amplified by multiplex PCR. Shrimp alkaline phosphatase treatment was performed to inactivate surplus nucleotides. A primer extension reaction (iPLEX Pro) was performed with mass-modified terminator nucleotides, and the product was spotted on SpectroCHIP. Wild-type and mutant alleles were then discriminated by molecular weights determined by MassARRAY analyzer.
Allele calls were performed with MassARRAY Typer Analyzer software (Typer 4.0.26). Additionally, at least two investigators independently reviewed the mass spectra to further confirm the automated calls by the software. To estimate mutant allele frequency, the heights of raw spectral peaks corresponding to the wild-type and mutation allele were quantified. Mutation allele frequency was estimated by calculating mutant peak /(mutant peak + wild type peak).
Targeted NGS and data analysis
Genomic DNA was fragmented to an average size of 200 to 500 bp using Bioruptor Pico sonication device (Diagenode, Denville, NJ, USA). Sequencing libraries were prepared using the KAPA LTP Library Preparation Kit for Illumina (Kapa Biosystems, Wilmington, MA, USA) according to manufacturer’s suggestions. Hybridization-based target enrichment was carried out with xGen Pan-Cancer Panel v2.4 (532 cancer-relevant genes), and xGen Lockdown Hybridization and Wash Reagents Kit (Integrated DNA Technologies). Captured libraries were amplified and purified using Agencourt AMPure XP Beads (Beckman Coulter, Atlanta, GA, USA). Concentrations and qualities of DNA libraries were analyzed by the Agilent High Sensitivity DNA Kit (Agilent Technologies, Santa Clara, CA, USA).
The libraries were paired-end sequenced on the Illumina HiSeq X platform (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions. Sequencing adapters and low quality bases were trimmed from raw sequencing reads using Trim Galore (v0.4.1; https:// github.com/FelixKrueger/TrimGalore). The filtered reads were then mapped to the reference Human Genome (hg38) using BWA-MEM (v0.7.12; https://github.com/lh3/bwa/tree/master/bwakit) with the default settings. The GATK (v188.8.131.52; https://software.broadinstitute.org/gatk/) was used for single nucleotide variation (SNV) identifications. SNV calls with at least 2.5% variant allele frequency (VAF) were retained, followed by annotation using ANNOVAR .
Mutation validation by capillary sequencing
A subset of samples were selected for Sanger sequencing to validate MS results. The capillary sequencing was performed with the BigDye Terminator Cycle Sequencing Kit and ABI 3730 Genetic Analyzer (Thermo Fisher Scientific).
Statistical analyses were performed with IBM SPSS Statistics 20.0. The χ2 test or Fisher exact test were used to compare baseline categorical variables, and the Kruskal-Wallis test was used to analyze the association between mutation and tumor size. The Cox proportional-hazards model of multivariate analysis was performed to analyze covariables, such as age, gender, clinical stage, location, differentiation, metastasis, treatment after surgery, and tumor size.
A total of 300 patients diagnosed with colorectal cancer between July 2015 and June 2019 were recruited. Seventy one patients were excluded (see exclusion criteria in the methods section). Frozen tissues (paired tumor and adjacent normal tissues) from 136 patients, frozen tumor tissues from 47 patients, and formalin fixed paraffin embedded (FFPE) tissues (paired tumor and adjacent normal tissues) from 46 patients with CRC were analyzed. Table 1 summarizes the patients’ baseline characteristics, tumor stage, location, differentiation, size, follow-up treatment after surgery and tumor marker (CEA) level.
Selection of somatic mutations and multiplex assay design
Somatic mutations for CRC were chosen based on their prevalence, and biological/clinical relevance. The final mutation panel consists of 299 mutations from 109 genes (Additional file 1: Supplementary Table 1).
We designed an 8-well multiple assay for a panel of 299 mutations (36-plex in well 1 and 2, 34-plex in well 3 and 5, 30-plex in well 4, 29-plex in well 6, 26-plex in well 7, and 12-plex in well 8) for MALDI-TOF MS analysis. The sequences of PCR primers and extension primers are provided in Additional file 1: Supplementary Table 2. Wild-type and mutant alleles were then discriminated by molecular weights of the extension products measured by MALDI-TOF MS (example results shown in Fig. 1).
To determine the sensitivity of the assays, we constructed plasmids for 35 different mutations. Plasmid mutations were confirmed by Sanger sequencing. Mixtures of plasmids containing wild-type and mutant sequences to mimic samples with 5 and 10% mutation frequencies were prepared to analyze the sensitivity and accuracy of the assay. Among them, 19 assays achieved a 5% sensitivity for mutations, and 16 assays achieved a 10% sensitivity (Fig. 2 and Supplement figure 1).
Mutation profiles of 229 patients from frozen and FFPE tissues
Using the multiplex assay, we analyzed 411 samples from 229 patients (319 freshly frozen tissues from 183 patients and 92 FFPE samples from 46 patients) (Fig. 3). In the frozen tissue cohort, we identified and quantified 52 different somatic mutations in 107 patients. The most frequently mutated genes were KRAS (52 patients, 28.4%), TP53 (48 patients, 26.2%) and APC (32 patients, 17.5%). We chose 20 mutations from 31 patients for Sanger sequencing validation. Sanger Sequencing showed complete concordance with MALDI-TOF MS (Additional file 2: Supplementary Figure 2).
In the cohort with 46 paired FFPE tissues (tumor and adjacent normal tissues), we first performed histopathological analysis with H&E staining (Supplement figure 3). With the multiplex MS assay, we detected a total of 39 mutations in 35 samples. Similarly, APC, TP53, and KRAS were most frequently mutated. We chose 15 different mutations (19 total mutations) from 13 patients for Sanger sequencing validation (Additional file 2: Supplementary Figure 4). Eighteen of the 19 mutations (95%) were validated successfully. Samples chosen for validation, as well as sequences of PCR primers and sequencing primers, are provided in Additional file 1: Supplementary Table 3.
Both frozen and FFPE tissues are commonly used for molecular pathology analysis. Tissue sampling differences may affect the somatic mutation frequencies in individual samples. For the three genes (TP53, APC, and KRAS) that are most commonly mutated in CRC, we found two genes (TP53 and APC) showed higher mutation frequencies in the frozen tissues than the FFPE tissues. The average TP53 mutation frequencies were 50.9 and 38.3% for the frozen tissues and FFPE tissues, respectively (p = 0.006, Mann Whitney U test). The average APC mutation frequencies were 47.5 and 26.8% for the frozen tissues and FFPE tissues, respectively (p = 0.005, Mann Whitney U test).
We next examined the association between the molecular alterations and patient characteristics (Table 2). Among the 142 patients with at least one somatic mutation identified, we found that the 63 patients without a RAS mutation exhibited different clinical stage distribution as compared with the 79 patients with RAS mutation (p = 0.019, χ2 test). TP53 mutations were more frequently found in male patients by univariate analysis (p = 0.002, χ2 test).
Comparison between targeted NGS and MS assays
We chose 13 patient samples and performed targeted NGS sequencing using xGen Pan-Cancer Panel V2.4. The NGS panel covers 532 genes. Deep sequencing was performed to have an average sequencing depth of 1400. Among these 13 patients, the MS method identified 11 somatic mutations. Targeted NGS detected 14 mutations, including the 11 mutations detected by MS. Three mutations were identified only by NGS. These three discordant mutations were further analyzed by DNA cloning and capillary sequencing. Two mutations (KRAS_c.G35A, APC_c.C2626T) were verified by cloning and sequencing. One mutation (KRAS_c.G34A) was present at 3% by targeted NGS. We sequenced 37 clones and did not observe any mutant.
We developed an 8-well, multiplexed assay based on automated MALDI-TOF mass spectrometry to detect 299 CRC related mutations in 229 southern China patients.
The selected 299 mutations cover the National Comprehensive Cancer Network (NCCN) guideline of colon cancer recommended gene list relevant to treatment and prognosis, such as KRAS, NRAS, BRAF V600E mutations [15,16,17]. The MALDI-TOF MS assay is sensitive and quantitative in analyzing somatic mutations, with extensive validation by capillary sequencing, as well as a head-to-head comparison with a targeted NGS panel covering 532 genes. We analyzed both frozen tissues and FFPE samples as both sample types are commonly used for molecular pathology analysis. In our cohorts, we found the frozen tissues may contain higher tumor content as evidenced by higher mutation frequencies for TP53 and APC, which may be due to sampling differences.
NGS-based methods can detect known and unknown mutations while the MALDI-TOF MS assay is more suitable to detect pre-selected, functionally relevant known mutations. The cost of the multiplex MS assay (about USD 60/sample) is substantially lower than targeted NGS (about USD 150 ~ 200/sample). The MS approach is highly automated and suitable for high throughput analysis of larger sample size, making it suitable as a routine testing platform. The overall time needed for the MS assay is about 9 h, with hands-on time of about 60 min. The MS assay is also flexible in panel expansion to add more important mutations.
We have developed and validated a highly multiplexed assay for the quantification of 299 somatic mutations in colorectal cancer tissues, offering a tool for studying the biological and clinical significance of somatic mutations with large numbers of cancer tissues. The multiplex assay may also be useful in clinical management of colorectal cancer patients.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request. All sequencing files are available from the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA662695).
- MALDI-TOF MS:
Matrix-assisted laser desorption/ionization-time of flight mass spectrometry
Polymerase chain reaction
- KRAS :
KRAS proto-oncogene, GTPase
- TP53 :
Tumor protein p53
- APC :
APC, WNT signaling pathway regulator
- EGFR :
Epidermal growth factor receptor
Tyrosine kinase inhibitor
- NRAS :
NRAS proto-oncogene, GTPase
The Cancer Genome Atlas
American Joint Committee on Cancer
International Cancer Genome Consortium
Single Nucleotide Variation
Single Nucleotide Polymorphisms
Variant Allele Frequency
- TCF7L2 :
Transcription factor 7 like 2
National Comprehensive Cancer Network
- BRAF :
B-Raf Proto-Oncogene, Serine/Threonine Kinase
Stratton MR, Campbell PJ, Futreal PA. The cancer genome. Nature. 2009;458(7239):719–24.
Santos R, Ursu O, Gaulton A, Bento AP, Donadi RS, Bologa CG, Karlsson A, Al-Lazikani B, Hersey A, Oprea TI, Overington JP. A comprehensive map of molecular drug targets. Nat Rev Drug Discov. 2017;16(1):19–34.
Soria JC, Ohe Y, Vansteenkiste J, Reungwetwattana T, Chewaskulyong B, Lee KH, Dechaphunkul A, Imamura F, Nogami N, Kurata T, Okamoto I, Zhou C, Cho BC, Cheng Y, Cho EK, Voon PJ, Planchard D, Su WC, Gray JE, Lee SM, Hodge R, Marotti M, Rukazenkov Y, Ramalingam SS, Investigators F. Osimertinib in untreated EGFR-mutated advanced non-small-cell lung Cancer. N Engl J Med. 2018;378(2):113–25.
Lievre A, Bachet JB, Boige V, Cayre A, Le Corre D, Buc E, Ychou M, Bouche O, Landi B, Louvet C, Andre T, Bibeau F, Diebold MD, Rougier P, Ducreux M, Tomasic G, Emile JF, Penault-Llorca F, Laurent-Puig P. KRAS mutations as an independent prognostic factor in patients with advanced colorectal cancer treated with cetuximab. J Clin Oncol. 2008;26(3):374–9.
Amado RG, Wolf M, Peeters M, Van Cutsem E, Siena S, Freeman DJ, Juan T, Sikorski R, Suggs S, Radinsky R, Patterson SD, Chang DD. Wild-type KRAS is required for panitumumab efficacy in patients with metastatic colorectal cancer. J Clin Oncol. 2008;26(10):1626–34.
Metzker ML. Sequencing technologies - the next generation. Nat Rev Genet. 2010;11(1):31–46.
Hutter C, Zenklusen JC. The Cancer Genome Atlas: creating lasting value beyond its data. Cell. 2018;173(2):283–5.
Zehir A, Benayed R, Shah RH, Syed A, Middha S, Kim HR, Srinivasan P, Gao J, Chakravarty D, Devlin SM, Hellmann MD, Barron DA, Schram AM, Hameed M, Dogan S, Ross DS, Hechtman JF, DF DL, Yao J, Mandelker DL, Cheng DT, Chandramohan R, Mohanty AS, Ptashkin RN, Jayakumaran G, Prasad M, Syed MH, Rema AB, Liu ZY, Nafa K, Borsu L, Sadowska J, Casanova J, Bacares R, Kiecka IJ, Razumova A, Son JB, Stewart L, Baldi T, Mullaney KA, Al-Ahmadie H, Vakiani E, Abeshouse AA, Penson AV, Jonsson P, Camacho N, Chang MT, Won HH, Gross BE, Kundra R, Heins ZJ, Chen HW, Phillips S, Zhang H, Wang J, Ochoa A, Wills J, Eubank M, Thomas SB, Gardos SM, Reales DN, Galle J, Durany R, Cambria R, Abida W, Cercek A, Feldman DR, Gounder MM, Hakimi AA, Harding JJ, Iyer G, Janjigian YY, Jordan EJ, Kelly CM, Lowery MA, LGT M, Omuro AM, Raj N, Razavi P, Shoushtari AN, Shukla N, Soumerai TE, Varghese AM, Yaeger R, Coleman J, Bochner B, Riely GJ, Saltz LB, Scher HI, Sabbatini PJ, Robson ME, Klimstra DS, Taylor BS, Baselga J, Schultz N, Hyman DM, Arcila ME, Solit DB, Ladanyi M, Berger MF. Erratum: Mutational landscape of metastatic cancer revealed from prospective clinical sequencing of 10,000 patients. Nat Med. 2017;23(8):1004.
Malapelle U, Sirera R, Jantus-Lewintre E, Reclusa P, Calabuig-Farinas S, Blasco A, Pisapia P, Rolfo C, Camps C. Profile of the Roche cobas(R) EGFR mutation test v2 for non-small cell lung cancer. Expert Rev Mol Diagn. 2017;17(3):209–15.
Milosevic D, Mills JR, Campion MB, Vidal-Folch N, Voss JS, Halling KC, Highsmith WE, Liu MC, Kipp BR, Grebe SKG. Applying standard clinical chemistry assay validation to droplet digital PCR quantitative liquid biopsy testing. Clin Chem. 2018;64(12):1732–42.
Fleitas T, Ibarrola-Villava M, Ribas G, Cervantes A. MassARRAY determination of somatic oncogenic mutations in solid tumors: moving forward to personalized medicine. Cancer Treat Rev. 2016;49:57–64.
Cancer Genome Atlas N. Comprehensive molecular characterization of human colon and rectal cancer. Nature. 2012;487(7407):330–7.
Froyen G, Broekmans A, Hillen F, Pat K, Achten R, Mebis J, Rummens JL, Willemse J, Maes B. Validation and application of a custom-designed targeted next-generation sequencing panel for the diagnostic mutational profiling of solid tumors. PLoS One. 2016;11(4):e0154038.
Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164.
Di Nicolantonio F, Martini M, Molinari F, Sartore-Bianchi A, Arena S, Saletti P, De Dosso S, Mazzucchelli L, Frattini M, Siena S, Bardelli A. Wild-type BRAF is required for response to panitumumab or cetuximab in metastatic colorectal cancer. J Clin Oncol. 2008;26(35):5705–12.
Kopetz S, McDonough SL, Morris VK, Lenz H-J, Magliocco AM, Atreya CE, Diaz LA, Allegra CJ, Wang SE, Lieu CH, Eckhardt SG, Semrad TJ, Kaberle K, Guthrie KA, Hochster HS. Randomized trial of irinotecan and cetuximab with or without vemurafenib in BRAF-mutant metastatic colorectal cancer (SWOG 1406). J Clin Oncol. 2017;35(4_suppl):520.
Van Cutsem E, Kohne CH, Lang I, Folprecht G, Nowacki MP, Cascinu S, Shchepotin I, Maurel J, Cunningham D, Tejpar S, Schlichting M, Zubel A, Celik I, Rougier P, Ciardiello F. Cetuximab plus irinotecan, fluorouracil, and leucovorin as first-line treatment for metastatic colorectal cancer: updated analysis of overall survival according to tumor KRAS and BRAF mutation status. J Clin Oncol. 2011;29(15):2011–9.
We would like to thank all the patients who participated in the study.
We would like to thank the support from the doctors and nurses from Department of Colorectal and Anal Surgery, The First and Second Affiliated Hospital of Wenzhou Medical University. We also thank Mr. Jin Weijiang for his help in the final revision of the manuscript.
This work was supported by the National Natural Sciences Foundation of China (No. 81672922), Medical Health Science and Technology Project of Zhejiang Provincial Health Commission (No. 2017174160), Wenzhou Municipal Science and Technology Bureau (No. Y20160044), Innovation Discipline of Zhejiang Province in Nucleic Acid Molecular Diagnostics (No. 437201702G) and Key Discipline of Zhejiang Province in Medical Technology (First Class, Category A) (No.437601607).
Ethics approval and consent to participate
The study was approved by the Clinical Research Ethics Committee of the First Affiliated Hospital of Wenzhou Medical University and Wenzhou Medical University (ref. no. 2016–161). Written informed consent was obtained from participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Mutation list of CRC panel. Table S2. Primer sequences for mutation detection by CRC panel. Table S3. Primer sequences for mutation validation by Sanger sequencing. (XLS 145 kb)
. MS assays could detect 5 and 10% mutations mixed by wild-type and mutant plasmids. Figure S2. Genomic alterations of frozen tissues were detected by MALDI-TOF MS and confirmed by Sanger sequencing (A-F). Figure S3. Diagnosis and histopathological determination of Formalin Fixation and Paraffin Embedding tumor tissues (A-U) were conducted by analysis of H&E staining. Figure S4. Genomic alterations of FFTE tissues were detected by MALDI-TOF MS and confirmed by Sanger sequencing (A-F).
About this article
Cite this article
Xu, C., Peng, D., Li, J. et al. Highly multiplexed quantifications of 299 somatic mutations in colorectal cancer patients by automated MALDI-TOF mass spectrometry. BMC Med Genomics 13, 143 (2020). https://doi.org/10.1186/s12920-020-00804-y
- Somatic mutation
- Colorectal cancer
- MALDI-TOF mass spectrometry
- Multiplex detection