de novo variant calling identifies cancer mutation signatures in the 1000 Genomes Project.

Detection of de novo variants (DNVs) is critical for studies of disease-related variation and mutation rates. To accelerate DNV calling, we developed a graphics processing units-based workflow. We applied our workflow to whole-genome sequencing data from three parent-child sequenced cohorts including the Simons Simplex Collection (SSC), Simons Foundation Powering Autism Research (SPARK), and the 1000 Genomes Project (1000G) that were sequenced using DNA from blood, saliva, and lymphoblastoid cell lines (LCLs), respectively. The SSC and SPARK DNV callsets were within expectations for number of DNVs, percent at CpG sites, phasing to the paternal chromosome of origin, and average allele balance. However, the 1000G DNV callset was not within expectations and contained excessive DNVs that are likely cell line artifacts. Mutation signature analysis revealed 30% of 1000G DNV signatures matched B-cell lymphoma. Furthermore, we found variants in DNA repair genes and at Clinvar pathogenic or likely-pathogenic sites and significant excess of protein-coding DNVs in IGLL5; a gene known to be involved in B-cell lymphomas. Our study provides a new rapid DNV caller for the field and elucidates important implications of using sequencing data from LCLs for reference building and disease-related projects.

Standard operating procedure for somatic variant refinement of sequencing data with paired tumor and normal samples.

PURPOSE: Following automated variant calling, manual review of aligned read sequences is required to identify a high-quality list of somatic variants. Despite widespread use in analyzing sequence data, methods to standardize manual review have not been described, resulting in high inter- and intralab variability. METHODS: This manual review standard operating procedure (SOP) consists of methods to annotate variants with four different calls and 19 tags. The calls indicate a reviewer's confidence in each variant and the tags indicate commonly observed sequencing patterns and artifacts that inform the manual review call. Four individuals were asked to classify variants prior to, and after, reading the SOP and accuracy was assessed by comparing reviewer calls with orthogonal validation sequencing. RESULTS: After reading the SOP, average accuracy in somatic variant identification increased by 16.7% (p value = 0.0298) and average interreviewer agreement increased by 12.7% (p value < 0.001). Manual review conducted after reading the SOP did not significantly increase reviewer time. CONCLUSION: This SOP supports and enhances manual somatic variant detection by improving reviewer accuracy while reducing the interreviewer variability for variant calling and annotation.

Recurrent somatic mutations affecting B-cell receptor signaling pathway genes in follicular lymphoma.

Follicular lymphoma (FL) is the most common form of indolent non-Hodgkin lymphoma, yet it remains only partially characterized at the genomic level. To improve our understanding of the genetic underpinnings of this incurable and clinically heterogeneous disease, whole-exome sequencing was performed on tumor/normal pairs from a discovery cohort of 24 patients with FL. Using these data and mutations identified in other B-cell malignancies, 1716 genes were sequenced in 113 FL tumor samples from 105 primarily treatment-naive individuals. We identified 39 genes that were mutated significantly above background mutation rates. CREBBP mutations were associated with inferior PFS. In contrast, mutations in previously unreported HVCN1, a voltage-gated proton channel-encoding gene and B-cell receptor signaling modulator, were associated with improved PFS. In total, 47 (44.8%) patients harbor mutations in the interconnected B-cell receptor (BCR) and CXCR4 signaling pathways. Histone gene mutations were more frequent than previously reported (identified in 43.8% of patients) and often co-occurred (17.1% of patients). A novel, recurrent hotspot was identified at a posttranslationally modified residue in the histone H2B family. This study expands the number of mutated genes described in several known signaling pathways and complexes involved in lymphoma pathogenesis (BCR, Notch, SWitch/sucrose nonfermentable (SWI/SNF), vacuolar ATPases) and identified novel recurrent mutations (EGR1/2, POU2AF1, BTK, ZNF608, HVCN1) that require further investigation in the context of FL biology, prognosis, and treatment.

Mutations associated with progression in follicular lymphoma predict inferior outcomes at diagnosis: Alliance A151303.

Follicular lymphoma (FL) is clinically heterogeneous, with select patients tolerating extended watch-and-wait, whereas others require prompt treatment, suffer progression of disease within 24 months of treatment (POD24), and/or experience aggressive histologic transformation (t-FL). Because our understanding of the relationship between genetic alterations in FL and patient outcomes remains limited, we conducted a clinicogenomic analysis of 370 patients with FL or t-FL (from Cancer and Leukemia Group B/Alliance trials 50402/50701/50803, or real-world cohorts from Washington University School of Medicine, Cleveland Clinic, or University of Miami). FL subsets by grade, stage, watch-and-wait, or POD24 status did not differ by mutation burden, whereas mutation burden was significantly higher in relapsed/refractory (rel/ref) FL and t-FL than in newly diagnosed (dx) FL. Nonetheless, mutation burden in dx FL was not associated with frontline progression-free survival (PFS). CREBBP was the only gene more commonly mutated in FL than in t-FL yet mutated CREBBP was associated with shorter frontline PFS in FL. Mutations in 20 genes were more common in rel/ref FL or t-FL than in dx FL, including 6 significantly mutated genes (SMGs): STAT6, TP53, IGLL5, B2M, SOCS1, and MYD88. We defined a mutations associated with progression (MAP) signature as ≥2 mutations in these 7 genes (6 rel/ref FL or t-FL SMGs plus CREBBP). Patients with dx FL possessing a MAP signature had shorter frontline PFS, revealing a 7-gene set offering insight into FL progression risk potentially more generalizable than the m7-Follicular Lymphoma International Prognostic Index (m7-FLIPI), which had modest prognostic value in our cohort. Future studies are warranted to validate the poor prognosis associated with a MAP signature in dx FL, potentially facilitating novel trials specifically in this high-risk subset of patients.

