Antigen presentation profiling reveals recognition of lymphoma immunoglobulin neoantigens.

Cancer somatic mutations can generate neoantigens that distinguish malignant from normal cells. However, the personalized identification and validation of neoantigens remains a major challenge. Here we discover neoantigens in human mantle-cell lymphomas by using an integrated genomic and proteomic strategy that interrogates tumour antigen peptides presented by major histocompatibility complex (MHC) class I and class II molecules. We applied this approach to systematically characterize MHC ligands from 17 patients. Remarkably, all discovered neoantigenic peptides were exclusively derived from the lymphoma immunoglobulin heavy- or light-chain variable regions. Although we identified MHC presentation of private polymorphic germline alleles, no mutated peptides were recovered from non-immunoglobulin somatically mutated genes. Somatic mutations within the immunoglobulin variable region were almost exclusively presented by MHC class II. We isolated circulating CD4+ T cells specific for immunoglobulin-derived neoantigens and found these cells could mediate killing of autologous lymphoma cells. These results demonstrate that an integrative approach combining MHC isolation, peptide identification, and exome sequencing is an effective platform to uncover tumour neoantigens. Application of this strategy to human lymphoma implicates immunoglobulin neoantigens as targets for lymphoma immunotherapy.

Minimal residual disease negativity using deep sequencing is a major prognostic factor in multiple myeloma.

The introduction of novel agents has led to major improvements in clinical outcomes for patients with multiple myeloma. To shorten evaluation times for new treatments, health agencies are currently examining minimal residual disease (MRD) as a surrogate end point in clinical trials. We assessed the prognostic value of MRD, measured during maintenance therapy by next-generation sequencing (NGS). MRD negativity was defined as the absence of tumor plasma cell within 1 000 000 bone marrow cells (<10-6). Data were analyzed from a recent clinical trial that evaluated the role of transplantation in newly diagnosed myeloma patients treated with lenalidomide, bortezomib, and dexamethasone (RVD). MRD negativity was achieved at least once during maintenance in 127 patients (25%). At the start of maintenance therapy, MRD was a strong prognostic factor for both progression-free survival (adjusted hazard ratio, 0.22; 95% confidence interval, 0.15-0.34; P < .001) and overall survival (adjusted hazard ratio, 0.24; 95% confidence interval, 0.11-0.54; P = .001). Patients who were MRD negative had a higher probability of prolonged progression-free survival than patients with detectable residual disease, regardless of treatment group (RVD vs transplant), cytogenetic risk profile, or International Staging System disease stage at diagnosis. These results were similar after completion of maintenance therapy. Our findings confirm the value of MRD status, as determined by NGS, as a prognostic biomarker in multiple myeloma, and suggest that this approach could be used to adapt treatment strategies in future clinical trials.

Next-generation sequencing-based detection of circulating tumour DNA After allogeneic stem cell transplantation for lymphoma.

Next-generation sequencing (NGS)-based circulating tumour DNA (ctDNA) detection is a promising monitoring tool for lymphoid malignancies. We evaluated whether the presence of ctDNA was associated with outcome after allogeneic haematopoietic stem cell transplantation (HSCT) in lymphoma patients. We studied 88 patients drawn from a phase 3 clinical trial of reduced-intensity conditioning HSCT in lymphoma. Conventional restaging and collection of peripheral blood samples occurred at pre-specified time points before and after HSCT and were assayed for ctDNA by sequencing of the immunoglobulin or T-cell receptor genes. Tumour clonotypes were identified in 87% of patients with adequate tumour samples. Sixteen of 19 (84%) patients with disease progression after HSCT had detectable ctDNA prior to progression at a median of 3·7 months prior to relapse/progression. Patients with detectable ctDNA 3 months after HSCT had inferior progression-free survival (PFS) (2-year PFS 58% vs. 84% in ctDNA-negative patients, P = 0·033). In multivariate models, detectable ctDNA was associated with increased risk of progression/death (Hazard ratio 3·9, P = 0·003) and increased risk of relapse/progression (Hazard ratio 10·8, P = 0·0006). Detectable ctDNA is associated with an increased risk of relapse/progression, but further validation studies are necessary to confirm these findings and determine the clinical utility of NGS-based minimal residual disease monitoring in lymphoma patients after HSCT.

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.

Detection of classical Hodgkin lymphoma specific sequence in peripheral blood using a next-generation sequencing approach.

