The development of noninvasive methods to detect and monitor tumors continues to be a major challenge in oncology. We used digital polymerase chain reaction–based technologies to evaluate the ability of circulating tumor DNA (ctDNA) to detect tumors in 640 patients with various cancer types. We found that ctDNA was detectable in >75% of patients with advanced pancreatic, ovarian, colorectal, bladder, gastroesophageal, breast, melanoma, hepatocellular, and head and neck cancers, but in less than 50% of primary brain, renal, prostate, or thyroid cancers. In patients with localized tumors, ctDNA was detected in 73, 57, 48, and 50% of patients with colorectal cancer, gastroesophageal cancer, pancreatic cancer, and breast adenocarcinoma, respectively. ctDNA was often present in patients without detectable circulating tumor cells, suggesting that these two biomarkers are distinct entities. In a separate panel of 206 patients with metastatic colorectal cancers, we showed that the sensitivity of ctDNA for detection of clinically relevant KRAS gene mutations was 87.2% and its specificity was 99.2%. Finally, we assessed whether ctDNA could provide clues into the mechanisms underlying resistance to epidermal growth factor receptor blockade in 24 patients who objectively responded to therapy but subsequently relapsed. Twenty-three (96%) of these patients developed one or more mutations in genes involved in the mitogen-activated protein kinase pathway. Together, these data suggest that ctDNA is a broadly applicable, sensitive, and specific biomarker that can be used for a variety of clinical and research purposes in patients with multiple different types of cancer.

Colorectal tumors that are wild-type (WT) for KRAS are often sensitive to EGFR blockade, but almost always develop resistance within several months of initiating therapy1,2. The mechanisms underlying this acquired resistance to anti-EGFR antibodies are largely unknown. This situation stands in marked contrast to that of small molecule targeted agents, such as inhibitors of ABL, EGFR, BRAF, and MEK, in which mutations in the genes encoding the protein targets render the tumors resistant to the effects of the drugs3–6. The simplest hypothesis to account for the development of resistance to EGFR blockade are that rare cells with KRAS mutations pre-exist at low levels in tumors with ostensibly WT KRAS genes. Though this hypothesis would seem readily testable, there is no evidence in pre-clinical models to support it, nor is there data from patients. To test this hypothesis, we determined whether mutant KRAS DNA could be detected in the circulation of 28 patients receiving monotherapy with panitumumab, a therapeutic anti-EGFR antibody. We found that nine of 24 (38%) patients whose tumors were initially KRAS WT developed detectable mutations in KRAS in their sera, three of which developed multiple different KRAS mutations. The appearance of these mutations was very consistent, generally occurring between five to six months following treatment. Mathematical modeling indicated that the mutations were present in expanded subclones prior to the initiation of panitumumab. These results suggest that the emergence of KRAS mutations is a mediator of acquired resistance to EGFR blockade and that these mutations can be detected in a non-invasive manner. Moreover, they explain why solid tumors develop resistance to targeted therapies in a highly reproducible fashion.

The identification of mutations that are present in a small fraction of DNA templates is essential for progress in several areas of biomedical research. Although massively parallel sequencing instruments are in principle well suited to this task, the error rates in such instruments are generally too high to allow confident identification of rare variants. We here describe an approach that can substantially increase the sensitivity of massively parallel sequencing instruments for this purpose. The keys to this approach, called the Safe-Sequencing System ("Safe-SeqS"), are (i) assignment of a unique identifier (UID) to each template molecule, (ii) amplification of each uniquely tagged template molecule to create UID families, and (iii) redundant sequencing of the amplification products. PCR fragments with the same UID are considered mutant ("supermutants") only if ≥95% of them contain the identical mutation. We illustrate the utility of this approach for determining the fidelity of a polymerase, the accuracy of oligonucleotides synthesized in vitro, and the prevalence of mutations in the nuclear and mitochondrial genomes of normal cells.diagnostics | early diagnosis | biomarkers | genetics | cancer

More than 2% of the adult U.S. population harbors a pancreatic cyst. These often pose a difficult management problem because conventional criteria cannot always distinguish cysts with malignant potential from those that are innocuous. One of the most common cystic neoplasms of the pancreas, and a bona fide precursor to invasive adenocarcinoma, is called intraductal papillary mucinous neoplasm (IPMN). To help reveal the pathogenesis of these lesions, we purified the DNA from IPMN cyst fluids from 19 patients and searched for mutations in 169 genes commonly altered in human cancers. In addition to the expected KRAS mutations, we identified recurrent mutations at codon 201 of GNAS. A larger number (113) of additional IPMNs were then analyzed to determine the prevalence of KRAS and GNAS mutations. In total, we found that GNAS mutations were present in 66% of IPMNs and that either KRAS or GNAS mutations could be identified in 96%. In eight cases, we could investigate invasive adenocarcinomas that developed in association with IPMNs containing GNAS mutations. In seven of these eight cases, the GNAS mutations present in the IPMNs were also found in the invasive lesion. GNAS mutations were not found in other types of cystic neoplasms of the pancreas or in invasive adenocarcinomas not associated with IPMNs. In addition to defining a new pathway for pancreatic neoplasia, these data suggest that GNAS mutations can inform the diagnosis and management of patients with cystic pancreatic lesions.

