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.

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.

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.

Whole-genome fetal and maternal DNA methylation analysis using MeDIP-NGS for the identification of differentially methylated regions.

DNA methylation is an epigenetic marker that has been shown to vary significantly across different tissues. Taking advantage of the methylation differences between placenta-derived cell-free DNA and maternal blood, several groups employed different approaches for the discovery of fetal-specific biomarkers. The aim of this study was to analyse whole-genome fetal and maternal methylomes in order to identify and confirm the presence of differentially methylated regions (DMRs). We have initially utilized methylated DNA immunoprecipitation (MeDIP) and next-generation sequencing (NGS) to identify genome-wide DMRs between chorionic villus sampling (CVS) and female non-pregnant plasma (PL) and peripheral blood (WBF) samples. Next, using specific criteria, 331 fetal-specific DMRs were selected and confirmed in eight CVS, eight WBF and eight PL samples by combining MeDIP and in-solution targeted enrichment followed by NGS. Results showed higher enrichment in CVS samples as compared to both WBF and PL samples, confirming the distinct methylation levels between fetal and maternal DNA for the selected DMRs. We have successfully implemented a novel approach for the discovery and confirmation of a significant number of fetal-specific DMRs by combining for the first time MeDIP and in-solution targeted enrichment followed by NGS. The implementation of this double-enrichment approach is highly efficient and enables the detailed analysis of multiple DMRs by targeted NGS. Also, this is, to our knowledge, the first reported application of MeDIP on plasma samples, which leverages the implementation of our enrichment methodology in the detection of fetal abnormalities in maternal plasma.

FoSTeS, MMBIR and NAHR at the human proximal Xp region and the mechanisms of human Xq isochromosome formation.

The recently described DNA replication-based mechanisms of fork stalling and template switching (FoSTeS) and microhomology-mediated break-induced replication (MMBIR) were previously shown to catalyze complex exonic, genic and genomic rearrangements. By analyzing a large number of isochromosomes of the long arm of chromosome X (i(Xq)), using whole-genome tiling path array comparative genomic hybridization (aCGH), ultra-high resolution targeted aCGH and sequencing, we provide evidence that the FoSTeS and MMBIR mechanisms can generate large-scale gross chromosomal rearrangements leading to the deletion and duplication of entire chromosome arms, thus suggesting an important role for DNA replication-based mechanisms in both the development of genomic disorders and cancer. Furthermore, we elucidate the mechanisms of dicentric i(Xq) (idic(Xq)) formation and show that most idic(Xq) chromosomes result from non-allelic homologous recombination between palindromic low copy repeats and highly homologous palindromic LINE elements. We also show that non-recurrent-breakpoint idic(Xq) chromosomes have microhomology-associated breakpoint junctions and are likely catalyzed by microhomology-mediated replication-dependent recombination mechanisms such as FoSTeS and MMBIR. Finally, we stress the role of the proximal Xp region as a chromosomal rearrangement hotspot.

KISSPLICE: de-novo calling alternative splicing events from RNA-seq data.

BACKGROUND: In this paper, we address the problem of identifying and quantifying polymorphisms in RNA-seq data when no reference genome is available, without assembling the full transcripts. Based on the fundamental idea that each polymorphism corresponds to a recognisable pattern in a De Bruijn graph constructed from the RNA-seq reads, we propose a general model for all polymorphisms in such graphs. We then introduce an exact algorithm, called KISSPLICE, to extract alternative splicing events. RESULTS: We show that KISSPLICE enables to identify more correct events than general purpose transcriptome assemblers. Additionally, on a 71 M reads dataset from human brain and liver tissues, KISSPLICE identified 3497 alternative splicing events, out of which 56% are not present in the annotations, which confirms recent estimates showing that the complexity of alternative splicing has been largely underestimated so far. CONCLUSIONS: We propose new models and algorithms for the detection of polymorphism in RNA-seq data. This opens the way to a new kind of studies on large HTS RNA-seq datasets, where the focus is not the global reconstruction of full-length transcripts, but local assembly of polymorphic regions. KISSPLICE is available for download at http://alcovna.genouest.org/kissplice/.

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.

