High resolution melting analysis for rapid and sensitive EGFR and KRAS mutation detection in formalin fixed paraffin embedded biopsies

doi:10.1186/1471-2407-8-142

BMC Cancer. 2008; 8: 142.

Published online 2008 May 21. doi: 10.1186/1471-2407-8-142

PMCID: PMC2408599

High resolution melting analysis for rapid and sensitive EGFR and KRAS mutation detection in formalin fixed paraffin embedded biopsies

Hongdo Do,^1,² Michael Krypuy,¹ Paul L Mitchell,³ Stephen B Fox,^1,² and Alexander Dobrovic^1,²

¹Molecular Pathology Research and Development Laboratory, Department of Pathology, Peter MacCallum Cancer Centre, Locked Bag 1, A'Beckett St, Melbourne, Victoria 8006, Australia

²Department of Pathology, University of Melbourne, Parkville, Victoria, 3010, Australia

³Ludwig Medical Oncology Department, Austin Hospital, Heidelberg, Victoria, 3084, Australia

Corresponding author.

Hongdo Do: hongdo.do/at/petermac.org; Michael Krypuy: mkrypuy/at/gmail.com; Paul L Mitchell: paul.mitchell/at/austin.org.au; Stephen B Fox: stephen.fox/at/petermac.org; Alexander Dobrovic: alexander.dobrovic/at/petermac.org

Received February 6, 2008; Accepted May 21, 2008.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

Abstract

Background

Epithelial growth factor receptor (EGFR) and KRAS mutation status have been reported as predictive markers of tumour response to EGFR inhibitors. High resolution melting (HRM) analysis is an attractive screening method for the detection of both known and unknown mutations as it is rapid to set up and inexpensive to operate. However, up to now it has not been fully validated for clinical samples when formalin-fixed paraffin-embedded (FFPE) sections are the only material available for analysis as is often the case.

Methods

We developed HRM assays, optimised for the analysis of FFPE tissues, to detect somatic mutations in EGFR exons 18 to 21. We performed HRM analysis for EGFR and KRAS on DNA isolated from a panel of 200 non-small cell lung cancer (NSCLC) samples derived from FFPE tissues.

Results

All 73 samples that harboured EGFR mutations previously identified by sequencing were correctly identified by HRM, giving 100% sensitivity with 90% specificity. Twenty five samples were positive by HRM for KRAS exon 2 mutations. Sequencing of these 25 samples confirmed the presence of codon 12 or 13 mutations. EGFR and KRAS mutations were mutually exclusive.

Conclusion

This is the first extensive validation of HRM on FFPE samples using the detection of EGFR exons 18 to 21 mutations and KRAS exon 2 mutations. Our results demonstrate the utility of HRM analysis for the detection of somatic EGFR and KRAS mutations in clinical samples and for screening of samples prior to sequencing. We estimate that by using HRM as a screening method, the number of sequencing reactions needed for EGFR and KRAS mutation detection can be reduced by up to 80% and thus result in substantial time and cost savings.

Background

Lung cancer is the leading cause of cancer-related death, accounting for one third of all cancer mortality worldwide due to high incidence, advanced stage at diagnosis and aggressive tumour behaviour [1]. Non-small cell lung cancer (NSCLC), comprising 80% of all lung cancer cases, has a poor prognosis if diagnosed at an advanced stage. There is a low median survival of less than one year after diagnosis when treated by conventional chemotherapy [2].

The epidermal growth factor receptor (EGFR) is a member of the ErbB receptor tyrosine kinase family. EGFR has been found to be over-expressed in a variety of human malignancies [3,4]. Activation of EGFR results in the initiation of a diverse range of cellular signalling pathways, including cell proliferation and protection of the cell from apoptosis [5,6]. Activating mutations in the tyrosine kinase domain of the EGFR gene (EGFR) have been shown to be associated with a dramatic response to the tyrosine kinase inhibitors (TKIs), such as gefitinib and erlotinib [7-10]. These mutations are located in exons 18 to 21 and are more common in females, non-smokers, tumours with a histological diagnosis of adenocarcinoma, and individuals of Asian descent [11]. In cell line studies, these mutations have been shown to induce oncogenic transformation of fibroblasts and lung epithelial cells [7,12-15].

KRAS is a member of the Ras gene family which encode small G proteins with intrinsic GTPase activity. The KRAS protein plays a key role in Ras/MAPK signalling which is involved in multiple pathways including proliferation, differentiation and apoptosis [16]. KRAS mutations which are found in 33% of NSCLC are restricted to specific codons; more commonly codons 12 and 13 in exon 2 and rarely codons 59 and 61 in exon 3 [16,17]. These mutations alter the conformation of KRAS causing impaired GTPase activity resulting in the protein being constitutively active. KRAS mutation testing is an important adjunct to EGFR testing because KRAS mutations are significantly associated with absence of responsiveness to EGFR inhibitors and are mutually exclusive to EGFR mutations [18-20]. Thus, the mutational status of EGFR and KRAS can provide important information for stratification of NSCLC patients to receive molecularly targeted treatment with tyrosine kinase inhibitors.

Currently, the most widely used method for EGFR and KRAS mutation detection is direct sequencing. To be successful, the sequencing methodology requires a sufficient amount of tumour material of relatively good quality, which is difficult to obtain from cancer patients with inoperable tumours. In addition, the high cost, limited sensitivity and time consuming nature of sequencing has prompted the development of alternative methods that are more cost effective, faster, easier to perform, and more sensitive. The relatively low sensitivity of sequencing in somatic mutation detection [21] is a particular problem in lung cancer where biopsies are often small and often contain only a small proportion of neoplastic cells. Studies using denaturing high performance liquid chromatography (DHPLC) have found additional EGFR mutations, which were undetected by sequencing [22-24]. However DHPLC, requiring extra sample handling after PCR amplification and expensive instrumentation, is relatively slow as the samples can only be analysed sequentially.

High resolution melting (HRM) analysis is a recently developed methodology that has enormous potential for the detection of DNA sequence changes [25]. New instruments combined with DNA intercalating dyes that can be used at saturating concentrations allow the discrimination of sequence changes in PCR amplicons without manual handling of PCR products. The recent application of HRM to mutation scanning and SNP genotyping as well as DNA methylation studies have been demonstrated [26-31]. In particular, the HRM methodology has shown great promise for the detection of heterozygous germline mutations as well as somatic mutations e.g. in the KIT, BRAF and TP53 genes [26,32,33].

