Optimized whole-genome sequencing workflow for tumor diagnostics in routine pathology practice.

Two decades after the genomics revolution, oncology is rapidly transforming into a genome-driven discipline, yet routine cancer diagnostics is still mainly microscopy based, except for tumor type-specific predictive molecular tests. Pathology laboratories struggle to quickly validate and adopt biomarkers identified by genomics studies of new targeted therapies. Consequently, clinical implementation of newly approved biomarkers suffers substantial delays, leading to unequal patient access to these therapies. Whole-genome sequencing (WGS) can successfully address these challenges by providing a stable molecular diagnostic platform that allows detection of a multitude of genomic alterations in a single cost-efficient assay and facilitating rapid implementation, as well as by the development of new genomic biomarkers. Recently, the Whole-genome sequencing Implementation in standard Diagnostics for Every cancer patient (WIDE) study demonstrated that WGS is a feasible and clinically valid technique in routine clinical practice with a turnaround time of 11 workdays. As a result, WGS was successfully implemented at the Netherlands Cancer Institute as part of routine diagnostics in January 2021. The success of implementing WGS has relied on adhering to a comprehensive protocol including recording patient information, sample collection, shipment and storage logistics, sequencing data interpretation and reporting, integration into clinical decision-making and data usage. This protocol describes the use of fresh-frozen samples that are necessary for WGS but can be challenging to implement in pathology laboratories accustomed to using formalin-fixed paraffin-embedded samples. In addition, the protocol outlines key considerations to guide uptake of WGS in routine clinical care in hospitals worldwide.

Feasibility of whole-genome sequencing-based tumor diagnostics in routine pathology practice.

The current increase in number and diversity of targeted anticancer agents poses challenges to the logistics and timeliness of molecular diagnostics (MolDx), resulting in underdiagnosis and treatment. Whole-genome sequencing (WGS) may provide a sustainable solution for addressing current as well as future diagnostic challenges. The present study therefore aimed to prospectively assess feasibility, validity, and value of WGS in routine clinical practice. WGS was conducted independently of, and in parallel with, standard of care (SOC) diagnostics on routinely obtained tumor samples from 1,200 consecutive patients with metastatic cancer. Results from both tests were compared and discussed in a dedicated tumor board. From 1,200 patients, 1,302 samples were obtained, of which 1,216 contained tumor cells. WGS was successful in 70% (854/1,216) of samples with a median turnaround time of 11 days. Low tumor purity (<20%) was the main reason for not completing WGS. WGS identified 99.2% and SOC MolDx 99.7% of the total of 896 biomarkers found in genomic regions covered by both tests. Actionable biomarkers were found in 603/848 patients (71%). Of the 936 associated therapy options identified by WGS, 343 were identified with SOC MolDx (36.6%). Biomarker-based therapy was started in 147 patients. WGS revealed 49 not previously identified pathogenic germline variants. Fresh-frozen, instead of formalin-fixed and paraffin-embedded, sample logistics were easily adopted as experienced by the professionals involved. WGS for patients with metastatic cancer is well feasible in routine clinical practice, successfully yielding comprehensive genomic profiling for the vast majority of patients. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Analytical demands to use whole-genome sequencing in precision oncology.

Interrogating the tumor genome in its entirety by whole-genome sequencing (WGS) offers an unprecedented insight into the biology and pathogenesis of cancer, with potential impact on diagnostics, prognostication and therapy selection. WGS is able to detect sequence as well as structural variants and thereby combines central domains of cytogenetics and molecular genetics. Given the potential of WGS in directing targeted therapeutics and clinical decision-making, we envision a gradual transition of the method from research to clinical routine. This review is one out of three within this issue aimed at facilitating this effort, by discussing in-depth analytical validation, clinical interpretation and clinical utility of WGS. The review highlights the requirements for implementing, validating and maintaining a clinical WGS pipeline to obtain high-quality patient-specific data in accordance with the local regulatory landscape. Every step of the WGS pipeline, which includes DNA extraction, library preparation, sequencing, bioinformatics analysis, and data storage, is considered with respect to its logistics, necessities, potential pitfalls, and the required quality management. WGS is likely to drive clinical diagnostics and patient care forward, if requirements and challenges of the technique are recognized and met.

