Detection Of Structural Variants Research Articles

Non-invasive monitoring of neoadjuvant radiation therapy response in soft tissue sarcomas by multiparametric MRI and quantification of circulating tumor DNA-A study protocol.

Wide resection remains the cornerstone of localized soft-tissue sarcomas (STS) treatment. Neoadjuvant radiation therapy (NRT) may decrease the risk of local recurrences; however, its effectiveness for different histological STS subtypes has not been systematically investigated. The proposed prospective study evaluates the NRT response in STS using liquid biopsies and the correlation of multiparametric magnetic resonance imaging (mpMRI) with histopathology and immunohistochemistry. Patients with localized high-grade STS, who qualify for NRT, are included in this study. Quantification of circulating tumor DNA (ctDNA) in patient blood samples is performed by targeted next-generation sequencing. Soft-tissue sarcoma subtype-specific panel sequencing in combination with patient-specific exome sequencing allows the detection of individual structural variants and point mutations. Circulating free DNA is isolated from peritherapeutically collected patient plasma samples and ctDNA quantified therein. Identification of breakpoints is carried out using FACTERA. Bioinformatic analysis is performed using samtools, picard, fgbio, and the MIRACUM Pipeline. Combination of conventional MRI sequences with diffusion-weighted imaging, intravoxel-incoherent motion, and dynamic contrast enhancement. Multiparametric MRI is performed before, during, and after NRT. We aim to correlate mpMRI data with the resected specimen's macroscopical, histological, and immunohistochemical findings. Preliminary data support the notion that quantification of ctDNA in combination with tumor mass characterization through co-registration of mpMRI and histopathology can predict NRT response of STS. The methods presented in this prospective study are necessary to assess therapy response in heterogeneous tumors and lay the foundation of future patient- and tumor-specific therapy concepts. These methods can be applied to various tumor entities. Thus, the participation and support of a wider group of oncologic surgeons are needed to validate these findings on a larger patient cohort.

MosaiCatcher v2: a single-cell structural variations detection and analysis reference framework based on Strand-seq.

Single-cell DNA template strand sequencing (Strand-seq) allows a range of various genomic analysis including chromosome length haplotype phasing and structural variation (SV) calling in individual cells. Here, we present MosaiCatcher v2, a standardized workflow and reference framework for single-cell SV detection using Strand-seq. This framework introduces a range of functionalities, including: an automated upstream Quality Control (QC) and assembly sub-workflow that relies on multiple genome assemblies and incorporates a multistep normalization module, integration of the single-cell nucleosome occupancy and genetic variation analysis SV functional characterization and of the ArbiGent SV genotyping modules, platform portability, as well as a user-friendly and shareable web report. These new features of MosaiCatcher v2 enable reproducible computational processing of Strand-seq data, which are increasingly used in human genetics and single-cell genomics, toward production environments. MosaiCatcher v2 is compatible with both container and conda environments, ensuring reproducibility and robustness and positioning the framework as a cornerstone in computational processing of Strand-seq data. MosaiCatcher v2 is a standardized workflow, implemented using the Snakemake workflow management system. The pipeline is available on GitHub: https://github.com/friendsofstrandseq/mosaicatcher-pipeline/ and on the snakemake-workflow-catalog: https://snakemake.github.io/snakemake-workflow-catalog/?usage=friendsofstrandseq/mosaicatcher-pipeline. Strand-seq example input data used in the publication can be found in the Data availability statement. Additionally, a lightweight dataset for test purposes can be found on the GitHub repository.

Design of optimal labeling patterns for optical genome mapping via information theory.

Optical genome mapping (OGM) is a technique that extracts partial genomic information from optically imaged and linearized DNA fragments containing fluorescently labeled short sequence patterns. This information can be used for various genomic analyses and applications, such as the detection of structural variations and copy-number variations, epigenomic profiling, and microbial species identification. Currently, the choice of labeled patterns is based on the available biochemical methods and is not necessarily optimized for the application. In this work, we develop a model of OGM based on information theory, which enables the design of optimal labeling patterns for specific applications and target organism genomes. We validated the model through experimental OGM on human DNA and simulations on bacterial DNA. Our model predicts up to 10-fold improved accuracy by optimal choice of labeling patterns, which may guide future development of OGM biochemical labeling methods and significantly improve its accuracy and yield for applications such as epigenomic profiling and cultivation-free pathogen identification in clinical samples. https://github.com/yevgenin/PatternCode.