Formalin fixation followed by paraffin embedding of tissue specimens is a widely used preservation method because it helps to maintain morphological features of the tissue specimen [34]. However, irreversible damage to DNA can occur during this process or subsequent prolonged storage, resulting in an adverse effect on DNA quality [35].

As FFPE tissues are the most common clinical specimens used for detection of the EGFR and KRAS mutations, the validation of HRM using DNA samples extracted from FFPE tissues is essential for its application as a screening method for mutation detection in these genes. In this study, we developed HRM assays to evaluate the efficacy of this methodology for screening EGFR mutations in exons 18 to 21 using a panel of 200 NSCLC FFPE biopsies. KRAS exon 2 mutations at codon 12 and 13 were also screened in this sample set using a previously described HRM assay.

Methods

Samples

A total of 200 lung cancer biopsy specimens were analysed from specimens sent to the Peter MacCallum Cancer Centre, Melbourne, Australia for EGFR mutation detection by direct sequencing. These patients were often referred for EGFR mutation testing by clinicians since their clinical features and history were suggestive of those previously reported for patients with EGFR-associated mutations; histological diagnosis of adenocarcinoma, female gender, non-smoker status and Asian ethnicity. Of the 200 samples, 141 were adenocarcinomas, 24 were large cell carcinomas, 10 were squamous cell carcinomas, and 25 were of other or unknown histologies. The samples were independently assessed by pathologists from both the referring hospital, and from the Peter MacCallum Cancer Centre. Tumour-rich areas on a hematoxylin and eosin slide were marked by the pathologist to ensure that the maximum amount of tumour material was collected for the genetic testing. Written informed consent was obtained from all patients or families prior to testing. This study was approved by the Ethics of Human Research Committee at the Peter MacCallum Cancer Centre (project number 03/90).

DNA extraction from formalin-fixed paraffin-embedded (FFPE) tissue

Tissue sections of 5 μm thickness were obtained from FFPE tissues and stained with Methyl green. The tumour-rich areas were micro-dissected using a 21G needle and the samples underwent proteinase K digestion in a rotating incubator at 56°C for 3 days. Genomic DNA was extracted using the DNeasy Tissue kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol and was kept at 4°C before use.

Design of HRM primers

Several factors were taken into account during the design of primers for this study. As FFPE DNA was used, primers that would yield reproducible amplification from degraded templates, or DNA in which the quality has been compromised were designed. As a result, primers giving rise to shorter amplicons were chosen as they are more likely to result in satisfactory amplification from degraded FFPE DNA. Furthermore, existing primer pairs that were used for the sequencing of EGFR exons 18 to 21 (Table (Table1)1) were found to be inappropriate for HRM as, with these longer amplicons (233, 294, 356 and 295 bp for exons 18 to 21 respectively), samples harbouring mutations were not efficiently distinguished from wild type samples in melting analysis. We thus chose to design primers that flanked the exons as closely as possible. As single nucleotide polymorphisms (SNPs) cannot be readily distinguished from mutations by HRM analysis, designing shorter amplicons in this way also minimises the inclusion of these confounding sequence variants. Most intronic SNPs can thus be excluded from the amplicons. If a SNP is close to the exon boundary, the primer can be placed over the SNP and a mismatched base with no allelic preference can be introduced at the SNP position [26]. Thus, the base 'A' was incorporated at the c.2184+19G>A position within the EGFR exon 18 reverse primer. The primers were designed to have annealing temperatures of around 60°C as predicted by Primer Express 1.5 software (Applied Biosystems, Foster City, CA). Each amplicon was analysed to ensure that it contained only a single melting domain using the Poland algorithm [36].

Table 1

EGFR HRM and sequencing primer sequences

HRM assays

The EGFR HRM primer sequences are shown in Table Table1.1. EGFR exon 18, 19 and 21 HRM primers were designed to span the entire exon with product sizes of 199, 250 and 210 bp respectively. EGFR exon 20 was amplified in two fragments (ex20a and ex20b) with product sizes of 121 bp and 146 bp. The KRAS exon 2 primers which amplified a 92 bp product were described previously [37].

PCR for HRM analysis was performed in 0.1 ml tubes on the Rotor-Gene 6000™ (Corbett Research, Sydney, Australia) in the presence of the fluorescent DNA intercalating dye, SYTO 9 (Invitrogen, Carlsbad, CA). The reaction mixture in a 20 μl final volume contained; 1× PCR buffer, 2.5 mM MgCl₂, 200–400 nM forward primer, 200–400 nM reverse primer, 5 ng of genomic DNA, 200 μM of dNTPs, 5 μM of SYTO 9, 0.5 U of HotStarTaq (Qiagen) polymerase and PCR grade water. The cycling and melting conditions for EGFR exons 18 to 21 were as follows; one cycle of 95°C for 15 min; 45–50 cycles of 95°C for 10 s, 65°C for 10 s with an initial 10 cycles of touchdown (1°C/cycle), 72°C for 30 s; one cycle of 97°C for 1 min and a melt from 70°C to 95°C rising 0.2°C per second. The cycling and melting conditions for KRAS exon 2 were as follows; one cycle of 95°C for 15 min; 50 cycles of 95°C for 10 s, 67.5°C for 5 s with an initial 10 cycles of touchdown (1°C/cycle), 72°C for 20 s; one cycle of 97°C for one min and a melt from 70°C to 95°C rising 0.2°C per second. The genomic DNA was diluted to 2.5 ng/μl (5 ng tested) to provide a consistent testing condition. All samples were tested in duplicate.

HRM analysis

High resolution melting analysis was performed on the Rotor-Gene 6000 Software (v1.7) and analysed by two scientists who were blinded to the sequencing results. The normalised graph and the difference graph were used to analyse the data. The normalised graph was generated by the monitoring of dissociation of the fluorescent dye from double-stranded DNA as the temperature increased. The dye (SYTO 9) used in the current study can only fluoresce when it is intercalated into double-strand DNA. The normalised graph shows the degree of reduction in fluorescence over a temperature range (typically 70°C to 95°C). All samples including the wild-type were plotted according to their melting profiles. In the difference graph, the melting profiles of each sample were compared to that of the wild-type which was converted to a horizontal line. Significant deviations from the horizontal line (relative to the spread of the wild type controls) were indicative of sequence changes within the amplicon analysed. Samples with aberrant melting curves were recorded as HRM mutation positive. HRM results were compared with sequencing results for the validation of HRM analysis in the detection of EGFR mutations.

