This article has Open Peer Review reports available.
Development of a semi-conductor sequencing-based panel for genotyping of colon and lung cancer by the Onconetwork consortium
- Bastiaan BJ Tops1Email author,
- Nicola Normanno2, 3Email author,
- Henriette Kurth4,
- Eliana Amato5,
- Andrea Mafficini5,
- Nora Rieber6,
- Delphine Le Corre7,
- Anna Maria Rachiglio2,
- Anne Reiman8,
- Orla Sheils9,
- Christoph Noppen4,
- Ludovic Lacroix10,
- Ian A Cree8,
- Aldo Scarpa5, 11,
- Marjolijn JL Ligtenberg1, 12 and
- Pierre Laurent-Puig7
© Tops et al.; licensee BioMed Central. 2015
Received: 1 July 2014
Accepted: 8 January 2015
Published: 31 January 2015
The number of predictive biomarkers that will be necessary to assess in clinical practice will increase with the availability of drugs that target specific molecular alterations. Therefore, diagnostic laboratories are confronted with new challenges: costs, turn-around-time and the amount of material required for testing will increase with the number of tests performed on a sample. Our consortium of European clinical research laboratories set out to test if semi-conductor sequencing provides a solution for these challenges.
We designed a multiplex PCR targeting 87 hotspot regions in 22 genes that are of clinical interest for lung and/or colorectal cancer. The gene-panel was tested by 7 different labs in their own clinical setting using ion-semiconductor sequencing.
We analyzed 155 samples containing 112 previously identified mutations in the KRAS, EGFR en BRAF genes. Only 1 sample failed analysis due to poor quality of the DNA. All other samples were correctly genotyped for the known mutations, even as low as 2%, but also revealed other mutations. Optimization of the primers used in the multiplex PCR resulted in a uniform coverage distribution over the amplicons that allows for efficient pooling of samples in a sequencing run.
We show that a semi-conductor based sequencing approach to stratify colon and lung cancer patients is feasible in a clinical setting.
Although molecular mechanisms underlying carcinogenesis have been the subject of medical research for decades, only in the last few years has this information translated into the clinical setting through ‘smart’ drugs targeting specific molecular aberrations of neoplastic cells . Target based agents were initially explored in unselected cohorts of cancer patients, and this approach led to failure of the majority of clinical trials with these new agents. More recently, the value of selecting patients to be treated with targeted agents on the basis of specific molecular alterations of neoplastic cells has been acknowledged [2,3]. This new therapeutic approach has led to the approval of novel drugs for molecularly selected patient populations. These agents have revolutionized the treatment of cancer patients, improving treatment outcome, while minimizing side effects.
The number of predictive biomarkers that are assessed in clinical practice is rapidly increasing with the availability of drugs that target specific molecular alterations . In some tumor types, such as non-small-cell lung carcinoma (NSCLC) and colorectal carcinoma (CRC), a new classification based on the identification of specific predictive and/or prognostic molecular alterations has emerged in the past decade [5,6]. However, molecular diagnostics is faced with its own specific challenges. Costs per test, turn-around-time and the amount of material required for testing increase with the number of tests performed on a sample. In particular, the availability of tissue is limiting in NSCLC, which is often diagnosed on small histological or cytological samples yielding limited amounts of DNA .
Comprehensive molecular characterization of tumor tissue in clinical practice will rely on the development of high throughput technologies allowing this process to be accomplished in a cost-effective and timely manner and by using a limited amount of tissue. Two different approaches, genotyping assays and next generation sequencing (NGS) have been proposed for molecular screening [4,8]. NGS has the advantage to providing information on known and novel molecular alterations. Indeed, by using NGS-based techniques multiple genes can be sequenced simultaneously, and benchtop sequencers, like the Ion PGM™ (Life Technologies), are fast and relatively inexpensive. In this context, the ‘OncoNetwork consortium’, a European collaborative effort developed and tested the performance of a next-generation sequencing approach in a clinical setting. The advantages of a multi-lab collaboration are that methods can be quickly validated, facilitating availability, and that tests are evaluated in different laboratory settings, assessing robustness. As a proof-of-principle we chose to focus on CRC and NSCLC, two frequent tumor types for which there are known mutations associated with treatment decisions. We decided to design an Ion AmpliSeq custom gene-panel as a single multiplex PCR requiring only 10 ng of input DNA for the detection of the most frequent hot spot mutations in NSCLC and CRC. We opted for a relatively small, focused gene-panel, instead of a broader panel (like the AmpliSeq Cancer panel). Since a smaller gene-panel requires less sequence capacity more samples can be pooled in a single sequence run resulting in lower costs per sample and an increased throughput compatible with a diagnostic activity. Finally, we assessed the performance of the panel by semiconductor sequencing of 155 CRC and NSCLC formalin-fixed, paraffin-embedded (FFPE) cancer samples containing mutations previously identified by alternative genotyping methods.
Tumor specimens and DNA isolation
FFPE tumor samples were retrieved from the archives of the collaborating institutes and anonymised. Given the study was interpreted by all institutions (Radboud university medical centre, Istituto di ricovero e cura a carattere scientifico (IRCCS), VIOLLIER, Applied Research On Cancer Centre (ARC-NET), Institut national de la santé et de la recherche médicale (INSERM), Trinity College, Institut Gustave-Roussy and University Hospitals Coventry and Warwickshire) as a service improvement, it did not require specific research ethics committee approval as stated in the EU Clinical Trials Directive (2001/20/EC). Specific approval was obtained from the SJH/AMNCH Joint Research Ethics Committee of Trinity College, Dublin. Approval was procured for the current study as part of a wider study involving molecular analysis of a wider cohort and for research purposes that does require approval from an ethics committee. A cover letter to this effect is available.
Tissue areas for DNA extraction were micro-dissected and quantified for the percentage of neoplastic cells by a pathologist. More details regarding sample characteristics are depicted in Additional file 1: Figure S1. DNA extraction differed per institute. Most laboratories (6/7) isolated DNA using the QIAamp DNA FFPE Tissue Kit, while one lab purified DNA by overnight digestion in a proteinase K/Chelex-100 solution.
Ion AmpliSeq™ primer design
Ion AmpliSeq™ primer pools were designed by Life Technologies (Additional file 2: Table S1 and Additional file 3: Table S2). To minimize the risk of allelic drop-out during the PCR due to primers that are not able to prime, no confirmed SNPs with a frequency >0.5% were allowed in the 5 most 3’ nucleotides. For the same reason low-frequency SNPs at other positions were kept as low as possible, with a maximum of four. Primer designs were confirmed by SNPCheck v3 analysis (www.snpcheck.net/, April 2012). The resulting oncopanel primer pool contained 87 primer pairs with an average amplicon length of 203 bp (Additional file 2: Table S1).
AmpliSeq enrichment and Ion Torrent sequencing
Library generation was performed according to the manufacturer’s protocol. In short, 10 ng of DNA per pool was amplified in 21 cycles by PCR using the Ion AmpliSeq™ mastermix, followed by barcode and adapter ligation. Amplified products were purified with Agencourt AMPure XP beads (Beckman Coulter Genomics, High Wycombe, UK). The library was diluted to 20 pM. Emulsion PCR was performed using the Ion OneTouch™ 200 Template kit following the protocol of the Ion OneTouch™ System. Next, Ion Sphere Particles (ISPs) were recovered and enriched for template positive ISPs using Dynabeads MyOne Streptavidin C1 beads (Life Technologies) in the Ion OneTouch™ ES instrument (Life Technologies). ISP enrichment was quantified using the Qubit 2.0 fluorometer (Life Technologies). Sequencing primer and polymerase were added to the final enriched spheres before loading onto an Ion 316 chip according to the Ion PGM™ 200 sequencing kit protocol.
Sequence data were directly uploaded from the Ion Torrent Suite Software to the IonReporter™ Software version 1.2 and data were analyzed using a pre-configured workflow specific for detecting low frequency somatic variants. The allele frequency cut-off was set at 0.04. In addition, the data were analyzed in parallel using a ‘hotspot bed file’ containing hotspot mutations in the regions of our gene-panel derived from the COSMIC database. The allele frequency cut-off for hotspot analysis was set at 0.02. These cutoffs were chosen following a preliminary phase in which, starting from samples with known sequence, the best bioinformatic algorithm ensuring the highest sensitivity and specificity was identified. The workflow includes mapping statistics, measurements of the quality of all reads, on and off target statistics, SNP and INDEL calling and comprehensive annotation using public databases (dbSNP version 137, COSMIC version 60, OMIM version 08012012, Gene Ontology version 1.218). The variant caller algorithm was used and variants that were identified for each sample per lab in each phase were further filtered. We applied quality filtering by discarding reads with a quality score <40 and checking against strand bias. We filtered against variants with a variant allele count less than 10, variants that showed up in every sample of one lab with the same allele frequency, the synonymous and non-coding variants and variants that were known as polymorphisms or with a Minor Allele Frequency (MAF) greater than 35%. As a result, a variant knowledgebase was generated with the variants detected in those samples. A specific workflow for this panel was configured in the online Ion Reporter™ Software and was used to perform analysis for data generated using the Colon/Lung panel.
Panel design and performance
Mutation frequency of genes for which hotspot regions are present in the gene-panel
% Mutated CRC#
% Mutated lung cancer#
Dilution experiments with an artificial control sample containing 10 hotspot mutations in the BRAF, EGFR, KRAS, NRAS and PIK3CA genes (Horizon Diagnostics) demonstrated that hotspot mutations could be confidentially identified as low as 2% mutant alleles provided that the coverage was >500× (results not shown). Preliminary testing of the lung/colon cancer primer pool showed that up to 5 samples could be pooled on a 316 chip (Ion PGM™ Sequencer) with a minimal average read-depth of 500× (results not shown).
Variants identified in the 5 control samples
Mutations identified by NGS‡
0.44 - 0.54
0.66 - 0.67
0.46 - 0.49
0.74 - 0.75
0.50 - 0.56
0.43 - 0.46
0.47 - 0.56
0.50 - 0.66
0.48 - 0.54
0.65 - 0.75
0.40 - 0.51
0.54 - 0.69
0.44 - 0.46
0.46 - 0.55
0.32 - 0.35
By using the panel, all 7 labs also identified 6 new variants in the two AcroMetrix® control samples (Table 2). To confirm that these new variants are not sequencing artifacts, we re-analyzed these positions by conventional Sanger sequencing. All 6 new variants were confirmed suggesting the specificity of our workflow (data not shown). The two xenograft samples were not analyzed for new variants since off-target amplification of mouse DNA hampered data analysis.
Known mutations present in the 60 samples that were analyzed in the blind during phase 2 of the panel validation
Thus, combining the data of all experiments, we analyzed 155 unique samples containing 112 previously identified mutations. All specimens were correctly genotyped, except one that could not be analyzed due to poor quality of the DNA, demonstrating excellent performance of the gene panel.
Sequencing efficiency and panel re-design
To check the performance of the panel, the 5 control samples of phase 1 were-resequenced by all 7 laboratories, and the expected variants as well as the new variants identified in the two Accrometrix control samples were detected using panel v2. In addition, the new panel was used by 5 labs to re-sequence their 15 in-house samples (75 samples in total) previously tested with panel v1, and 6 labs also tested 5 samples not previously analyzed (not all labs were able to participate). All mutations previously detected with panel v1 were again identified using panel v2 (data not shown), demonstrating that this panel performs as well as v1 for these mutations. However, since 8 samples can be multiplexed on a 316 chip with panel v2, instead of 5 samples using v1, the sequencing costs per sample are significantly reduced.
Specification of the molecular class of NSCLC and CRC is mandatory to optimize personalized medicine for these diseases. However, a comprehensive molecular characterization of the tumor is faced with a number of challenges, including cost per test, turn-around-time and the limited amount of material available for genotyping. NGS-based approaches can overcome these challenges by providing in a single analysis information on tumor inter- and intra-heterogeneity that can be relevant in a clinical setting. In this respect, a number of NGS panels for tumor profiling have become commercially available. However, these panels lack of validation in clinical samples and their inter-laboratory reproducibility has not been shown. In this study we propose a novel model of collaboration between clinical research laboratories and NGS-companies that led to the development of a tool that is ready for application in a clinical setting. Indeed, the validation of panels through a consortium of academic institutions represents a novel approach to speed the development of new diagnostic tools, being able to provide information on a significant number of clinical specimens and on inter-laboratory reproducibility.
We must acknowledge that there are still several issues that need to be addressed to implement NGS panels in the clinic. Our panel included 22 genes, although the approved predictive biomarkers in colon and lung cancer are quite few (KRAS, NRAS, EGFR, ALK). Nevertheless, in many cancer comprehensive centers clinical trials with new drugs are open and molecular pathologists are requested to screen for a number of biomarkers that are included in the Colon and Lung cancer panel that we developed. For example, among the 61 difficult NSCLC samples analyzed in the third phase, we identified 18 EGFR, 7 PIK3CA and 3 BRAF mutants. Since drugs targeting these latter molecular alterations are in clinical trials, the availability of this information might significantly improve the possibility of these patients to receive a personalized therapy . It is possible that novel biomarkers will be discovered in the future. In this respect, one advantage of this panel is that it can be easily adapted by adding new amplicons/genes of interest. Of course, a re-validation on a minimum number of samples is recommended for any modification of the panel.
The performance of a new test is a key element in evaluating the possibility to introduce it in clinical diagnostics. In this regard, we chose to analyze our data with an allele frequency cut-off of 4% for all variants and 2% for known hotspot mutations corresponding to 5-10% neoplastic cells carrying heterozygous mutations. This estimate is in agreement with a previous study that assessed somatic mutations in FFPE material by using semiconductor-based massive parallel sequencing . The allele frequency cut-off required for the detection of somatic mutations is debatable, but samples with <10% neoplastic cells are also difficult to diagnose by classical histomorphological methods. Using this allele-frequency cut-off we tested 155 clinical FFPE samples that were previously analyzed by conventional methods. A substantial number of these specimens were specifically selected based on the fact that they were difficult to analyze with routine techniques. Even with this bias in sample selection we were still able to correctly analyze 154 of the 155 samples. Although our sample set was selected for certain mutations, our specimens contained variants that represent the full spectrum of mutations that are identified in NSCLC and CRC specimens. These include a variety of indel and missense mutations of which also hotspot mutations present in only 2% of the alleles were readily detected. Other studies have previously suggested that semiconductor-based sequencing of FFPE tissue is indeed feasible [18-20]. However, this is the first study that evaluated a so large collection of FFPE samples.
Other mutation-detection techniques, like quantitative PCR (qPCR) or high resolution melting (HRM), can also reach this level of sensitivity, but lack information regarding the allele-frequency of the mutation. For example, EGFR mutations showed allelic frequencies ranging between 20% and 210% (amplification), after normalization for the percentage of neoplastic cells. This information could be clinically relevant in the near future as suggested by a report showing that the duration of the response to target based agents might be related to the relative content of mutant alleles in the tumor . A quantitative assessment of the different mutant clones in polyclonal tumors will also be useful for treatment decision, as recently shown for EGFR sensitizing and resistance mutations in NSCLC . In addition, genotyping techniques can only provide information on known mutations whereas sequencing-based methods can also detect novel and rare mutations. Some of the new variants we identified might not be true somatic mutations. Variants with high allelic frequencies could also represent infrequent SNPs, while especially non-hotspot variants with low allelic frequencies could be false positive calls. These are difficult to confirm however with standard methods like Sanger sequencing. It must be emphasized that during this project we used IonReporter Software version 1.2, while currently version 1.6 is available with improved mutation calling, especially in homopolymer regions.
Conventional methods, like Sanger sequencing or qPCR, are singleplex tests and, as the number of targets increases, the amount of DNA required for multiple tests is a limiting factor. This can be problematic especially in lung cancer, since molecular analyses are often performed on fine needle aspiration (FNA) biopsies or cytology samples yielding limited amounts of DNA. Our single multiplex PCR, consisting of 87 amplicons targeting 22 genes and requiring as little as 10 ng of input DNA, overcomes this problem. Analyzing a similar number of amplicons by conventional methods would require 500–1000 ng of DNA, clearly demonstrating the advantage of NGS.
Turn-around-time (TAT) is an important variable to take into account in clinical diagnostics. We opted for the IonPGM in our institutes since benchtop sequencers, like the Ion PGM or the MiSeq only need hours to generate sequence data. In combination with a sample preparation that is fast, like the Ion AmpliSeq™ gene-panel presented here, the TAT from DNA isolation to results is between 48 and 72 hours. While this may not be as fast as some conventional methods, it is adequate for most routine clinical applications.
An additional important consideration is cost. NGS in general is cheap if costs are calculated per sequenced base. However, clinical diagnostic laboratories do not need to sequence hundreds of genes (yet), and neither can they afford to batch numerous samples because of the short TAT needed for clinical specimens. It is therefore important that the sequencing capacity of the NGS platform is flexible and can be adjusted based on the requirements per experiment (e.g. the number of clinical specimens may vary from day to day). A second consideration is that the available sequence capacity on the chip is efficiently used. Since our custom panel specifically investigates genetic aberrations of interest for CRC and NSCLC, the required sequencing capacity for our panel is smaller compared to commercial panels that have been designed for a broad spectrum of cancers. We estimate the bill of materials costs between €130 and €175 per sample (excluding sales tax, depreciation and overhead costs) using the re-designed panel. While this may be higher compared to a single conventional diagnostic test, it is actually a lot cheaper considering that the NGS approach replaces multiple conventional tests (Additional file 6: Table S5 and ).
The rapidly changing landscape in the field of tumor molecular characterization requires the analysis of an increasing amount of targets. With the availability of fast bench top sequencers such as the Ion Torrent PGM, massive parallel sequencing becomes available for clinical laboratories. Here we showed that the development and validation of new tests is feasible with a multi-lab effort, and believe it will facilitate the introduction of these tests into clinical settings. We also demonstrated that implementation of massive parallel sequencing using dedicated gene-panels can be fast, cost-effective and accurate. More importantly, this novel approach might facilitate the adoption of a molecular classification of lung and colon carcinoma in clinical research, thus improving the possibility to identify tumors carrying actionable molecular alterations. Finally, we want to highlight that organization of consortia similar to the one that we established in this study might represent a novel approach to have a rapid, un-biased validation among academic Institutions of new testing methods in a clinical scenario.
We would like to thank Chrysanthi Ainali, Alain Rico, Mark Andersen and Rosella Petraroli from Life Technologies for their expertise and helpful discussion. This work is supported by: AIRC 5xmille (grant n. 12182) and the Italian Cancer Genome Project (FIRB RBAP10AHJB) to AS. AIRC also supports the laboratory of NN (Grant number: IG12118).
- Mendelsohn J. Personalizing oncology: perspectives and prospects. J Clin Oncol. 2013;31(15):1904–11.View ArticlePubMedGoogle Scholar
- Douillard J-Y, Oliner KS, Siena S, Tabernero J, Burkes R, Barugel M, et al. Panitumumab-FOLFOX4 treatment and RAS mutations in colorectal cancer. N Engl J Med. 2013;369(11):1023–34.View ArticlePubMedGoogle Scholar
- Mok TS, Wu YL, Thongprasert S, Yang CH, Chu DT, Saijo N, et al. Gefitinib or carboplatin-paclitaxel in pulmonary adenocarcinoma. N Engl J Med. 2009;361(10):947–57.View ArticlePubMedGoogle Scholar
- Normanno N, Rachiglio AM, Roma C, Fenizia F, Esposito C, Pasquale R, et al. Molecular diagnostics and personalized medicine in oncology: Challenges and opportunities. J Cell Biochem. 2013;114(3):514–24.View ArticlePubMedGoogle Scholar
- Garraway LA. Genomics-Driven Oncology: Framework for an Emerging Paradigm. J Clin Oncol. 2013;31(15):1806–14.View ArticlePubMedGoogle Scholar
- Li TH, Kung HJ, Mack PC, Gandara DR. Genotyping and Genomic Profiling of Non-Small-Cell Lung Cancer: Implications for Current and Future Therapies. J Clin Oncol. 2013;31(8):1039–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Coghlin CL, Smith LJ, Bakar S, Stewart KN, Devereux GS, Nicolson MC, et al. Quantitative analysis of tumour in bronchial biopsy specimens. J Thorac Oncol. 2010;5(4):448–52.View ArticlePubMedGoogle Scholar
- Dienstmann R, Rodon J, Barretina J, Tabernero J. Genomic medicine frontier in human solid tumors: prospects and challenges. J Clin Oncol. 2013;31(15):1874–84.View ArticlePubMedGoogle Scholar
- Tian S, Simon I, Moreno V, Roepman P, Tabernero J, Snel M, et al. A combined oncogenic pathway signature of BRAF, KRAS and PI3KCA mutation improves colorectal cancer classification and cetuximab treatment prediction. Gut. 2013;62(4):540–9.View ArticlePubMedGoogle Scholar
- Yokota T, Ura T, Shibata N, Takahari D, Shitara K, Nomura M, et al. BRAF mutation is a powerful prognostic factor in advanced and recurrent colorectal cancer. Br J Cancer. 2011;104(5):856–62.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Day E, Waters B, Spiegel K, Alnadaf T, Manley PW, Buchdunger E, et al. Inhibition of collagen-induced discoidin domain receptor 1 and 2 activation by imatinib, nilotinib and dasatinib. Eur J Pharmacol. 2008;599(1-3):44–53.View ArticlePubMedGoogle Scholar
- Hammerman PS, Sos ML, Ramos AH, Xu C, Dutt A, Zhou W, et al. Mutations in the DDR2 kinase gene identify a novel therapeutic target in squamous cell lung cancer. Canc Discov. 2011;1(1):78–89.View ArticleGoogle Scholar
- Arcila ME, Chaft JE, Nafa K, Roy-Chowdhuri S, Lau C, Zaidinski M, et al. Prevalence, clinicopathologic associations, and molecular spectrum of ERBB2 (HER2) tyrosine kinase mutations in lung adenocarcinomas. Clin Canc Res. 2012;18(18):4910–8.View ArticleGoogle Scholar
- De Greve J, Teugels E, Geers C, Decoster L, Galdermans D, De Mey J, et al. Clinical activity of afatinib (BIBW 2992) in patients with lung adenocarcinoma with mutations in the kinase domain of HER2/neu. Lung Cancer. 2012;76(1):123–7.View ArticlePubMedGoogle Scholar
- Mazieres J, Peters S, Lepage B, Cortot AB, Barlesi F, Beau-Faller M, et al. Lung cancer that harbors an HER2 mutation: epidemiologic characteristics and therapeutic perspectives. J Clin Oncol. 2013;31(16):1997–2003.View ArticlePubMedGoogle Scholar
- Williams C, Ponten F, Moberg C, Soderkvist P, Uhlen M, Ponten J, et al. A high frequency of sequence alterations is due to formalin fixation of archival specimens. Am J Pathol. 1999;155(5):1467–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Endris V, Penzel R, Warth A, Muckenhuber A, Schirmacher P, Stenzinger A, et al. Molecular diagnostic profiling of lung cancer specimens with a semiconductor-based massive parallel sequencing approach: feasibility, costs, and performance compared with conventional sequencing. The Journal of molecular diagnostics. 2013;15(6):765–75.View ArticlePubMedGoogle Scholar
- Beadling C, Neff TL, Heinrich MC, Rhodes K, Thornton M, Leamon J, et al. Combining highly multiplexed PCR with semiconductor-based sequencing for rapid cancer genotyping. The Journal of molecular diagnostics: JMD. 2013;15(2):171–6.View ArticlePubMedGoogle Scholar
- Hadd AG, Houghton J, Choudhary A, Sah S, Chen L, Marko AC, et al. Targeted, high-depth, next-generation sequencing of cancer genes in formalin-fixed, paraffin-embedded and fine-needle aspiration tumor specimens. The Journal of molecular diagnostics: JMD. 2013;15(2):234–47.View ArticlePubMedGoogle Scholar
- Zhou Q, Zhang XC, Chen ZH, Yin XL, Yang JJ, Xu CR, et al. Relative abundance of EGFR mutations predicts benefit from gefitinib treatment for advanced non-small-cell lung cancer. J Clin Oncol. 2011;29(24):3316–21.View ArticlePubMedGoogle Scholar
- Morabito A, Costanzo R, Rachiglio AM, Pasquale R, Sandomenico C, Franco R, et al. Activity of Gefitinib in a Non-Small-Cell Lung Cancer Patient with Both Activating and Resistance EGFR Mutations. J Thorac Oncol. 2013;8(7):e59–60.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.