Optical Genome Mapping as a Potential Routine Clinical Diagnostic Method.

Chromosome analysis (CA) and chromosomal microarray analysis (CMA) have been successfully used to diagnose genetic disorders. However, many conditions remain undiagnosed due to limitations in resolution (CA) and detection of only unbalanced events (CMA). Optical genome mapping (OGM) has the potential to address these limitations by capturing both structural variants (SVs) resulting in copy number changes and balanced rearrangements with high resolution. In this study, we investigated OGM's concordance using 87 SVs previously identified by CA, CMA, or Southern blot. Overall, OGM was 98% concordant with only three discordant cases: (1) uncalled translocation with one breakpoint in a centromere; (2) uncalled duplication with breakpoints in the pseudoautosomal region 1; and (3) uncalled mosaic triplication originating from a marker chromosome. OGM provided diagnosis for three previously unsolved cases: (1) disruption of the SON gene due to a balanced reciprocal translocation; (2) disruption of the NBEA gene due to an inverted insertion; (3) disruption of the TSC2 gene due to a mosaic deletion. We show that OGM is a valid method for the detection of many types of SVs in a single assay and is highly concordant with legacy cytogenomic methods; however, it has limited SV detection capabilities in centromeric and pseudoautosomal regions.

Closing the Gap - Detection of 5q-Spinal Muscular Atrophy by Short-Read Next-Generation Sequencing and Unexpected Results in a Diagnostic Patient Cohort.

BACKGROUND: The importance of early diagnosis of 5q-Spinal muscular atrophy (5q-SMA) has heightened as early intervention can significantly improve clinical outcomes. In 96% of cases, 5q-SMA is caused by a homozygous deletion of SMN1. Around 4 % of patients carry a SMN1 deletion and a single-nucleotide variant (SNV) on the other allele. Traditionally, diagnosis is based on multiplex ligation probe amplification (MLPA) to detect homozygous or heterozygous exon 7 deletions in SMN1. Due to high homologies within the SMN1/SMN2 locus, sequence analysis to identify SNVs of the SMN1 gene is unreliable by standard Sanger or short-read next-generation sequencing (srNGS) methods. OBJECTIVE: The objective was to overcome the limitations in high-throughput srNGS with the aim of providing SMA patients with a fast and reliable diagnosis to enable their timely therapy. METHODS: A bioinformatics workflow to detect homozygous SMN1 deletions and SMN1 SNVs on srNGS analysis was applied to diagnostic whole exome and panel testing for suggested neuromuscular disorders (1684 patients) and to fetal samples in prenatal diagnostics (260 patients). SNVs were detected by aligning sequencing reads from SMN1 and SMN2 to an SMN1 reference sequence. Homozygous SMN1 deletions were identified by filtering sequence reads for the ,, gene-determining variant" (GDV). RESULTS: 10 patients were diagnosed with 5q-SMA based on (i) SMN1 deletion and hemizygous SNV (2 patients), (ii) homozygous SMN1 deletion (6 patients), and (iii) compound heterozygous SNVs in SMN1 (2 patients). CONCLUSIONS: Applying our workflow in srNGS-based panel and whole exome sequencing (WES) is crucial in a clinical laboratory, as otherwise patients with an atypical clinical presentation initially not suspected to suffer from SMA remain undiagnosed.

Transcript capture and ultradeep long-read RNA sequencing (CAPLRseq) to diagnose HNPCC/Lynch syndrome.

PURPOSE: Whereas most human genes encode multiple mRNA isoforms with distinct function, clinical workflows for assessing this heterogeneity are not readily available. This is a substantial shortcoming, considering that up to 25% of disease-causing gene variants are suspected of disrupting mRNA splicing or mRNA abundance. Long-read sequencing can readily portray mRNA isoform diversity, but its sensitivity is relatively low due to insufficient transcriptome penetration. METHODS: We developed and applied capture-based target enrichment from patient RNA samples combined with Oxford Nanopore long-read sequencing for the analysis of 123 hereditary cancer transcripts (capture and ultradeep long-read RNA sequencing (CAPLRseq)). RESULTS: Validating CAPLRseq, we confirmed 17 cases of hereditary non-polyposis colorectal cancer/Lynch syndrome based on the demonstration of splicing defects and loss of allele expression of mismatch repair genes MLH1, PMS2, MSH2 and MSH6. Using CAPLRseq, we reclassified two variants of uncertain significance in MSH6 and PMS2 as either likely pathogenic or benign. CONCLUSION: Our data show that CAPLRseq is an automatable and adaptable workflow for effective transcriptome-based identification of disease variants in a clinical diagnostic setting.

