Structural and non-coding variants increase the diagnostic yield of clinical whole genome sequencing for rare diseases.

BACKGROUND: Whole genome sequencing is increasingly being used for the diagnosis of patients with rare diseases. However, the diagnostic yields of many studies, particularly those conducted in a healthcare setting, are often disappointingly low, at 25-30%. This is in part because although entire genomes are sequenced, analysis is often confined to in silico gene panels or coding regions of the genome. METHODS: We undertook WGS on a cohort of 122 unrelated rare disease patients and their relatives (300 genomes) who had been pre-screened by gene panels or arrays. Patients were recruited from a broad spectrum of clinical specialties. We applied a bioinformatics pipeline that would allow comprehensive analysis of all variant types. We combined established bioinformatics tools for phenotypic and genomic analysis with our novel algorithms (SVRare, ALTSPLICE and GREEN-DB) to detect and annotate structural, splice site and non-coding variants. RESULTS: Our diagnostic yield was 43/122 cases (35%), although 47/122 cases (39%) were considered solved when considering novel candidate genes with supporting functional data into account. Structural, splice site and deep intronic variants contributed to 20/47 (43%) of our solved cases. Five genes that are novel, or were novel at the time of discovery, were identified, whilst a further three genes are putative novel disease genes with evidence of causality. We identified variants of uncertain significance in a further fourteen candidate genes. The phenotypic spectrum associated with RMND1 was expanded to include polymicrogyria. Two patients with secondary findings in FBN1 and KCNQ1 were confirmed to have previously unidentified Marfan and long QT syndromes, respectively, and were referred for further clinical interventions. Clinical diagnoses were changed in six patients and treatment adjustments made for eight individuals, which for five patients was considered life-saving. CONCLUSIONS: Genome sequencing is increasingly being considered as a first-line genetic test in routine clinical settings and can make a substantial contribution to rapidly identifying a causal aetiology for many patients, shortening their diagnostic odyssey. We have demonstrated that structural, splice site and intronic variants make a significant contribution to diagnostic yield and that comprehensive analysis of the entire genome is essential to maximise the value of clinical genome sequencing.

Whole-genome sequencing of chronic lymphocytic leukemia identifies subgroups with distinct biological and clinical features.

The value of genome-wide over targeted driver analyses for predicting clinical outcomes of cancer patients is debated. Here, we report the whole-genome sequencing of 485 chronic lymphocytic leukemia patients enrolled in clinical trials as part of the United Kingdom's 100,000 Genomes Project. We identify an extended catalog of recurrent coding and noncoding genetic mutations that represents a source for future studies and provide the most complete high-resolution map of structural variants, copy number changes and global genome features including telomere length, mutational signatures and genomic complexity. We demonstrate the relationship of these features with clinical outcome and show that integration of 186 distinct recurrent genomic alterations defines five genomic subgroups that associate with response to therapy, refining conventional outcome prediction. While requiring independent validation, our findings highlight the potential of whole-genome sequencing to inform future risk stratification in chronic lymphocytic leukemia.

Clinical utility of whole-genome sequencing in precision oncology.

Precision diagnostics is one of the two pillars of precision medicine. Sequencing efforts in the past decade have firmly established cancer as a primarily genetically driven disease. This concept is supported by therapeutic successes aimed at particular pathways that are perturbed by specific driver mutations in protein-coding domains and reflected in three recent FDA tissue agnostic cancer drug approvals. In addition, there is increasing evidence from studies that interrogate the entire genome by whole-genome sequencing that acquired global and complex genomic aberrations including those in non-coding regions of the genome might also reflect clinical outcome. After addressing technical, logistical, financial and ethical challenges, national initiatives now aim to introduce clinical whole-genome sequencing into real-world diagnostics as a rational and potentially cost-effective tool for response prediction in cancer and to identify patients who would benefit most from 'expensive' targeted therapies and recruitment into clinical trials. However, so far, this has not been accompanied by a systematic and prospective evaluation of the clinical utility of whole-genome sequencing within clinical trials of uniformly treated patients of defined clinical outcome. This approach would also greatly facilitate novel predictive biomarker discovery and validation, ultimately reducing size and duration of clinical trials and cost of drug development. This manuscript is the third in a series of three to review and critically appraise the potential and challenges of clinical whole-genome sequencing in solid tumors and hematological malignancies.

