Jianbiao Zheng

Diverse somatic mutation patterns and pathway alterations in human cancers.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

Deep-sequencing approach for minimal residual disease detection in acute lymphoblastic leukemia

The persistence of minimal residual disease (MRD) during therapy is the strongest adverse prognostic factor in acute lymphoblastic leukemia (ALL). We developed a high-throughput sequencing method that universally amplifies antigen-receptor gene segments and identifies all clonal gene rearrangements (ie, leukemia-specific sequences) at diagnosis, allowing monitoring of disease progression and clonal evolution during therapy. In the present study, the assay specifically detected 1 leukemic cell among greater than 1 million leukocytes in spike-in experiments. We compared this method with the gold-standard MRD assays multiparameter flow cytometry and allele-specific oligonucleotide polymerase chain reaction (ASO-PCR) using diagnostic and follow-up samples from 106 patients with ALL. Sequencing detected MRD in all 28 samples shown to be positive by flow cytometry and in 35 of the 36 shown to be positive by ASO-PCR and revealed MRD in 10 and 3 additional samples that were negative by flow cytometry and ASO-PCR, respectively. We conclude that this new method allows monitoring of treatment response in ALL and other lymphoid malignancies with great sensitivity and precision. The www.clinicaltrials.gov identifier number for the Total XV study is NCT00137111.

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.

Minimal residual disease quantification using consensus primers and high-throughput IGH sequencing predicts post-transplant relapse in chronic lymphocytic leukemia

Quantification of minimal residual disease (MRD) following allogeneic hematopoietic cell transplantation (allo-HCT) predicts post-transplant relapse in patients with chronic lymphocytic leukemia (CLL). We utilized an MRD-quantification method that amplifies immunoglobulin heavy chain (IGH) loci using consensus V and J segment primers followed by high-throughput sequencing (HTS), enabling quantification with a detection limit of one CLL cell per million mononuclear cells. Using this IGH–HTS approach, we analyzed MRD patterns in over 400 samples from 40 CLL patients who underwent reduced-intensity allo-HCT. Nine patients relapsed within 12 months post-HCT. Of the 31 patients in remission at 12 months post-HCT, disease-free survival was 86% in patients with MRD <10−4 and 20% in those with MRD ⩾10−4 (relapse hazard ratio (HR) 9.0; 95% confidence interval (CI) 2.5–32; P<0.0001), with median follow-up of 36 months. Additionally, MRD predicted relapse at other time points, including 9, 18 and 24 months post-HCT. MRD doubling time <12 months with disease burden ⩾10−5 was associated with relapse within 12 months of MRD assessment in 50% of patients, and within 24 months in 90% of patients. This IGH–HTS method may facilitate routine MRD quantification in clinical trials.

A comprehensive assay for targeted multiplex amplification of human DNA sequences

We developed a robust and reproducible methodology to amplify human sequences in parallel for use in downstream multiplexed sequence analyses. We call the methodology SMART (Spacer Multiplex Amplification Reaction), and it is based, in part, on padlock probe technology. As a proof of principle, we used SMART technology to simultaneously amplify 485 human exons ranging from 100 to 500 bp from human genomic DNA. In multiple repetitions, >90% of the targets were successfully amplified with a high degree of uniformity, with 70% of targets falling within a 10-fold range and all products falling within a 100-fold range of each other in abundance. We used long padlock probes (LPPs) >300 bases in length for the assay, and the increased length of these probes allowed for the capture of human sequences up to 500 bp in length, which is optimal for capturing most human exons. To engineer the LPPs, we developed a method that generates ssDNA molecules with precise ends, using an appropriately designed dsDNA template. The template has appropriate restriction sites engineered into it that can be digested to generate nucleotide overhangs that are suitable for lambda exonuclease digestion, producing a single-stranded probe from dsDNA. The SMART technology is flexible and can be easily adapted to multiplex tens of thousands of target sequences in a single reaction.

High-throughput, high-accuracy array-based resequencing