HQAlign: aligning nanopore reads for SV detection using current-level modeling.

Detection of structural variants (SVs) from the alignment of sample DNA reads to the reference genome is an important problem in understanding human diseases. Long reads that can span repeat regions, along with an accurate alignment of these long reads play an important role in identifying novel SVs. Long-read sequencers, such as nanopore sequencing, can address this problem by providing very long reads but with high error rates, making accurate alignment challenging. Many errors induced by nanopore sequencing have a bias because of the physics of the sequencing process and proper utilization of these error characteristics can play an important role in designing a robust aligner for SV detection problems. In this article, we design and evaluate HQAlign, an aligner for SV detection using nanopore sequenced reads. The key ideas of HQAlign include (i) using base-called nanopore reads along with the nanopore physics to improve alignments for SVs, (ii) incorporating SV-specific changes to the alignment pipeline, and (iii) adapting these into existing state-of-the-art long-read aligner pipeline, minimap2 (v2.24), for efficient alignments. We show that HQAlign captures about 4%-6% complementary SVs across different datasets, which are missed by minimap2 alignments while having a standalone performance at par with minimap2 for real nanopore reads data. For the common SV calls between HQAlign and minimap2, HQAlign improves the start and the end breakpoint accuracy by about 10%-50% for SVs across different datasets. Moreover, HQAlign improves the alignment rate to 89.35% from minimap2 85.64% for nanopore reads alignment to recent telomere-to-telomere CHM13 assembly, and it improves to 86.65% from 83.48% for nanopore reads alignment to GRCh37 human genome. https://github.com/joshidhaivat/HQAlign.git.

Long-read whole-genome analysis of human single cells

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

AML with complex karyotype: extreme genomic complexity revealed by combined long-read sequencing and Hi-C technology

AbstractAcute myeloid leukemia with complex karyotype (CK-AML) is associated with poor prognosis, which is only in part explained by underlying TP53 mutations. Especially in the presence of complex chromosomal rearrangements, such as chromothripsis, the outcome of CK-AML is dismal. However, this degree of complexity of genomic rearrangements contributes to the leukemogenic phenotype and treatment resistance of CK-AML remains largely unknown. Applying an integrative workflow for the detection of structural variants (SVs) based on Oxford Nanopore (ONT) genomic DNA long-read sequencing (gDNA-LRS) and high-throughput chromosome confirmation capture (Hi-C) in a well-defined cohort of CK-AML identified regions with an extreme density of SVs. These rearrangements consisted to a large degree of focal amplifications enriched in the proximity of mammalian-wide interspersed repeat elements, which often result in oncogenic fusion transcripts, such as USP7::MVD, or the deregulation of oncogenic driver genes as confirmed by RNA-seq and ONT direct complementary DNA sequencing. We termed this novel phenomenon chromocataclysm. Thus, our integrative SV detection workflow combing gDNA-LRS and Hi-C enables to unravel complex genomic rearrangements at a very high resolution in regions hard to analyze by conventional sequencing technology, thereby providing an important tool to identify novel important drivers underlying cancer with complex karyotypic changes.

Resolving complex structural variants via nanopore sequencing.

