Skip to main content

Advertisement

Towards standardization of next-generation sequencing of FFPE samples for clinical oncology: intrinsic obstacles and possible solutions

Abstract

Background

Next generation sequencing has a potential to revolutionize the management of cancer patients within the framework of precision oncology. Nevertheless, lack of standardization decelerated entering of the technology into the clinical testing space. Here we dissected a number of common problems of NGS diagnostics in oncology and introduced ways they can be resolved.

Methods

DNA was extracted from 26 formalin fixed paraffin embedded (FFPE) specimens and processed with the TrueSeq Amplicon Cancer Panel (Illumina Inc, San Diego, California) targeting 48 cancer-related genes and sequenced in single run. Sequencing data were comparatively analyzed by several bioinformatics pipelines.

Results

Libraries yielded sufficient coverage to detect even low prevalent mutations. We found that the number of FFPE sequence artifacts significantly correlates with pre-normalization concentration of libraries (rank correlation −0.81; p < 1e−10), thus, contributing to sample-specific variant detection cut-offs. Surprisingly, extensive validation of EGFR mutation calls by a combination of aligners and variant callers resulted in identification of two false negatives and one false positive that were due to complexity of underlying genomic change, confirmed by Sanger sequencing. Additionally, the study of the non-EGFR amplicons revealed 33 confirmed unique mutations in 17 genes, with TP53 being the most frequently mutated. Clinical relevance of these finding is discussed.

Conclusions

Reporting of entire mutational spectrum revealed by targeted sequencing is questionable, at least until the clinically-driven guidelines on reporting of somatic mutations are established. The standardization of sequencing protocols, especially their data analysis components, requires assay-, disease-, and, in many cases, even sample-specific customization that could be performed only in cooperation with clinicians.

Background

The advent of next-generation sequencing and other genomic technologies makes it routinely possible today to investigate DNA specimens extracted from human tumors for gene mutations, chromosomal aberrations, differential expression of mRNA and epigenetic alterations. Such detected molecular changes are extensively used for integrative analyses aimed at the identification of key dysregulated pathways and importantly, the establishment of molecular classifiers indispensable for targeted therapy and personalized management of malignant disorders [1].

Due to relatively high prevalence of the mutations in the epidermal growth factor receptor (EGFR) that has led to the development of EGFR tyrosine kinase inhibitors (EGFR-TKIs), non-small cell lung cancer (NSCLC) has become a proving ground for the development of novel approaches for molecular typing [2]. Unfortunately, it has been established that EGFR-responsive patients with activating EGFR mutations eventually develop resistance to the treatment and progress either due to commonly acquired T790M mutations or to a variety of other molecular changes [3]. On the other hand, initial success in exploiting EGFR-TKI sensitivity paved the way for the development of other small molecular inhibitors as well as widespread introduction of the mutation analyses into the molecular subtyping of other tumors. Examples of actionable mutations include RAS (KRAS and NRAS), BRAF and PI3 K in colorectal cancer [4], BRAF or NRAS as well as KIT in malignant melanoma [5], and others. These and other genetic alterations continue to gain importance in companion diagnostics associated with recently marketed targeted therapies and investigational drugs [68]. As a consequence of efforts invested in this area, the number of clinically actionable variations that contribute to sensitivity or resistance of tumors to each of these drugs is dramatically expanding alongside available therapeutic solutions, especially in cases of metastatic and refractory diseases. Moreover, the numerous case reports demonstrating clinical response to the drugs prescribed off-label based on matching molecular evidence ensure the spread of a molecular-informed therapeutic decision-concept into oncologist office routine [912].

This avalanche of changes was precipitated by the replacement of conventional real-time polymerase chain reaction (PCR) or direct sequencing by Sanger with Next Generation Sequencing (NGS) which has recently revolutionized the field of molecular diagnostics. In addition to higher sensitivity and remarkable throughput, NGS allows gaining non-traditional kinds of information about extracted DNA specimen, including the prevalence of shorter fragments [13] and the degree of its sequence heterogeneity within each locus [14]. Nevertheless, a majority of NGS applications are currently limited to research use only, while its implementation in clinical practice awaits rigorous validation requiring establishing Clinical Laboratory Improvement Amendments (CLIA) and College of American Pathologists (CAP)-compliant performance characteristics. Thus, rapid entering of NGS into the space of clinical diagnostics is complicated by a number of loosely defined unknowns, including a coverage required for confident detection, data quality control, benchmarking bioinformatics pipelines currently producing discordant variant calls and overall validation of the robustness of NGS techniques.

Besides the technical obstacles outlined above, there are basic issues with existing “gold standard”: a number of novel, potentially very important insights into the tumorigenesis were gained by massive parallel sequencing for the first time, as they originated from observations that cannot be performed earlier due to intrinsic shortcomings of Sanger sequencing and conventional PCR. As example we can point to the clinical significance of intratumor heterogeneity with therapy resistance mutations being present with low prevalence already at the pretreatment assessment [15]. Currently, there is no consensus on how to differentiate these mutations from the FFPE artifacts [16]. Another set of important issues are incidental findings that may be revealed by NGS reads even if they were targeted for detecting of the mutation spectrum defined beforehand [17]. Examples of such incidental findings may include novel, not yet characterize mutations within target gene, common variants that may influence risks of other diseases or even germline mutations [18]. Such findings are generally hard to interpret in terms of their clinical relevance and lack a consensus on whether they should be communicated to physician or the patient at all. Finally, NGS allows one to detect mutations, which remain unseen if not being targeted by analysis. For instance, owing to implementation of the NGS into clinical practice, EGFR kinase domain duplication mutation has recently emerged as novel EGFR TKI sensitizing variant [19].

Despite recent publications of a number of clinical NGS guidelines clinical test development by a variety of organizations including American College of Medical Genetics [20] and College of American Pathologist [21], overcoming the obstacles mentioned above remains challenging, especially when the detected mutations are somatic, as the majority of guidelines are established only for the reporting of germline mutations. It is likely that both general recommendations and relevant protocols for sequencing, data analysis and clinical interpretation will continue to evolve in the process of sharing the outcomes of the tailoring of commercially available tests to the needs of particular practice or country, or panel customization. In this paper we describe particular obstacles we encountered while analyzing the clinical NGS dataset obtained using the TruSeq Amplicon—Cancer Panel (TSACP) and possible solutions to the problems presented.

Methods

Sample collection

Thirteen archived clinical tumor specimens from twelve lung cancer patients treated at Blokhin Russian Cancer Research Centre (RCRC) in 2014–2015 were randomly selected from respective existing registry. Another set of 13 samples was retrospectively randomly selected from a collection of Russian Scientific Center of Roentgenology and Radiology (RSCRR). Specimens of the latter set had already been screened for the presence of EGFR mutations; hence, this RSCRR set of samples was enriched with EGFR positive patients by design. Overall, we studied 26 tumor specimens including lung adenocarcinoma (n = 11, 42%), squamous cell carcinoma (n = 11, 42%) and the tumors of mixed or unknown histology (n = 4, 16%). The study protocol was approved by Atlas Biomed Internal Review Board. Written informed consent was provided by all patients at inception of the study, all analyses were based on archival data and stored in database with no connections to the patient identifiers.

DNA extraction and sample quality control

