Gramtools enables multiscale variation analysis with genome graphs.
Ontology highlight
ABSTRACT:
SUBMITTER: Letcher B
PROVIDER: S-EPMC8420074 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Analysis of sequence diversity in the human genome is fundamental for genetic studies. Structural variants (SVs) are frequently omitted in sequence analysis studies, although each has a relatively large impact on the genome. Here, we present GraphTyper2, which uses pangenome graphs to genotype SVs and small variants using short-reads. Comparison to the syndip benchmark dataset shows that our SV genotyping is sensitive and variant segregation in families demonstrates the accuracy of our approach. We demonstrate that incorporating public assembly data into our pipeline greatly improves sensitivity, particularly for large insertions. We validate 6,812 SVs on average per genome using long-read data of 41 Icelanders. We show that GraphTyper2 can simultaneously genotype tens of thousands of whole-genomes by characterizing 60 million small variants and half a million SVs in 49,962 Icelanders, including 80 thousand SVs with high-confidence.
| S-EPMC6881350 | biostudies-literature
Project description:BackgroundGenotyping of sequence variants typically involves, as a first step, the alignment of sequencing reads to a linear reference genome. Because a linear reference genome represents only a small fraction of all the DNA sequence variation within a species, reference allele bias may occur at highly polymorphic or divergent regions of the genome. Graph-based methods facilitate the comparison of sequencing reads to a variation-aware genome graph, which incorporates a collection of non-redundant DNA sequences that segregate within a species. We compared the accuracy and sensitivity of graph-based sequence variant genotyping using the Graphtyper software to two widely-used methods, i.e., GATK and SAMtools, which rely on linear reference genomes using whole-genome sequencing data from 49 Original Braunvieh cattle.ResultsWe discovered 21,140,196, 20,262,913, and 20,668,459 polymorphic sites using GATK, Graphtyper, and SAMtools, respectively. Comparisons between sequence variant genotypes and microarray-derived genotypes showed that Graphtyper outperformed both GATK and SAMtools in terms of genotype concordance, non-reference sensitivity, and non-reference discrepancy. The sequence variant genotypes that were obtained using Graphtyper had the smallest number of Mendelian inconsistencies between sequence-derived single nucleotide polymorphisms and indels in nine sire-son pairs. Genotype phasing and imputation using the Beagle software improved the quality of the sequence variant genotypes for all the tools evaluated, particularly for animals that were sequenced at low coverage. Following imputation, the concordance between sequence- and microarray-derived genotypes was almost identical for the three methods evaluated, i.e., 99.32, 99.46, and 99.24% for GATK, Graphtyper, and SAMtools, respectively. Variant filtration based on commonly used criteria improved genotype concordance slightly but it also decreased sensitivity. Graphtyper required considerably more computing resources than SAMtools but less than GATK.ConclusionsSequence variant genotyping using Graphtyper is accurate, sensitive and computationally feasible in cattle. Graph-based methods enable sequence variant genotyping from variation-aware reference genomes that may incorporate cohort-specific sequence variants, which is not possible with the current implementation of state-of-the-art methods that rely on linear reference genomes.
| S-EPMC6521551 | biostudies-literature
Project description:Volume electron microscopy enables the ultrastructural analysis of biological tissue. Currently, the techniques involving ultramicrotomy (ATUM, ssTEM) allow large fields of view but afford only limited z-resolution, whereas ion beam-milling approaches (FIB-SEM) yield isotropic voxels but are restricted in volume size. Now we present a hybrid method, named ATUM-FIB, which combines the advantages of both approaches. ATUM-FIB is based on serial sectioning of tissue into "semithick" (2-10 ?m) sections collected onto tape. Serial light and electron microscopy allows the identification of regions of interest that are then directly accessible for targeted FIB-SEM. The set of semithick sections thus represents a tissue "library" which provides three-dimensional context information that can be probed "on demand" by local high-resolution analysis. We demonstrate the potential of this technique to reveal the ultrastructure of rare but pathologically important events by identifying microglia contact sites with amyloid plaques in a mouse model of familial Alzheimer's disease.
| S-EPMC7334410 | biostudies-literature
Project description:MotivationThe de Bruijn graph has become a ubiquitous graph model for biological data ever since its initial introduction in the late 1990s. It has been used for a variety of purposes including genome assembly (Zerbino and Birney, 2008; Bankevich et al., 2012; Peng et al., 2012), variant detection (Alipanahi et al., 2020b; Iqbal et al., 2012), and storage of assembled genomes (Chikhi et al., 2016). For this reason, there have been over a dozen methods for building and representing the de Bruijn graph and its variants in a space and time efficient manner.ResultsWith the exception of a few data structures (Muggli et al., 2019; Holley and Melsted, 2020; Crawford et al.,2018), compressed and compact de Bruijn graphs do not allow for the graph to be efficiently updated, meaning that data can be added or deleted. The most recent compressed dynamic de Bruijn graph (Alipanahi et al., 2020a), relies on dynamic bit vectors which are slow in theory and practice. To address this shortcoming, we present a compressed dynamic de Bruijn graph that removes the necessity of dynamic bit vectors by buffering data that should be added or removed from the graph. We implement our method, which we refer to as BufBOSS, and compare its performance to Bifrost, DynamicBOSS, and FDBG. Our experiments demonstrate that BufBOSS achieves attractive trade-offs compared to other tools in terms of time, memory and disk, and has the best deletion performance by an order of magnitude.
| S-EPMC8326735 | biostudies-literature
Project description:Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis.
| S-EPMC7912537 | biostudies-literature
Project description:A graph grammar with parallel replacement of subgraphs, based on the single-pushout approach in graph rewriting, was designed which constructs Cayley graphs of monoids of transformations of a finite set, with permutation groups as a special case. As input, graph-based representations of a finite number of generating transformations have to be specified; they will then correspond to the edge types of the Cayley graph which is the final result of the rewriting process. The grammar has
| S-EPMC7314706 | biostudies-literature
Project description:Despite their widespread use in research, there has not yet been a systematic genomic analysis of human embryonic stem cell (hESC) lines at a single-nucleotide resolution. We therefore performed whole-genome sequencing (WGS) of 143 hESC lines and annotated their single-nucleotide and structural genetic variants. We found that while a substantial fraction of hESC lines contained large deleterious structural variants, finer-scale structural and single-nucleotide variants (SNVs) that are ascertainable only through WGS analyses were present in hESC genomes and human blood-derived genomes at similar frequencies. Moreover, WGS allowed us to identify SNVs associated with cancer and other diseases that could alter cellular phenotypes and compromise the safety of hESC-derived cellular products transplanted into humans. As a resource to enable reproducible hESC research and safer translation, we provide a user-friendly WGS data portal and a data-driven scheme for cell line maintenance and selection.
| S-EPMC8900618 | biostudies-literature
Project description:There is an acute need for better tools to extract knowledge from the growing flood of sequence data. For example, thousands of complete genomes have been sequenced, and their metabolic networks inferred. Such data should enable a better understanding of evolution. However, most existing network analysis methods are based on pair-wise comparisons, and these do not scale to thousands of genomes. Here we propose the use of weighted graphs as a data structure to enable large-scale phylogenetic analysis of networks. We have developed three types of weighted graph for enzymes: taxonomic (these summarize phylogenetic importance), isoenzymatic (these summarize enzymatic variety/redundancy), and sequence-similarity (these summarize sequence conservation); and we applied these types of weighted graph to survey prokaryotic metabolism. To demonstrate the utility of this approach we have compared and contrasted the large-scale evolution of metabolism in Archaea and Eubacteria. Our results provide evidence for limits to the contingency of evolution.
| S-EPMC3949676 | biostudies-literature
Project description:We introduce, and provide examples of, the application of the bond graph formalism to explicitly represent biophysical processes between and within modular biological compartments in ApiNATOMY. In particular, we focus on modelling scenarios from acid-base physiology to link distinct process modalities as bond graphs over an ApiNATOMY circuit of multiscale compartments. The embedding of bond graphs onto ApiNATOMY compartments provides a semantically and mathematically explicit basis for the coherent representation, integration and visualisation of multiscale physiology processes together with the compartmental topology of those biological structures that convey these processes.
| S-EPMC5740224 | biostudies-literature
Project description:BackgroundTo realize the potential of genome engineering therapeutics, tractable strategies must be identified that balance personalized therapy with the need for off-the-shelf availability. We hypothesized that regional clustering of pathogenic variants can inform the design of rational prime editing therapeutics to treat the majority of genetic cardiovascular diseases with a limited number of reagents.MethodsWe collated 2435 high-confidence pathogenic/likely pathogenic (P/LP) variants in 82 cardiovascular disease genes from ClinVar. We assessed the regional density of these variants by defining a regional clustering index. We then combined a highly active base editor with prime editing to demonstrate the feasibility of a P/LP hotspot-directed genome engineering therapeutic strategy in vitro.ResultsP/LP variants in cardiovascular disease genes display higher regional density than rare variants found in the general population. P/LP missense variants displayed higher average regional density than P/LP truncating variants. Following hypermutagenesis at a pathogenic hotspot, mean prime editing efficiency across introduced variants was 57±27%.ConclusionsDesigning therapeutics that target pathogenic hotspots will not only address known missense P/LP variants but also novel P/LP variants identified in these hotspots as well. Moreover, the clustering of P/LP missense rather than truncating variants in these hotspots suggests that prime editing technology is particularly valuable for dominant negative disease. Although prime editing technology in relation to cardiac health continues to improve, this study presents an approach to targeting the most impactful regions of the genome for inherited cardiovascular disease.
| S-EPMC11020015 | biostudies-literature