Parallel in-depth analysis of repeat expansions in ataxia patients by long-read sequencing.

Instability of simple DNA repeats has been known as a common cause of hereditary ataxias for over 20 years. Routine genetic diagnostics of these phenotypically similar diseases still rely on an iterative workflow for quantification of repeat units by PCR-based methods of limited precision. We established and validated clinical nanopore Cas9-targeted sequencing, an amplification-free method for simultaneous analysis of 10 repeat loci associated with clinically overlapping hereditary ataxias. The method combines target enrichment by CRISPR-Cas9, Oxford Nanopore long-read sequencing and a bioinformatics pipeline using the tools STRique and Megalodon for parallel detection of length, sequence, methylation and composition of the repeat loci. Clinical nanopore Cas9-targeted sequencing allowed for the precise and parallel analysis of 10 repeat loci associated with adult-onset ataxia and revealed additional parameter such as FMR1 promotor methylation and repeat sequence required for diagnosis at the same time. Using clinical nanopore Cas9-targeted sequencing we analysed 100 clinical samples of undiagnosed ataxia patients and identified causative repeat expansions in 28 patients. Parallel repeat analysis enabled a molecular diagnosis of ataxias independent of preconceptions on the basis of clinical presentation. Biallelic expansions within RFC1 were identified as the most frequent cause of ataxia. We characterized the RFC1 repeat composition of all patients and identified a novel repeat motif, AGGGG. Our results highlight the power of clinical nanopore Cas9-targeted sequencing as a readily expandable workflow for the in-depth analysis and diagnosis of phenotypically overlapping repeat expansion disorders.

Methylation of the 4q35 D4Z4 repeat defines disease status in facioscapulohumeral muscular dystrophy.

Genetic diagnosis of facioscapulohumeral muscular dystrophy (FSHD) remains a challenge in clinical practice as it cannot be detected by standard sequencing methods despite being the third most common muscular dystrophy. The conventional diagnostic strategy addresses the known genetic parameters of FSHD: the required presence of a permissive haplotype, a size reduction of the D4Z4 repeat of chromosome 4q35 (defining FSHD1) or a pathogenic variant in an epigenetic suppressor gene (consistent with FSHD2). Incomplete penetrance and epistatic effects of the underlying genetic parameters as well as epigenetic parameters (D4Z4 methylation) pose challenges to diagnostic accuracy and hinder prediction of clinical severity. In order to circumvent the known limitations of conventional diagnostics and to complement genetic parameters with epigenetic ones, we developed and validated a multistage diagnostic workflow that consists of a haplotype analysis and a high-throughput methylation profile analysis (FSHD-MPA). FSHD-MPA determines the average global methylation level of the D4Z4 repeat array as well as the regional methylation of the most distal repeat unit by combining bisulphite conversion with next-generation sequencing and a bioinformatics pipeline and uses these as diagnostic parameters. We applied the diagnostic workflow to a cohort of 148 patients and compared the epigenetic parameters based on FSHD-MPA to genetic parameters of conventional genetic testing. In addition, we studied the correlation of repeat length and methylation level within the most distal repeat unit with age-corrected clinical severity and age at disease onset in FSHD patients. The results of our study show that FSHD-MPA is a powerful tool to accurately determine the epigenetic parameters of FSHD, allowing discrimination between FSHD patients and healthy individuals, while simultaneously distinguishing FSHD1 and FSHD2. The strong correlation between methylation level and clinical severity indicates that the methylation level determined by FSHD-MPA accounts for differences in disease severity among individuals with similar genetic parameters. Thus, our findings further confirm that epigenetic parameters rather than genetic parameters represent FSHD disease status and may serve as a valuable biomarker for disease status.

Somatic copy number alteration and fragmentation analysis in circulating tumor DNA for cancer screening and treatment monitoring in colorectal cancer patients.