Ultra-Deep Sequencing Reveals the Mutational Landscape of Classical Hodgkin Lymphoma.

The malignant Hodgkin and Reed Sternberg (HRS) cells of classical Hodgkin lymphoma (cHL) are scarce in affected lymph nodes, creating a challenge to detect driver somatic mutations. As an alternative to cell purification techniques, we hypothesized that ultra-deep exome sequencing would allow genomic study of HRS cells, thereby streamlining analysis and avoiding technical pitfalls. To test this, 31 cHL tumor/normal pairs were exome sequenced to approximately 1,000× median depth of coverage. An orthogonal error-corrected sequencing approach verified >95% of the discovered mutations. We identified mutations in genes novel to cHL including: CDH5 and PCDH7, novel stop gain mutations in IL4R, and a novel pattern of recurrent mutations in pathways regulating Hippo signaling. As a further application of our exome sequencing, we attempted to identify expressed somatic single-nucleotide variants (SNV) in single-nuclei RNA sequencing (snRNA-seq) data generated from a patient in our cohort. Our snRNA analysis identified a clear cluster of cells containing a somatic SNV identified in our deep exome data. This cluster has differentially expressed genes that are consistent with genes known to be dysregulated in HRS cells (e.g., PIM1 and PIM3). The cluster also contains cells with an expanded B-cell clonotype further supporting a malignant phenotype. This study provides proof-of-principle that ultra-deep exome sequencing can be utilized to identify recurrent mutations in HRS cells and demonstrates the feasibility of snRNA-seq in the context of cHL. These studies provide the foundation for the further analysis of genomic variants in large cohorts of patients with cHL. SIGNIFICANCE: Our data demonstrate the utility of ultra-deep exome sequencing in uncovering somatic variants in Hodgkin lymphoma, creating new opportunities to define the genes that are recurrently mutated in this disease. We also show for the first time the successful application of snRNA-seq in Hodgkin lymphoma and describe the expression profile of a putative cluster of HRS cells in a single patient.

Recurrent WNT pathway alterations are frequent in relapsed small cell lung cancer.

Nearly all patients with small cell lung cancer (SCLC) eventually relapse with chemoresistant disease. The molecular mechanisms driving chemoresistance in SCLC remain un-characterized. Here, we describe whole-exome sequencing of paired SCLC tumor samples procured at diagnosis and relapse from 12 patients, and unpaired relapse samples from 18 additional patients. Multiple somatic copy number alterations, including gains in ABCC1 and deletions in MYCL, MSH2, and MSH6, are identifiable in relapsed samples. Relapse samples also exhibit recurrent mutations and loss of heterozygosity in regulators of WNT signaling, including CHD8 and APC. Analysis of RNA-sequencing data shows enrichment for an ASCL1-low expression subtype and WNT activation in relapse samples. Activation of WNT signaling in chemosensitive human SCLC cell lines through APC knockdown induces chemoresistance. Additionally, in vitro-derived chemoresistant cell lines demonstrate increased WNT activity. Overall, our results suggest WNT signaling activation as a mechanism of chemoresistance in relapsed SCLC.

Melorheostosis: Exome sequencing of an associated dermatosis implicates postzygotic mosaicism of mutated KRAS.

Melorheostosis (MEL) is the rare sporadic dysostosis characterized by monostotic or polyostotic osteosclerosis and hyperostosis often distributed in a sclerotomal pattern. The prevailing hypothesis for MEL invokes postzygotic mosaicism. Sometimes scleroderma-like skin changes, considered a representation of the pathogenetic process of MEL, overlie the bony changes, and sometimes MEL becomes malignant. Osteopoikilosis (OPK) is the autosomal dominant skeletal dysplasia that features symmetrically distributed punctate osteosclerosis due to heterozygous loss-of-function mutation within LEMD3. Rarely, radiographic findings of MEL occur in OPK. However, germline mutation of LEMD3 does not explain sporadic MEL. To explore if mosaicism underlies MEL, we studied a boy with polyostotic MEL and characteristic overlying scleroderma-like skin, a few bony lesions consistent with OPK, and a large epidermal nevus known to usually harbor a HRAS, FGFR3, or PIK3CA gene mutation. Exome sequencing was performed to ~100× average read depth for his two dermatoses, two areas of normal skin, and peripheral blood leukocytes. As expected for non-malignant tissues, the patient's mutation burden in his normal skin and leukocytes was low. He, his mother, and his maternal grandfather carried a heterozygous, germline, in-frame, 24-base-pair deletion in LEMD3. Radiographs of the patient and his mother revealed bony foci consistent with OPK, but she showed no MEL. For the patient, somatic variant analysis, using four algorithms to compare all 20 possible pairwise combinations of his five DNA samples, identified only one high-confidence mutation, heterozygous KRAS Q61H (NM_033360.3:c.183A>C, NP_203524.1:p.Gln61His), in both his dermatoses but absent in his normal skin and blood. Thus, sparing our patient biopsy of his MEL bone, we identified a heterozygous somatic KRAS mutation in his scleroderma-like dermatosis considered a surrogate for MEL. This implicates postzygotic mosaicism of mutated KRAS, perhaps facilitated by germline LEMD3 haploinsufficiency, causing his MEL.