We applied a highly sensitive next-generation sequencing method to identify lymphoma-specific immunoglobulin gene segments in classical Hodgkin lymphoma (CHL) at initial diagnosis or recurrence, and assessed the ability of detecting such lymphoma-specific sequences in peripheral blood (PB). Seventeen CHL cases were tested and lymphoma-specific sequences were identified in 12 of the primary tumour biopsies. In 11 of these patients whose paired PB samples were available, tumour-specific clonotypes were detected in PB in eight patients. This data demonstrates the feasibility of detecting circulating tumour-specific sequences, creating an unprecedented opportunity to optimize the future treatment and monitoring strategies for patients with CHL.

Immunoglobulin and T cell receptor gene high-throughput sequencing quantifies minimal residual disease in acute lymphoblastic leukemia and predicts post-transplantation relapse and survival.

Minimal residual disease (MRD) quantification is an important predictor of outcome after treatment for acute lymphoblastic leukemia (ALL). Bone marrow ALL burden ≥ 10(-4) after induction predicts subsequent relapse. Likewise, MRD ≥ 10(-4) in bone marrow before initiation of conditioning for allogeneic (allo) hematopoietic cell transplantation (HCT) predicts transplantation failure. Current methods for MRD quantification in ALL are not sufficiently sensitive for use with peripheral blood specimens and have not been broadly implemented in the management of adults with ALL. Consensus-primed immunoglobulin (Ig), T cell receptor (TCR) amplification and high-throughput sequencing (HTS) permit use of a standardized algorithm for all patients and can detect leukemia at 10(-6) or lower. We applied the LymphoSIGHT HTS platform (Sequenta Inc., South San Francisco, CA) to quantification of MRD in 237 samples from 29 adult B cell ALL patients before and after allo-HCT. Using primers for the IGH-VDJ, IGH-DJ, IGK, TCRB, TCRD, and TCRG loci, MRD could be quantified in 93% of patients. Leukemia-associated clonotypes at these loci were identified in 52%, 28%, 10%, 35%, 28%, and 41% of patients, respectively. MRD ≥ 10(-4) before HCT conditioning predicted post-HCT relapse (hazard ratio [HR], 7.7; 95% confidence interval [CI], 2.0 to 30; P = .003). In post-HCT blood samples, MRD ≥10(-6) had 100% positive predictive value for relapse with median lead time of 89 days (HR, 14; 95% CI, 4.7 to 44, P < .0001). The use of HTS-based MRD quantification in adults with ALL offers a standardized approach with sufficient sensitivity to quantify leukemia MRD in peripheral blood. Use of this approach may identify a window for clinical intervention before overt relapse.

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.

Massive evolution of the immunoglobulin heavy chain locus in children with B precursor acute lymphoblastic leukemia.

The ability to distinguish clonal B-cell populations based on the sequence of their rearranged immunoglobulin heavy chain (IgH) locus is an important tool for diagnosing B-cell neoplasms and monitoring treatment response. Leukemic precursor B cells may continue to undergo recombination of the IgH gene after malignant transformation; however, the magnitude of evolution at the IgH locus is currently unknown. We used next-generation sequencing to characterize the repertoire of IgH sequences in diagnostic samples of 51 children with B precursor acute lymphoblastic leukemia (B-ALL). We identified clonal IgH rearrangements in 43 of 51 (84%) cases and found that the number of evolved IgH sequences per patient ranged dramatically from 0 to 4024. We demonstrate that the evolved IgH sequences are not the result of amplification artifacts and are unique to leukemic precursor B cells. In addition, the evolution often follows an allelic exclusion pattern, where only 1 of 2 rearranged IgH loci exhibit ongoing recombination. Thus, precursor B-cell leukemias maintain evolution at the IgH locus at levels that were previously underappreciated. This finding sheds light on the mechanisms associated with leukemic clonal evolution and may fundamentally change approaches for monitoring minimal residual disease burden.

High-throughput method for analyzing methylation of CpGs in targeted genomic regions.

A unique microarray-based method for determining the extent of DNA methylation has been developed. It relies on a selective enrichment of the regions to be assayed by target amplification by capture and ligation (mTACL). The assay is quantitatively accurate, relatively precise, and lends itself to high-throughput determination using nanogram amounts of DNA. The measurements using mTACLs are highly reproducible and in excellent agreement with those obtained by sequencing (r = 0.94). In the present work, the methylation status of >145,000 CpGs from 5,472 promoters in 221 samples was measured. The methylation levels of nearby CpGs are correlated, but the correlation falls off dramatically over several hundred base pairs. In some instances, nearby CpGs have very different levels of methylation. Comparison of normal and tumor samples indicates that in tumors, the promoter regions of genes involved in differentiation and signaling are preferentially hypermethylated, whereas those of housekeeping genes remain hypomethylated. mTACL is a platform for profiling the state of methylation of a large number of CpG in many samples in a cost-effective fashion, and is capable of scaling to much larger numbers of CpGs than those collected here.

