Current genotyping approaches for single-nucleotide variations rely on short, accurate reads from second-generation sequencing devices. Presently, third-generation sequencing platforms are rapidly ...becoming more widespread, yet approaches for leveraging their long but error-prone reads for genotyping are lacking. Here, we introduce a novel statistical framework for the joint inference of haplotypes and genotypes from noisy long reads, which we term diplotyping. Our technique takes full advantage of linkage information provided by long reads. We validate hundreds of thousands of candidate variants that have not yet been included in the high-confidence reference set of the Genome-in-a-Bottle effort.
A draft human pangenome reference Liao, Wen-Wei; Asri, Mobin; Ebler, Jana ...
Nature,
05/2023, Letnik:
617, Številka:
7960
Journal Article
Recenzirano
Odprti dostop
Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse ...individuals
. These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample.
Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled ...haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.
Human genomes are typically assembled as consensus sequences that lack information on parental haplotypes. Here we describe a reference-free workflow for diploid de novo genome assembly that combines ...the chromosome-wide phasing and scaffolding capabilities of single-cell strand sequencing
with continuous long-read or high-fidelity
sequencing data. Employing this strategy, we produced a completely phased de novo genome assembly for each haplotype of an individual of Puerto Rican descent (HG00733) in the absence of parental data. The assemblies are accurate (quality value > 40) and highly contiguous (contig N50 > 23 Mbp) with low switch error rates (0.17%), providing fully phased single-nucleotide variants, indels and structural variants. A comparison of Oxford Nanopore Technologies and Pacific Biosciences phased assemblies identified 154 regions that are preferential sites of contig breaks, irrespective of sequencing technology or phasing algorithms.
Pangenome Graphs Eizenga, Jordan M; Novak, Adam M; Sibbesen, Jonas A ...
Annual review of genomics and human genetics,
08/2020, Letnik:
21, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Low-cost whole-genome assembly has enabled the collection of haplotype-resolved pangenomes for numerous organisms. In turn, this technological change is encouraging the development of methods that ...can precisely address the sequence and variation described in large collections of related genomes. These approaches often use graphical models of the pangenome to support algorithms for sequence alignment, visualization, functional genomics, and association studies. The additional information provided to these methods by the pangenome allows them to achieve superior performance on a variety of bioinformatic tasks, including read alignment, variant calling, and genotyping. Pangenome graphs stand to become a ubiquitous tool in genomics. Although it is unclear whether they will replace linearreference genomes, their ability to harmoniously relate multiple sequence and coordinate systems will make them useful irrespective of which pangenomic models become most common in the future.
Resolving genomes at haplotype level is crucial for understanding the evolutionary history of polyploid species and for designing advanced breeding strategies. Polyploid phasing still presents ...considerable challenges, especially in regions of collapsing haplotypes.We present WHATSHAP POLYPHASE, a novel two-stage approach that addresses these challenges by (i) clustering reads and (ii) threading the haplotypes through the clusters. Our method outperforms the state-of-the-art in terms of phasing quality. Using a real tetraploid potato dataset, we demonstrate how to assemble local genomic regions of interest at the haplotype level. Our algorithm is implemented as part of the widely used open source tool WhatsHap.
The telomere-to-telomere (T2T) complete human reference has significantly improved our ability to characterize genome structural variation. To understand its impact on inversion polymorphisms, we ...remapped data from 41 genomes against the T2T reference genome and compared it to the GRCh38 reference. We find a ~ 21% increase in sensitivity improving mapping of 63 inversions on the T2T reference. We identify 26 misorientations within GRCh38 and show that the T2T reference is three times more likely to represent the correct orientation of the major human allele. Analysis of 10 additional samples reveals novel rare inversions at chromosomes 15q25.2, 16p11.2, 16q22.1-23.1, and 22q11.21.
The mutational landscape of B‐cell precursor acute lymphoblastic leukemia (BCP‐ALL), the most common pediatric cancer, is not fully described partially because commonly applied short‐read next ...generation sequencing has a limited ability to identify structural variations. By combining comprehensive analysis of structural variants (SVs), single‐nucleotide variants (SNVs), and small insertions‐deletions, new subtype‐defining and therapeutic targets may be detected. We analyzed the landscape of somatic alterations in 60 pediatric patients diagnosed with the most common BCP‐ALL subtypes, ETV6::RUNX1+ and classical hyperdiploid (HD), using conventional cytogenetics, single nucleotide polymorphism (SNP) array, whole exome sequencing (WES), and the novel optical genome mapping (OGM) technique. Ninety‐five percent of SVs detected by cytogenetics and SNP‐array were verified by OGM. OGM detected an additional 677 SVs not identified using the conventional methods, including (subclonal) IKZF1 deletions. Based on OGM, ETV6::RUNX1+ BCP‐ALL harbored 2.7 times more SVs than HD BCP‐ALL, mainly focal deletions. Besides SVs in known leukemia development genes (ETV6, PAX5, BTG1, CDKN2A), we identified 19 novel recurrently altered regions (in n ≥ 3) including 9p21.3 (FOCAD/HACD4), 8p11.21 (IKBKB), 1p34.3 (ZMYM1), 4q24 (MANBA), 8p23.1 (MSRA), and 10p14 (SFMBT2), as well as ETV6::RUNX1+ subtype‐specific SVs (12p13.1 (GPRC5A), 12q24.21 (MED13L), 18q11.2 (MIB1), 20q11.22 (NCOA6)). We detected 3 novel fusion genes (SFMBT2::DGKD, PDS5B::STAG2, and TDRD5::LPCAT2), for which the sequence and expression were validated by long‐read and whole transcriptome sequencing, respectively. OGM and WES identified double hits of SVs and SNVs (ETV6, BTG1, STAG2, MANBA, TBL1XR1, NSD2) in the same patient demonstrating the power of the combined approach to define the landscape of genomic alterations in BCP‐ALL.
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.
Display omitted
•Whole-genome analyses using multiple state-of-the-art techniques•Renewable COLO829 tumor-normal cell line pairs•Generation of a carefully curated and validated somatic SV truth set•Publicly available benchmark dataset for sequencing technology and analysis software
Structural variation (SV) is part of the cancer genome mutational landscape driving tumor characteristics. Due to its nature, detection of somatic SVs is complicated and unresolved. DNA sequencing technology and data analysis tools are developing rapidly. Valle-Inclan et al. describe the generation of a carefully curated and validated somatic SV truth set based on the COLO829 cell line that can be used as a reference for benchmarking existing and novel somatic SV detection techniques.
Abstract
Motivation
With the fast development of sequencing technology, accurate de novo genome assembly is now possible even for larger genomes. Graph-based representations of genomes arise both as ...part of the assembly process, but also in the context of pangenomes representing a population. In both cases, polymorphic loci lead to bubble structures in such graphs. Detecting bubbles is hence an important task when working with genomic variants in the context of genome graphs.
Results
Here, we present a fast general-purpose tool, called BubbleGun, for detecting bubbles and superbubbles in genome graphs. Furthermore, BubbleGun detects and outputs runs of linearly connected bubbles and superbubbles, which we call bubble chains. We showcase its utility on de Bruijn graphs and compare our results to vg’s snarl detection. We show that BubbleGun is considerably faster than vg especially in bigger graphs, where it reports all bubbles in less than 30 min on a human sample de Bruijn graph of around 2 million nodes.
Availability and implementation
BubbleGun is available and documented as a Python3 package at https://github.com/fawaz-dabbaghieh/bubble_gun under MIT license.
Supplementary information
Supplementary data are available at Bioinformatics online.