The recent development of high-throughput sequencing platforms provided impressive insights into the field of human genetics and contributed to considering structural variants (SVs) as the hallmark of genome instability, leading to the establishment of several pathologic conditions, including neoplasia and neurodegenerative and cognitive disorders. While SV detection is addressed by next-generation sequencing (NGS) technologies, the introduction of more recent long-read sequencing technologies have already been proven to be invaluable in overcoming the inaccuracy and limitations of NGS technologies when applied to resolve wide and structurally complex SVs due to the short length (100-500bp) of the sequencing read utilized. Among the long-read sequencing technologies, Oxford Nanopore Technologies developed a sequencing platform based on a protein nanopore that allows the sequencing of "native" long DNA molecules of virtually unlimited length (typical range 1-100Kb). In this review, we focus on the bioinformatics methods that improve the identification and genotyping of known and novel SVs to investigate human pathological conditions, discussing the possibility of introducing nanopore sequencing technology into routine diagnostics.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches.

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Linked-read based analysis of the medulloblastoma genome.

Medulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies. Here, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements. We demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it. We detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

Structural variants integration and visualization: A comprehensive R package for integration of somatic structural variations from multiple callers and visualization

Whole genome sequencing (WGS) emerges as a powerful tool for detecting structural variations (SVs) in genomes. However, different SV callers can produce variable results due to the distinct rationale and sensitivity of pipelines, highlighting the need for effective tools to compare and merge results from multiple callers. Here, we developed an R package, structural variants integration and visualization, to facilitate the integration, classification, and visualization of SV results from multiple callers, allowing for accurate identification of the most reliable SVs. Our package relies on a complex translocation projection and clustering method, enabling the projection of each translation to a point in a Cartesian coordinate system and visualization of SVs at both whole-genome and individual chromosome levels. Thus, our approach provides a valuable framework for analyzing SVs from WGS data, improving the accuracy and efficiency of SV detection, and enhancing the potential of WGS for clinical and research applications.

A survey of algorithms for the detection of genomic structural variants from long-read sequencing data.

As long-read sequencing technologies are becoming increasingly popular, a number of methods have been developed for the discovery and analysis of structural variants (SVs) from long reads. Long reads enable detection of SVs that could not be previously detected from short-read sequencing, but computational methods must adapt to the unique challenges and opportunities presented by long-read sequencing. Here, we summarize over 50 long-read-based methods for SV detection, genotyping and visualization, and discuss how new telomere-to-telomere genome assemblies and pangenome efforts can improve the accuracy and drive the development of SV callers in the future.

A powerful random forest-based pipeline for accurately identifying the structural variants on panel sequencing data.

e18898 Background: Numerous studies have reported that the identification of structural variants (SVs) can provide ideal targets for cancer targeting therapy. However, on the one hand, the current algorithms exhibited different performances in detecting different types and sizes of SVs; on the other hand, the whole-genome/exome sequencing for detecting SVs is costly, and only a portion of SVs can be used to guide the clinical treatment of cancer patients. Therefore, developing a powerful, highly accurate, and cost-effective pipeline to detect SVs for clinical application is urgently required. Methods: Based on simulated sequencing data, four common SV detection tools (Delly, Lumpy, SvABA, and Manta) were used to select the best performance one. This tool will be used for detecting SVs in panel sequencing data. Then, Integrative Genomics Viewer (IGV) was used to annotate these SVs as true positive (TP) and false positive (FP) SVs. These annotated TP and FP SVs were used as input in the random forest classifier (RFC) to train a model for predicting the true or false positive SVs. Two independent testing cohorts and standards were used to validate this constructed pipeline. In addition, IHC/FISH experiments were performed to further validate these predicted TP SVs. Results: In the simulation data, the detection tool with the highest SVs identification sensitivity was Delly. By Delly, a total of 1,303 SVs in 384 tumor samples were detected. Based on annotated TP and FP SVs, an RFC model was constructed in the training cohort, which predicted 334 and 560 TP and FP SVs, respectively. The accuracy rates (AR) were 99.85% and 100% of TP and FP SVs, respectively. The predicted TP and FP SVs by constructed pipeline also showed a high accuracy rate in cohort 1 (98.96% and 99.24% of TP and FP SVs, respectively), cohort 2 (98.48% and 100% of TP and FP SVs, respectively), and standards (the mean accuracy rate was 90%). Finally, based on these TP SVs predicted by the constructed pipeline, FISH/IHC experiment results further verified the robustness of this pipeline (the man accuracy rate was 95.83%). Conclusions: The constructed random forest-based pipeline can robustly and accurately identify SVs based on panel sequencing data, which will be helpful for clinical decision-making.