BACKGROUND: Analysis of circulating free DNA (cfDNA) is a promising tool for personalized management of colorectal cancer (CRC) patients. Untargeted cfDNA analysis using whole-genome sequencing (WGS) does not need a priori knowledge of the patient´s mutation profile. METHODS: Here we established LIquid biopsy Fragmentation, Epigenetic signature and Copy Number Alteration analysis (LIFE-CNA) using WGS with ~ 6× coverage for detection of circulating tumor DNA (ctDNA) in CRC patients as a marker for CRC detection and monitoring. RESULTS: We describe the analytical validity and a clinical proof-of-concept of LIFE-CNA using a total of 259 plasma samples collected from 50 patients with stage I-IV CRC and 61 healthy controls. To reliably distinguish CRC patients from healthy controls, we determined cutoffs for the detection of ctDNA based on global and regional cfDNA fragmentation patterns, transcriptionally active chromatin sites, and somatic copy number alterations. We further combined global and regional fragmentation pattern into a machine learning (ML) classifier to accurately predict ctDNA for cancer detection. By following individual patients throughout their course of disease, we show that LIFE-CNA enables the reliable prediction of response or resistance to treatment up to 3.5 months before commonly used CEA. CONCLUSION: In summary, we developed and validated a sensitive and cost-effective method for untargeted ctDNA detection at diagnosis as well as for treatment monitoring of all CRC patients based on genetic as well as non-genetic tumor-specific cfDNA features. Thus, once sensitivity and specificity have been externally validated, LIFE-CNA has the potential to be implemented into clinical practice. To the best of our knowledge, this is the first study to consider multiple genetic and non-genetic cfDNA features in combination with ML classifiers and to evaluate their potential in both cancer detection and treatment monitoring. Trial registration DRKS00012890.

Constitutional chromothripsis of the APC locus as a cause of genetic predisposition to colon cancer.

PURPOSE: Approximately 20% of patients with clinical familial adenomatous polyposis (FAP) remain unsolved after molecular genetic analysis of the APC and other polyposis genes, suggesting additional pathomechanisms. METHODS: We applied multidimensional genomic analysis employing chromosomal microarray profiling, optical mapping, long-read genome and RNA sequencing combined with FISH and standard PCR of genomic and complementary DNA to decode a patient with an attenuated FAP that had remained unsolved by Sanger sequencing and multigene panel next-generation sequencing for years. RESULTS: We identified a complex 3.9 Mb rearrangement involving 14 fragments from chromosome 5q22.1q22.3 of which three were lost, 1 reinserted into chromosome 5 and 10 inserted into chromosome 10q21.3 in a seemingly random order and orientation thus fulfilling the major criteria of chromothripsis. The rearrangement separates APC promoter 1B from the coding ORF (open reading frame) thus leading to allele-specific downregulation of APC mRNA. The rearrangement also involves three additional genes implicated in the APC-Axin-GSK3B-ß-catenin signalling pathway. CONCLUSIONS: Based on comprehensive genomic analysis, we propose that constitutional chromothripsis dampening APC expression, possibly modified by additional APC-Axin-GSK3B-ß-catenin pathway disruptions, underlies the patient's clinical phenotype. The combinatorial approach we deployed provides a powerful tool set for deciphering unsolved familial polyposis and potentially other tumour syndromes and monogenic diseases.

Highly sensitive liquid biopsy Duplex sequencing complements tissue biopsy to enhance detection of clinically relevant genetic variants.

Background: Liquid biopsy (LB) is a promising complement to tissue biopsy for detection of clinically relevant genetic variants in cancer and mosaic diseases. A combined workflow to enable parallel tissue and LB analysis is required to maximize diagnostic yield for patients. Methods: We developed and validated a cost-efficient combined next-generation sequencing (NGS) workflow for both tissue and LB samples, and applied Duplex sequencing technology for highly accurate detection of low frequency variants in plasma. Clinically relevant cutoffs for variant reporting and quantification were established. Results: We investigated assay performance characteristics for very low amounts of clinically relevant variants. In plasma, the assay achieved 100% sensitivity and 92.3% positive predictive value (PPV) for single nucleotide variants (SNVs) and 91.7% sensitivity and 100% PPV for insertions and deletions (InDel) in clinically relevant hotspots with 0.5-5% variant allele frequencies (VAFs). We further established a cutoff for reporting variants (i.e. Limit of Blank, LOB) at 0.25% VAF and a cutoff for quantification (i.e. Limit of Quantification, LOQ) at 5% VAF in plasma for accurate clinical interpretation of analysis results. With our LB approach, we were able to identify the molecular cause of a clinically confirmed asymmetric overgrowth syndrome in a 10-year old child that would have remained undetected with tissue analysis as well as other molecular diagnostic approaches. Conclusion: Our flexible and cost-efficient workflow allows analysis of both tissue and LB samples and provides clinically relevant cutoffs for variant reporting and precise quantification. Complementing tissue analysis by LB is likely to increase diagnostic yield for patients with molecular diseases.