Although genomewide association studies have successfully identified associations of many common single-nucleotide polymorphisms (SNPs) with common diseases, the SNPs implicated so far account for only a small proportion of the genetic variability of tested diseases. It has been suggested that common diseases may often be caused by rare alleles missed by genomewide association studies. To identify these rare alleles we need high-throughput, high-accuracy resequencing technologies. Although array-based genotyping has allowed genomewide association studies of common SNPs in tens of thousands of samples, array-based resequencing has been limited for 2 main reasons: the lack of a fully multiplexed pipeline for high-throughput sample processing, and failure to achieve sufficient performance. We have recently solved both of these problems and created a fully multiplexed high-throughput pipeline that results in high-quality data. The pipeline consists of target amplification from genomic DNA, followed by allele enrichment to generate pools of purified variant (or nonvariant) DNA and ends with interrogation of purified DNA on resequencing arrays. We have used this pipeline to resequence ≈5 Mb of DNA (on 3 arrays) corresponding to the exons of 1,500 genes in >473 samples; in total >2,350 Mb were sequenced. In the context of this large-scale study we obtained a false positive rate of ≈1 in 500,000 bp and a false negative rate of ≈10%.

Combining Next-Generation Sequencing and Immune Assays: A Novel Method for Identification of Antigen- Specific T Cells

In this study, we combined a novel sequencing method, which can identify individual clonotypes based on their unique T cell receptor (TCR) rearrangement, with existing immune assays to characterize antigen-specific T cell responses. We validated this approach using three types of assays routinely used to measure antigen-specific responses: pentamers which enable identification of T cells bearing specific TCRs, activation marker expression following antigen stimulation and antigen-induced proliferation to identify cytomegalovirus (CMV) specific clonotypes. In one individual, 8 clonotypes were identified using a pentamer reagent derived from the CMV pp65 protein. The same 8 clonotypes were also identified following sequencing of cells that upregulated an activation marker following incubation with an identical peptide derived from pp65. These 8 and an additional 8 clonotypes were identified using a more sensitive CFSE-based proliferation assay. We found clear sequence homology among some of the clonotypes identified, and the CDR3 region in one clonotype was identical to a previously published pp65-specific clonotype sequence. Many of these CMV-specific clonotypes were present at frequencies below 10−5 which are undetectable using standard flow-cytometric methods. These studies suggest that an immune response is comprised of a diverse set of clones, many of which are present at very low frequencies. Thus, the combination of immune assays and sequencing depicts the richness and diversity of an immune response at a level that is not possible using standard immune assays alone. The methods articulated in this work provide an enhanced understanding of T cell-mediated immune responses at the clonal level.

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. Hum Mutat 29(3), 441–450, 2008.

Multiplex Identification of Antigen-Specific T Cell Receptors Using a Combination of Immune Assays and Immune Receptor Sequencing

Monitoring antigen-specific T cells is critical for the study of immune responses and development of biomarkers and immunotherapeutics. We developed a novel multiplex assay that combines conventional immune monitoring techniques and immune receptor repertoire sequencing to enable identification of T cells specific to large numbers of antigens simultaneously. We multiplexed 30 different antigens and identified 427 antigen-specific clonotypes from 5 individuals with frequencies as low as 1 per million T cells. The clonotypes identified were validated several ways including repeatability, concordance with published clonotypes, and high correlation with ELISPOT. Applying this technology we have shown that the vast majority of shared antigen-specific clonotypes identified in different individuals display the same specificity. We also showed that shared antigen-specific clonotypes are simpler sequences and are present at higher frequencies compared to non-shared clonotypes specific to the same antigen. In conclusion this technology enables sensitive and quantitative monitoring of T cells specific for hundreds or thousands of antigens simultaneously allowing the study of T cell responses with an unprecedented resolution and scale.

Discovery of T Cell Receptor β Motifs Specific to HLA–B27–Positive Ankylosing Spondylitis by Deep Repertoire Sequence Analysis

Ankylosing spondylitis (AS), a chronic inflammatory disorder, has a notable association with HLA–B27. One hypothesis suggests that a common antigen that binds to HLA–B27 is important for AS disease pathogenesis. This study was undertaken to determine sequences and motifs that are shared among HLA–B27–positive AS patients, using T cell repertoire next‐generation sequencing.