A fast and efficient algorithm for mapping short sequences to a reference genome.

Novel high-throughput (Deep) sequencing technology methods have redefined the way genome sequencing is performed. They are able to produce tens of millions of short sequences (reads) in a single experiment and with a much lower cost than previous sequencing methods. In this paper, we present a new algorithm for addressing the problem of efficiently mapping millions of short reads to a reference genome. In particular, we define and solve the Massive Approximate Pattern Matching problem for mapping short sequences to a reference genome.

Clinical significance of TP53, BIRC3, ATM and MAPK-ERK genes in chronic lymphocytic leukaemia: data from the randomised UK LRF CLL4 trial.

Despite advances in chronic lymphocytic leukaemia (CLL) treatment, globally chemotherapy remains a central treatment modality, with chemotherapy trials representing an invaluable resource to explore disease-related/genetic features contributing to long-term outcomes. In 499 LRF CLL4 cases, a trial with >12 years follow-up, we employed targeted resequencing of 22 genes, identifying 623 mutations. After background mutation rate correction, 11/22 genes were recurrently mutated at frequencies between 3.6% (NFKBIE) and 24% (SF3B1). Mutations beyond Sanger resolution (<12% VAF) were observed in all genes, with KRAS mutations principally composed of these low VAF variants. Firstly, employing orthogonal approaches to confirm <12% VAF TP53 mutations, we assessed the clinical impact of TP53 clonal architecture. Whilst ≥ 12% VAF TP53mut cases were associated with reduced PFS and OS, we could not demonstrate a difference between <12% VAF TP53 mutations and either wild type or ≥12% VAF TP53mut cases. Secondly, we identified biallelic BIRC3 lesions (mutation and deletion) as an independent marker of inferior PFS and OS. Finally, we observed that mutated MAPK-ERK genes were independent markers of poor OS in multivariate survival analysis. In conclusion, our study supports using targeted resequencing of expanded gene panels to elucidate the prognostic impact of gene mutations.

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.

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.

Bisulfite Conversion of DNA: Performance Comparison of Different Kits and Methylation Quantitation of Epigenetic Biomarkers that Have the Potential to Be Used in Non-Invasive Prenatal Testing.

INTRODUCTION: Epigenetic alterations, including DNA methylation, play an important role in the regulation of gene expression. Several methods exist for evaluating DNA methylation, but bisulfite sequencing remains the gold standard by which base-pair resolution of CpG methylation is achieved. The challenge of the method is that the desired outcome (conversion of unmethylated cytosines) positively correlates with the undesired side effects (DNA degradation and inappropriate conversion), thus several commercial kits try to adjust a balance between the two. The aim of this study was to compare the performance of four bisulfite conversion kits [Premium Bisulfite kit (Diagenode), EpiTect Bisulfite kit (Qiagen), MethylEdge Bisulfite Conversion System (Promega) and BisulFlash DNA Modification kit (Epigentek)] regarding conversion efficiency, DNA degradation and conversion specificity. METHODS: Performance was tested by combining fully methylated and fully unmethylated λ-DNA controls in a series of spikes by means of Sanger sequencing (0%, 25%, 50% and 100% methylated spikes) and Next-Generation Sequencing (0%, 3%, 5%, 7%, 10%, 25%, 50% and 100% methylated spikes). We also studied the methylation status of two of our previously published differentially methylated regions (DMRs) at base resolution by using spikes of chorionic villus sample in whole blood. RESULTS: The kits studied showed different but comparable results regarding DNA degradation, conversion efficiency and conversion specificity. However, the best performance was observed with the MethylEdge Bisulfite Conversion System (Promega) followed by the Premium Bisulfite kit (Diagenode). The DMRs, EP6 and EP10, were confirmed to be hypermethylated in the CVS and hypomethylated in whole blood. CONCLUSION: Our findings indicate that the MethylEdge Bisulfite Conversion System (Promega) was shown to have the best performance among the kits. In addition, the methylation level of two of our DMRs, EP6 and EP10, was confirmed. Finally, we showed that bisulfite amplicon sequencing is a suitable approach for methylation analysis of targeted regions.