Complex inheritance of familial hypercholanemia with associated mutations in TJP2 and BAAT.

Familial hypercholanemia (FHC) is characterized by elevated serum bile acid concentrations, itching, and fat malabsorption. We show here that FHC in Amish individuals is associated with mutations in tight junction protein 2 (encoded by TJP2, also known as ZO-2) and bile acid Coenzyme A: amino acid N-acyltransferase (encoded by BAAT). The mutation of TJP2, which occurs in the first PDZ domain, reduces domain stability and ligand binding in vitro. We noted a morphological change in hepatic tight junctions. The mutation of BAAT, a bile acid-conjugating enzyme, abrogates enzyme activity; serum of individuals homozygous with respect to this mutation contains only unconjugated bile acids. Mutations in both TJP2 and BAAT may disrupt bile acid transport and circulation. Inheritance seems to be oligogenic, with genotype at BAAT modifying penetrance in individuals homozygous with respect to the mutation in TJP2.

A Systematic Method for Detecting Abnormal mRNA Splicing and Assessing Its Clinical Impact in Individuals Undergoing Genetic Testing for Hereditary Cancer Syndromes.

Nearly 14% of disease-causing germline variants result from the disruption of mRNA splicing. Most (67%) DNA variants predicted in silico to disrupt splicing are classified as variants of uncertain significance. An analytic workflow-splice effect event resolver (SPEER)-was developed and validated to use mRNA sequencing to reveal significant deviations in splicing, pinpoint the DNA variants potentially involved, and measure the deleterious effects of the altered splicing on mRNA transcripts, providing evidence for assessing the pathogenicity of the variant. SPEER was used to analyze leukocyte RNA encoding 63 hereditary cancer syndrome-related genes in 20,317 patients. Among 3563 patients (17.5%) with at least one DNA variant predicted to affect splicing, 971 (4.8%) had altered splicing with a deleterious effect on the transcript, and 40 had altered splicing due to a DNA variant located outside of the reportable range of the test. Integrating SPEER results into the interpretation of variants allowed variants of uncertain significance to be reclassified as pathogenic or likely pathogenic in 0.4%, and as benign or likely benign in 5.9%, of the 20,317 patients. SPEER-based evidence was associated with a significantly greater effect on classifications of pathogenic or likely pathogenic and benign or likely benign in nonwhite versus non-Hispanic white patients, illustrating that evidence derived from mRNA splicing analysis may help to reduce ethnic/ancestral disparities in genetic testing.

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 approximately 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 approximately 1 in 500,000 bp and a false negative rate of approximately 10%.

Prediction of epigenetically regulated genes in breast cancer cell lines.

BACKGROUND: Methylation of CpG islands within the DNA promoter regions is one mechanism that leads to aberrant gene expression in cancer. In particular, the abnormal methylation of CpG islands may silence associated genes. Therefore, using high-throughput microarrays to measure CpG island methylation will lead to better understanding of tumor pathobiology and progression, while revealing potentially new biomarkers. We have examined a recently developed high-throughput technology for measuring genome-wide methylation patterns called mTACL. Here, we propose a computational pipeline for integrating gene expression and CpG island methylation profiles to identify epigenetically regulated genes for a panel of 45 breast cancer cell lines, which is widely used in the Integrative Cancer Biology Program (ICBP). The pipeline (i) reduces the dimensionality of the methylation data, (ii) associates the reduced methylation data with gene expression data, and (iii) ranks methylation-expression associations according to their epigenetic regulation. Dimensionality reduction is performed in two steps: (i) methylation sites are grouped across the genome to identify regions of interest, and (ii) methylation profiles are clustered within each region. Associations between the clustered methylation and the gene expression data sets generate candidate matches within a fixed neighborhood around each gene. Finally, the methylation-expression associations are ranked through a logistic regression, and their significance is quantified through permutation analysis. RESULTS: Our two-step dimensionality reduction compressed 90% of the original data, reducing 137,688 methylation sites to 14,505 clusters. Methylation-expression associations produced 18,312 correspondences, which were used to further analyze epigenetic regulation. Logistic regression was used to identify 58 genes from these correspondences that showed a statistically significant negative correlation between methylation profiles and gene expression in the panel of breast cancer cell lines. Subnetwork enrichment of these genes has identified 35 common regulators with 6 or more predicted markers. In addition to identifying epigenetically regulated genes, we show evidence of differentially expressed methylation patterns between the basal and luminal subtypes. CONCLUSIONS: Our results indicate that the proposed computational protocol is a viable platform for identifying epigenetically regulated genes. Our protocol has generated a list of predictors including COL1A2, TOP2A, TFF1, and VAV3, genes whose key roles in epigenetic regulation is documented in the literature. Subnetwork enrichment of these predicted markers further suggests that epigenetic regulation of individual genes occurs in a coordinated fashion and through common regulators.

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.