Provocation Testing and Therapeutic Response in a Newly Described Channelopathy: RyR2 Calcium Release Deficiency Syndrome.

BACKGROUND: A novel familial arrhythmia syndrome, cardiac ryanodine receptor (RyR2) calcium release deficiency syndrome (CRDS), has recently been described. We evaluated a large and well characterized family to assess provocation testing, risk factor stratification and response to therapy in CRDS. METHODS: We present a family with multiple unheralded sudden cardiac deaths and aborted cardiac arrests, primarily in children and young adults, with no clear phenotype on standard clinical testing. RESULTS: Genetic analysis, including whole genome sequencing, firmly established that a missense mutation in RYR2, Ala4142Thr, was the underlying cause of disease in the family. Functional study of the variant in a cell model showed RyR2 loss-of-function, indicating that the family was affected by CRDS. EPS (Electrophysiological Study) was undertaken in 9 subjects known to carry the mutation, including a survivor of aborted sudden cardiac death, and the effects of flecainide alone and in combination with metoprolol were tested. There was a clear gradation in inducibility of nonsustained and sustained ventricular arrhythmia between subjects at EPS, with the survivor of aborted sudden cardiac death being the most inducible subject. Administration of flecainide substantially reduced arrhythmia inducibility in this subject and abolished arrhythmia in all others. Finally, the effects of additional metoprolol were tested; it increased inducibility in 4/9 subjects. CONCLUSIONS: The Ala4142Thr mutation of RYR2 causes the novel heritable arrhythmia syndrome CRDS, which is characterized by familial sudden death in the absence of prior symptoms or a recognizable phenotype on ambulatory monitoring or exercise stress testing. We increase the experience of a specific EPS protocol in human subjects and show that it is helpful in establishing the clinical status of gene carriers, with potential utility for risk stratification. Our data provide evidence that flecainide is protective in human subjects with CRDS, consistent with the effect previously shown in a mouse model.

FOXN1 forms higher-order nuclear condensates displaced by mutations causing immunodeficiency.

The transcription factor FOXN1 is a master regulator of thymic epithelial cell (TEC) development and function. Here, we demonstrate that FOXN1 expression is differentially regulated during organogenesis and participates in multimolecular nuclear condensates essential for the factor's transcriptional activity. FOXN1's C-terminal sequence regulates the diffusion velocity within these aggregates and modulates the binding to proximal gene regulatory regions. These dynamics are altered in a patient with a mutant FOXN1 that is modified in its C-terminal sequence. This mutant is transcriptionally inactive and acts as a dominant negative factor displacing wild-type FOXN1 from condensates and causing athymia and severe lymphopenia in heterozygotes. Expression of the mutated mouse ortholog selectively impairs mouse TEC differentiation, revealing a gene dose dependency for individual TEC subtypes. We have therefore identified the cause for a primary immunodeficiency disease and determined the mechanism by which this FOXN1 gain-of-function mutant mediates its dominant negative effect.

Short and long-read genome sequencing methodologies for somatic variant detection; genomic analysis of a patient with diffuse large B-cell lymphoma.