More than 2% of adults harbor a pancreatic cyst, a subset of which progresses to invasive lesions with lethal consequences. To assess the genomic landscapes of neoplastic cysts of the pancreas, we determined the exomic sequences of DNA from the neoplastic epithelium of eight surgically resected cysts of each of the major neoplastic cyst types: serous cystadenomas (SCAs), intraductal papillary mucinous neoplasms (IPMNs), mucinous cystic neoplasms (MCNs), and solid pseudopapillary neoplasms (SPNs). SPNs are low-grade malignancies, and IPMNs and MCNs, but not SCAs, have the capacity to progress to cancer. We found that SCAs, IPMNs, MCNs, and SPNs contained 10 ± 4.6, 27 ± 12, 16 ± 7.6, and 2.9 ± 2.1 somatic mutations per tumor, respectively. Among the mutations identified, E3 ubiquitin ligase components were of particular note. Four of the eight SCAs contained mutations of the von Hippel-Lindau gene (VHL), a key component of the VHL ubiquitin ligase complex that has previously been associated with renal cell carcinomas, SCAs, and other neoplasms. Six of the eight IPMNs and three of the eight MCNs harbored mutations of RNF43, a gene coding for a protein with intrinsic E3 ubiquitin ligase activity that has not previously been found to be genetically altered in any human cancer. The preponderance of inactivating mutations in RNF43 unequivocally establish it as a suppressor of both IPMNs and MCNs. SPNs contained remarkably few genetic alterations but always contained mutations of CTNNB1, previously demonstrated to inhibit degradation of the encoded protein (β-catenin) by E3 ubiquitin ligases. These results highlight the essential role of ubiquitin ligases in these neoplasms and have important implications for the diagnosis and treatment of patients with cystic tumors.

The presence of hundreds of copies of mitochondrial (mt) DNA in each human cell poses a challenge for complete characterization of mtDNA genomes by conventional sequencing technologies1. Here, we describe digital sequencing of mtDNA genomes using massively parallel sequencing-by-synthesis. Though the mtDNA of human cells is considered to be homogeneous, we found widespread heterogeneity (heteroplasmy) in the mtDNA of normal human cells. Moreover, the frequency of heteroplasmic variants among different tissues of the same individual varied considerably. In addition to the variants identified in normal tissues, cancer cells harbored additional homoplasmic and heteroplasmic mutations that could also be detected in patient plasma. These studies provide new insights into the nature and variability of mtDNA sequences and have intriguing implications for mitochondrial processes during embryogenesis, cancer biomarker development, and forensic analysis. In particular, they demonstrate that individual humans are characterized by a complex mixture of related mitochondrial genotypes rather than a single genotype.

More information is needed about genetic factors that initiate development of pancreatic intraepithelial neoplasms (PanINs)—the most common precursors of pancreatic ductal adenocarcinoma. We show that >99% of the earliest stage, lowest-grade, PanIN-1 lesions contain mutations in KRAS, p16/CDKN2A, GNAS, or BRAF. These findings could improve our understanding the development and progression of these pre-malignant lesions.

Neuroblastomas are tumors of peripheral sympathetic neurons and are the most common solid tumor in children. To determine the genetic basis for neuroblastoma we performed whole-genome sequencing (6 cases), exome sequencing (16 cases), genome-wide rearrangement analyses (32 cases), and targeted analyses of specific genomic loci (40 cases) using massively parallel sequencing. On average each tumor had 19 somatic alterations in coding genes (range, 3–70). Among genes not previously known to be involved in neuroblastoma, chromosomal deletions and sequence alterations of chromatin remodeling genes, ARID1A and ARID1B, were identified in 8 of 71 tumors (11%) and were associated with early treatment failure and decreased survival. Using tumor-specific structural alterations, we developed an approach to identify rearranged DNA fragments in sera, providing personalized biomarkers for minimal residual disease detection and monitoring. These results highlight dysregulation of chromatin remodeling in pediatric tumorigenesis and provide new approaches for the management of neuroblastoma patients.