Nonsynonymous coding single-nucleotide polymorphisms spanning the genome in relation to glioblastoma survival and age at diagnosis.

PURPOSE: Our aim was to discover possible inherited factors associated with glioblastoma age at diagnosis and survival. Although new genotyping technologies allow greatly expanded exploration of such factors, they pose many challenges. EXPERIMENTAL DESIGN: In this pilot study, we (a) genotyped 112 newly diagnosed glioblastoma patients ascertained through a population-based study (group 1) with the ParAllele assay panel of approximately 10,000 nonsynonymous coding single-nucleotide polymorphisms (SNP), (b) used several statistical and bioinformatic techniques to identify 17 SNPs potentially related to either glioblastoma age at diagnosis or survival, and (c) genotyped 16 of these SNPs using conventional PCR methods in an independent group of 195 glioblastoma patients (group 2). RESULTS: In group 2, only one of the 16 SNPs, rs8057643 (located on 16p13.2), was significantly associated with glioblastoma age at diagnosis (nominal P = 0.0017; Bonferroni corrected P = 0.054). Median ages at diagnosis for those with 0, 1, or 2 T alleles were 66, 57, and 59 years in group 1 and 64, 57, and 55 years in group 2 (combined P = 0.001). Furthermore, Cox regression analyses of time to death with number of T alleles adjusted for gender and patient group yielded a hazard ratio of 0.82 (95% confidence interval, 0.68-0.98; P = 0.03). CONCLUSIONS: Although limited by a relatively small sample size, this pilot study, using well-characterized, unambiguous disease characteristics, illustrates the necessity of independent replication owing to the likelihood of false positives. Several other challenges are discussed, including attempts to incorporate information on the potential functional importance of SNPs in genome-disease association studies.

Functional single nucleotide polymorphism-based association studies.

Association studies hold great promise for the elucidation of the genetic basis of diseases. Studies based on functional single nucleotide polymorphisms (SNPs) or on linkage disequilibrium (LD) represent two main types of designs. LD-based association studies can be comprehensive for common causative variants, but they perform poorly for rare alleles. Conversely, functional SNP-based studies are efficient because they focus on the SNPs with the highest a priori chance of being associated. Our poor ability to predict the functional effect of SNPs, however, hampers attempts to make these studies comprehensive. Recent progress in comparative genomics, and evidence that functional elements tend to lie in conserved regions, promises to change the landscape, permitting functional SNP association studies to be carried out that comprehensively assess common and rare alleles. SNP genotyping technologies are already sufficient for such studies, but studies will require continued genomic sequencing of multiple species, research on the functional role of conserved sequences and additional SNP discovery and validation efforts (including targeted SNP discovery to identify the rare alleles in functional regions). With these resources, we expect that comprehensive functional SNP association studies will soon be possible.

The prognostic value of clonal heterogeneity and quantitative assessment of plasma circulating clonal IG-VDJ sequences at diagnosis in patients with follicular lymphoma.

Recent advances in next-generation sequencing (NGS) have enabled the quantitation of circulating tumour DNA (ctDNA) encoding the clonal rearranged V(D)J immunoglobulin locus. We aimed to evaluate the clonal heterogeneity of follicular lymphoma (FL) in the tumour and the plasma at diagnosis and to assess the prognostic value of the ctDNA level. Plasma samples at diagnosis were available for 34 patients registered in the PRIMA trial (NCT00140582). One tumour clonotype or more could be detected for 29 (85%) and 25 (74%) patients, respectively, in the tumour or plasma samples. In 18 patients, several subclones were detected in the tumour (2 to 71 subclones/cases) and/or in the plasma (2 to 20 subclones/cases). In more than half of the cases, the distribution of subclones differed between the tumour and plasma samples, reflecting high clonal heterogeneity and diversity in lymphoma subclone dissemination. In multivariate analysis, a high level of ctDNA was the only independent factor associated with patients' progression-free survival (HR 4, IC 95 (1.1-37), p=.039). In conclusion, an NGS-based immunosequencing method reveals the marked clonal heterogeneity of follicular lymphoma in patients with FL, and quantification of ctDNA at diagnosis represents a potential powerful prognostic biomarker that needs to be investigated in larger cohorts.