Recent advances in throughput and accuracy mean that the Oxford Nanopore Technologies PromethION platform is a now a viable solution for genome sequencing. Much of the validation of bioinformatic tools for this long-read data has focussed on calling germline variants (including structural variants). Somatic variants are outnumbered many-fold by germline variants and their detection is further complicated by the effects of tumour purity/subclonality. Here, we evaluate the extent to which Nanopore sequencing enables detection and analysis of somatic variation. We do this through sequencing tumour and germline genomes for a patient with diffuse B-cell lymphoma and comparing results with 150 bp short-read sequencing of the same samples. Calling germline single nucleotide variants (SNVs) from specific chromosomes of the long-read data achieved good specificity and sensitivity. However, results of somatic SNV calling highlight the need for the development of specialised joint calling algorithms. We find the comparative genome-wide performance of different tools varies significantly between structural variant types, and suggest long reads are especially advantageous for calling large somatic deletions and duplications. Finally, we highlight the utility of long reads for phasing clinically relevant variants, confirming that a somatic 1.6 Mb deletion and a p.(Arg249Met) mutation involving TP53 are oriented in trans.

Bi-allelic MCM10 variants associated with immune dysfunction and cardiomyopathy cause telomere shortening.

Minichromosome maintenance protein 10 (MCM10) is essential for eukaryotic DNA replication. Here, we describe compound heterozygous MCM10 variants in patients with distinctive, but overlapping, clinical phenotypes: natural killer (NK) cell deficiency (NKD) and restrictive cardiomyopathy (RCM) with hypoplasia of the spleen and thymus. To understand the mechanism of MCM10-associated disease, we modeled these variants in human cell lines. MCM10 deficiency causes chronic replication stress that reduces cell viability due to increased genomic instability and telomere erosion. Our data suggest that loss of MCM10 function constrains telomerase activity by accumulating abnormal replication fork structures enriched with single-stranded DNA. Terminally-arrested replication forks in MCM10-deficient cells require endonucleolytic processing by MUS81, as MCM10:MUS81 double mutants display decreased viability and accelerated telomere shortening. We propose that these bi-allelic variants in MCM10 predispose specific cardiac and immune cell lineages to prematurely arrest during differentiation, causing the clinical phenotypes observed in both NKD and RCM patients.

Validation of clinical-grade whole genome sequencing reproduces cytogenetic analysis and identifies mutational landscape in newly-diagnosed multiple myeloma patients: A pilot study from the 100,000 Genomes Project.

Multiple myeloma is characterized by chromosomal abnormalities and genetic variation, which may inform prognosis and guide treatment. This pilot study sought to examine the feasibility of incorporating Whole Genome Sequencing (WGS) alongside the routine laboratory evaluation of 14 patients with newly diagnosed multiple myeloma who had enrolled in the 100,000 Genomes Project. In all 14 cases, WGS data could be obtained in a timely fashion within existing clinical frameworks in a tertiary hospital setting. The data not only replicated standard-of-care FISH analysis of chromosomal abnormalities but also provided further chromosomal and molecular genetic insights that may influence patient management.

Genomic and transcriptomic correlates of Richter transformation in chronic lymphocytic leukemia.

The transformation of chronic lymphocytic leukemia (CLL) to high-grade B-cell lymphoma is known as Richter syndrome (RS), a rare event with dismal prognosis. In this study, we conducted whole-genome sequencing (WGS) of paired circulating CLL (PB-CLL) and RS biopsies (tissue-RS) from 17 patients recruited into a clinical trial (CHOP-O). We found that tissue-RS was enriched for mutations in poor-risk CLL drivers and genes in the DNA damage response (DDR) pathway. In addition, we identified genomic aberrations not previously implicated in RS, including the protein tyrosine phosphatase receptor (PTPRD) and tumor necrosis factor receptor-associated factor 3 (TRAF3). In the noncoding genome, we discovered activation-induced cytidine deaminase-related and unrelated kataegis in tissue-RS affecting regulatory regions of key immune-regulatory genes. These include BTG2, CXCR4, NFATC1, PAX5, NOTCH-1, SLC44A5, FCRL3, SELL, TNIP2, and TRIM13. Furthermore, differences between the global mutation signatures of pairs of PB-CLL and tissue-RS samples implicate DDR as the dominant mechanism driving transformation. Pathway-based clonal deconvolution analysis showed that genes in the MAPK and DDR pathways demonstrate high clonal-expansion probability. Direct comparison of nodal-CLL and tissue-RS pairs from an independent cohort confirmed differential expression of the same pathways by RNA expression profiling. Our integrated analysis of WGS and RNA expression data significantly extends previous targeted approaches, which were limited by the lack of germline samples, and it facilitates the identification of novel genomic correlates implicated in RS transformation, which could be targeted therapeutically. Our results inform the future selection of investigative agents for a UK clinical platform study. This trial was registered at www.clinicaltrials.gov as #NCT03899337.