Genomic DNA was extracted from formalin fixed, paraffin embedded (FFPE) tissues with the QIAamp DNA FFPE kit (Qiagen, USA) according to the manufacturer’s instructions and was eluted in a 25 μL volume. The extracted DNA specimens were further quantified using the Qubit dsDNA HS assay kit (Life Technologies/Fisher Scientific, USA) and samples with DNA concentration lower 50 ng/μl were concentrated to this value. DNA quality and quantity were further assessed using the Illumina FFPE QC Kit according to the manufacturer’s instructions.

Library preparation and quality control

Sequencing libraries were prepared with the TSACP (Illumina, San Diego, USA), according to manufacturer’s protocol. Briefly, an oligo pool was hybridized to each genomic DNA sample. Following the removal of unbound oligos, target regions of interest flanked by sequences required for amplification were generated by extension and ligation and libraries were further PCR amplified. Library quality was assessed on a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, California). Prior sequencing, the libraries were normalized following the manufacturer protocol and equal volumes were pooled to generate the final sequencing library. Primer panel designed to generate 212 amplicons within 48 cancer-related genes: ABL1, AKT1, ALK, APC, ATM, BRAF, CDH1, CDKN2A, CSF1R, CTNNB1, EGFR, ERBB2, ERBB4, FBXW7, FGFR1, FGFR2, FGFR3, FLT3, GNA11, GNAQ, GNAS, HNF1A, HRAS, IDH1, JAK2, JAK3, KDR, KIT, KRAS, MET, MLH1, MPL, NOTCH1, NPM1, NRAS, PDGFRA, PIK3CA, PTEN, PTPN11, RB1, RET, SMAD4, SMARCB1, SMO, SRC, STK11, TP53, and VHL.

Sequencing and data analysis

Pooled libraries were sequenced using MiSeqDx (llumina) with a 2 × 150 paired-end sequencing design. Image processing and fastq file generation were further performed with CASAVA version 1.8.2 and RTA version 1.17.28 (Illumina). Reads were preprocessed for a removal of low-quality and too short nucleotide sequences using the Prinseq-lite program [22]. Minimum mean read quality score was set to Q30, and minimum length to 75 base pairs. Remaining paired-end reads were mapped to the GRCh37.p13 human genome employing Bowtie2 [23] software with varying parameters. After alignment, an exclusion of primers was performed employing in-house software. SAMtools (version 1.2) [24] was applied for the calling of germline variants, while somatic mutations were identified by Strelka (version 1.014) [25], Varscan (version 2.3.9) [26], ScalPel (version 0.5.3) [27] and Illumina Somatic Variant Caller (SVC) (version 3.1.6.4). GATK version 3.6 [28] was additionally used for indel realignment and other analysis. Germline polymorphisms were discriminated from the somatic based on their frequencies in human populations and presence in dbSNP [29] and COSMIC [30] databases. Recurrent artifact variant calls were discarded from the analysis employing in-house software. Protein variant annotation was performed using ANNOVAR [31], with non-coding or synonymous mutations discarded from further analysis. Copy number variations were detected using CNVPanelizer [32].

Mutation verification

Following NGS sequencing, EGFR and KRAS (including codons 12, 13) mutations in samples from the RCRC set as well as KRAS mutations from RSCRR set were validated either by Sanger sequencing or Real-Time PCR. Samples from RSCRR set were pre-screened for EGFR mutations and, thus, not required further validation with orthogonal methods. EGFR mutations (including exons 18–21) were examined by Sanger Sequencing. Primer pairs were 5′-CTGAGGTGACCCTTGTCTCTG-3′ and 5′-CCAAACACTCAGTGAAAC-3′, 5′-TGCCAGTTAACGTCTTCCTT-3′ and 5′-CAGGGTCTAGAGCAGAGCAG-3′, 5′-CATTCATGCGTCTTCACCTG-3′ and 5′-TTATCTCCCCTCCCCGTATC-3′, and 5′-TGATCTGTCCCTCACAGCAG-3′ and 5′-GGCTGACCTAAAGCCACCTC-3′ for exons 18, 19, 20 and 21 respectively. Thermal cycling conditions included 5 min at 95 °C followed by 35 cycles of 95 °C for 30 s, 60 °C for 30 s, 72 °C for 1 min and one cycle of 72 °C for 7 min. The PCR products were further purified with USB ExoSapit (GE Healthcare, Uppsala, Sweden) followed by cycle sequencing with the BigDye Terminator version 3.1 cycle sequencing kit (Applied Biosystems) according to the manufacturer’s protocol and resolved on an ABI 3500xL sequencer (Applied Biosystems). Sequence chromatograms were analyzed by Sequencher software (Gene Codes Corp, Ann Arbor, MI), followed by manual review. KRAS codon 12, 13 mutations were examined by pyrosequencing as described previously [33].

Results

NGS sequencing results

We performed high-throughput sequencing employing amplicon library construction with primer panel targeting 212 regions of 48 oncogenes and tumor suppressor genes in 26 FFPE tumor samples. Approximately 91% of all reads from the Illumina MiSeq sequencer were successfully mapped to the reference genome, with the coverage depth of 720× for 82% of the target bases and 220× for 93% of the target bases. Across the samples, the overall mean coverage was 2084× and median—2016×. The frequency of amplicon drop-out was at 0.5%. These results indicate that a coverage resolution was high enough to identify somatic point mutations and short indels as well as the copy number variations.

In contrast to high coverage observed across entire panel, one or more amplicons within the genes CDKN2A, FGFR3, GNA11, HRAS, MPL, SRC, STK11, VHL and SMO displayed relatively low depth of the coverage in substantial proportion of DNA specimens. For many samples, amplicon-specific coverages were seen to drop as low as 140×. An extremely low coverage at less than 4% of the average amounts of reads across all amplicons and samples with the mean of 51× was detected in MPL gene, which, to the best of our knowledge, was not previously noted as deviating from average. Additionally, the mining of existing literature highlights that the sequencing of some regions within RB1, HNF1A, NOTCH1 and RET might also result in consistently low coverage [34, 35]. In our study, uniformly high coverage of these genes was observed across all samples. It is possible that amplicon-specific coverage may vary form run to run, or depend on the batch of the primers or the library kit, or clinic-specific protocol of FFPE preparation.

For the mutations informing the treatment strategy specifically in the lung cancer, the coverage was consistently high, with the lower border of the 99.9% CI of coverage depth 1.5 times higher than the average across all the genes comprising TCASP panel. This quality of the coverage allowed us to proceed with other types of analysis aimed at highlighting problematic aspects of NGS testing in clinical oncology.

DNA degradation artifacts may limit the accuracy of the test

The quantities and the quality of DNA extracted from tumor biopsies might vary from sample to sample. In particular, the tissue fixation process might result in extensive DNA degradation and, therefore, to introduction of artifactual bases. According to previous observations made using FFPE samples, false positive findings are predominantly represented by C:G > T:A base substitutions with frequencies in range of 0–10% of the covering reads [36, 37]. Typical assessment by Sanger sequencing is insensitive to the variants present at frequencies lower than 20% [38]. In contrast, the diagnostic with NGS overcomes this obstacle by increase in the coverage proportional to the rarity of the variant. Given that clinical guidelines recommend detection of the variants that are present in 5 or 2.5% of the covering reads for mutations associated with the response and the resistance, respectively, poses the problem of discrimination between verifiable mutations and artifacts of the FFPE processing.