HPO-driven virtual gene panel: a new efficient approach in molecular autopsy of sudden unexplained death.

BACKGROUND: Molecular autopsy represents an efficient tool to save the diagnosis in up to one-third of sudden unexplained death (SUD). A defined gene panel is usually used for the examination. Alternatively, it is possible to carry out a comprehensive genetic assessment (whole exome sequencing, WES), which also identifies rare, previously unknown variants. The disadvantage is that a dramatic number of variants must be assessed to identify the causal variant. To improve the evaluation of WES, the human phenotype ontology (HPO) annotation is used internationally for deep phenotyping in the field of rare disease. However, a HPO-based evaluation of WES in SUD has not been described before. METHODS: We performed WES in tissue samples from 16 people after SUD. Instead of a fixed gene panel, we defined a set of HPO terms and thus created a flexible "virtual gene panel", with the advantage, that recently identified genes are automatically associated by HPO terms in the HPO database. RESULTS: We obtained a mean value of 68,947 variants per sample. Stringent filtering ended up in a mean value of 276 variants per sample. Using the HPO-driven virtual gene panel we developed an algorithm that prioritized 1.4% of the variants. Variant interpretation resulted in eleven potentially causative variants in 16 individuals. CONCLUSION: Our data introduce an effective diagnostic procedure in molecular autopsy of SUD with a non-specific clinical phenotype.

Critical assessment of secondary findings in genes linked to primary arrhythmia syndromes.

As comprehensive sequencing technologies gain widespread use, questions about so-called secondary findings (SF) require urgent consideration. The American College of Medical Genetics and Genomics has recommended to report SF in 59 genes (ACMG SF v2.0) including four actionable genes associated with inherited primary arrhythmia syndromes (IPAS) such as catecholaminergic polymorphic ventricular tachycardia, long QT syndrome, and Brugada syndrome. Databases provide conflicting results for the purpose of identifying pathogenic variants in SF associated with IPAS at a level of sufficient evidence for clinical return. As IPAS account for a significant proportion of sudden cardiac deaths (SCD) in young and apparently healthy individuals, variant interpretation has a great impact on diagnosis and prevention of disease. Of 6381 individuals, 0.4% carry pathogenic variants in one of the four actionable genes related to IPAS: RYR2, KCNQ1, KCNH2, and SCN5A. Comparison of the databases ClinVar, Leiden Open-source Variant Database, and Human Gene Mutation Database showed impactful differences (0.2% to 1.3%) in variant interpretation improvable by expert-curation depending on database and classification system used. These data further highlight the need for international consensus regarding the variant interpretation, and subsequently management of SF in particular with regard to treatable arrhythmic disorders with increased risk of SCD.

Mitochondrial and nuclear disease panel (Mito-aND-Panel): Combined sequencing of mitochondrial and nuclear DNA by a cost-effective and sensitive NGS-based method.

BACKGROUND: The diagnosis of mitochondrial disorders is challenging because of the clinical variability and genetic heterogeneity of these conditions. Next-Generation Sequencing (NGS) technology offers a robust high-throughput platform for nuclear and mitochondrial DNA (mtDNA) analyses. METHOD: We developed a custom Agilent SureSelect Mitochondrial and Nuclear Disease Panel (Mito-aND-Panel) capture kit that allows parallel enrichment for subsequent NGS-based sequence analysis of nuclear mitochondrial disease-related genes and the complete mtDNA genome. Sequencing of enriched mtDNA simultaneously with nuclear genes was compared with the separated sequencing of the mitochondrial genome and whole exome sequencing (WES). RESULTS: The Mito-aND-Panel permits accurate detection of low-level mtDNA heteroplasmy due to a very high sequencing depth compared to standard diagnostic procedures using Sanger sequencing/SNaPshot and WES which is crucial to identify maternally inherited mitochondrial disorders. CONCLUSION: We established a NGS-based method with combined sequencing of the complete mtDNA and nuclear genes which enables a more sensitive heteroplasmy detection of mtDNA mutations compared to traditional methods. Because the method promotes the analysis of mtDNA variants in large cohorts, it is cost-effective and simple to setup, we anticipate this is a highly relevant method for sequence-based genetic diagnosis in clinical diagnostic applications.