De Novo Structural Variations of Escherichia coli Detected by Nanopore Long-Read Sequencing.

Spontaneous mutations power evolution, whereas large-scale structural variations (SVs) remain poorly studied, primarily because of the lack of long-read sequencing techniques and powerful analytical tools. Here, we explore the SVs of Escherichia coli by running 67 wild-type (WT) and 37 mismatch repair (MMR)-deficient (ΔmutS) mutation accumulation lines, each experiencing more than 4,000 cell divisions, by applying Nanopore long-read sequencing and Illumina PE150 sequencing and verifying the results by Sanger sequencing. In addition to precisely repeating previous mutation rates of base-pair substitutions and insertion and deletion (indel) mutation rates, we do find significant improvement in insertion and deletion detection using long-read sequencing. The long-read sequencing and corresponding software can particularly detect bacterial SVs in both simulated and real data sets with high accuracy. These lead to SV rates of 2.77 × 10-4 (WT) and 5.26 × 10-4 (MMR-deficient) per cell division per genome, which is comparable with previous reports. This study provides the SV rates of E. coli by applying long-read sequencing and SV detection programs, revealing a broader and more accurate picture of spontaneous mutations in bacteria.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30–40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

Novel NUP98::ASH1L Gene Fusion in Acute Myeloid Leukemia Detected by Optical Genome Mapping.

Optical genome mapping (OGM) recently has demonstrated the potential to improve genetic diagnostics in acute myeloid leukemia (AML). In this study, OGM was utilized as a tool for the detection of genome-wide structural variants and disease monitoring. A previously unrecognized NUP98::ASH1L fusion was detected in an adult patient with secondary AML. OGM identified the fusion of NUP98 to Absent, Small, or Homeotic-Like Histone Lysine Methyltransferase (ASH1L) as result of a complex structural rearrangement between chromosomes 1 and 11. A pipeline for the measurement of rare structural variants (Rare Variant Pipeline, Bionano Genomics, San Diego, CA, USA) was used for detection. As NUP98 and other fusions are relevant for disease classification, this demonstrates the necessity for methods such as OGM for cytogenetic diagnostics in AML. Furthermore, other structural variants showed discordant variant allele frequencies at different time points over the course of the disease and treatment pressure, indicating clonal evolution. These results support OGM to be a valuable tool for primary diagnostics in AML as well as longitudinal testing for disease monitoring and deepening our understanding of genetically heterogenous diseases.

The AnnotSV webserver in 2023: updated visualization and ranking.