Across all samples studied, we identified a bit over 3300 mutations with the mean C > T base substitution proportion of 44%. Among these, approximately 3200 mutations were detected with the allele frequency ranges of 0–10%. Yet, four out of 26 FFPE specimens displayed high mutation rate with average allele frequencies of more than 10%, with eight, seventeen, forty and fifty-six of these highly prevalent mutations detected in same DNA sample, respectively (Fig. 1a). Importantly, a majority of these highly frequent mutations were C > T substitutions. In four specimens of interest, they were at 34, 50, 63 and 65%, respectively. In each sample with aberrant count of high prevalent substitution, prevalence of the C > T substitutions in allele frequency range from 0 to 10% tended to be lower as compared with the other samples (with p value ranging from 0.015 to 2e−7), thus, demonstrating a peculiar trend.

Fig. 1
figure1

The description of potential artifacts detected in FFPE samples. a Mutation counts (solid line) and C/T base substitution rates (bar histogram) at two allele frequency ranges across all samples were sorted according to the total mutation counts for given specimen. There was a trend of increasing C/T substitution rate along with total mutation count. Four samples did not follow this trend. Only these samples and none of the others showed very high mutation counts at allele frequencies of more than 10% with the accordingly high C/T substitution rate in this range. Pre-normalization library concentrations were low for all four samples in question. b Per-base Phred quality scores for each position across all samples. Four samples with high mutation counts at allele frequencies of more than 10% are highlighted by bold lines. c Overall sequencing results for all samples. Four samples of interest are respectively highlighted

In addition, for three out of four samples with high mutation loads, the coverage was significantly reduced due to low total read counts (p value ranging from 0.026 to 5e−3) (Fig. 1c). Similarly, the percentages of reads passing base quality filters for these three samples were also smaller (Fig. 1b). For one out of four samples with high mutation loads, the coverage was trending along with the majority of samples (p-value = 0.61).

An insight on initial DNA quality may be derived by assessing sequencing library concentrations before the normalization. In specimens of interest, these concentrations were at 1.0, 1.0, 1.3 and 1.3 ng/μl, while for the rest the lowest sequencing library concentration was at 1.9 ng/μl (p < 1e−5, by Wilcoxon-Mann–Whitney Test). Furthermore, the DNA specimens converted to the libraries with high pre-normalization concentrations displayed lowest total mutation counts (Fig. 2). It seems that pre-normalization concentrations of libraries are indicative of the degree of DNA degradation in a given specimen and could be useful for calculations of sample-specific expected rates of artifactual findings.

Fig. 2
figure2

The mutations counts correlate with the pre-normalization library concentrations. Samples with highest pre-normalization library concentrations were found to harbor fewer mutations, most probably due to lower amounts of DNA degradation artifacts, while samples with lowest pre-normalization library concentrations harbored a diverse spectrum of the detectable variants (Spearman’s rank correlation −0.81; p < 1e−10)

DNA specimens may yield multiple somatic mutations associated with potentially important incidental findings

Owing to advances in the development of TKIs that target EGFR, the screening for point mutations within EGFR already became an integral part of routine examination of lung tumor specimens [39]. However, recent studies indicated that the mutations in the genes of MAPK and PI3K-AKT pathways may influence the response to EGFR TKI and, therefore, are useful as both predictive and prognostic markers. These observations prompted us to perform retrospective validation of the presence of NGS-detected somatic variations in PIK3CA, BRAF, KRAS, NRAS, HER2, ALK and MET genes employing Sanger sequencing or quantitative Real-Time PCR.

In a study of 26 FFPE tumor specimens, activating EGFR mutations associated with sensitivity to the EGFR TKI therapy were detected in eight cases. Among these, three tumor specimens were found to harbor coding frame point mutation p.L858R (or p.Leu858Arg), four tumors had the deletion of exon 19 and a single sample had an insertion of exon 19. The EGFR TKI resistance associated mutation p.T790M (or p.Thr790Met) was detected in a single sample with p.L8585R variant of EGFR.

As TSACP panel covers 48 oncogenes and tumor suppressor genes, we were able to extend the analysis of 26 FFPE lung cancer specimens beyond EGFR gene. Assuming that non-EGFR mutations gain their therapeutic significance only when present in a relatively major clone, the allele frequency cutoffs were redefined at 10%. After filtering out all non-recurrent mutations detected exclusively in samples with low pre-normalization library concentrations, we identified 33 unique non-EGFR point mutations including 25 single nucleotide variations and 8 small insertions/deletions (Table 1). Of these, 28 variants were likely of somatic origin, and further 5 were likely germline. In each particular case, the determination was made after taking into account the presence of the mutation in dbSNP and COSMIC databases, and the reported prevalence in human populations or in cancer specimens. The highest rate of somatic mutations was observed in TP53 gene with eleven alterations detected. Overall, previously described non-EGFR mutations known to modify the response to drug therapy were identified in 8 patients, and including the mutations in TP53, KRAS, PIK3CA, HRAS, AKT1 and CTNNB1. Of particular interest, the mutations of KRAS were detected in 5 samples including one sample with simultaneous presence of two KRAS mutations (p.Gly12Cys and p.Ala146Thr), though only two of the KRAS positive samples harbored high prevalent mutations, while the rest, thus, could not be validated with orthogonal methods. Additionally, PIK3CA mutation p.Glu545Lys was detected in single case.

Table 1 Identified exonic mutations in non-EGFR genes

Standard analytic pipelines as a source of possible false positive and false negative findings

Importantly, two Sanger-confirmed EGFR mutations were identified by NGS sequencing, but missed by SVC that is built into the MiSeq Reporter v1.3 + software available in Ilumina maintained BaseSpace™. In particular, inframe deletion in exon 19 that was compounded with single nucleotide variant, a complex mutation p.Glu746_Ser752delinsAlaPhe was mislabeled as two frameshift mutations p.Gly746fs and p.Ser752fs. Though these two molecular changes were found within the same haplotype and, therefore, should be designed as complex mutation, SVC provided no information on this haplotype. Since SVC-reported frameshift deletion lead to total inactivation of the protein rather than activated EGFR that is generated by Sanger-confirmed inframe deletion in exon 19, incorrected calling resulted in a false-negative result. In clinical setting, this false positive would lead to missed opportunity for EGFR TKI treatment. On the other hand, this mutation was successfully detected using Bowtie2 in conjunction with Strelka, Varscan2 and Scalpel pipelines.

Another example of the false negative result was a previously described [52] insertion in exon 19 (p.Ile744_Lys745insLysIleProValAlaIle) that is located 138 bp away from the end of the amplicon. Therefore, at final cycles, reverse reads end within the insertion sequence. Moreover, 3′ end of insertion sequence is similar to the reference sequence. Thus, when BWA [53] or Bowtie2 are run with the default parameters, the reverse reads got misaligned and incorrectly confirm another, previously undescribed shorter insertion. The presence of two insertions in both strand was not reported due to strand bias, and, as a consequence, this mutation was not reported by any tested variant caller, including Strelka, Varscan2 and Scalpel and also Illumina Somatic Variant Caller. However, after adjusting alignment gap open penalty to high and gap extension penalty to low, which resulted in proper alignment, this mutation was successfully identified. Interestingly, realignment with GATK and its standard parameters also allowed successful detection of this mutation, indicating that this aligner should be recommended either as a primary alignment tool or as a backup for independent validation of the findings.