A multi-platform reference for somatic structural variation detection.

Accurate detection of somatic structural variation (SV) in cancer genomes remains a challenging problem. This is in part due to the lack of high-quality, gold-standard datasets that enable the benchmarking of experimental approaches and bioinformatic analysis pipelines. Here, we performed somatic SV analysis of the paired melanoma and normal lymphoblastoid COLO829 cell lines using four different sequencing technologies. Based on the evidence from multiple technologies combined with extensive experimental validation, we compiled a comprehensive set of carefully curated and validated somatic SVs, comprising all SV types. We demonstrate the utility of this resource by determining the SV detection performance as a function of tumor purity and sequence depth, highlighting the importance of assessing these parameters in cancer genomics projects. The truth somatic SV dataset as well as the underlying raw multi-platform sequencing data are freely available and are an important resource for community somatic benchmarking efforts.

Clinical Validation of Whole Genome Sequencing for Cancer Diagnostics.

Whole genome sequencing (WGS) using fresh-frozen tissue and matched blood samples from cancer patients may become the most complete genetic tumor test. With the increasing availability of small biopsies and the need to screen more number of biomarkers, the use of a single all-inclusive test is preferable over multiple consecutive assays. To meet high-quality diagnostics standards, we optimized and clinically validated WGS sample and data processing procedures, resulting in a technical success rate of 95.6% for fresh-frozen samples with sufficient (≥20%) tumor content. Independent validation of identified biomarkers against commonly used diagnostic assays showed a high sensitivity (recall; 98.5%) and precision (positive predictive value; 97.8%) for detection of somatic single-nucleotide variants and insertions and deletions (across 22 genes), and high concordance for detection of gene amplification (97.0%; EGFR and MET) as well as somatic complete loss (100%; CDKN2A/p16). Gene fusion analysis showed a concordance of 91.3% between DNA-based WGS and an orthogonal RNA-based gene fusion assay. Microsatellite (in)stability assessment showed a sensitivity of 100% with a precision of 94%, and virus detection (human papillomavirus), an accuracy of 100% compared with standard testing. In conclusion, whole genome sequencing has a >95% sensitivity and precision compared with routinely used DNA techniques in diagnostics, and all relevant mutation types can be detected reliably in a single assay.

Micro-costing diagnostics in oncology: from single-gene testing to whole- genome sequencing.

Purpose: Predictive diagnostics play an increasingly important role in personalized medicine for cancer treatment. Whole-genome sequencing (WGS)-based treatment selection is expected to rapidly increase worldwide. This study aimed to calculate and compare the total cost of currently used diagnostic techniques and of WGS in treatment of non-small cell lung carcinoma (NSCLC), melanoma, colorectal cancer (CRC), and gastrointestinal stromal tumor (GIST) in the Netherlands.Methods: The activity-based costing (ABC) method was conducted to calculate total cost of included diagnostic techniques based on data provided by Dutch pathology laboratories and the Dutch-centralized cancer WGS facility. Costs were allocated to four categories: capital costs, maintenance costs, software costs, and operational costs.Results: The total cost per cancer patient per technique varied from € 58 (Sanger sequencing, three amplicons) to € 2925 (paired tumor-normal WGS). The operational costs accounted for the vast majority (over 90%) of the total per cancer patient technique costs.Conclusion: This study outlined in detail all costing aspects and cost prices of current and new diagnostic modalities used in treatment of NSCLC, melanoma, CRC, and GIST in the Netherlands. Detailed cost differences and value comparisons between these diagnostic techniques enable future economic evaluations to support decision-making.

Study protocol: Whole genome sequencing Implementation in standard Diagnostics for Every cancer patient (WIDE).

