No abstract

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

BackgroundSignificant clinical and research applications are driving large scale adoption of individualized tumor sequencing in cancer in order to identify tumors-specific mutations. When a matched germline sample is available, somatic mutations may be identified using comparative callers. However, matched germline samples are frequently not available such as with archival tissues, which makes it difficult to distinguish somatic from germline variants. While population databases may be used to filter out known germline variants, recent studies have shown private germline variants result in an inflated false positive rate in unmatched tumor samples, and the number germline false positives in an individual may be related to ancestry.MethodsFirst, we examined the relationship between the germline false positives and ancestry. Then we developed and implemented a tumor only caller (LumosVar) that leverages differences in allelic frequency between somatic and germline variants in impure tumors. We used simulated data to systematically examine how copy number alterations, tumor purity, and sequencing depth should affect the sensitivity of our caller. Finally, we evaluated the caller on real data.ResultsWe find the germline false-positive rate is significantly higher for individuals of non-European Ancestry largely due to the limited diversity in public polymorphism databases and due to population-specific characteristics such as admixture or recent expansions. Our Bayesian tumor only caller (LumosVar) is able to greatly reduce false positives from private germline variants, and our sensitivity is similar to predictions based on simulated data.ConclusionsTaken together, our results suggest that studies of individuals of non-European ancestry would most benefit from our approach. However, high sensitivity requires sufficiently impure tumors and adequate sequencing depth. Even in impure tumors, there are copy number alterations that result in germline and somatic variants having similar allele frequencies, limiting the sensitivity of the approach. We believe our approach could greatly improve the analysis of archival samples in a research setting where the normal is not available.Electronic supplementary materialThe online version of this article (10.1186/s12920-017-0296-8) contains supplementary material, which is available to authorized users.

Current understanding of the mutation spectrum of relapsed/refractory (RR) tumors is limited. We performed whole exome sequencing (WES) on 47 diffuse large B cell lymphoma (DLBCL) tumors that persisted after R-CHOP treatment, 8 matched to primary biopsies. We compared genomic alterations from the RR cohort against two treatment-naïve DLBCL cohorts (n=112). While the overall number and types of mutations did not differ significantly, we identified frequency changes in DLBCL driver genes. The overall frequency of MYD88 mutant samples increased (12% to 19%), but we noted a decrease in p.L265P (8% to 4%) and increase in p.S219C mutations (2% to 6%). CARD11 p.D230N, PIM1 p.K115N and CD79B p.Y196C mutations were not observed in the RR cohort, although these mutations were prominent in the primary DLBCL samples. We observed an increase in BCL2 mutations (21% to 38% of samples), BCL2 amplifications (3% to 6% of samples) and CREBBP mutations (31% to 42% of samples) in the RR cohort, supported by acquisition of mutations in these genes in relapsed compared to diagnostic biopsies from the same patient. These increases may reflect the genetic characteristics of R-CHOP RR tumors expected to be enriched for during clinical trial enrollment. These findings hold significance for a number of emerging targeted therapies aligned to genetic targets and biomarkers in DLBCL, reinforcing the importance of time-of-treatment biomarker screening during DLBCL therapy selection.

BackgroundThe clinical significance of tumor-specific genomic alterations in metastatic renal cell carcinoma (mRCC) is emerging, with several studies suggesting an association between PBRM1 mutations and response with immunotherapy (IO). We sought to determine genomic predictors of differential response to vascular endothelial growth factor–tyrosine kinase inhibitors (VEGF-TKIs) and IO.MethodsConsecutive patients who underwent genomic profiling were identified; patients receiving either VEGF-TKIs or IO were included. Clinical tumor-normal whole exome sequencing and tumor whole transcriptome sequencing test were performed using a Clinical Laboratory Improvement Amendments (CLIA)-certified assay (Ashion Analytics; Phoenix, Arizona, USA). Genomic findings were compared between patients with clinical benefit (CB; complete/partial response or stable disease for >6 months) and no clinical benefit (NCB) in VEGF-TKI-treated patient cohort and IO-treated patient cohort.Results91 patients received genomic profiling and 58 patients received VEGF-TKI and/or IO therapy. 17 received sequenced treatment involving both VEGF-TKI and IO, resulting in 32 patients in the IO cohort and 43 patients in the VEGF-TKI cohort. The most commonly used IO and VEGF-TKIs were nivolumab (66%) and sunitinib (40%). The most frequently detected alterations in the overall cohort were in VHL (64%), PBRM1 (38%), SETD2 (24%), KDM5C (17%) and TERT (12%). TERT promoter mutations were associated with NCB in the IO cohort (p=0.038); transcriptomic analysis revealed multiple differentially regulated pathways downstream of TERT. TERT promoter mutations and PBRM1 mutations were found to be mutually exclusive. While PBRM1 mutations were more prevalent in patients with CB with IO and VEGF-TKIs, no statistically significant association was found.ConclusionsOur analysis found that TERT promoter mutations may be a negative predictor of outcome with IO and are mutually exclusive with PBRM1 loss-of-function mutations.