In one patient, the sequencing of EGFR amplicons resulted in calling of both p.G719V (or p.Gly719Val) mutation and p.Leu718fs deletion that were present within the same copy of the gene. According to NGS results, all the sequencing reads with mutation p.G719V also had mutation p.Leu718fs, while all other reads from the same amplicon were mutation-free. While the mutation in codon 719 confers the response to EGFR TKI, the frameshift deletion in the preceding codon results in altered protein sequence and premature stop. Therefore, even though p.G719V mutation is present in the genome, it does not express at the level of the protein, and, therefore, the therapy with EGFR TKIs would not be effective. Thus, the detection of p.G719V is, indeed, a false positive, despite its proper identification. Manual analysis of the reads allowed the phasing of p.Leu718fs into the same haplotype, and the correct assessment of the combined effect of both molecular changes.

The Detection of Somatic Copy Number Variations

Even in absence of the sequenced control samples, the Copy Number Variations (CNVs) may be identified by employing bootstrapping that generates control sets by randomly subsampling experimental ones with replacement. Using this approach, in two specimens an amplification of EGFR was detected, with respective increase in its copy number by 5.7 and 17.6 folds. In both cases, this molecular event was detected across all 8 TSACP amplicons. Though EGFR amplification does not yield consistent predictive power as previously been shown [5456], it still may serve as important prognostic marker [57]. In another patient, a significant increase in MET gene coverage across its five TSACP amplicons was detected. Typically, the amplification of MET is defined by Fluorescent In Situ Hybridization (FISH) as MET:CEP7 ratio that is greater than 2.2. This and higher degrees of amplifications are associated with either primary or acquired resistance to the therapy with EGFR TKIs as well as with enhanced tumorigenesis, invasion and metastasis [5861]. Moreover, high levels of MET amplification defined as MET:CEP7 > 5 were previously identified as actionable driver mutations predicting sensitivity of the tumor to the inhibitors of kinase MET, including crizotinib, currently approved for the treatment of ALK or ROS1-rearranged lung cancers [62]. To investigate the levels of MET amplification in specimen of interest, we normalized the coverage by combining amplicons covering three other genes located on chromosome 7, EGFR, BRAF and SMO. The level of presumable amplification that we observed were low, at 1.65× before normalization and at 1.6× after normalization, thus, indicating that that, in specimen of study, amplification of MET had not reached clinical significance [63].

Discussion

Advances in target therapy development have introduced an opportunity of informed management of the malignant disorders using the guidance of their molecular profiles. Most conveniently, NGS techniques are capable of the detection of multiple genetic alterations at once, while retaining both an accuracy, informativeness and cost-efficiency. Despite the advantages that NGS techniques provides over conventional methods, their implementation often meets difficulties and should be accompanied equivalency of its analytical performance. For numerous targeted panels, a number of published NGS protocols have already demonstrated satisfactory results. Nevertheless, some of the obstacles remain.

In this study, TSACP technology was tested in the retrospective analysis of somatic mutations in 26 FFPE specimens collected from NSCLC patients. All hotspot alterations detected in clinically relevant genes, including EGFR, PIK3CA, BRAF, KRAS, NRAS, ERBB2, ALK and MET, were successfully validated by Sanger sequencing or Real-Time PCR. EGFR mutations were detected in eight patients out of 25, while high-prevalent mutations of KRAS were detected in two, and low prevalent in additional three samples. Some findings were less than typical. For example, one of the specimens harbored both activating EGFR p.L858R mutation and PIK3CA p.Glu545Lys mutation that was present at allele frequency of 15%. Another specimen with the deletion of EGFR exon 19 also harbored KRAS p.Gly12Val mutation that was detected at mutant allele frequency of 6%. Though cases of concomitant presence of EGFR and KRAS or PIK3CA mutations have already been described, they are very rare. In both cases, compound mutations were below the limits for the detection by Sanger sequencing [38]. Importantly, the shortcomings of Sanger sequencing and other conventional methods influence historic accumulation of the data on cancer related variants with lower prevalence and, therefore, provide a ground for underrating of these, potentially important, findings. Widespread use of high-sensitive mutation detection techniques should be accompanied by appropriate data analysis and their deposition into panel-specific databases that may facilitate further insights on co-occurrences of various mutations, their allelic frequencies and clinical significance of each finding. The latter aspect remains the most challenging due to common underreporting of the minor alleles.

Clear clinical need of identification of the low-prevalent mutations remains, at least in part, unmet due to use of FFPE for the tumor specimen long-term preservation. FFPE is convenient, cost-effective, and efficient for immunohistochemical staining and morphology analyses. In the meantime, it is unsuitable for the high resolution mutation analysis due to DNA damage artifacts introduced by fixation and paraffin-embedding procedures. In this study, the tissues were fixed using uniform protocol implemented at two independent laboratories. These FFPE samples are routinely produced by clinical workflow.

With the average coverage of 2084×, at least in theory, we were capable to detect mutations at allele frequencies of 1%. According to TCGA data, both lung adenocarcinoma and squamous cell lung cancer specimens harbor approximately 9 mutations per Mb, or, on average, 0.36 mutations per 40 kb, or the size of the TSACP panel [64, 65]. Considering that the sequencing efforts are concentrated in areas intrinsically important for tumorigenesis and, therefore, undergoing the selective pressure, one may expect the mutation rates a bit higher. In our experiment, we detected a bit over 500 variants per 40 kb of sequenced DNA in each sample, on average. In both histological types of lung cancer samples described in TCGA, the proportion of C > T substitutions was at approximately 20%, while in experimentally assessed DNA specimens extracted from FFPE blocks the percentage of C > T base substitutions was substantially higher (p < 0.01 for both lung adenocarcinoma or squamous cell carcinoma, by t-test), in one sample reaching 74%. Moreover, a majority of detected variants were present at frequencies of 0–10% (Fig. 1). Hence, we have to conclude that these changes are more likely to represent DNA degradation artifacts rather than to be a consequence of intratumor heterogeneity [66].

As it was already shown previously, in a majority of FFPE specimens, DNA degradation artifacts are seen in the allele frequency range of 0–10% and represented predominantly by C:G > T:A nucleotide substitutions. It is tempting to dismiss these artifacts by applying simple allele frequency threshold. Nevertheless, in our study we observed that some of the specimens produced a number of artifactual findings presenting as high frequency variants. We noted that the processing of these specimens resulted in remarkably low pre-normalization concentrations of libraries. One of these specimens harbored hotspot AKT1 p.E17K (or p.Glu17Lys) variant at allele frequency of 4%, which significance was overshadowed by highly prevalent artifacts. Therefore, limiting the analysis only to known hotspots while assuming that probability of artifactual detection of the molecular lesion within known hotspot is low would not resolve the reporting dilemma.

In attempt to determine the applicability of using simple mutation allele frequency cut-offs for sorting out FFPE derived artefacts, the relationship of allele frequencies and mutation counts at these frequencies was studied. Figure 3 illustrates that the DNA specimens that resulted in the production of libraries with high pre-normalization concentrations demonstrate a steady decrease in respective mutation counts along with growing frequencies of individual alleles. In contrast, the sequencing of DNA specimens that resulted in the production of libraries with low pre-normalization concentrations identified subset of called variants, with allele frequencies ranging between 9 and 20% (p < 4e−3). This subset could not be clearly filtered by recommended reportable allelic frequency ranges for an analysis of FFPE derived DNA specimens. Thus, it is more difficult to deal with, as it could be identified only when taking into account the range of pre-normalization concentrations for libraries entered into the sequencing run.