BACKGROUND: 'Precision oncology' can ensure the best suitable treatment at the right time by tailoring treatment towards individual patient and comprehensive tumour characteristics. In current molecular pathology, diagnostic tests which are part of the standard of care (SOC) only cover a limited part of the spectrum of genomic changes, and often are performed in an iterative way. This occurs at the expense of valuable patient time, available tissue sample, and interferes with 'first time right' treatment decisions. Whole Genome Sequencing (WGS) captures a near complete view of genomic characteristics of a tumour in a single test. Moreover, WGS facilitates faster implementation of new treatment relevant biomarkers. At present, WGS mainly has been applied in study settings, but its performance in a routine diagnostic setting remains to be evaluated. The WIDE study aims to investigate the feasibility and validity of WGS-based diagnostics in clinical practice. METHODS: 1200 consecutive patients in a single comprehensive cancer centre with (suspicion of) a metastasized solid tumour will be enrolled with the intention to analyse tumour tissue with WGS, in parallel to SOC diagnostics. Primary endpoints are (1) feasibility of implementation of WGS-based diagnostics into routine clinical care and (2) clinical validation of WGS by comparing identification of treatment-relevant variants between WGS and SOC molecular diagnostics. Secondary endpoints entail (1) added clinical value in terms of additional treatment options and (2) cost-effectiveness of WGS compared to SOC diagnostics through a Health Technology Assessment (HTA) analysis. Furthermore, the (3) perceived impact of WGS-based diagnostics on clinical decision making will be evaluated through questionnaires. The number of patients included in (experimental) therapies initiated based on SOC or WGS diagnostics will be reported with at least 3 months follow-up. The clinical efficacy is beyond the scope of WIDE. Key performance indicators will be evaluated after every 200 patients enrolled, and procedures optimized accordingly, to continuously improve the diagnostic performance of WGS in a routine clinical setting. DISCUSSION: WIDE will yield the optimal conditions under which WGS can be implemented in a routine molecular diagnostics setting and establish the position of WGS compared to SOC diagnostics in routine clinical care.

Pan-cancer whole-genome analyses of metastatic solid tumours.

Metastatic cancer is a major cause of death and is associated with poor treatment efficacy. A better understanding of the characteristics of late-stage cancer is required to help adapt personalized treatments, reduce overtreatment and improve outcomes. Here we describe the largest, to our knowledge, pan-cancer study of metastatic solid tumour genomes, including whole-genome sequencing data for 2,520 pairs of tumour and normal tissue, analysed at median depths of 106× and 38×, respectively, and surveying more than 70 million somatic variants. The characteristic mutations of metastatic lesions varied widely, with mutations that reflect those of the primary tumour types, and with high rates of whole-genome duplication events (56%). Individual metastatic lesions were relatively homogeneous, with the vast majority (96%) of driver mutations being clonal and up to 80% of tumour-suppressor genes being inactivated bi-allelically by different mutational mechanisms. Although metastatic tumour genomes showed similar mutational landscape and driver genes to primary tumours, we find characteristics that could contribute to responsiveness to therapy or resistance in individual patients. We implement an approach for the review of clinically relevant associations and their potential for actionability. For 62% of patients, we identify genetic variants that may be used to stratify patients towards therapies that either have been approved or are in clinical trials. This demonstrates the importance of comprehensive genomic tumour profiling for precision medicine in cancer.

Mapping and phasing of structural variation in patient genomes using nanopore sequencing.

Despite improvements in genomics technology, the detection of structural variants (SVs) from short-read sequencing still poses challenges, particularly for complex variation. Here we analyse the genomes of two patients with congenital abnormalities using the MinION nanopore sequencer and a novel computational pipeline-NanoSV. We demonstrate that nanopore long reads are superior to short reads with regard to detection of de novo chromothripsis rearrangements. The long reads also enable efficient phasing of genetic variations, which we leveraged to determine the parental origin of all de novo chromothripsis breakpoints and to resolve the structure of these complex rearrangements. Additionally, genome-wide surveillance of inherited SVs reveals novel variants, missed in short-read data sets, a large proportion of which are retrotransposon insertions. We provide a first exploration of patient genome sequencing with a nanopore sequencer and demonstrate the value of long-read sequencing in mapping and phasing of SVs for both clinical and research applications.

Sensitive Monogenic Noninvasive Prenatal Diagnosis by Targeted Haplotyping.

