Whole genome sequencing provides comprehensive genetic testing in childhood B-cell acute lymphoblastic leukaemia.

Childhood B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by recurrent genetic abnormalities that drive risk-directed treatment strategies. Using current techniques, accurate detection of such aberrations can be challenging, due to the rapidly expanding list of key genetic abnormalities. Whole genome sequencing (WGS) has the potential to improve genetic testing, but requires comprehensive validation. We performed WGS on 210 childhood B-ALL samples annotated with clinical and genetic data. We devised a molecular classification system to subtype these patients based on identification of key genetic changes in tumour-normal and tumour-only analyses. This approach detected 294 subtype-defining genetic abnormalities in 96% (202/210) patients. Novel genetic variants, including fusions involving genes in the MAP kinase pathway, were identified. WGS results were concordant with standard-of-care methods and whole transcriptome sequencing (WTS). We expanded the catalogue of genetic profiles that reliably classify PAX5alt and ETV6::RUNX1-like subtypes. Our novel bioinformatic pipeline improved detection of DUX4 rearrangements (DUX4-r): a good-risk B-ALL subtype with high survival rates. Overall, we have validated that WGS provides a standalone, reliable genetic test to detect all subtype-defining genetic abnormalities in B-ALL, accurately classifying patients for the risk-directed treatment stratification, while simultaneously performing as a research tool to identify novel disease biomarkers.

Functional Analysis Identifies Damaging CHEK2 Missense Variants Associated with Increased Cancer Risk.

Heterozygous carriers of germline loss-of-function variants in the tumor suppressor gene checkpoint kinase 2 (CHEK2) are at an increased risk for developing breast and other cancers. While truncating variants in CHEK2 are known to be pathogenic, the interpretation of missense variants of uncertain significance (VUS) is challenging. Consequently, many VUS remain unclassified both functionally and clinically. Here we describe a mouse embryonic stem (mES) cell-based system to quantitatively determine the functional impact of 50 missense VUS in human CHEK2. By assessing the activity of human CHK2 to phosphorylate one of its main targets, Kap1, in Chek2 knockout mES cells, 31 missense VUS in CHEK2 were found to impair protein function to a similar extent as truncating variants, while 9 CHEK2 missense VUS resulted in intermediate functional defects. Mechanistically, most VUS impaired CHK2 kinase function by causing protein instability or by impairing activation through (auto)phosphorylation. Quantitative results showed that the degree of CHK2 kinase dysfunction correlates with an increased risk for breast cancer. Both damaging CHEK2 variants as a group [OR 2.23; 95% confidence interval (CI), 1.62-3.07; P < 0.0001] and intermediate variants (OR 1.63; 95% CI, 1.21-2.20; P = 0.0014) were associated with an increased breast cancer risk, while functional variants did not show this association (OR 1.13; 95% CI, 0.87-1.46; P = 0.378). Finally, a damaging VUS in CHEK2, c.486A>G/p.D162G, was also identified, which cosegregated with familial prostate cancer. Altogether, these functional assays efficiently and reliably identified VUS in CHEK2 that associate with cancer. SIGNIFICANCE: Quantitative assessment of the functional consequences of CHEK2 variants of uncertain significance identifies damaging variants associated with increased cancer risk, which may aid in the clinical management of patients and carriers.

Correction: Rare germline variants in DNA repair genes and the angiogenesis pathway predispose prostate cancer patients to develop metastatic disease.

This article was originally published under the standard License to Publish, but has now been made available under a CC BY 4.0 license. The PDF and HTML versions of the paper have been modified accordingly.

Fine-mapping of prostate cancer susceptibility loci in a large meta-analysis identifies candidate causal variants.

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling.

Rare germline variants in DNA repair genes and the angiogenesis pathway predispose prostate cancer patients to develop metastatic disease.