Fig. 3
figure3

The counts of mutations detected in various ranges of allele frequencies. Samples with pre-normalization library concentrations as low as 1.5 ng/μl were shown to yield substantial amounts of presumably artifactual variants with high prevalence in the sequencing reads. These mutations could not be filtered out by using preset frequency thresholds

In this situation, an assessment of pre-normalization concentrations of libraries may be part for the sample-specific calculation of the detection cut-offs. Moreover, these cut-offs may be even mutation-specific: if the frequency of allele in question is below 10%, but its nature is different from C > T, the probabilities that the variant originated in damaged DNA base would be relatively low. On the other hand, there is some evidence that the DNA degradation artifacts are not exclusively represented by substitutions C:G > T:A. For example, when we applied Fisher’s exact test to mutation counts for all nucleotide substitutions detected in our samples and their frequency ranges, we observed that A > T, T > C, G > T and G > C substitutions may also be associated with increased probability of being FFPE derived artefacts, while for G > C and T > C substitutions these probabilities were significantly decreased (data not shown).

This means that instead of simple empirical rules based on allele frequency cut-offs or nucleotide substitution types, this task may require supervised machine learning approaches executed on paired tumor samples obtained from the same patient and validated on the similar datasets. In any case, the necessity of further research in this field is clear since formalin fixation remains the most common technique for long-term preservation of cancer specimens. An alternative is in introduction of additional library preparation steps that may reduce the counts of FFPE-related artifacts, for example, by treating samples with uracil-DNA glycosylase [16]. It should be also noted that an introduction of additional sample processing steps may influence sequencing outputs, the interpretation of these outputs, and, therefore, clinical decision. In particular, the treatment of samples with uracil-DNA glycosylase may significantly impact mutation allele frequencies [67], which, in turn, are paramount for evaluating potential benefits from the targeted therapy [68]. Therefore, an introduction of additional enzymatic processing requires a comparative study and a through validation before its adoption into the routine.

Despite relatively small size of the study, a single false positive and two false negative calls of the EGFR mutations were uncovered. Importantly, using of two or more combinations of aligning software and further indel realignment with variant callers may help to solve these discrepancies. It is important to note that, for the specimens of question, using of unified and, therefore, simplified analysis protocols would definitely result in non-optimal management of the disease. Aiming at eventual replacement of conventional molecular diagnostics with high-throughput NGS-based methods of NGS, we should remember that commonly used analytic pipelines may remain inefficient when dealing with insertions and deletions, thus, justifying the need for in-house pipelines and scripts.

In addition to the mutation in known hotspots, a total of 24 unique somatic mutations were detected in 26 studied samples. TP53 was the most mutated gene with a total of 11 alterations detected in 11 patients. Assuming that these non-hotspot mutations were never prospectively validated in NSCLC, and, therefore, never received a designation of actionable item, all of these would have to be classified as incidental findings. Among these were the molecular changes validated as actionable in other types of the tumors, or the mutations that influence the prognosis but not yet targeted by approved medicines and non-hotspot mutations associated with treatment resistance. It seems that the validation of these mutations as actionable may be hastened by deposition of these, not-yet-actionable findings into panel-specific databases that may mined for further insights into their clinical significance.

In addition, a total of five presumably germline coding frame variants were identified, including KDR p.Q472H (or p.Gln472His), KIT p.M541L (or p.Met541Leu), TP53 p.P72R (or p.Pro72Arg), HNF1A p.G226A (or p.Gly226Ala) and ATM p.F858L (or p.Phe858Leu). For a majority of these variants, the frequencies in human populations were high, thus indicating that these variant were likely inherited form the parent. Nevertheless, it may be possible that the presence of one or more of these variants in patient’ genome may have an impact on tumorigenesis or on the therapy outcomes. For instance, the presence of KDR p.Q472H variant was previously shown to alter tumor angiogenesis and vascularization [69]. Moreover, in vitro studies showed that the sensitivity to variant-harboring cells to VEGFR2 inhibition is higher than that in the cells with wild-type genotype [69, 70]. As compared to the most common TP53-72P variant, the presence of TP53-72R protein augments apoptosis up to 15-folds, thus, to some degree protecting individual from neoplastic development and, possibly, modifying the response to conventional chemotherapy [71]. Though prognostic value of this polymorphism have already been shown in several case series and prospective studies in several tumor types [72, 73], its predictive role remains not well understood [7476]. Furthermore, KIT p.M541L variant was reported to confer an enhanced proliferative response to low levels of stem cell factor, though evidences supporting its predictive effect are controversial [7783]. Finally, despite relatively low frequency of the ATM p.F858L variant in human populations described within the 1000 Genomes Project (0.5%), according to previous publications, its nature is likely germline. Interestingly, in large population based cohorts, this variant was associated with an increase in risks for prostate [84], breast [85] and colorectal [86] carcinomas as well as chronic lymphocytic leukemia [87]. Additionally, in patients with childhood T-lineage acute lymphoblastic leukemia, carrying ATM p.F858L was associated with a variety of negative predictors, and worsened outcomes [88].

It is clear that accumulation of the knowledge on the prevalence of these variants in tumor specimens may lead to eventual recognition of these incidental findings as relevant to either personalized selection of therapeutic strategy, or to the risks of the cancer development in the relatives of the proband. No database currently accumulates these kinds of findings in depersonalized form. Therefore, in a majority of the cases, these germline variants remain unreported and, therefore, missed for further analysis. In our opinion, creation of panel-specific incidental finding databases is warranted, as it may hasten overall understanding of tumorigenesis.

Conclusion

Reporting of entire mutational spectrum revealed by targeted sequencing is questionable, at least until the clinically-driven guidelines on reporting of somatic mutations are established. The need for the development of panel-specific databases allowing analysis of co-occurrence and relevance of somatic mutations, CNVs and germline variants in de-identified form is evident. Further standardization of sequencing protocols, especially their data analysis components, may require assay-, disease-, and, in many cases, even sample-specific customization that could be performed only in cooperation with clinicians.

Abbreviations

TKI:

tyrosine kinase inhibitor

NSCLC:

non-small cell lung cancer

PCR:

polymerase chain reaction

NGS:

Next Generation Sequencing

CLIA:

Clinical Laboratory Improvement Amendments

CAP:

College of American Pathologists

TSACP:

TruSeq Amplicon—Cancer Panel

RCRC:

Russian Cancer Research Centre

RSCRR:

Russian Scientific Center of Roentgenology and Radiology

FFPE:

Formalin Fixed, Parffine Embedded

ASCO:

American Society of Clinical Oncology

CDF:

cumulative distribution function

SVC:

Illumina Somatic Variant Caller (Illumina Inc, San Diego, California)

CNV:

Copy Number Variations

FISH:

Fluorescent In Situ Hybridization