During pregnancy, cell-free DNA (cfDNA) in maternal blood encompasses a small percentage of cell-free fetal DNA (cffDNA), an easily accessible source for determination of fetal disease status in risk families through non-invasive procedures. In case of monogenic heritable disease, background maternal cfDNA prohibits direct observation of the maternally inherited allele. Non-invasive prenatal diagnostics (NIPD) of monogenic diseases therefore relies on parental haplotyping and statistical assessment of inherited alleles from cffDNA, techniques currently unavailable for routine clinical practice. Here, we present monogenic NIPD (MG-NIPD), which requires a blood sample from both parents, for targeted locus amplification (TLA)-based phasing of heterozygous variants selectively at a gene of interest. Capture probes-based targeted sequencing of cfDNA from the pregnant mother and a tailored statistical analysis enables predicting fetal gene inheritance. MG-NIPD was validated for 18 pregnancies, focusing on CFTR, CYP21A2, and HBB. In all cases we could predict the inherited alleles with >98% confidence, even at relatively early stages (8 weeks) of pregnancy. This prediction and the accuracy of parental haplotyping was confirmed by sequencing of fetal material obtained by parallel invasive procedures. MG-NIPD is a robust method that requires standard instrumentation and can be implemented in any clinic to provide families carrying a severe monogenic disease with a prenatal diagnostic test based on a simple blood draw.

The PLETHORA Gene Regulatory Network Guides Growth and Cell Differentiation in Arabidopsis Roots.

Organ formation in animals and plants relies on precise control of cell state transitions to turn stem cell daughters into fully differentiated cells. In plants, cells cannot rearrange due to shared cell walls. Thus, differentiation progression and the accompanying cell expansion must be tightly coordinated across tissues. PLETHORA (PLT) transcription factor gradients are unique in their ability to guide the progression of cell differentiation at different positions in the growing Arabidopsis thaliana root, which contrasts with well-described transcription factor gradients in animals specifying distinct cell fates within an essentially static context. To understand the output of the PLT gradient, we studied the gene set transcriptionally controlled by PLTs. Our work reveals how the PLT gradient can regulate cell state by region-specific induction of cell proliferation genes and repression of differentiation. Moreover, PLT targets include major patterning genes and autoregulatory feedback components, enforcing their role as master regulators of organ development.

Epigenomic annotation of gene regulatory alterations during evolution of the primate brain.

Although genome sequencing has identified numerous noncoding alterations between primate species, which of those are regulatory and potentially relevant to the evolution of the human brain is unclear. Here we annotated cis-regulatory elements (CREs) in the human, rhesus macaque and chimpanzee genomes using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in different anatomical regions of the adult brain. We found high similarity in the genomic positioning of rhesus macaque and human CREs, suggesting that the majority of these elements were already present in a common ancestor 25 million years ago. Most of the observed regulatory changes between humans and rhesus macaques occurred before the ancestral separation of humans and chimpanzees, leaving a modest set of regulatory elements with predicted human specificity. Our data refine previous predictions and hypotheses on the consequences of genomic changes between primate species and allow the identification of regulatory alterations relevant to the evolution of the brain.

Genomic and functional overlap between somatic and germline chromosomal rearrangements.

Genomic rearrangements are a common cause of human congenital abnormalities. However, their origin and consequences are poorly understood. We performed molecular analysis of two patients with congenital disease who carried de novo genomic rearrangements. We found that the rearrangements in both patients hit genes that are recurrently rearranged in cancer (ETV1, FOXP1, and microRNA cluster C19MC) and drive formation of fusion genes similar to those described in cancer. Subsequent analysis of a large set of 552 de novo germline genomic rearrangements underlying congenital disorders revealed enrichment for genes rearranged in cancer and overlap with somatic cancer breakpoints. Breakpoints of common (inherited) germline structural variations also overlap with cancer breakpoints but are depleted for cancer genes. We propose that the same genomic positions are prone to genomic rearrangements in germline and soma but that timing and context of breakage determines whether developmental defects or cancer are promoted.

Large-scale identification of coregulated enhancer networks in the adult human brain.