The complete costs of genome sequencing: a microcosting study in cancer and rare diseases from a single center in the United Kingdom.

PURPOSE: The translation of genome sequencing into routine health care has been slow, partly because of concerns about affordability. The aspirational cost of sequencing a genome is $1000, but there is little evidence to support this estimate. We estimate the cost of using genome sequencing in routine clinical care in patients with cancer or rare diseases. METHODS: We performed a microcosting study of Illumina-based genome sequencing in a UK National Health Service laboratory processing 399 samples/year. Cost data were collected for all steps in the sequencing pathway, including bioinformatics analysis and reporting of results. Sensitivity analysis identified key cost drivers. RESULTS: Genome sequencing costs £6841 per cancer case (comprising matched tumor and germline samples) and £7050 per rare disease case (three samples). The consumables used during sequencing are the most expensive component of testing (68-72% of the total cost). Equipment costs are higher for rare disease cases, whereas consumable and staff costs are slightly higher for cancer cases. CONCLUSION: The cost of genome sequencing is underestimated if only sequencing costs are considered, and likely surpasses $1000/genome in a single laboratory. This aspirational sequencing cost will likely only be achieved if consumable costs are considerably reduced and sequencing is performed at scale.

Identification of a new VHL exon and complex splicing alterations in familial erythrocytosis or von Hippel-Lindau disease.