DNA sequencing

All samples were sequenced for the detection of EGFR mutations in exons 18 to 21 using M13 tagged primers (Table (Table1).1). The reaction mixture in a total of 25 μl contained the following; 1× PCR buffer, 2.5 mM MgCl₂, 200 nM of each primer, 50 ng of genomic DNA (if possible), 200 μM of dNTPs, 0.5 U of HotStarTaq polymerase and PCR grade water. The PCR reaction was performed using the following conditions; initial denaturation at 95°C for 15 min; 40 cycles of 94°C for 45 s, 65°C for 45 s, 72°C for 45 s; one cycle of 72°C for 10 min. 6 μl of the PCR products were purified with ExoSapIT (GE Healthcare, Little Chalfont, England) followed by a sequencing reaction with Big Dye Terminator v3.1 (Applied Biosystems, Foster City, CA) according to the manufacturer's protocol. The sequencing products were ethanol precipitated before running on a 3100 Genetic Analyser (Applied Biosystems). The sequencing data was visualized using Sequencher 4.6 (Gene Codes Corporation, Ann Arbor, MI) and was independently analysed by two scientists. Each mutation was confirmed by sequencing a second independent PCR reaction.

Samples were sequenced for KRAS exon 2 mutations using an 189 bp amplicon. The sequencing primers were described previously [37]. The reaction mixture in a total of 20 μl was made using HotStarTaq (Qiagen) and contained the following; 1× PCR buffer, 2.5 mM MgCl₂, 200 nM of each primer, 5 ng of genomic DNA, 200 μM of dNTPs, 0.5 U of Taq polymerase and PCR grade water PCR reactions were performed using the following conditions; initial denaturation at 95°C for 15 min; 50 cycles of 95°C for 10 s, 67.5°C for 10 s, 72°C for 20 s with an initial 10 cycles of touchdown (1°C/cycle); one cycle of 72°C for 10 min. Sequencing reactions were performed as above except that an annealing temperature of 60.7°C and an extension temperature of 72°C were used.

Results

EGFR mutation detection by direct sequencing

EGFR exons 18 to 21 were sequenced using genomic DNA extracted from FFPE tissues from the 200 NSCLC patients. In total, 73 samples were identified with EGFR mutations by sequencing (Table (Table2).2). The EGFR mutations were most common in exon 19, comprising 61% (46/75) of all mutations found, followed by exon 21 missense mutations in 18% (14/75). Ten mutations in exon 20 (13%) and five mutations in exon 18 (6%) were also detected.

Table 2

Summary of EGFR mutations detected by sequencing from the 200 NSCLC samples