Understanding the complexity of the human brain and its functional diversity remain a major challenge. Distinct anatomical regions are involved in an array of processes, including organismal homeostasis, cognitive functions, and susceptibility to neurological pathologies, many of which define our species. Distal enhancers have emerged as key regulatory elements that acquire histone modifications in a cell- and species-specific manner, thus enforcing specific gene expression programs. Here, we survey the epigenomic landscape of promoters and cis-regulatory elements in 136 regions of the adult human brain. We identify a total of 83,553 promoter-distal H3K27ac-enriched regions showing global characteristics of brain enhancers. We use coregulation of enhancer elements across many distinct regions of the brain to uncover functionally distinct networks at high resolution and link these networks to specific neuroglial functions. Furthermore, we use these data to understand the relevance of noncoding genomic variations previously linked to Parkinson's disease incidence.

Extensive localization of long noncoding RNAs to the cytosol and mono- and polyribosomal complexes.

BACKGROUND: Long noncoding RNAs (lncRNAs) form an abundant class of transcripts, but the function of the majority of them remains elusive. While it has been shown that some lncRNAs are bound by ribosomes, it has also been convincingly demonstrated that these transcripts do not code for proteins. To obtain a comprehensive understanding of the extent to which lncRNAs bind ribosomes, we performed systematic RNA sequencing on ribosome-associated RNA pools obtained through ribosomal fractionation and compared the RNA content with nuclear and (non-ribosome bound) cytosolic RNA pools. RESULTS: The RNA composition of the subcellular fractions differs significantly from each other, but lncRNAs are found in all locations. A subset of specific lncRNAs is enriched in the nucleus but surprisingly the majority is enriched in the cytosol and in ribosomal fractions. The ribosomal enriched lncRNAs include H19 and TUG1. CONCLUSIONS: Most studies on lncRNAs have focused on the regulatory function of these transcripts in the nucleus. We demonstrate that only a minority of all lncRNAs are nuclear enriched. Our findings suggest that many lncRNAs may have a function in cytoplasmic processes, and in particular in ribosome complexes.

Identification of Srp9 as a febrile seizure susceptibility gene.

OBJECTIVE: Febrile seizures (FS) are the most common seizure type in young children. Complex FS are a risk factor for mesial temporal lobe epilepsy (mTLE). To identify new FS susceptibility genes we used a forward genetic strategy in mice and subsequently analyzed candidate genes in humans. METHODS: We mapped a quantitative trait locus (QTL1) for hyperthermia-induced FS on mouse chromosome 1, containing the signal recognition particle 9 (Srp9) gene. Effects of differential Srp9 expression were assessed in vivo and in vitro. Hippocampal SRP9 expression and genetic association were analyzed in FS and mTLE patients. RESULTS: Srp9 was differentially expressed between parental strains C57BL/6J and A/J. Chromosome substitution strain 1 (CSS1) mice exhibited lower FS susceptibility and Srp9 expression than C57BL/6J mice. In vivo knockdown of brain Srp9 reduced FS susceptibility. Mice with reduced Srp9 expression and FS susceptibility, exhibited reduced hippocampal AMPA and NMDA currents. Downregulation of neuronal Srp9 reduced surface expression of AMPA receptor subunit GluA1. mTLE patients with antecedent FS had higher SRP9 expression than patients without. SRP9 promoter SNP rs12403575(G/A) was genetically associated with FS and mTLE. INTERPRETATION: Our findings identify SRP9 as a novel FS susceptibility gene and indicate that SRP9 conveys its effects through endoplasmic reticulum (ER)-dependent synthesis and trafficking of membrane proteins, such as glutamate receptors. Discovery of this new FS gene and mechanism may provide new leads for early diagnosis and treatment of children with complex FS at risk for mTLE.

Combined sequence-based and genetic mapping analysis of complex traits in outbred rats.