BACKGROUND: Prostate cancer (PrCa) demonstrates a heterogeneous clinical presentation ranging from largely indolent to lethal. We sought to identify a signature of rare inherited variants that distinguishes between these two extreme phenotypes. METHODS: We sequenced germline whole exomes from 139 aggressive (metastatic, age of diagnosis < 60) and 141 non-aggressive (low clinical grade, age of diagnosis ≥60) PrCa cases. We conducted rare variant association analyses at gene and gene set levels using SKAT and Bayesian risk index techniques. GO term enrichment analysis was performed for genes with the highest differential burden of rare disruptive variants. RESULTS: Protein truncating variants (PTVs) in specific DNA repair genes were significantly overrepresented among patients with the aggressive phenotype, with BRCA2, ATM and NBN the most frequently mutated genes. Differential burden of rare variants was identified between metastatic and non-aggressive cases for several genes implicated in angiogenesis, conferring both deleterious and protective effects. CONCLUSIONS: Inherited PTVs in several DNA repair genes distinguish aggressive from non-aggressive PrCa cases. Furthermore, inherited variants in genes with roles in angiogenesis may be potential predictors for risk of metastases. If validated in a larger dataset, these findings have potential for future clinical application.

Off-Target V(D)J Recombination Drives Lymphomagenesis and Is Escalated by Loss of the Rag2 C Terminus.

Genome-wide analysis of thymic lymphomas from Tp53(-/-) mice with wild-type or C-terminally truncated Rag2 revealed numerous off-target, RAG-mediated DNA rearrangements. A significantly higher fraction of these errors mutated known and suspected oncogenes/tumor suppressor genes than did sporadic rearrangements (p < 0.0001). This tractable mouse model recapitulates recent findings in human pre-B ALL and allows comparison of wild-type and mutant RAG2. Recurrent, RAG-mediated deletions affected Notch1, Pten, Ikzf1, Jak1, Phlda1, Trat1, and Agpat9. Rag2 truncation substantially increased the frequency of off-target V(D)J recombination. The data suggest that interactions between Rag2 and a specific chromatin modification, H3K4me3, support V(D)J recombination fidelity. Oncogenic effects of off-target rearrangements created by this highly regulated recombinase may need to be considered in design of site-specific nucleases engineered for genome modification.

RAG2 mutants alter DSB repair pathway choice in vivo and illuminate the nature of 'alternative NHEJ'.

DNA double-stranded breaks (DSBs) can be repaired by several mechanisms, including classical NHEJ (c-NHEJ) and a poorly defined, error-prone process termed alternative NHEJ (a-NHEJ). How cells choose between these alternatives to join physiologic DSBs remains unknown. Here, we show that deletion of RAG2's C-terminus allows a-NHEJ to repair RAG-mediated DSBs in developing lymphocytes from both c-NHEJ-proficient and c-NHEJ-deficient mice, demonstrating that the V(D)J recombinase influences repair pathway choice in vivo. Analysis of V(D)J junctions revealed that, contrary to expectation, junctional characteristics alone do not reliably distinguish between a-NHEJ and c-NHEJ. These data suggest that a-NHEJ is not necessarily mutagenic, and may be more prevalent than previously appreciated. Whole genome sequencing of a lymphoma arising in a p53(-/-) mouse bearing a C-terminal RAG2 truncation reveals evidence of a-NHEJ and also of aberrant recognition of DNA sequences resembling RAG recognition sites.

A streamlined method for detecting structural variants in cancer genomes by short read paired-end sequencing.

Defining the architecture of a specific cancer genome, including its structural variants, is essential for understanding tumor biology, mechanisms of oncogenesis, and for designing effective personalized therapies. Short read paired-end sequencing is currently the most sensitive method for detecting somatic mutations that arise during tumor development. However, mapping structural variants using this method leads to a large number of false positive calls, mostly due to the repetitive nature of the genome and the difficulty of assigning correct mapping positions to short reads. This study describes a method to efficiently identify large tumor-specific deletions, inversions, duplications and translocations from low coverage data using SVDetect or BreakDancer software and a set of novel filtering procedures designed to reduce false positive calls. Applying our method to a spontaneous T cell lymphoma arising in a core RAG2/p53-deficient mouse, we identified 40 validated tumor-specific structural rearrangements supported by as few as 2 independent read pairs.