There were three cases of rare mutations. A 4 bp deletion and 1 bp insertion (p.E709_T710delinsD) and a 9 bp deletion (p.A767_V769del) mutations were detected in exons 18 and 20 respectively. In exon 19, there was a novel in-frame 18 bp insertion mutation that replaced the glutamic acid at p.746 by valine and inserted six amino acids (PVAIKE) afterwards (Table (Table22).

We also found two samples harbouring double mutations in exons 18 and 20. The mutation in exon 20 in both individuals was the same missense substitution (p.S768I). The exon 18 mutations were p.G719S and p.G724S. The combination of p.G719S and p.S768I mutations has been reported previously in NSCLC samples [38].

Two novel synonymous mutations, p.A763A (c.2289C>T) and p.N771N (c.2313C>T), were found in exon 20. A novel intronic variant with C>T at c.2469+21 was found in a single patient. These may represent passenger mutations or rare SNPs although the germline DNA was not able to be tested to distinguish the two possibilities

HRM assay with different amounts of DNA template

In many NSCLC cases, the FFPE specimens available for analysis are often small in size. To investigate whether HRM assay can be performed on a small amount of DNA template, DNA amounts ranging from 1 to 100 ng (1, 5, 10, 25, 50 and 100 ng) were tested using a sample with an exon 19 p.E746_A750 deletion mutation present at 50% mutant allele frequency.

Distinctive heteroduplex melting patterns were detected across the entire range of DNA template amounts, featuring an earlier melting of amplicon at the initial melting stage and becoming more stable at later stages compared with the wild-type (Figure (Figure1.1. Panel A). In particular, the melting curve from 1 ng of mutant template was sufficiently different from wild-type to identify the mutation, and this distinct melting profile was consistently seen across all the other template amounts investigated. This illustrates the sensitivity of the technique as mutations can be detected in samples with as little as 1 ng of template DNA by HRM analysis.

Figure 1

Difference plots showing mutation detection with various amounts of starting template. In difference plots, the melting profile of a wild type control is chosen as a horizontal base line and the relative differences in the melting of all the other samples (more ...)

After confirming that HRM can be performed with a small amount of template DNA, we further tested the practicality of HRM using a panel of samples with different DNA concentrations, ranging from 8 ng to 651 ng. All samples harboured the same EGFR exon 19 p.E746_A750 deletion mutation. Again, all mutant samples exhibited distinct melting characteristics from wild-type, and could easily be identified as mutants (Figure (Figure1.1. Panel B).

These analyses demonstrated that reliable results could also be obtained when more than 1 ng of template DNA was tested by HRM. For the samples that amplified well, the variation in the amount of template within samples had minimal influence on the HRM analysis, indicating that adjusting the amount of template might normally be unnecessary.

EGFR and KRAS mutation detection by HRM analysis

All 200 NSCLC samples were tested by HRM for the detection of mutations in EGFR exons 18 to 21 and KRAS exon 2. The EGFR HRM assays gave 23, 67, 23 and 41 results that were scored as HRM positive in exons 18 to 21 respectively (Table (Table3).3). All mutations identified by sequencing were correctly identified by HRM assays giving a sensitivity of 100%. The difference plots for EGFR exons 18 to 21 and sequencing traces of positive samples are shown in Figures Figures22 and and33.

Table 3

Summary of EGFR mutation testing by sequencing and HRM

Figure 2

Difference plots and sequence traces for EGFR exons 18 and 19. Panel A: The difference plot of EGFR exon 18 shows melting profiles for three positive samples (TX36 in purple, TX53 in orange and TX290 in red). Panel B: Sequencing chromatograms show a p.E709A (more ...)

Figure 3

Difference plots and sequence traces for EGFR exons 20 and 21. Panel A: The difference plot of EGFR HRM exon 20a shows different melting profiles of three positive samples (TX53 in red, TX383 in brown and TX440 in orange). Panel B: A p.S768I in TX53, (more ...)

However, HRM indicated more positive samples than sequencing for all of the EGFR exons. A total of 45 samples were positive only by HRM. There were 18, 21, 13 and 18 apparently false positive results from exons 18 to 21 respectively. Significantly, seven samples were positive in three or four EGFR HRM assays, and eight samples were positive in two assays. Most of these are likely to be true false positives due to degraded DNA from the FFPE specimen. Thirty samples were positive in a single assay only. Compared to other assays, the number of samples with discrepant results was lower in the exon 20 assays where relatively shorter amplicons were tested, further highlighting the importance of using shorter HRM amplicons to enhance the reliability of the results.

KRAS exon 2 mutations at codons 12 to 13 were screened by our previously described HRM assay [28] that generated a 92 bp PCR amplicon (Figure (Figure4).4). Twenty-five samples of the 200 tested (12.5%) were scored as HRM positive. These samples were negative for EGFR mutations by both sequencing and HRM.

Figure 4

Difference plots and sequence traces for KRAS exon 2 mutations. In the difference plot for KRAS exon 2, the melting curves from the four samples (TX23 in purple, TX41 in green, TX55 in orange and TX64 in brown) and a positive mutation control (RPMI8226 (more ...)

All 25 KRAS exon 2 HRM positive samples proved to have a mutation at codons 12 or 13 when tested by sequencing (Table (Table4).4). Seven different KRAS mutations were observed amongst the 25 mutation positive samples, comprising 5 mutations at codon 12 (8 p.G12C, 7 p.G12D, 5 p.G12V, 1 p.G12F and 1 p.G12A) and 2 mutations at codon 13 (2 p.G13C and 1 p.G13V). In addition, 25 samples that were both EGFR mutation negative and KRAS mutation negative were sequenced for KRAS exon 2. All showed normal KRAS sequence.

Table 4

Summary of KRAS exon 2 mutations detected by HRM and sequencing

Discussion

The advent of molecularly targeted therapy is shifting the paradigm of management of cancer patients from generalised chemotherapy and/or radiotherapy to personalised treatments with better efficacy and lower toxicity. An example of this is the development of EGFR tyrosine kinase inhibitors in the management of NSCLC patients where the presence of predictive markers such as EGFR and KRAS mutations in their tumours stratifies patients to receive the appropriate treatment.

Currently, direct sequencing is the standard method for EGFR mutation detection, but its limited sensitivity, high cost and long turnaround time have prompted the development of alternative methods for routine clinical testing which have greater diagnostic practicality for somatic mutation detection.

HRM has recently been introduced as a screening method for mutation detection. It is an in-tube method which can be performed in a fast, cheap, and robust manner. It has been applied to the detection of both germline and somatic mutations. For heterozygous germline mutations, HRM has sensitivities approaching 100% [32,39]. For somatic mutations in tumours, detection can be compromised by a low proportion of tumour cells in the biopsy. In practice, it has been shown previously that mutant alleles at levels as low as 5 to 10% can be detected by HRM for KRAS exon 2 mutations [29,37].

HRM has previously been used for detection of EGFR mutations. In lung cancer, the two most common types of EGFR mutation, exon 19 deletions and exon 21 p.L858R, were screened by HRM, with a reported 92% sensitivity compared with direct sequencing [29,40]. In head and neck cancer, EGFR exons 18 to 21 were screened by HRM, resulting in the detection of two EGFR mutations in 24 squamous cell carcinomas [28].

We designed this study to determine whether HRM can be a diagnostically useful screening method for EGFR and KRAS mutations in clinical FFPE specimens by validating HRM against sequencing in a large sample cohort containing various types of EGFR and KRAS mutations.

Each of the previously described HRM mutation screening assays for EGFR have limitations preventing them from being practical for one or more of the following reasons: they do not cover all of exons 18 to 21 [29,40], they do not discriminate against common SNPs leading to unnecessary sequencing [28], they do not detect all mutations that are detectable by sequencing [29,40] and they have not been validated against a large panel of clinical samples with previously determined mutations [28].

In one assay, a primer is located over the c.2184+19 common SNP potentially resulting in non-amplification of the mutant allele [28]. Allele dropout due to sequence variation at the primer binding site has been previously demonstrated [41]. If the allele that has not amplified contains the mutation, the mutation will not be detected.

In this study, HRM assays were designed to maximise the benefits of HRM screening and to minimise SNP interference. EGFR exon 18 to 21 assays covering the entire coding regions were developed. The incorporation of mismatched bases at SNP loci within the primer sequences made it possible to exclude two SNPs, c.2361G>A and c.2184+19G>A. Under the conditions used in our assays, this did not preclude efficient amplification of both alleles, allowing us to detect all 15 mutations from exons 18 and 20.

Due to its high frequency, the exonic SNP, c.2361G>A, would normally necessitate nearly 50% of samples being sequenced for exon 20 because the heterozygous melting pattern given by the SNP can not readily be distinguished from mutation by HRM. Fortunately, its position, in the middle of a large exon, allowed us to divide exon 20 into two fragments using two overlapping amplicons with PCR product sizes of 121 bp and 146 bp. Primers which overlaid the SNP and contained a mismatched residue at the SNP location were used for both amplicons to exclude the SNP from being detected by HRM. A 'G' and an 'A' respectively (mismatched to both alleles) were introduced at c.2361 into the exon 20a reverse and exon 20b forward primers. This strategy led to a short region of 15 bp (c.2353_2367) flanking the SNP for which mutations could not be detected. This is not a serious limitation as no mutations have so far been reported in that region [42].

The other common exonic SNP, c.2508C>T in exon 21, is present at a much lower frequency and thus necessitated only a comparatively small increase in the amount of sequencing that would be required. The SNP was detected in nine of the 200 samples.

It is important to consider the amplification information when interpreting HRM results. Instruments allowing the real time monitoring of amplification are advantageous in this regard. We observed that melting curves from samples with insufficient amplification tended to be shifted to the right relative to the wild-type curves in the normalised plots. This implies that the amount of the amplifiable (functional) templates varies in each sample depending on the degree of DNA degradation even though they are all adjusted to the same concentration (2.5 ng/μl). For those samples, we diluted the wild-type control DNA to get similar amplification to the samples with insufficient amplification. We also increased the amplification cycle number to 60 to allow sufficient amplification for melting analysis. As shown in Figure Figure5,5, with sufficient amplification, the right shifting of melting curves was corrected and thus the patient DNA could be reliably compared to the wild-type.

Figure 5

Influence of amplification status on melting analysis. Panel A: The amplification plot of EGFR exon 19 shows insufficient amplification of several samples relative to the wild-type (blue). Panel B: The subsequent melting analysis of the same samples on (more ...)

Accurate identification of mutations is a crucial aspect of all mutation screening methods. This current study demonstrates the accuracy of HRM in the detection of EGFR and KRAS mutations in a panel of 200 NSCLC samples. Seventy-five EGFR mutations and 25 KRAS mutations were identified by HRM analysis in concordance with sequencing results. The mutation types for each exon were similar to those reported in previous studies, in-frame deletions in exon 19, insertions in exon 20 and missense mutations in both exon 18 and 21 [7-9]. However, the overall mutation rate (37.5%) in this Australian study was much higher than in North America and Western Europe (10%) but similar to East Asian populations (30–50%) [43-47]. Although most Australians are of European descent, the higher EGFR mutation rate is likely to be due to selection bias of patients by referring oncologists for features associated with EGFR mutation. The pre-selection criteria included tumours that were histologically adenocarcinomas (70% of samples tested) and patients that were female (60% of patients tested) and/or of Asian ethnicity as many of the patients had a name that was consistent with Asian ancestry. The KRAS mutation rate (12.5%) which is lower than that reported in previous studies [16,17] also supports the notion of pre-selection bias.

Although all positive sequencing results were detected by HRM, some samples were considered positive by HRM but were negative by sequencing (Table (Table3).3). There are several possible explanations. One possible explanation is that the adverse effects of formalin fixation on DNA can cause PCR artefacts during amplification. At least four chemical reactions occur between formaldehyde and DNA; methylol formation, methylene bridge formation, apurinic and apyrimidinic site formation, and hydrolysis of the phosphodiester bonds [48]. Compared to the DNA extracted from frozen tissues, a higher frequency of non-reproducible sequence alterations have been reported with DNA isolated from the formalin fixed tissues [49,50]. Therefore, the cumulative effects of the PCR artefacts either from Taq polymerase error and errors attributed from chemical reactions of formalin on DNA influence the melting profile of the amplicon depending on the degree of DNA damages.

Seven samples gave positive results in more than three EGFR HRM assays, supporting the hypothesis that the quality of DNA is one of the factors causing aberrant variation of melting in HRM analysis. It has been observed that the wild-type variation in melting analysis is much greater with FFPE DNA than DNA from frozen tissues or peripheral blood (data not shown). The amount of false positives decreased with decreasing amplicon length, with EGFR exon 19 (amplicon size 250 bp) giving the most false positives and KRAS exon 2 (92 bp) giving the least.

Another possibility is that some samples contained levels of mutation below the sensitivity of sequencing detection as a result of a low percentage of tumour in the sample or genetic heterogeneity within the tumour. Where sequencing requires the mutation to be present at a level of 20% of the sample, HRM can detect heterozygous genetic changes down to 10% or below [26,29,37]. We are now adopting digital techniques to confirm that some specimens have true mutations present at low levels. Other HRM false positives can arise from PCR errors due to amplification from very low levels of template. True mutations can be distinguished from artefacts by confirming the identical sequence variations from independent amplification (Do and Dobrovic, manuscript submitted).

HRM is a suitable methodology to test FFPE samples as well as samples with a very low quantity of DNA. It has been reported that ten percent buffered formalin, an aqueous dilution of formaldehyde, can interact with DNA and initiate irreversible DNA degradation resulting in an adverse effect on DNA quality [35]. In our HRM assays, all the samples were successfully amplified and analysed. Our results show that HRM can be performed with as little as 1 ng template in the EGFR exon 19 HRM assay. This level of sensitivity of HRM, together with the possibility of sequencing of the HRM product, will extend our ability to screen even clinical samples with extremely low DNA quantity such as samples taken from patients with inoperable tumours. HRM analysis can now provide a genetic testing option for these patients, in which the results might prove useful in directing treatment and may ultimately improve outcomes.

EGFR and KRAS mutations are predominantly mutually exclusive with very rare tumours containing both genes mutated [44,51]. The coexistence of mutations in both EGFR and KRAS has only been reported in two patients [52]. In the current study, all 25 KRAS positive samples were wild-type for EGFR, supporting the general mutual exclusiveness of the two mutations.

Conclusion

We have established a fast, efficient and reproducible screening method for EGFR and KRAS mutation detection by which degraded DNA from FFPE tissues can be tested. This is the first study in which the HRM method has been properly validated as a scanning method for EGFR and KRAS mutations in a large sample set showing accurate mutation detection with 100% sensitivity. It is estimated that up to 80% of sequencing reactions can be eliminated for these two genes if samples are screened by HRM. These results further demonstrate the utility of HRM for the detection of somatic mutations in clinical samples and for screening of samples prior to sequencing.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HD participated in the development of the assays carried out the HRM studies, prepared the figures and wrote the manuscript. MK assisted in development of the assays and in data analysis. SBF contributed to the manuscript and provided access to the sequencing data and to the specimens. This study was initiated as part of a grant to PLM and AD. AD initiated the study, was responsible for primer design, participated in the data analysis, and co-wrote the manuscript. All authors read and approved the final manuscript.

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2407/8/142/prepub

Acknowledgements

We thank Sergey Kovalenko and the Molecular Pathology Diagnostic Lab, Peter MacCallum Cancer Centre, for DNA preparation, sequencing, and sequencing analysis for this study. Also we would like to acknowledge the Molecular Pathology Research and Development Lab, for invaluable discussion regarding HRM assay development and analysis. Chelsee Hewitt, Renato Salemi and Lasse Kristensen critically read the manuscript. This project was supported by funding from the National Health and Medical Research Council to PLM and AD.

References

Sasaki H, Endo K, Konishi A, Takada M, Kawahara M, Iuchi K, Matsumura A, Okumura M, Tanaka H, Kawaguchi T, Shimizu T, Takeuchi H, Yano M, Fukai I, Fujii Y. EGFR Mutation status in Japanese lung cancer patients: genotyping analysis using LightCycler. Clin Cancer Res. 2005;11:2924–2929. doi: 10.1158/1078-0432.CCR-04-1904. [PubMed] [Cross Ref]
Schiller JH, Harrington D, Belani CP, Langer C, Sandler A, Krook J, Zhu J, Johnson DH. Comparison of four chemotherapy regimens for advanced non-small-cell lung cancer. N Engl J Med. 2002;346:92–98. doi: 10.1056/NEJMoa011954. [PubMed] [Cross Ref]
Nicholson RI, Gee JM, Harper ME. EGFR and cancer prognosis. Eur J Cancer. 2001;37 Suppl 4:S9–15. doi: 10.1016/S0959-8049(01)00231-3. [PubMed] [Cross Ref]
Ohsaki Y, Tanno S, Fujita Y, Toyoshima E, Fujiuchi S, Nishigaki Y, Ishida S, Nagase A, Miyokawa N, Hirata S, Kikuchi K. Epidermal growth factor receptor expression correlates with poor prognosis in non-small cell lung cancer patients with p53 overexpression. Oncol Rep. 2000;7:603–607. [PubMed]
Mendelsohn J, Baselga J. The EGF receptor family as targets for cancer therapy. Oncogene. 2000;19:6550–6565. doi: 10.1038/sj.onc.1204082. [PubMed] [Cross Ref]
Scaltriti M, Baselga J. The epidermal growth factor receptor pathway: a model for targeted therapy. Clin Cancer Res. 2006;12:5268–5272. doi: 10.1158/1078-0432.CCR-05-1554. [PubMed] [Cross Ref]
Lynch TJ, Bell DW, Sordella R, Gurubhagavatula S, Okimoto RA, Brannigan BW, Harris PL, Haserlat SM, Supko JG, Haluska FG, Louis DN, Christiani DC, Settleman J, Haber DA. Activating mutations in the epidermal growth factor receptor underlying responsiveness of non-small-cell lung cancer to gefitinib. N Engl J Med. 2004;350:2129–2139. doi: 10.1056/NEJMoa040938. [PubMed] [Cross Ref]
Pao W, Miller V, Zakowski M, Doherty J, Politi K, Sarkaria I, Singh B, Heelan R, Rusch V, Fulton L, Mardis E, Kupfer D, Wilson R, Kris M, Varmus H. EGF receptor gene mutations are common in lung cancers from "never smokers" and are associated with sensitivity of tumors to gefitinib and erlotinib. Proc Natl Acad Sci U S A. 2004;101:13306–13311. doi: 10.1073/pnas.0405220101. [PubMed] [Cross Ref]
Paez JG, Janne PA, Lee JC, Tracy S, Greulich H, Gabriel S, Herman P, Kaye FJ, Lindeman N, Boggon TJ, Naoki K, Sasaki H, Fujii Y, Eck MJ, Sellers WR, Johnson BE, Meyerson M. EGFR mutations in lung cancer: correlation with clinical response to gefitinib therapy. Science. 2004;304:1497–1500. doi: 10.1126/science.1099314. [PubMed] [Cross Ref]
Han SW, Kim TY, Hwang PG, Jeong S, Kim J, Choi IS, Oh DY, Kim JH, Kim DW, Chung DH, Im SA, Kim YT, Lee JS, Heo DS, Bang YJ, Kim NK. Predictive and prognostic impact of epidermal growth factor receptor mutation in non-small-cell lung cancer patients treated with gefitinib. J Clin Oncol. 2005;23:2493–2501. doi: 10.1200/JCO.2005.01.388. [PubMed] [Cross Ref]
Shigematsu H, Lin L, Takahashi T, Nomura M, Suzuki M, Wistuba. Fong KM, Lee H, Toyooka S, Shimizu N, Fujisawa T, Feng Z, Roth JA, Herz J, Minna JD, Gazdar AF. Clinical and biological features associated with epidermal growth factor receptor gene mutations in lung cancers. J Natl Cancer Inst. 2005;97:339–346. [PubMed]
Amann J, Kalyankrishna S, Massion PP, Ohm JE, Girard L, Shigematsu H, Peyton M, Juroske D, Huang Y, Stuart Salmon J, Kim YH, Pollack JR, Yanagisawa K, Gazdar A, Minna JD, Kurie JM, Carbone DP. Aberrant epidermal growth factor receptor signaling and enhanced sensitivity to EGFR inhibitors in lung cancer. Cancer Res. 2005;65:226–235. [PubMed]
Engelman JA, Janne PA, Mermel C, Pearlberg J, Mukohara T, Fleet C, Cichowski K, Johnson BE, Cantley LC. ErbB-3 mediates phosphoinositide 3-kinase activity in gefitinib-sensitive non-small cell lung cancer cell lines. Proc Natl Acad Sci U S A. 2005;102:3788–3793. doi: 10.1073/pnas.0409773102. [PubMed] [Cross Ref]
Greulich H, Chen TH, Feng W, Janne PA, Alvarez JV, Zappaterra M, Bulmer SE, Frank DA, Hahn WC, Sellers WR, Meyerson M. Oncogenic transformation by inhibitor-sensitive and -resistant EGFR mutants. PLoS Med. 2005;2:e313. doi: 10.1371/journal.pmed.0020313. [PMC free article] [PubMed] [Cross Ref]
Sordella R, Bell DW, Haber DA, Settleman J. Gefitinib-sensitizing EGFR mutations in lung cancer activate anti-apoptotic pathways. Science. 2004;305:1163–1167. doi: 10.1126/science.1101637. [PubMed] [Cross Ref]
Adjei AA. Blocking oncogenic Ras signaling for cancer therapy. J Natl Cancer Inst. 2001;93:1062–1074. doi: 10.1093/jnci/93.14.1062. [PubMed] [Cross Ref]
Bos JL. ras oncogenes in human cancer: a review. Cancer Res. 1989;49:4682–4689. [PubMed]
van Zandwijk N, Mathy A, Boerrigter L, Ruijter H, Tielen I, de Jong D, Baas P, Burgers S, Nederlof P. EGFR and KRAS mutations as criteria for treatment with tyrosine kinase inhibitors: retro- and prospective observations in non-small-cell lung cancer. Ann Oncol. 2007;18:99–103. doi: 10.1093/annonc/mdl323. [PubMed] [Cross Ref]
Uchida A, Hirano S, Kitao H, Ogino A, Rai K, Toyooka S, Takigawa N, Tabata M, Takata M, Kiura K, Tanimoto M. Activation of downstream epidermal growth factor receptor (EGFR) signaling provides gefitinib-resistance in cells carrying EGFR mutation. Cancer Sci. 2007;98:357–363. doi: 10.1111/j.1349-7006.2007.00387.x. [PubMed] [Cross Ref]
Kim YT, Kim TY, Lee DS, Park SJ, Park JY, Seo SJ, Choi HS, Kang HJ, Hahn S, Kang CH, Sung SW, Kim JH. Molecular changes of epidermal growth factor receptor (EGFR) and KRAS and their impact on the clinical outcomes in surgically resected adenocarcinoma of the lung. Lung Cancer. 2007.
Tsao MS, Sakurada A, Cutz JC, Zhu CQ, Kamel-Reid S, Squire J, Lorimer I, Zhang T, Liu N, Daneshmand M, Marrano P, da Cunha Santos G, Lagarde A, Richardson F, Seymour L, Whitehead M, Ding K, Pater J, Shepherd FA. Erlotinib in lung cancer - molecular and clinical predictors of outcome. N Engl J Med. 2005;353:133–144. doi: 10.1056/NEJMoa050736. [PubMed] [Cross Ref]
Chin TM, Anuar D, Soo R, Salto-Tellez M, Li WQ, Ahmad B, Lee SC, Goh BC, Kawakami K, Segal A, Iacopetta B, Soong R. Detection of epidermal growth factor receptor variations by partially denaturing HPLC. Clin Chem. 2007;53:62–70. doi: 10.1373/clinchem.2006.074831. [PubMed] [Cross Ref]
Janne PA, Borras AM, Kuang Y, Rogers AM, Joshi VA, Liyanage H, Lindeman N, Lee JC, Halmos B, Maher EA, Distel RJ, Meyerson M, Johnson BE. A rapid and sensitive enzymatic method for epidermal growth factor receptor mutation screening. Clin Cancer Res. 2006;12:751–758. doi: 10.1158/1078-0432.CCR-05-2047. [PubMed] [Cross Ref]
Cohen V, Agulnik JS, Jarry J, Batist G, Small D, Kreisman H, Tejada NA, Miller WH, Jr., Chong G. Evaluation of denaturing high-performance liquid chromatography as a rapid detection method for identification of epidermal growth factor receptor mutations in nonsmall-cell lung cancer. Cancer. 2006;107:2858–2865. doi: 10.1002/cncr.22331. [PubMed] [Cross Ref]
Reed GH, Kent JO, Wittwer CT. High-resolution DNA melting analysis for simple and efficient molecular diagnostics. Pharmacogenomics. 2007;8:597–608. doi: 10.2217/14622416.8.6.597. [PubMed] [Cross Ref]
Krypuy M, Ahmed AA, Etemadmoghadam D, Hyland SJ, Group AO, Brenton JD, Fox SB, Defazio A, Bowtell DD, Dobrovic A. High resolution melting for mutation scanning of TP53 exons 5-8. BMC Cancer. 2007;7:168. doi: 10.1186/1471-2407-7-168. [PMC free article] [PubMed] [Cross Ref]
Holden JA, Willmore-Payne C, Coppola D, Garrett CR, Layfield LJ. High-resolution melting amplicon analysis as a method to detect c-kit and platelet-derived growth factor receptor alpha activating mutations in gastrointestinal stromal tumors. Am J Clin Pathol. 2007;128:230–238. doi: 10.1309/7TEH56K6WWXENNQY. [PubMed] [Cross Ref]
Willmore-Payne C, Holden JA, Layfield LJ. Detection of epidermal growth factor receptor and human epidermal growth factor receptor 2 activating mutations in lung adenocarcinoma by high-resolution melting amplicon analysis: correlation with gene copy number, protein expression, and hormone receptor expression. Hum Pathol. 2006;37:755–763. doi: 10.1016/j.humpath.2006.02.004. [PubMed] [Cross Ref]
Nomoto K, Tsuta K, Takano T, Fukui T, Fukui T, Yokozawa K, Sakamoto H, Yoshida T, Maeshima AM, Shibata T, Furuta K, Ohe Y, Matsuno Y. Detection of EGFR mutations in archived cytologic specimens of non-small cell lung cancer using high-resolution melting analysis. Am J Clin Pathol. 2006;126:608–615. doi: 10.1309/N5PQNGW2QKMX09X7. [PubMed] [Cross Ref]
Wojdacz TK, Dobrovic A. Methylation-sensitive high resolution melting (MS-HRM): a new approach for sensitive and high-throughput assessment of methylation. Nucleic Acids Res. 2007;35:e41. doi: 10.1093/nar/gkm013. [PMC free article] [PubMed] [Cross Ref]
White HE, Hall VJ, Cross NC. Methylation-Sensitive High-Resolution Melting-Curve Analysis of the SNRPN Gene as a Diagnostic Screen for Prader-Willi and Angelman Syndromes. Clin Chem. 2007.
Chou LS, Lyon E, Wittwer CT. A comparison of high-resolution melting analysis with denaturing high-performance liquid chromatography for mutation scanning: cystic fibrosis transmembrane conductance regulator gene as a model. Am J Clin Pathol. 2005;124:330–338. doi: 10.1309/BF3M-LJN8-J527-MWQY. [PubMed] [Cross Ref]
Willmore-Payne C, Holden JA, Hirschowitz S, Layfield LJ. BRAF and c-kit gene copy number in mutation-positive malignant melanoma. Hum Pathol. 2006;37:520–527. doi: 10.1016/j.humpath.2006.01.003. [PubMed] [Cross Ref]
von Ahlfen S, Missel A, Bendrat K, Schlumpberger M. Determinants of RNA Quality from FFPE Samples. PLoS ONE. 2007;2:e1261. doi: 10.1371/journal.pone.0001261. [PMC free article] [PubMed] [Cross Ref]
Koshiba M, Ogawa K, Hamazaki S, Sugiyama T, Ogawa O, Kitajima T. The effect of formalin fixation on DNA and the extraction of high-molecular-weight DNA from fixed and embedded tissues. Pathol Res Pract. 1993;189:66–72. [PubMed]
The thermal denaturation profile by Poland's algorithm http://www.biophys.uni-duesseldorf.de/local/POLAND/poland.html
Krypuy M, Newnham GM, Thomas DM, Conron M, Dobrovic A. High resolution melting analysis for the rapid and sensitive detection of mutations in clinical samples: KRAS codon 12 and 13 mutations in non-small cell lung cancer. BMC Cancer. 2006;6:295. doi: 10.1186/1471-2407-6-295. [PMC free article] [PubMed] [Cross Ref]
Yokoyama T, Kondo M, Goto Y, Fukui T, Yoshioka H, Yokoi K, Osada H, Imaizumi K, Hasegawa Y, Shimokata K, Sekido Y. EGFR point mutation in non-small cell lung cancer is occasionally accompanied by a second mutation or amplification. Cancer Sci. 2006;97:753–759. doi: 10.1111/j.1349-7006.2006.00233.x. [PubMed] [Cross Ref]
Takano EA, Mitchell G, Fox SB, Dobrovic A. Rapid detection of carriers with BRCA1 and BRCA2 mutations using high resolution melting analysis. BMC Cancer. 2008;8:59. doi: 10.1186/1471-2407-8-59. [PMC free article] [PubMed] [Cross Ref]
Takano T, Ohe Y, Tsuta K, Fukui T, Sakamoto H, Yoshida T, Tateishi U, Nokihara H, Yamamoto N, Sekine I, Kunitoh H, Matsuno Y, Furuta K, Tamura T. Epidermal growth factor receptor mutation detection using high-resolution melting analysis predicts outcomes in patients with advanced non small cell lung cancer treated with gefitinib. Clin Cancer Res. 2007;13:5385–5390. doi: 10.1158/1078-0432.CCR-07-0627. [PubMed] [Cross Ref]
Mullins FM, Dietz L, Lay M, Zehnder JL, Ford J, Chun N, Schrijver I. Identification of an intronic single nucleotide polymorphism leading to allele dropout during validation of a CDH1 sequencing assay: implications for designing polymerase chain reaction-based assays. Genet Med. 2007;9:752–760. [PubMed]
EGFR mutation batabase http://www.cityofhope.org/cmdl/egfr_db
Haneda H, Sasaki H, Lindeman N, Kawano O, Endo K, Suzuki E, Shimizu S, Yukiue H, Kobayashi Y, Yano M, Fujii Y. A correlation between EGFR gene mutation status and bronchioloalveolar carcinoma features in Japanese patients with adenocarcinoma. Jpn J Clin Oncol. 2006;36:69–75. doi: 10.1093/jjco/hyi228. [PubMed] [Cross Ref]
Marchetti A, Martella C, Felicioni L, Barassi F, Salvatore S, Chella A, Camplese PP, Iarussi T, Mucilli F, Mezzetti A, Cuccurullo F, Sacco R, Buttitta F. EGFR mutations in non-small-cell lung cancer: analysis of a large series of cases and development of a rapid and sensitive method for diagnostic screening with potential implications on pharmacologic treatment. J Clin Oncol. 2005;23:857–865. doi: 10.1200/JCO.2005.08.043. [PubMed] [Cross Ref]
Janne PA, Engelman JA, Johnson BE. Epidermal growth factor receptor mutations in non-small-cell lung cancer: implications for treatment and tumor biology. J Clin Oncol. 2005;23:3227–3234. doi: 10.1200/JCO.2005.09.985. [PubMed] [Cross Ref]
Cappuzzo F, Hirsch FR, Rossi E, Bartolini S, Ceresoli GL, Bemis L, Haney J, Witta S, Danenberg K, Domenichini I, Ludovini V, Magrini E, Gregorc V, Doglioni C, Sidoni A, Tonato M, Franklin WA, Crino L, Bunn PA, Jr., Varella-Garcia M. Epidermal growth factor receptor gene and protein and gefitinib sensitivity in non-small-cell lung cancer. J Natl Cancer Inst. 2005;97:643–655. [PubMed]
Kosaka T, Yatabe Y, Endoh H, Kuwano H, Takahashi T, Mitsudomi T. Mutations of the epidermal growth factor receptor gene in lung cancer: biological and clinical implications. Cancer Res. 2004;64:8919–8923. doi: 10.1158/0008-5472.CAN-04-2818. [PubMed] [Cross Ref]
Srinivasan M, Sedmak D, Jewell S. Effect of fixatives and tissue processing on the content and integrity of nucleic acids. Am J Pathol. 2002;161:1961–1971. [PubMed]
Quach N, Goodman MF, Shibata D. In vitro mutation artifacts after formalin fixation and error prone translesion synthesis during PCR. BMC Clin Pathol. 2004;4:1. doi: 10.1186/1472-6890-4-1. [PMC free article] [PubMed] [Cross Ref]
Williams C, Ponten F, Moberg C, Soderkvist P, Uhlen M, Ponten J, Sitbon G, Lundeberg J. A high frequency of sequence alterations is due to formalin fixation of archival specimens. Am J Pathol. 1999;155:1467–1471. [PubMed]
Gazdar AF, Shigematsu H, Herz J, Minna JD. Mutations and addiction to EGFR: the Achilles 'heal' of lung cancers? Trends Mol Med. 2004;10:481–486. doi: 10.1016/j.molmed.2004.08.008. [PubMed] [Cross Ref]
Han SW, Kim TY, Jeon YK, Hwang PG, Im SA, Lee KH, Kim JH, Kim DW, Heo DS, Kim NK, Chung DH, Bang YJ. Optimization of patient selection for gefitinib in non-small cell lung cancer by combined analysis of epidermal growth factor receptor mutation, K-ras mutation, and Akt phosphorylation. Clin Cancer Res. 2006;12:2538–2544. doi: 10.1158/1078-0432.CCR-05-2845. [PubMed] [Cross Ref]

Articles from BMC Cancer are provided here courtesy of
BioMed Central