Archival tumor samples represent a rich resource of annotated specimens for translational genomics research. However, standard variant calling approaches require a matched normal sample from the same individual, which is often not available in the retrospective setting, making it difficult to distinguish between true somatic variants and individual-specific germline variants. Archival sections often contain adjacent normal tissue, but this tissue can include infiltrating tumor cells. As existing comparative somatic variant callers are designed to exclude variants present in the normal sample, a novel approach is required to leverage adjacent normal tissue with infiltrating tumor cells for somatic variant calling. Here we present lumosVar 2.0, a software package designed to jointly analyze multiple samples from the same patient, built upon our previous single sample tumor only variant caller lumosVar 1.0. The approach assumes that the allelic fraction of somatic variants and germline variants follow different patterns as tumor content and copy number state change. lumosVar 2.0 estimates allele specific copy number and tumor sample fractions from the data, and uses a to model to determine expected allelic fractions for somatic and germline variants and to classify variants accordingly. To evaluate the utility of lumosVar 2.0 to jointly call somatic variants with tumor and adjacent normal samples, we used a glioblastoma dataset with matched high and low tumor content and germline whole exome sequencing data (for true somatic variants) available for each patient. Both sensitivity and positive predictive value were improved when analyzing the high tumor and low tumor samples jointly compared to analyzing the samples individually or in-silico pooling of the two samples. Finally, we applied this approach to a set of breast and prostate archival tumor samples for which tumor blocks containing adjacent normal tissue were available for sequencing. Joint analysis using lumosVar 2.0 detected several variants, including known cancer hotspot mutations that were not detected by standard somatic variant calling tools using the adjacent tissue as presumed normal reference. Together, these results demonstrate the utility of leveraging paired tissue samples to improve somatic variant calling when a constitutional sample is not available.

PURPOSE Adult T-cell leukemia/lymphoma (ATLL) is an aggressive disease caused by the human T-cell leukemia virus type 1. Real-world data of ATLL in Latin America are lacking. PATIENTS AND METHODS We analyzed patients with ATLL (acute, lymphomatous, chronic, and smoldering) encountered in 11 Latin American countries between 1995 and 2019. Treatment response was assessed according to the 2009 consensus report. Survival curves were estimated using the Kaplan-Meier method and log-rank test. RESULTS We identified 253 patients; 226 (lymphomatous: n = 122, acute: n = 73, chronic: n = 26, and smoldering: n = 5) had sufficient data for analysis (median age 57 years). Most patients with ATLL were from Peru (63%), Chile (17%), Argentina (8%), and Colombia (7%). Hypercalcemia was positively associated with acute type (57% v lymphomatous 27%, P = .014). The median survival times (months) were 4.3, 7.9, 21.1, and not reached for acute, lymphomatous, chronic, and smoldering forms, with 4-year survival rates of 8%, 22%, 40%, and 80%, respectively. First-line zidovudine (AZT)-interferon alfa (IFN) resulted in an overall response rate of 63% (complete response [CR] 24%) for acute. First-line chemotherapy yielded an overall response rate of 41% (CR 29%) for lymphomatous. CR rate was 42% for etoposide, cyclophosphamide, vincristine, doxorubicin, and prednisone versus 12% for cyclophosphamide, vincristine, doxorubicin, and prednisone–like regimen ( P < .001). Progression-free survival at 1 year for acute type patients treated with AZT-IFN was 67%, whereas 2-year progression-free survival in lymphomatous type patients who achieved CR after chemotherapy was 77%. CONCLUSION This study confirms Latin American ATLL presents at a younger age and has a high incidence of lymphomatous type, low incidence of indolent subtypes, and worse survival rates as compared with Japanese patients. In aggressive ATLL, chemotherapy remains the preferred choice for lymphomatous favoring etoposide-based regimen (etoposide, cyclophosphamide, vincristine, doxorubicin, and prednisone), whereas AZT-IFN remains a good first-line option for acute subtype.

Chordoma is a rare, orphan cancer arising from embryonal precursors of bone. Surgery and radiotherapy (RT) provide excellent local control, often at the price of significant morbidity because of the structures involved and the need for relatively high doses of RT; however, recurrence remains high. Although our understanding of the genetic changes that occur in chordoma is evolving rapidly, this knowledge has yet to translate into treatments. We performed comprehensive DNA (paired tumor/normal whole-exome and shallow whole-genome) and RNA (tumor whole-transcriptome) next-generation sequencing analyses of archival sacral and clivus chordoma specimens. Incorporation of transcriptomic data enabled the identification of gene overexpression and expressed DNA alterations, thus providing additional support for potential therapeutic targets. In three patients, we identified alterations that may be amenable to off-label FDA-approved treatments for other tumor types. These alterations include FGFR1 overexpression (ponatinib, pazopanib) and copy-number duplication of CDK4 (palbociclib) and ERBB3 (gefitinib). In a third patient, germline DNA demonstrated predicted pathogenic changes in CHEK2 and ATM, which may have predisposed the patient to developing chordoma at a young age and may also be associated with potential sensitivity to PARP inhibitors because of homologous recombination repair deficiency. Last, in the fourth patient, a missense mutation in IGF1R was identified, suggesting potential activity for investigational anti-IGF1R strategies. Our findings demonstrate that chordoma patients present with aberrations in overlapping pathways. These results provide support for targeting the IGF1R/FGFR/EGFR and CDK4/6 pathways as treatment strategies for chordoma patients. This study underscores the value of comprehensive genomic and transcriptomic analysis in the development of rational, individualized treatment plans for chordoma.