Genetic mapping on fully sequenced individuals is transforming understanding of the relationship between molecular variation and variation in complex traits. Here we report a combined sequence and genetic mapping analysis in outbred rats that maps 355 quantitative trait loci for 122 phenotypes. We identify 35 causal genes involved in 31 phenotypes, implicating new genes in models of anxiety, heart disease and multiple sclerosis. The relationship between sequence and genetic variation is unexpectedly complex: at approximately 40% of quantitative trait loci, a single sequence variant cannot account for the phenotypic effect. Using comparable sequence and mapping data from mice, we show that the extent and spatial pattern of variation in inbred rats differ substantially from those of inbred mice and that the genetic variants in orthologous genes rarely contribute to the same phenotype in both species.

Systematic biases in DNA copy number originate from isolation procedures.

BACKGROUND: The ability to accurately detect DNA copy number variation in both a sensitive and quantitative manner is important in many research areas. However, genome-wide DNA copy number analyses are complicated by variations in detection signal. RESULTS: While GC content has been used to correct for this, here we show that coverage biases are tissue-specific and independent of the detection method as demonstrated by next-generation sequencing and array CGH. Moreover, we show that DNA isolation stringency affects the degree of equimolar coverage and that the observed biases coincide with chromatin characteristics like gene expression, genomic isochores, and replication timing. CONCLUSION: These results indicate that chromatin organization is a main determinant for differential DNA retrieval. These findings are highly relevant for germline and somatic DNA copy number variation analyses.

A systematic genome-wide analysis of zebrafish protein-coding gene function.

Since the publication of the human reference genome, the identities of specific genes associated with human diseases are being discovered at a rapid rate. A central problem is that the biological activity of these genes is often unclear. Detailed investigations in model vertebrate organisms, typically mice, have been essential for understanding the activities of many orthologues of these disease-associated genes. Although gene-targeting approaches and phenotype analysis have led to a detailed understanding of nearly 6,000 protein-coding genes, this number falls considerably short of the more than 22,000 mouse protein-coding genes. Similarly, in zebrafish genetics, one-by-one gene studies using positional cloning, insertional mutagenesis, antisense morpholino oligonucleotides, targeted re-sequencing, and zinc finger and TAL endonucleases have made substantial contributions to our understanding of the biological activity of vertebrate genes, but again the number of genes studied falls well short of the more than 26,000 zebrafish protein-coding genes. Importantly, for both mice and zebrafish, none of these strategies are particularly suited to the rapid generation of knockouts in thousands of genes and the assessment of their biological activity. Here we describe an active project that aims to identify and phenotype the disruptive mutations in every zebrafish protein-coding gene, using a well-annotated zebrafish reference genome sequence, high-throughput sequencing and efficient chemical mutagenesis. So far we have identified potentially disruptive mutations in more than 38% of all known zebrafish protein-coding genes. We have developed a multi-allelic phenotyping scheme to efficiently assess the effects of each allele during embryogenesis and have analysed the phenotypic consequences of over 1,000 alleles. All mutant alleles and data are available to the community and our phenotyping scheme is adaptable to phenotypic analysis beyond embryogenesis.

Genetic basis of transcriptome differences between the founder strains of the rat HXB/BXH recombinant inbred panel.

BACKGROUND: With the advent of next generation sequencing it has become possible to detect genomic variation on a large scale. However, predicting which genomic variants are damaging to gene function remains a challenge, as knowledge of the effects of genomic variation on gene expression is still limited. Recombinant inbred panels are powerful tools to study the cis and trans effects of genetic variation on molecular phenotypes such as gene expression. RESULTS: We generated a comprehensive inventory of genomic differences between the two founder strains of the rat HXB/BXH recombinant inbred panel: SHR/OlaIpcv and BN-Lx/Cub. We identified 3.2 million single nucleotide variants, 425,924 small insertions and deletions, 907 copy number changes and 1,094 large structural genetic variants. RNA-sequencing analyses on liver tissue of the two strains identified 532 differentially expressed genes and 40 alterations in transcript structure. We identified both coding and non-coding variants that correlate with differential expression and alternative splicing. Furthermore, structural variants, in particular gene duplications, show a strong correlation with transcriptome alterations. CONCLUSIONS: We show that the panel is a good model for assessing the genetic basis of phenotypic heterogeneity and for providing insights into possible underlying molecular mechanisms. Our results reveal a high diversity and complexity underlying quantitative and qualitative transcriptional differences.