References

  1. 1.

    Yu KH, Snyder M. Omics profiling in precision oncology. Mol Cell Proteomics. 2016. doi:10.1074/mcp.O116.059253.

  2. 2.

    Thomas A, Rajan A, Lopez-Chavez A, Wang Y, Giaccone G. From targets to targeted therapies and molecular profiling in non-small cell lung carcinoma. Ann Oncol. 2013;24(3):577–85.

  3. 3.

    Stewart EL, Tan SZ, Liu G, Tsao MS. Known and putative mechanisms of resistance to EGFR targeted therapies in NSCLC patients with EGFR mutations-a review. Transl Lung Cancer Res. 2015;4(1):67–81.

  4. 4.

    Therkildsen C, Bergmann TK, Henrichsen-Schnack T, Ladelund S, Nilbert M. The predictive value of KRAS, NRAS, BRAF, PIK3CA and PTEN for anti-EGFR treatment in metastatic colorectal cancer: a systematic review and meta-analysis. Acta Oncol. 2014;53(7):852–64.

  5. 5.

    Foth M, Wouters J, de Chaumont C, Dynoodt P, Gallagher WM. Prognostic and predictive biomarkers in melanoma: an update. Expert Rev Mol Diagn. 2016;16(2):223–37.

  6. 6.

    Hainsworth JD, Meric-Bernstam F, Swanton C, Hurwitz H, Spigel DR, Sweeney C, et al. Targeted therapy for advanced solid tumors based on molecular profiles: early results from MyPathway, an open-label, phase IIa multiple basket study. J Clin Oncol. 2016;34 Suppl abstract LBA11511.

  7. 7.

    Kim ES, Herbst RS, Wistuba II, Lee JJ, Blumenschein GR Jr, Tsao A, et al. The BATTLE trial: personalizing therapy for lung cancer. Cancer Discov. 2011;1(1):44–53.

  8. 8.

    Hyman DM, Puzanov I, Subbiah V, Faris JE, Chau I, Blay JY, et al. Vemurafenib in multiple nonmelanoma cancers with BRAF V600 mutations. N Engl J Med. 2015;373(8):726–36.

  9. 9.

    Iyer G, Hanrahan AJ, Milowsky MI, Al-Ahmadie H, Scott SN, Janakiraman M, et al. Genome sequencing identifies a basis for everolimus sensitivity. Science. 2012;338(6104):221.

  10. 10.

    Pricl S, Cortelazzi B, Dal Col V, Marson D, Laurini E, Fermeglia M, et al. Smoothened (SMO) receptor mutations dictate resistance to vismodegib in basal cell carcinoma. Mol Oncol. 2015;9(2):389–97.

  11. 11.

    Gao J, Adams RP, Swain SM. Does CDKN2A loss predict palbociclib benefit? Curr Oncol. 2015;22(6):e498–501.

  12. 12.

    Villaruz LC, Socinski MA. Temsirolimus therapy in a patient with lung adenocarcinoma harboring an FBXW7 mutation. Lung Cancer. 2014;83(2):300–1.

  13. 13.

    Valsesia A, Macé A, Jacquemont S, Beckmann JS, Kutalik Z. The growing importance of CNVs: new insights for detection and clinical interpretation. Front Genet. 2013;4:92.

  14. 14.

    Aparicio S, Mardis E. Tumor heterogeneity: next-generation sequencing enhances the view from the pathologist’s microscope. Genome Biol. 2014;15(9):463.

  15. 15.

    Watanabe M, Kawaguchi T, Isa S, Ando M, Tamiya A, Kubo A, et al. Ultra-sensitive detection of the pretreatment EGFR T790 M mutation in non-small cell lung cancer patients with an EGFR-activating mutation using droplet digital PCR. Clin Cancer Res. 2015;21(15):3552–60.

  16. 16.

    Do H, Wong SQ, Li J, Dobrovic A. Reducing sequence artifacts in amplicon-based massively parallel sequencing of formalin-fixed paraffin-embedded DNA by enzymatic depletion of uracil-containing templates. Clin Chem. 2013;59(9):1376–83.

  17. 17.

    Blackburn HL, Schroeder B, Turner C, Shriver CD, Ellsworth DL, Ellsworth RE. Management of incidental findings in the era of next-generation sequencing. Curr Genomics. 2015;16(3):159–74.

  18. 18.

    Gagan J, Van Allen EM. Next-generation sequencing to guide cancer therapy. Genome Med. 2015;7(1):80.

  19. 19.

    Gallant JN, Sheehan JH, Shaver TM, Bailey M, Lipson D, Chandramohan R, et al. EGFR kinase domain duplication (EGFR-KDD) is a novel oncogenic driver in lung cancer that is clinically responsive to afatinib. Cancer Discov. 2015;5(11):1155–63.

  20. 20.

    Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.

  21. 21.

    Aziz N, Zhao Q, Bry L, Driscoll DK, Funke B, Gibson JS, et al. College of American Pathologists’ laboratory standards for next-generation sequencing clinical tests. Arch Pathol Lab Med. 2015;139(4):481–93.

  22. 22.

    Schmieder R, Edwards R. Fast identification and removal of sequence contamination from genomic and metagenomic datasets. PLoS One. 2011;6(3):e17288. doi:10.1371/journal.pone.0017288.

  23. 23.

    Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.

  24. 24.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.

  25. 25.

    Saunders CT, Wong WS, Swamy S, Becq J, Murray LJ, Cheetham RK. Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs. Bioinformatics. 2012;28(14):1811–7.

  26. 26.

    Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.

  27. 27.

    Narzisi G, O’Rawe JA, Iossifov I, Fang H, Lee YH, Wang Z, et al. Accurate de novo and transmitted indel detection in exome-capture data using microassembly. Nat Methods. 2014;11(10):1033–6.

  28. 28.

    McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.

  29. 29.

    Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11.

  30. 30.

    Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world’s knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43(Database issue):D805–11.

  31. 31.

    Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010;38(16):e164. doi:10.1093/nar/gkq603.

  32. 32.

    Oliveira C, Wolf T. CNVPanelizer: reliable CNV detection in targeted sequencing applications. R package version 1.0.0. 2015.

  33. 33.

    Tsiatis AC, Norris-Kirby A, Rich RG, Hafez MJ, Gocke CD, Eshleman JR, et al. Comparison of Sanger sequencing, pyrosequencing, and melting curve analysis for the detection of KRAS mutations: diagnostic and clinical implications. J Mol Diagn. 2010;12(4):425–32.

  34. 34.

    Nummela P, Saarinen L, Thiel A, Järvinen P, Lehtonen R, Lepistö A, et al. Genomic profile of pseudomyxoma peritonei analyzed using next-generation sequencing and immunohistochemistry. Int J Cancer. 2015;136(5):E282–9.

  35. 35.

    Simen BB, Yin L, Goswami CP, Davis KO, Bajaj R, Gong JZ, et al. Validation of a next-generation-sequencing cancer panel for use in the clinical laboratory. Arch Pathol Lab Med. 2015;139(4):508–17.

  36. 36.

    Wong SQ, Li J, Tan AY, Vedururu R, Pang JM, Do H, et al. Sequence artefacts in a prospective series of formalin-fixed tumours tested for mutations in hotspot regions by massively parallel sequencing. BMC Med Genomics. 2014;7:23.

  37. 37.

    Do H, Dobrovic A. Dramatic reduction of sequence artefacts from DNA isolated from formalin-fixed cancer biopsies by treatment with uracil-DNA glycosylase. Oncotarget. 2012;3(5):546–58.

  38. 38.

    Arsenic R, Treue D, Lehmann A, Hummel M, Dietel M, Denkert C, et al. Comparison of targeted next-generation sequencing and sanger sequencing for the detection of PIK3CA mutations in breast cancer. BMC Clin Pathol. 2015;15:20.

  39. 39.

    D’Addario G, Früh M, Reck M, Baumann P, Klepetko W, et al. Metastatic non-small-cell lung cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2010;21(Suppl 5):v116–9.

  40. 40.

    Kato S, Han SY, Liu W, Otsuka K, Shibata H, Kanamaru R, et al. Understanding the function-structure and function-mutation relationships of p53 tumor suppressor protein by high-resolution missense mutation analysis. Proc Natl Acad Sci USA. 2003;100(14):8424–9.

  41. 41.

    Takahashi M, Shimodaira H, Andreutti-Zaugg C, Iggo R, Kolodner RD, Ishioka C. Functional analysis of human MLH1 variants using yeast and in vitro mismatch repair assays. Cancer Res. 2007;67(10):4595–604.

  42. 42.

    Fukuchi T, Sakamoto M, Tsuda H, Maruyama K, Nozawa S, Hirohashi S. Beta-catenin mutation in carcinoma of the uterine endometrium. Cancer Res. 1998;58(16):3526–8.

  43. 43.

    Sulis ML, Williams O, Palomero T, Tosello V, Pallikuppam S, Real PJ, et al. NOTCH1 extracellular juxtamembrane expansion mutations in T-ALL. Blood. 2008;112(3):733–40.

  44. 44.

    Rajalingam K, Schreck R, Rapp UR, Albert S. Ras oncogenes and their downstream targets. Biochim Biophys Acta. 2007;1773(8):1177–95.

  45. 45.

    Carpten JD, Faber AL, Horn C, Donoho GP, Briggs SL, Robbins CM, et al. A transforming mutation in the pleckstrin homology domain of AKT1 in cancer. Nature. 2007;448(7152):439–44.

  46. 46.

    Cai Q, Robertson ES. Ubiquitin/SUMO modification regulates VHL protein stability and nucleocytoplasmic localization. PLoS ONE. 2010;5(9):e12636.

  47. 47.

    Cai Q, Verma SC, Kumar P, Ma M, Robertson ES. Hypoxia inactivates the VHL tumor suppressor through PIASy-mediated SUMO modification. PLoS ONE. 2010;5(3):e9720.

  48. 48.

    Wang H, Bao W, Jiang F, Che Q, Chen Z, Wang F, et al. Mutant p53 (p53-R248Q) functions as an oncogene in promoting endometrial cancer by up-regulating REGγ. Cancer Lett. 2015;360(2):269–79.

  49. 49.

    Muller PA, Vousden KH. Mutant p53 in cancer: new functions and therapeutic opportunities. Cancer Cell. 2014;25(3):304–17.

  50. 50.

    Zahari MS, Wu X, Blair BG, Pinto SM, Nirujogi RS, Jelinek CA, et al. Activating mutations in PIK3CA Lead to widespread modulation of the tyrosine phosphoproteome. J Proteome Res. 2015;14(9):3882–91.

  51. 51.

    Wu X, Renuse S, Sahasrabuddhe NA, Zahari MS, Chaerkady R, Kim MS, et al. Activation of diverse signalling pathways by oncogenic PIK3CA mutations. Nat Commun. 2014;5:4961.

  52. 52.

    Satya RV, DiCarlo J. Edge effects in calling variants from targeted amplicon sequencing. BMC Genom. 2014;15:1073.

  53. 53.

    Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60.

  54. 54.

    Kim ES, Hirsh V, Mok T, Socinski MA, Gervais R, Wu YL, et al. Gefitinib versus docetaxel in previously treated non-small-cell lung cancer (INTEREST): a randomised phase III trial. Lancet. 2008;372(9652):1809–18.

  55. 55.

    Crinò L, Cappuzzo F, Zatloukal P, Reck M, Pesek M, Thompson JC, et al. Gefitinib versus vinorelbine in chemotherapy-naive elderly patients with advanced non-small-cell lung cancer (INVITE): a randomized, phase II study. J Clin Oncol. 2008;26(26):4253–60.

  56. 56.

    Fukuoka M, Wu Y, Thongprasert S, Yang C, Chu D, Saijo N, et al. Biomarker analyses from a phase III, randomized, open-label, first-line study of gefitinib (G) versus carboplatin/paclitaxel (C/P) in clinically selected patients (pts) with advanced non-small cell lung cancer (NSCLC) in Asia (IPASS). J Clin Oncol. 2009;27:15 Suppl abstract 8006.

  57. 57.

    Jia XF, Li J, Zhao HB, Liu J, Liu JJ. Correlation of EGFR gene amplification with invasion and metastasis of non-small cell lung cancer. Genet Mol Res. 2015;14(3):11006–12.

  58. 58.

    Benedettini E, Sholl LM, Peyton M, Reilly J, Ware C, Davis L, et al. Met activation in non-small cell lung cancer is associated with de novo resistance to EGFR inhibitors and the development of brain metastasis. Am J Pathol. 2010;177(1):415–23.

  59. 59.

    Bean J, Brennan C, Shih JY, Riely G, Viale A, Wang L, et al. MET amplification occurs with or without T790 M mutations in EGFR mutant lung tumors with acquired resistance to gefitinib or erlotinib. Proc Natl Acad Sci USA. 2007;104(52):20932–7.

  60. 60.

    Chen HJ, Mok TS, Chen ZH, Guo AL, Zhang XC, Su J, et al. Clinicopathologic and molecular features of epidermal growth factor receptor T790M mutation and c-MET amplification in tyrosine kinase inhibitor-resistant Chinese non-small cell lung cancer. Pathol Oncol Res. 2009;15(4):651–8.

  61. 61.

    Iommelli F, De Rosa V, Gargiulo S, Panico M, Monti M, Greco A, et al. Monitoring reversal of MET-mediated resistance to EGFR tyrosine kinase inhibitors in non-small cell lung cancer using 3′-deoxy-3′-[18F]-fluorothymidine positron emission tomography. Clin Cancer Res. 2014;20(18):4806–15.

  62. 62.

    Ou SH, Kwak EL, Siwak-Tapp C, Dy J, Bergethon K, Clark JW, et al. Activity of crizotinib (PF02341066), a dual mesenchymal-epithelial transition (MET) and anaplastic lymphoma kinase (ALK) inhibitor, in a non-small cell lung cancer patient with de novo MET amplification. J Thorac Oncol. 2011;6(5):942–6.

  63. 63.

    Noonan SA, Berry L, Lu X, Gao D, Barón AE, Chesnut P, et al. Identifying the appropriate FISH criteria for defining MET copy number-driven lung adenocarcinoma through oncogene overlap analysis. J Thorac Oncol. 2016;S1556–0864(16):30487.

  64. 64.

    Cancer Genome Atlas Research Network. Comprehensive genomic characterization of squamous cell lung cancers. Nature. 2012;489(7417):519–25.

  65. 65.

    Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511(7511):543–50.

  66. 66.

    Gerlinger M, Rowan AJ, Horswell S, Larkin J, Endesfelder D, Gronroos E, et al. Intratumor heterogeneity and branched evolution revealed by multiregion sequencing. N Engl J Med. 2012;366(10):883–92.

  67. 67.

    Serizawa M, Yokota T, Hosokawa A, Kusafuka K, Sugiyama T, Tsubosa Y, et al. The efficacy of uracil DNA glycosylase pretreatment in amplicon-based massively parallel sequencing with DNA extracted from archived formalin-fixed paraffin-embedded esophageal cancer tissues. Cancer Genet. 2015;208(9):415–27.

  68. 68.

    Atreya CE, Corcoran RB, Kopetz S. Expanded RAS: refining the patient population. J Clin Oncol. 2015;33(7):682–5.

  69. 69.

    Silva IP, Salhi A, Giles KM, Vogelsang M, Han SW, Ismaili N, et al. Identification of a novel pathogenic germline KDR variant in melanoma. Clin Cancer Res. 2016;22(10):2377–85.

  70. 70.

    Glubb DM, Cerri E, Giese A, Zhang W, Mirza O, Thompson EE, et al. Novel functional germline variants in the VEGF receptor 2 gene and their effect on gene expression and microvessel density in lung cancer. Clin Cancer Res. 2011;17(16):5257–67.

  71. 71.

    Dumont P, Leu JI, Della Pietra AC 3rd, George DL, Murphy M. The codon 72 polymorphic variants of p53 have markedly different apoptotic potential. Nat Genet. 2003;33(3):357–65.

  72. 72.

    Zeichner SB, Alghamdi S, Elhammady G, Poppiti RJ. Prognostic significance of TP53 mutations and single nucleotide polymorphisms in acute myeloid leukemia: a case series and literature review. Asian Pac J Cancer Prev. 2014;15(4):1603–9.

  73. 73.

    Lo Iacono M, Buttigliero C, Monica V, Bollito E, Garrou D, Cappia S, et al. Retrospective study testing next generation sequencing of selected cancer-associated genes in resected prostate cancer. Oncotarget. 2016;7(12):14394–404.

  74. 74.

    Masago K, Fujita S, Muraki M, Hata A, Okuda C, Otsuka K, et al. Next-generation sequencing of tyrosine kinase inhibitor-resistant non-small-cell lung cancers in patients harboring epidermal growth factor-activating mutations. BMC Cancer. 2015;15:908.

  75. 75.

    Sullivan A, Syed N, Gasco M, Bergamaschi D, Trigiante G, Attard M, et al. Polymorphism in wild-type p53 modulates response to chemotherapy in vitro and in vivo. Oncogene. 2004;23(19):3328–37.

  76. 76.

    Pim D, Banks L. p53 polymorphic variants at codon 72 exert different effects on cell cycle progression. Int J Cancer. 2004;108(2):196–9.

  77. 77.

    Inokuchi K, Yamaguchi H, Tarusawa M, Futaki M, Hanawa H, Tanosaki S, et al. Abnormality of c-kit oncoprotein in certain patients with chronic myelogenous leukemia–potential clinical significance. Leukemia. 2002;16(2):170–7.

  78. 78.

    Foster R, Byrnes E, Meldrum C, Griffith R, Ross G, Upjohn E, et al. Association of paediatric mastocytosis with a polymorphism resulting in an amino acid substitution (M541L) in the transmembrane domain of c-KIT. Br J Dermatol. 2008;159(5):1160–9.

  79. 79.

    Dufresne A, Alberti L, Brahmi M, Kabani S, Philippon H, Pérol D, et al. Impact of KIT exon 10 M541L allelic variant on the response to imatinib in aggressive fibromatosis: analysis of the desminib series by competitive allele specific Taqman PCR technology. BMC Cancer. 2014;14:632.

  80. 80.

    Dufresne A, Bertucci F, Penel N, Le Cesne A, Bui B, Tubiana-Hulin M, et al. Identification of biological factors predictive of response to imatinib mesylate in aggressive fibromatosis. Br J Cancer. 2010;103(4):482–5.

  81. 81.

    Grabellus F, Worm K, Sheu SY, Siffert W, Schmid KW, Bachmann HS. The prevalence of the c-kit exon 10 variant, M541L, in aggressive fibromatosis does not differ from the general population. J Clin Pathol. 2011;64(11):1021–4.

  82. 82.

    Gonçalves A, Monges G, Yang Y, Palmerini F, Dubreuil P, Noguchi T, et al. Response of a KIT-positive extra-abdominal fibromatosis to imatinib mesylate and KIT genetic analysis. J Natl Cancer Inst. 2006;98(8):562–3.

  83. 83.

    Seinfeld J, Kleinschmidt-Demasters BK, Tayal S, Lillehei KO. Desmoid-type fibromatoses involving the brachial plexus: treatment options and assessment of c-KIT mutational status. J Neurosurg. 2006;104(5):749–56.

  84. 84.

    Meyer A, Wilhelm B, Dörk T, Bremer M, Baumann R, Karstens JH, et al. ATM missense variant P1054R predisposes to prostate cancer. Radiother Oncol. 2007;83(3):283–8.

  85. 85.

    Stredrick DL, Garcia-Closas M, Pineda MA, Bhatti P, Alexander BH, Doody MM, et al. The ATM missense mutation p.Ser49Cys (c.146C > G) and the risk of breast cancer. Hum Mutat. 2006;27(6):538–44.

  86. 86.

    Webb EL, Rudd MF, Sellick GS, El Galta R, Bethke L, Wood W, et al. Search for low penetrance alleles for colorectal cancer through a scan of 1467 non-synonymous SNPs in 2575 cases and 2707 controls with validation by kin-cohort analysis of 14,704 first-degree relatives. Hum Mol Genet. 2006;15(21):3263–71.

  87. 87.

    Rudd MF, Sellick GS, Webb EL, Catovsky D, Houlston RS. Variants in the ATM-BRCA2-CHEK2 axis predispose to chronic lymphocytic leukemia. Blood. 2006;108(2):638–44.

  88. 88.

    Meier M, den Boer ML, Hall AG, Irving JA, Passier M, Minto L, et al. Relation between genetic variants of the ataxia telangiectasia-mutated (ATM) gene, drug resistance, clinical outcome and predisposition to childhood T-lineage acute lymphoblastic leukaemia. Leukemia. 2005;19(11):1887–95.

