Circulating tumour DNA and CT monitoring in patients with untreated diffuse large B-cell lymphoma: a correlative biomarker study.

BACKGROUND: Diffuse large-B-cell lymphoma is curable, but when treatment fails, outcome is poor. Although imaging can help to identify patients at risk of treatment failure, they are often imprecise, and radiation exposure is a potential health risk. We aimed to assess whether circulating tumour DNA encoding the clonal immunoglobulin gene sequence could be detected in the serum of patients with diffuse large-B-cell lymphoma and used to predict clinical disease recurrence after frontline treatment. METHODS: We used next-generation DNA sequencing to retrospectively analyse cell-free circulating tumour DNA in patients assigned to one of three treatment protocols between May 8, 1993, and June 6, 2013. Eligible patients had diffuse large-B-cell lymphoma, no evidence of indolent lymphoma, and were previously untreated. We obtained serial serum samples and concurrent CT scans at specified times during most treatment cycles and up to 5 years of follow-up. VDJ gene segments of the rearranged immunoglobulin receptor genes were amplified and sequenced from pretreatment specimens and serum circulating tumour DNA encoding the VDJ rearrangements was quantitated. FINDINGS: Tumour clonotypes were identified in pretreatment specimens from 126 patients who were followed up for a median of 11 years (IQR 6·8-14·2). Interim monitoring of circulating tumour DNA at the end of two treatment cycles in 108 patients showed a 5-year time to progression of 41·7% (95% CI 22·2-60·1) in patients with detectable circulating tumour DNA and 80·2% (69·6-87·3) in those without detectable circulating tumour DNA (p<0·0001). Detectable interim circulating tumour DNA had a positive predictive value of 62·5% (95% CI 40·6-81·2) and a negative predictive value of 79·8% (69·6-87·8). Surveillance monitoring of circulating tumour DNA was done in 107 patients who achieved complete remission. A Cox proportional hazards model showed that the hazard ratio for clinical disease progression was 228 (95% CI 51-1022) for patients who developed detectable circulating tumour DNA during surveillance compared with patients with undetectable circulating tumour DNA (p<0·0001). Surveillance circulating tumour DNA had a positive predictive value of 88·2% (95% CI 63·6-98·5) and a negative predictive value of 97·8% (92·2-99·7) and identified risk of recurrence at a median of 3·5 months (range 0-200) before evidence of clinical disease. INTERPRETATION: Surveillance circulating tumour DNA identifies patients at risk of recurrence before clinical evidence of disease in most patients and results in a reduced disease burden at relapse. Interim circulating tumour DNA is a promising biomarker to identify patients at high risk of treatment failure. FUNDING: National Cancer Institute and Adaptive Biotechnologies.

Population structure, differential bias and genomic control in a large-scale, case-control association study.

The main problems in drawing causal inferences from epidemiological case-control studies are confounding by unmeasured extraneous factors, selection bias and differential misclassification of exposure. In genetics the first of these, in the form of population structure, has dominated recent debate. Population structure explained part of the significant +11.2% inflation of test statistics we observed in an analysis of 6,322 nonsynonymous SNPs in 816 cases of type 1 diabetes and 877 population-based controls from Great Britain. The remainder of the inflation resulted from differential bias in genotype scoring between case and control DNA samples, which originated from two laboratories, causing false-positive associations. To avoid excluding SNPs and losing valuable information, we extended the genomic control method by applying a variable downweighting to each SNP.

A second generation human haplotype map of over 3.1 million SNPs.

We describe the Phase II HapMap, which characterizes over 3.1 million human single nucleotide polymorphisms (SNPs) genotyped in 270 individuals from four geographically diverse populations and includes 25-35% of common SNP variation in the populations surveyed. The map is estimated to capture untyped common variation with an average maximum r2 of between 0.9 and 0.96 depending on population. We demonstrate that the current generation of commercial genome-wide genotyping products captures common Phase II SNPs with an average maximum r2 of up to 0.8 in African and up to 0.95 in non-African populations, and that potential gains in power in association studies can be obtained through imputation. These data also reveal novel aspects of the structure of linkage disequilibrium. We show that 10-30% of pairs of individuals within a population share at least one region of extended genetic identity arising from recent ancestry and that up to 1% of all common variants are untaggable, primarily because they lie within recombination hotspots. We show that recombination rates vary systematically around genes and between genes of different function. Finally, we demonstrate increased differentiation at non-synonymous, compared to synonymous, SNPs, resulting from systematic differences in the strength or efficacy of natural selection between populations.