A case of sigmoid endometriosis difficult to differentiate from colon cancer.

BACKGROUND: Although endometriosis with sigmoid serosal involvement is not uncommon in women of childbearing age, the mucosal involvement is rare and differential diagnosis from colon cancer may be difficult due to the lack of pathognomonic symptoms and the poor diagnostic yield of colonoscopy and colonic biopsies. CASE PRESENTATION: We present a case of a young woman with sigmoid endometriosis, in which the initial diagnostic workup suggested colon cancer. Histologic evidence, obtained from a second colonoscopy, along with pelvic ultrasound findings led to the final diagnosis of intestinal endometriosis which was confirmed by laparoscopy. CONCLUSION: Colonic endometriosis is often a diagnostic challenge and should be considered in young women with symptoms from the lower gastrointestinal tract.

Implementation of high resolution whole genome array CGH in the prenatal clinical setting: advantages, challenges, and review of the literature.

Array Comparative Genomic Hybridization analysis is replacing postnatal chromosomal analysis in cases of intellectual disabilities, and it has been postulated that it might also become the first-tier test in prenatal diagnosis. In this study, array CGH was applied in 64 prenatal samples with whole genome oligonucleotide arrays (BlueGnome, Ltd.) on DNA extracted from chorionic villi, amniotic fluid, foetal blood, and skin samples. Results were confirmed with Fluorescence In Situ Hybridization or Real-Time PCR. Fifty-three cases had normal karyotype and abnormal ultrasound findings, and seven samples had balanced rearrangements, five of which also had ultrasound findings. The value of array CGH in the characterization of previously known aberrations in five samples is also presented. Seventeen out of 64 samples carried copy number alterations giving a detection rate of 26.5%. Ten of these represent benign or variables of unknown significance, giving a diagnostic capacity of the method to be 10.9%. If karyotype is performed the additional diagnostic capacity of the method is 5.1% (3/59). This study indicates the ability of array CGH to identify chromosomal abnormalities which cannot be detected during routine prenatal cytogenetic analysis, therefore increasing the overall detection rate. In addition a thorough review of the literature is presented.

Algorithms for mapping short degenerate and weighted sequences to a reference genome.

Novel high-throughput (Deep) sequencing technologies have redefined the way genome sequencing is performed. They are able to produce millions of short sequences in a single experiment and with a much lower cost than previous methods. In this paper, we address the problem of efficiently mapping and classifying millions of short sequences to a reference genome, based on whether they occur exactly once in the genome or not, and by taking into consideration probability scores. In particular, we design algorithms for Massive Exact and Approximate Pattern Matching of short degenerate and weighted sequences, derived from Deep sequencing technologies, to a reference genome.

Hypercoagulable states in patients with hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) patients have an increased risk for venous thromboembolism, mainly portal venous thrombosis (PVT). The aim of this study was to assess the role of acquired and hereditary thrombotic risk factors in HCC patients. Thirty-one patients with HCC, 30 patients with cirrhosis but without HCC or PVT, and 48 matched healthy controls were studied. Mean levels of plasma protein C, protein S, antithrombin, and serum lipoprotein (a) were significantly lower in patients with HCC and in the cirrhotic group compared to the healthy controls. Mean serum homocysteine levels were significantly higher in patients with HCC compared to cirrhotics and healthy controls. The prevalence of activated protein C resistance, factor V Leiden mutation, prothrombin gene mutation G20210GA, and C677T methylenetetrahydrofolate reductase polymorphism was not significantly different among the three groups. In conclusion, thrombophilic defects are common in HCC patients and they might contribute to the observed thrombotic complications in this malignancy.