Chuvash polycythemia is an autosomal recessive form of erythrocytosis associated with a homozygous p.Arg200Trp mutation in the von Hippel-Lindau (VHL) gene. Since this discovery, additional VHL mutations have been identified in patients with congenital erythrocytosis, in a homozygous or compound-heterozygous state. VHL is a major tumor suppressor gene, mutations in which were first described in patients presenting with VHL disease, which is characterized by the development of highly vascularized tumors. Here, we identify a new VHL cryptic exon (termed E1') deep in intron 1 that is naturally expressed in many tissues. More importantly, we identify mutations in E1' in 7 families with erythrocytosis (1 homozygous case and 6 compound-heterozygous cases with a mutation in E1' in addition to a mutation in VHL coding sequences) and in 1 large family with typical VHL disease but without any alteration in the other VHL exons. In this study, we show that the mutations induced a dysregulation of VHL splicing with excessive retention of E1' and were associated with a downregulation of VHL protein expression. In addition, we demonstrate a pathogenic role for synonymous mutations in VHL exon 2 that altered splicing through E2-skipping in 5 families with erythrocytosis or VHL disease. In all the studied cases, the mutations differentially affected splicing, correlating with phenotype severity. This study demonstrates that cryptic exon retention and exon skipping are new VHL alterations and reveals a novel complex splicing regulation of the VHL gene. These findings open new avenues for diagnosis and research regarding the VHL-related hypoxia-signaling pathway.

The Light Chain IgLV3-21 Defines a New Poor Prognostic Subgroup in Chronic Lymphocytic Leukemia: Results of a Multicenter Study.

Purpose: Unmutated (UM) immunoglobulin heavy chain variable region (IgHV) status or IgHV3-21 gene usage is associated with poor prognosis in chronic lymphocytic leukemia (CLL) patients. Interestingly, IgHV3-21 is often co-expressed with light chain IgLV3-21, which is potentially able to trigger cell-autonomous BCR-mediated signaling. However, this light chain has never been characterized independently of the heavy chain IgHV3-21.Experimental Design: We performed total RNA sequencing in 32 patients and investigated IgLV3-21 prognostic impact in terms of treatment-free survival (TFS) and overall survival (OS) in 3 other independent cohorts for a total of 813 patients. IgLV3-21 presence was tested by real-time PCR and confirmed by Sanger sequencing.Results: Using total RNA sequencing to characterize 32 patients with high-risk CLL, we found a high frequency (28%) of IgLV3-21 rearrangements. Gene set enrichment analysis revealed that these patients express higher levels of genes responsible for ribosome biogenesis and translation initiation (P < 0.0001) as well as MYC target genes (P = 0.0003). Patients with IgLV3-21 rearrangements displayed a significantly shorter TFS and OS (P < 0.05), particularly those with IgHV mutation. In each of the three independent validation cohorts, we showed that IgLV3-21 rearrangements-similar to UM IgHV status-conferred poor prognosis compared with mutated IgHV (P < 0.0001). Importantly, we confirmed by multivariate analysis that this was independent of IgHV mutational status or subset #2 stereotyped receptor (P < 0.0001).Conclusions: We have demonstrated for the first time that a light chain can affect CLL prognosis and that IgLV3-21 light chain usage defines a new subgroup of CLL patients with poor prognosis. Clin Cancer Res; 24(20); 5048-57. ©2018 AACR.

Clinical-grade validation of whole genome sequencing reveals robust detection of low-frequency variants and copy number alterations in CLL.

The 100 000 Genome Project aims to develop a diagnostics platform by introducing whole genome sequencing (WGS) into clinical practice. Samples from patients with chronic lymphocytic leukaemia were subjected to WGS. WGS detection of single nucleotide variants and insertion/deletions were validated by targeted next generation sequencing showing high concordance (96·3%), also for detection of sub-clonal variants and low-frequency TP53 variants. Copy number alteration detection was verified by fluorescent in situ hybridisation and genome-wide single nucleotide polymorphism array (concordances of 86·7% and 92·9%, respectively), confirming adequate sensitivity by WGS. Our results confirm that WGS can provide comprehensive genomic characterisation for clinical trials, drug discovery and, ultimately, precision medicine.

Clinically actionable mutation profiles in patients with cancer identified by whole-genome sequencing.

Next-generation sequencing (NGS) efforts have established catalogs of mutations relevant to cancer development. However, the clinical utility of this information remains largely unexplored. Here, we present the results of the first eight patients recruited into a clinical whole-genome sequencing (WGS) program in the United Kingdom. We performed PCR-free WGS of fresh frozen tumors and germline DNA at 75× and 30×, respectively, using the HiSeq2500 HTv4. Subtracted tumor VCFs and paired germlines were subjected to comprehensive analysis of coding and noncoding regions, integration of germline with somatically acquired variants, and global mutation signatures and pathway analyses. Results were classified into tiers and presented to a multidisciplinary tumor board. WGS results helped to clarify an uncertain histopathological diagnosis in one case, led to informed or supported prognosis in two cases, leading to de-escalation of therapy in one, and indicated potential treatments in all eight. Overall 26 different tier 1 potentially clinically actionable findings were identified using WGS compared with six SNVs/indels using routine targeted NGS. These initial results demonstrate the potential of WGS to inform future diagnosis, prognosis, and treatment choice in cancer and justify the systematic evaluation of the clinical utility of WGS in larger cohorts of patients with cancer.

Clinical whole-genome sequencing from routine formalin-fixed, paraffin-embedded specimens: pilot study for the 100,000 Genomes Project.

PURPOSE: Fresh-frozen (FF) tissue is the optimal source of DNA for whole-genome sequencing (WGS) of cancer patients. However, it is not always available, limiting the widespread application of WGS in clinical practice. We explored the viability of using formalin-fixed, paraffin-embedded (FFPE) tissues, available routinely for cancer patients, as a source of DNA for clinical WGS. METHODS: We conducted a prospective study using DNAs from matched FF, FFPE, and peripheral blood germ-line specimens collected from 52 cancer patients (156 samples) following routine diagnostic protocols. We compared somatic variants detected in FFPE and matching FF samples. RESULTS: We found the single-nucleotide variant agreement reached 71% across the genome and somatic copy-number alterations (CNAs) detection from FFPE samples was suboptimal (0.44 median correlation with FF) due to nonuniform coverage. CNA detection was improved significantly with lower reverse crosslinking temperature in FFPE DNA extraction (80 °C or 65 °C depending on the methods). Our final data showed somatic variant detection from FFPE for clinical decision making is possible. We detected 98% of clinically actionable variants (including 30/31 CNAs). CONCLUSION: We present the first prospective WGS study of cancer patients using FFPE specimens collected in a routine clinical environment proving WGS can be applied in the clinic.

Characterisation of the changing genomic landscape of metastatic melanoma using cell free DNA.

Cancer is characterised by complex somatically acquired genetic aberrations that manifest as intra-tumour and inter-tumour genetic heterogeneity and can lead to treatment resistance. In this case study, we characterise the genome-wide somatic mutation dynamics in a metastatic melanoma patient during therapy using low-input (50 ng) PCR-free whole genome sequencing of cell-free DNA from pre-treatment and post-relapse blood samples. We identify de novo tumour-specific somatic mutations from cell-free DNA, while the sequence context of single nucleotide variants showed the characteristic UV-damage mutation signature of melanoma. To investigate the behaviour of individual somatic mutations during proto-oncogene B-Raf -targeted and immune checkpoint inhibition, amplicon-based deep sequencing was used to verify and track frequencies of 212 single nucleotide variants at 10 distinct time points over 13 months of treatment. Under checkpoint inhibition therapy, we observed an increase in mutant allele frequencies indicating progression on therapy 88 days before clinical determination of non-response positron emission tomogrophy-computed tomography. We also revealed mutations from whole genome sequencing of cell-free DNA that were not present in the tissue biopsy, but that later contributed to relapse. Our findings have potential clinical applications where high quality tumour-tissue derived DNA is not available.

Whole-genome sequencing identifies homozygous BRCA2 deletion guiding treatment in dedifferentiated prostate cancer.

Whole-genome sequencing (WGS) has transformed the understanding of the genetic drivers of cancer and is increasingly being used in cancer medicine to identify personalized therapies. Here we describe a case in which the application of WGS identified a tumoral BRCA2 deletion in a patient with aggressive dedifferentiated prostate cancer that was repeat-biopsied after disease progression. This would not have been detected by standard BRCA testing, and it led to additional treatment with a maintenance poly ADP ribose polymerase (PARP) inhibitor following platinum-based chemotherapy. This case demonstrates that repeat biopsy upon disease progression and application of WGS to tumor samples has meaningful clinical utility and the potential to transform outcomes in patients with cancer.

A novel 33-Gene targeted resequencing panel provides accurate, clinical-grade diagnosis and improves patient management for rare inherited anaemias.

Accurate diagnosis of rare inherited anaemias is challenging, requiring a series of complex and expensive laboratory tests. Targeted next-generation-sequencing (NGS) has been used to investigate these disorders, but the selection of genes on individual panels has been narrow and the validation strategies used have fallen short of the standards required for clinical use. Clinical-grade validation of negative results requires the test to distinguish between lack of adequate sequencing reads at the locations of known mutations and a real absence of mutations. To achieve a clinically-reliable diagnostic test and minimize false-negative results we developed an open-source tool (CoverMi) to accurately determine base-coverage and the 'discoverability' of known mutations for every sample. We validated our 33-gene panel using Sanger sequencing and microarray. Our panel demonstrated 100% specificity and 99·7% sensitivity. We then analysed 57 clinical samples: molecular diagnoses were made in 22/57 (38·6%), corresponding to 32 mutations of which 16 were new. In all cases, accurate molecular diagnosis had a positive impact on clinical management. Using a validated NGS-based platform for routine molecular diagnosis of previously undiagnosed congenital anaemias is feasible in a clinical diagnostic setting, improves precise diagnosis and enhances management and counselling of the patient and their family.