Genome-wide detection and characterization of positive selection in human populations.

With the advent of dense maps of human genetic variation, it is now possible to detect positive natural selection across the human genome. Here we report an analysis of over 3 million polymorphisms from the International HapMap Project Phase 2 (HapMap2). We used 'long-range haplotype' methods, which were developed to identify alleles segregating in a population that have undergone recent selection, and we also developed new methods that are based on cross-population comparisons to discover alleles that have swept to near-fixation within a population. The analysis reveals more than 300 strong candidate regions. Focusing on the strongest 22 regions, we develop a heuristic for scrutinizing these regions to identify candidate targets of selection. In a complementary analysis, we identify 26 non-synonymous, coding, single nucleotide polymorphisms showing regional evidence of positive selection. Examination of these candidates highlights three cases in which two genes in a common biological process have apparently undergone positive selection in the same population:LARGE and DMD, both related to infection by the Lassa virus, in West Africa;SLC24A5 and SLC45A2, both involved in skin pigmentation, in Europe; and EDAR and EDA2R, both involved in development of hair follicles, in Asia.

Rapid identification of somatic mutations in colorectal and breast cancer tissues using mismatch repair detection (MRD).

Mismatch repair detection (MRD) was used to screen 93 matched tumor-normal sample pairs and 22 cell lines for somatic mutations in 30 cancer relevant genes. Using a starting amount of only 150 ng of genomic DNA, we screened 102 kb of sequence for somatic mutations in colon and breast cancer. A total of 152 somatic mutations were discovered, encompassing previously reported mutations, such as BRAF V600E and KRAS G12S, G12V, and G13D, as well as novel mutations, including some in genes in which somatic mutations have not previously been reported, such as MAP2K1 and MAP2K2. The distribution of mutations ranged widely within and across tumor types. The functional significance of many of these mutations is not understood, with patterns of selection only evident in KRAS and BRAF in colon cancer. These results present a novel approach to high-throughput mutation screening using small amounts of starting material and reveal a mutation spectrum across 30 genes in a large cohort of breast and colorectal cancers.

Positive selection of a pre-expansion CAG repeat of the human SCA2 gene.

A region of approximately one megabase of human Chromosome 12 shows extensive linkage disequilibrium in Utah residents with ancestry from northern and western Europe. This strikingly large linkage disequilibrium block was analyzed with statistical and experimental methods to determine whether natural selection could be implicated in shaping the current genome structure. Extended Haplotype Homozygosity and Relative Extended Haplotype Homozygosity analyses on this region mapped a core region of the strongest conserved haplotype to the exon 1 of the Spinocerebellar ataxia type 2 gene (SCA2). Direct DNA sequencing of this region of the SCA2 gene revealed a significant association between a pre-expanded allele [(CAG)8CAA(CAG)4CAA(CAG)8] of CAG repeats within exon 1 and the selected haplotype of the SCA2 gene. A significantly negative Tajima's D value (-2.20, p < 0.01) on this site consistently suggested selection on the CAG repeat. This region was also investigated in the three other populations, none of which showed signs of selection. These results suggest that a recent positive selection of the pre-expansion SCA2 CAG repeat has occurred in Utah residents with European ancestry.

Multiplexed genotyping with sequence-tagged molecular inversion probes.

We report on the development of molecular inversion probe (MIP) genotyping, an efficient technology for large-scale single nucleotide polymorphism (SNP) analysis. This technique uses MIPs to produce inverted sequences, which undergo a unimolecular rearrangement and are then amplified by PCR using common primers and analyzed using universal sequence tag DNA microarrays, resulting in highly specific genotyping. With this technology, multiplex analysis of more than 1,000 probes in a single tube can be done using standard laboratory equipment. Genotypes are generated with a high call rate (95%) and high accuracy (>99%) as determined by independent sequencing.

Allele quantification using molecular inversion probes (MIP).

Detection of genomic copy number changes has been an important research area, especially in cancer. Several high-throughput technologies have been developed to detect these changes. Features that are important for the utility of technologies assessing copy number changes include the ability to interrogate regions of interest at the desired density as well as the ability to differentiate the two homologs. In addition, assessing formaldehyde fixed and paraffin embedded (FFPE) samples allows the utilization of the vast majority of cancer samples. To address these points we demonstrate the use of molecular inversion probe (MIP) technology to the study of copy number. MIP is a high-throughput genotyping technology capable of interrogating >20 000 single nucleotide polymorphisms in the same tube. We have shown the ability of MIP at this multiplex level to provide copy number measurements while obtaining the allele information. In addition we have demonstrated a proof of principle for copy number analysis in FFPE samples.