Download references

Authors’ contributions

Study concept and design: VM, AB, MI, SM, KL, VB, PC. Subjects recruitment and specimen collection: KL, VB, PC, EN, VM. DNA sequencing and supervision: EN, VM. Bioinformatic analysis of sequencing data: MI. Data analysis and interpretation: MI, VM. Drafting manuscript: MI, AB. Revising article critically for important intellectual content: VM, AB, KL, VB, PC, EN. All authors read and approved the final manuscript.

Acknowledgements

The authors wish to gratefully acknowledge technical support from Readsense, Ltd (Moscow, 142190, Russia).

Competing interests

Ivanov M, Musienko S and Mileyko V are employees of Atlas OncoDiagnostics, a company commercially involved in development and validation of various research and diagnostic kits.

Ethics approval and consent to participate

The study protocol was approved by Atlas Biomed Internal Review Board. Written informed consent was provided by all patients at inception, all analyses were based on archival data and stored in database with no connections to the patient identifiers.

Funding

This study was supported by Ministry of Science and Education, Russia (Project No. RFMEFI60714X0098).

Author information

Correspondence to Ancha Baranova.

Additional information

Ancha Baranova and Vladislav Mileyko contributed equally to this work

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • NGS
  • Non-small cell lung cancer
  • Cancer
  • Oncology
  • Precision oncology
  • Targeted therapy