Much of the human genetics variant repertoire is composed of single nucleotide variants (SNV) and small insertion/deletions (indel) but structural variants (SV) remain a major part of our modified DNA. SV detection has often been a complex question to answer either because of the necessity to use different technologies (array CGH, SNP array, Karyotype, Optical Genome Mapping…) to detect each category of SV or to get an appropriate resolution (Whole Genome Sequencing). Thanks to the deluge of pangenomic analysis, Human geneticists are accumulating SV and their interpretation remains time consuming and challenging. The AnnotSV webserver (https://www.lbgi.fr/AnnotSV/) aims at being an efficient tool to (i) annotate and interpret SV potential pathogenicity in the context of human diseases, (ii) recognize potential false positive variants from all the SV identified and (iii) visualize the patient variants repertoire. The most recent developments in the AnnotSV webserver are: (i) updated annotations sources and ranking, (ii) three novel output formats to allow diverse utilization (analysis, pipelines), as well as (iii) two novel user interfaces including an interactive circos view.

Faster detection of somatic structural variants.

Faster detection of somatic structural variants.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

Comparison and benchmark of structural variants detected from long read and long-read assembly.

Structural variant (SV) detection is essential for genomic studies, and long-read sequencing technologies have advanced our capacity to detect SVs directly from read or de novo assembly, also known as read-based and assembly-based strategy. However, to date, no independent studies have compared and benchmarked the two strategies. Here, on the basis of SVs detected by 20 read-based and eight assembly-based detection pipelines from six datasets of HG002 genome, we investigated the factors that influence the two strategies and assessed their performance with well-curated SVs. We found that up to 80% of the SVs could be detected by both strategies among different long-read datasets, whereas variant type, size, and breakpoint detected by read-based strategy were greatly affected by aligners. For the high-confident insertions and deletions at non-tandem repeat regions, a remarkable subset of them (82% in assembly-based calls and 93% in read-based calls), accounting for around 4000 SVs, could be captured by both reads and assemblies. However, discordance between two strategies was largely caused by complex SVs and inversions, which resulted from inconsistent alignment of reads and assemblies at these loci. Finally, benchmarking with SVs at medically relevant genes, the recall of read-based strategy reached 77% on 5X coverage data, whereas assembly-based strategy required 20X coverage data to achieve similar performance. Therefore, integrating SVs from read and assembly is suggested for general-purpose detection because of inconsistently detected complex SVs and inversions, whereas assembly-based strategy is optional for applications with limited resources.

Abstract LB080: SAVANA: a computational method to characterize structural variation in human cancer genomes using nanopore sequencing

Abstract Whole-genome sequencing (WGS) of human cancers has revealed that structural variation, which refers to the rearrangement of the genome leading to the deletion, amplification of reshuffling of DNA segments ranging from a few hundred bp to entire chromosomes, is a key mutational process in cancer evolution. Notably, pan-cancer analyses have revealed that both simple and complex forms of structural variation are pervasive across diverse human cancers, and often underpin drug resistance and metastasis. To date, the study of cancer genomes has relied on the analysis of short-read WGS on the dominant Illumina platform, which generates short, highly-accurate reads of 100-300bp that allow the study of point mutations at high resolution. However, detection of structural variants (SVs) using short reads is limited, as breakpoints falling in repetitive regions cannot be reliably mapped to the human genome. As a result, our understanding of the patterns and mechanisms underpinning structural variation in cancer genomes remains incomplete. In contrast to short-read sequencing, long-read sequencing technologies, such as Oxford Nanopore and PacBio, permit continuous reading of individual DNA molecules over 10 kilobases, thus providing unparalleled information to resolve SVs in repetitive regions and complex genome rearrangements. However, novel bioinformatics methods that account for the higher error rate of long-read methods are needed to take advantage of their capabilities for cancer genome analysis. Here, we present SAVANA, a novel structural variant caller for long-read sequencing data specifically designed for the analysis of cancer genomes. To identify both somatic and germline SVs, SAVANA takes as input long-read WGS data from a tumor and normal sample pair. SAVANA scans sequencing reads to detect split reads and gapped alignments, which are then clustered to define putative SVs. Next, SAVANA applies a machine learning-informed set of heuristics to remove false positives arising from mapping errors and sequencing artifacts. Extensively validated against a multi-platform truthset, we show that SAVANA identifies a range of somatic rearrangements with high recall and precision, outperforming existing tools while maintaining a lower execution time than competing methods. In patient samples, SAVANA identifies clinically relevant alterations, such as oncogenic gene fusions, with high accuracy. Additionally, SAVANA permits the reconstruction of double minutes, multi-chromosomal chromothripsis events, and SVs mapping to highly repetitive regions, including centromeres. In sum, SAVANA permits the characterization of complex structural variants and can uncover clinically relevant mutations across diverse cancer types with high accuracy. Citation Format: Hillary Elrick, Jose Espejo Valle-Inclan, Katherine E. Trevers, Francesc Muyas, Rita Cascão, Angela Afonso, Cláudia C. Faria, Adrienne M. Flanagan, Isidro Cortés-Ciriano. SAVANA: a computational method to characterize structural variation in human cancer genomes using nanopore sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 2 (Clinical Trials and Late-Breaking Research); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(8_Suppl):Abstract nr LB080.