A Phase I Trial of BKM120 (Buparlisib) in Combination with Fulvestrant in Postmenopausal Women with Estrogen Receptor-Positive Metastatic Breast Cancer.

PURPOSE: This trial was conducted to determine the maximum tolerated dose (MTD) and preliminary efficacy of buparlisib, an oral pan-class I PI3K inhibitor, plus fulvestrant in postmenopausal women with metastatic estrogen receptor positive (ER(+)) breast cancer. EXPERIMENTAL DESIGN: Phase IA employed a 3+3 design to determine the MTD of buparlisib daily plus fulvestrant. Subsequent cohorts (phase IB and cohort C) evaluated intermittent (5/7-day) and continuous dosing of buparlisib (100 mg daily). No more than 3 prior systemic treatments in the metastatic setting were allowed in these subsequent cohorts. RESULTS: Thirty-one patients were enrolled. MTD was defined as buparlisib 100 mg daily plus fulvestrant. Common adverse events (AE) included fatigue (38.7%), transaminases elevation (35.5%), rash (29%), and diarrhea (19.4%). C-peptide was significantly increased during treatment, consistent with on-target effect of buparlisib. Compared with intermittent dosing, daily buparlisib was associated with more frequent early onset AEs and higher buparlisib plasma concentrations. Among the 29 evaluable patients, the clinical benefit rate was 58.6% (95% CI, 40.7%-74.5%). Response was not associated with PIK3CA mutation or treatment cohort; however, loss of PTEN, progesterone receptor (PgR) expression, or mutation in TP53 was most common in resistant cases, and mutations inAKT1 and ESR1 did not exclude treatment response. CONCLUSIONS: Buparlisib plus fulvestrant is clinically active with manageable AEs in patients with metastatic ER(+)breast cancer. Weekend breaks in buparlisib dosing reduced toxicity. Patients with PgR negative and TP53 mutation did poorly, suggesting buparlisib plus fulvestrant may not be adequately effective against tumors with these poor prognostic molecular features.

Truncating Prolactin Receptor Mutations Promote Tumor Growth in Murine Estrogen Receptor-Alpha Mammary Carcinomas.

Estrogen receptor alpha-positive (ERα+) luminal tumors are the most frequent subtype of breast cancer. Stat1(-/-) mice develop mammary tumors that closely recapitulate the biological characteristics of this cancer subtype. To identify transforming events that contribute to tumorigenesis, we performed whole genome sequencing of Stat1(-/-) primary mammary tumors and matched normal tissues. This investigation identified somatic truncating mutations affecting the prolactin receptor (PRLR) in all tumor and no normal samples. Targeted sequencing confirmed the presence of these mutations in precancerous lesions, indicating that this is an early event in tumorigenesis. Functional evaluation of these heterozygous mutations in Stat1(-/-) mouse embryonic fibroblasts showed that co-expression of truncated and wild-type PRLR led to aberrant STAT3 and STAT5 activation downstream of the receptor, cellular transformation in vitro, and tumor formation in vivo. In conclusion, truncating mutations of PRLR promote tumor growth in a model of human ERα+ breast cancer and warrant further investigation.

Optimizing cancer genome sequencing and analysis.

Tumors are typically sequenced to depths of 75-100× (exome) or 30-50× (whole genome). We demonstrate that current sequencing paradigms are inadequate for tumors that are impure, aneuploid or clonally heterogeneous. To reassess optimal sequencing strategies, we performed ultra-deep (up to ~312×) whole genome sequencing (WGS) and exome capture (up to ~433×) of a primary acute myeloid leukemia, its subsequent relapse, and a matched normal skin sample. We tested multiple alignment and variant calling algorithms and validated ~200,000 putative SNVs by sequencing them to depths of ~1,000×. Additional targeted sequencing provided over 10,000× coverage and ddPCR assays provided up to ~250,000× sampling of selected sites. We evaluated the effects of different library generation approaches, depth of sequencing, and analysis strategies on the ability to effectively characterize a complex tumor. This dataset, representing the most comprehensively sequenced tumor described to date, will serve as an invaluable community resource (dbGaP accession id phs000159).