Optimal genotype determination in highly multiplexed SNP data.

High-throughput genotyping technologies that enable large association studies are already available. Tools for genotype determination starting from raw signal intensities need to be automated, robust, and flexible to provide optimal genotype determination given the specific requirements of a study. The key metrics describing the performance of a custom genotyping study are assay conversion, call rate, and genotype accuracy. These three metrics can be traded off against each other. Using the highly multiplexed Molecular Inversion Probe technology as an example, we describe a methodology for identifying the optimal trade-off. The methodology comprises: a robust clustering algorithm and assessment of a large number of data filter sets. The clustering algorithm allows for automatic genotype determination. Many different sets of filters are then applied to the clustered data, and performance metrics resulting from each filter set are calculated. These performance metrics relate to the power of a study and provide a framework to choose the most suitable filter set to the particular study.

Large-scale characterization of public database SNPs causing non-synonymous changes in three ethnic groups.

Single nucleotide polymorphisms (SNPs) that lead to non-synonymous changes in proteins may have functional effects and be subject to selection. Hence they are of particular interest in the study of genetic diseases. We have genotyped approximately 28,000 such SNPs in three ethnic populations (the HapMap plates) and ten primate species and analyzed these data for evidence of selection. We find SNPs predicted by PolyPhen to be damaging, have lower allele frequencies, and are particularly likely to be population-specific. We have also grouped SNPs by molecular function or biological process of the associated genes and find evidence that selection may be acting in concert on classes of genes.

Multiplexed variation scanning for 1,000 amplicons in hundreds of patients using mismatch repair detection (MRD) on tag arrays.

Identification of the genetic basis of common disease may require comprehensive sequence analysis of coding regions and regulatory elements in patients and controls to find genetic effects caused by rare or heterogeneous mutations. In this study, we demonstrate how mismatch repair detection on tag arrays can be applied in a case-control study. Mismatch repair detection allows >1,000 amplicons to be screened for variations in a single laboratory reaction. Variation scanning in 939 amplicons, mostly in coding regions within a linkage peak, was done for 372 patients and 404 controls. In total, >180 Mb of DNA was scanned. Several variants more prevalent in patients than in controls were identified. This study demonstrates an approach to the discovery of susceptibility genes for common disease: large-scale direct sequence comparison between patients and controls. We believe this approach can be scaled up to allow sequence comparison in the whole-genome coding regions among large sets of cases and controls at a reasonable cost in the near future.

Highly multiplexed molecular inversion probe genotyping: over 10,000 targeted SNPs genotyped in a single tube assay.

Large-scale genetic studies are highly dependent on efficient and scalable multiplex SNP assays. In this study, we report the development of Molecular Inversion Probe technology with four-color, single array detection, applied to large-scale genotyping of up to 12,000 SNPs per reaction. While generating 38,429 SNP assays using this technology in a population of 30 trios from the Centre d'Etude Polymorphisme Humain family panel as part of the International HapMap project, we established SNP conversion rates of approximately 90% with concordance rates >99.6% and completeness levels >98% for assays multiplexed up to 12,000plex levels. Furthermore, these individual metrics can be "traded off" and, by sacrificing a small fraction of the conversion rate, the accuracy can be increased to very high levels. No loss of performance is seen when scaling from 6,000plex to 12,000plex assays, strongly validating the ability of the technology to suppress cross-reactivity at high multiplex levels. The results of this study demonstrate the suitability of this technology for comprehensive association studies that use targeted SNPs in indirect linkage disequilibrium studies or that directly screen for causative mutations.

SNP discovery in pooled samples with mismatch repair detection.

A targeted discovery effort is required to identify low frequency single nucleotide polymorphisms (SNPs) in human coding and regulatory regions. We here describe combining mismatch repair detection (MRD) with dideoxy terminator sequencing to detect SNPs in pooled DNA samples. MRD enriches for variant alleles in the pooled sample, and sequencing determines the nature of the variants. By using a genomic DNA pool as a template, approximately 100 fragments were amplified and subsequently combined and subjected en masse to the MRD procedure. The variant-enriched pool from this one MRD reaction is enriched for the population variants of all the tested fragments. Each fragment was amplified from the variant-enriched pool and sequenced, allowing the discovery of alleles with frequencies as low as 1% in the initial population. Our results support that MRD-based SNP discovery can be used for large-scale discovery of SNPs at low frequencies in a population.