STRONG: metagenomics strain resolution on assembly graphs.
Ontology highlight
ABSTRACT:
SUBMITTER: Quince C
PROVIDER: S-EPMC8311964 | biostudies-literature | 2021 Jul
REPOSITORIES: biostudies-literature
Publications
STRONG: metagenomics strain resolution on assembly graphs.
Quince Christopher C Nurk Sergey S Raguideau Sebastien S James Robert R Soyer Orkun S OS Summers J Kimberly JK Limasset Antoine A Eren A Murat AM Chikhi Rayan R Darling Aaron E AE
Genome biology 20210726 1
We introduce STrain Resolution ON assembly Graphs (STRONG), which identifies strains de novo, from multiple metagenome samples. STRONG performs coassembly, and binning into metagenome assembled genomes (MAGs), and stores the coassembly graph prior to variant simplification. This enables the subgraphs and their unitig per-sample coverages, for individual single-copy core genes (SCGs) in each MAG, to be extracted. A Bayesian algorithm, BayesPaths, determines the number of strains present, their ha ...[more]
Similar Datasets
Project description:Dental implants are installed in an increasing number of patients. Mucositis and peri-implantitis are common microbial-biofilm-associated diseases affecting the tissues that surround the dental implant and are a major medical and socioeconomic burden. By metagenomic sequencing of the plaque microbiome in different peri-implant health and disease conditions (113 samples from 72 individuals), we found microbial signatures for peri-implantitis and mucositis and defined the peri-implantitis-related complex (PiRC) composed by the 7 most discriminative bacteria. The peri-implantitis microbiome is site specific as contralateral healthy sites resembled more the microbiome of healthy implants, while mucositis was specifically enriched for Fusobacterium nucleatum acting as a keystone colonizer. Microbiome-based machine learning showed high diagnostic and prognostic power for peri-implant diseases and strain-level profiling identified a previously uncharacterized subspecies of F. nucleatum to be particularly associated with disease. Altogether, we associated the plaque microbiome with peri-implant diseases and identified microbial signatures of disease severity.
| S-EPMC7603341 | biostudies-literature
Project description:Dysbiosis of airway microbiomes has been found in various respiratory diseases, but its molecular details in terms of taxonomic profile, metabolic characteristics, defensive function, and inhabit adaption are far from clear. Shotgun metagenome sequencing provides detailed information for microbes, whereas its application is rather limited in airways due to host DNA contaminants that overwhelm a minute amount of microbial content. Here, we describe a microfluidics-based enrichment device and an emulsion-based whole-genome amplification procedure (MEEA) for the preparation of DNA from sputa for shotgun sequencing in a metagenomics study. The two protocols coupled in MEEA are first separately assayed with mock samples and are both promising in efficiency and bias. The efficiency and consistency of MEEA are further evaluated in six clinical sputum samples against direct sequencing without enrichment, and MEEA enables 2 to 14 times enrichment for microbial reads, which take 14.68% to 33.52% of total reads. The dominant pathogens detected in MEEA are in excellent agreement with those from clinical etiological tests. Meanwhile, MEEA presents much more microbiome complexity and genome information at a strain level than direct sequencing, exhibiting high sensitivity for identifying prophages and DNA viruses. MEEA provides better microbiome profiling than direct sequencing without a preference for specific microorganisms. The more detailed functional and taxonomic characterization of their species constituents, including both bacterium and virus, facilitates metagenomics studies on the pathogenesis of respiratory microbiomes.IMPORTANCE The airway microbial community, which takes important pathogenic roles for respiratory diseases, is far from clear in terms of taxonomy and gene functions. One of the critical reasons is the heavy contamination of host cell/DNA in airway samples, which hinders the subsequent sequencing of the whole genomic contents of the microbial community-the metagenome. Here, we describe a protocol for airway sample preparation which couples a microbe enrichment microfluidic device and a DNA amplification method performed in numerous droplets. When evaluated with mock and clinical sputum samples, the microfluidics-based enrichment device and emulsion-based whole-genome amplification (MEEA) procedure efficiently removes host cells, amplifies the microbial genome, and shows no obvious bias among microbes. The efficiency of MEEA makes it a promising method in research of respiratory microbial communities and their roles in diseases.
| S-EPMC6589435 | biostudies-literature
Project description:MotivationAlthough many tools are available to study variation and its impact in single genomes, there is a lack of algorithms for finding such variation in metagenomes. This hampers the interpretation of metagenomics sequencing datasets, which are increasingly acquired in research on the (human) microbiome, in environmental studies and in the study of processes in the production of foods and beverages. Existing algorithms often depend on the use of reference genomes, which pose a problem when a metagenome of a priori unknown strain composition is studied. In this article, we develop a method to perform reference-free detection and visual exploration of genomic variation, both within a single metagenome and between metagenomes.ResultsWe present the MaryGold algorithm and its implementation, which efficiently detects bubble structures in contig graphs using graph decomposition. These bubbles represent variable genomic regions in closely related strains in metagenomic samples. The variation found is presented in a condensed Circos-based visualization, which allows for easy exploration and interpretation of the found variation. We validated the algorithm on two simulated datasets containing three respectively seven Escherichia coli genomes and showed that finding allelic variation in these genomes improves assemblies. Additionally, we applied MaryGold to publicly available real metagenomic datasets, enabling us to find within-sample genomic variation in the metagenomes of a kimchi fermentation process, the microbiome of a premature infant and in microbial communities living on acid mine drainage. Moreover, we used MaryGold for between-sample variation detection and exploration by comparing sequencing data sampled at different time points for both of these datasets.AvailabilityMaryGold has been written in C++ and Python and can be downloaded from http://bioinformatics.tudelft.nl/software
| S-EPMC3916741 | biostudies-literature
Project description:Metagenome assembly is a core yet methodologically challenging step for taxonomic classification and functional annotation of a microbiome. This study aims to generate the high-resolution human gut metagenome using both Illumina and Nanopore platforms. Assembly was achieved using four assemblers, including Flye (Nanopore), metaSPAdes (Illumina), hybridSPAdes (Illumina and Nanopore), and OPERA-MS (Illumina and Nanopore). Hybrid metagenome assembly was shown to generate contigs with almost same sizes comparable to those produced using Illumina reads alone, but was more contiguous, informative, and longer compared with those assembled with Illumina reads only. In addition, hybrid metagenome assembly enables us to obtain complete plasmid sequences and much more AMR gene-encoding contigs than the Illumina method. Most importantly, using our workflow, 58 novel high-quality metagenome bins were obtained from four assembly algorithms, particularly hybrid assembly (47/58), although metaSPAdes could provide 11 high-quality bins independently. Among them, 29 bins were currently uncultured bacterial metagenome-assembled genomes. These findings were highly consistent and supported by mock community data tested. In the analysis of biosynthetic gene clusters (BGCs), the number of BGCs in the contigs from hybridSPAdes (241) is higher than that of contigs from metaSPAdes (233). In conclusion, hybrid metagenome assembly could significantly enhance the efficiency of contig assembly, taxonomic binning, and genome construction compared with procedures using Illumina short-read data alone, indicating that nanopore long reads are highly useful in metagenomic applications. This technique could be used to create high-resolution references for future human metagenome studies.
| S-EPMC9134245 | biostudies-literature
Project description:Long-read sequencing and novel long-range assays have revolutionized de novo genome assembly by automating the reconstruction of reference-quality genomes. In particular, Hi-C sequencing is becoming an economical method for generating chromosome-scale scaffolds. Despite its increasing popularity, there are limited open-source tools available. Errors, particularly inversions and fusions across chromosomes, remain higher than alternate scaffolding technologies. We present a novel open-source Hi-C scaffolder that does not require an a priori estimate of chromosome number and minimizes errors by scaffolding with the assistance of an assembly graph. We demonstrate higher accuracy than the state-of-the-art methods across a variety of Hi-C library preparations and input assembly sizes. The Python and C++ code for our method is openly available at https://github.com/machinegun/SALSA.
| S-EPMC6719893 | biostudies-literature
Project description:Haplotype-resolved de novo assembly is the ultimate solution to the study of sequence variations in a genome. However, existing algorithms either collapse heterozygous alleles into one consensus copy or fail to cleanly separate the haplotypes to produce high-quality phased assemblies. Here we describe hifiasm, a de novo assembler that takes advantage of long high-fidelity sequence reads to faithfully represent the haplotype information in a phased assembly graph. Unlike other graph-based assemblers that only aim to maintain the contiguity of one haplotype, hifiasm strives to preserve the contiguity of all haplotypes. This feature enables the development of a graph trio binning algorithm that greatly advances over standard trio binning. On three human and five nonhuman datasets, including California redwood with a ~30-Gb hexaploid genome, we show that hifiasm frequently delivers better assemblies than existing tools and consistently outperforms others on haplotype-resolved assembly.
| S-EPMC7961889 | biostudies-literature
Project description:RNA-Seq is a powerful transcriptome profiling technology enabling transcript discovery and quantification. Whilst most commonly used for gene-level quantification, the data can be used for the analysis of transcript isoforms. However, when the underlying transcript assemblies are complex, current visualization approaches can be limiting, with splicing events a challenge to interpret. Here, we report on the development of a graph-based visualization method as a complementary approach to understanding transcript diversity from short-read RNA-Seq data. Following the mapping of reads to a reference genome, a read-to-read comparison is performed on all reads mapping to a given gene, producing a weighted similarity matrix between reads. This is used to produce an RNA assembly graph, where nodes represent reads and edges similarity scores between them. The resulting graphs are visualized in 3D space to better appreciate their sometimes large and complex topology, with other information being overlaid on to nodes, e.g. transcript models. Here we demonstrate the utility of this approach, including the unusual structure of these graphs and how they can be used to identify issues in assembly, repetitive sequences within transcripts and splice variants. We believe this approach has the potential to significantly improve our understanding of transcript complexity.
| S-EPMC6698738 | biostudies-literature
Project description:The "Graphical Fragment Assembly" (GFA) is an emerging format for the representation of sequence assembly graphs, which can be adopted by both de Bruijn graph- and string graph-based assemblers. Here we present RGFA, an implementation of the proposed GFA specification in Ruby. It allows the user to conveniently parse, edit and write GFA files. Complex operations such as the separation of the implicit instances of repeats and the merging of linear paths can be performed. A typical application of RGFA is the editing of a graph, to finish the assembly of a sequence, using information not available to the assembler. We illustrate a use case, in which the assembly of a repetitive metagenomic fosmid insert was completed using a script based on RGFA. Furthermore, we show how the API provided by RGFA can be employed to design complex graph editing algorithms. As an example, we developed a detection algorithm for CRISPRs in a de Bruijn graph. Finally, RGFA can be used for comparing assembly graphs, e.g., to document the changes in a graph after applying a GUI editor. A program, GFAdiff is provided, which compares the information in two graphs, and generate a report or a Ruby script documenting the transformation steps between the graphs.
| S-EPMC5103826 | biostudies-literature
Project description:MotivationThe assembly of contiguous sequence from metagenomic samples presents a particular challenge, due to the presence of multiple species, often closely related, at varying levels of abundance. Capturing diversity within species, for example, viral haplotypes, or bacterial strain-level diversity, is even more challenging.ResultsWe present MetaCortex, a metagenome assembler that captures intra-species diversity by searching for signatures of local variation along assembled sequences in the underlying assembly graph and outputting these sequences in sequence graph format. We show that MetaCortex produces accurate assemblies with higher genome coverage and contiguity than other popular metagenomic assemblers on mock viral communities with high levels of strain-level diversity and on simulated communities containing simulated strains.Availability and implementationSource code is freely available to download from https://github.com/SR-Martin/metacortex, is implemented in C and supported on MacOS and Linux. The version used for the results presented in this article is available at doi.org/10.5281/zenodo.7273627.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC9889960 | biostudies-literature
Project description:Synthetic long read sequencing techniques such as UST's TELL-Seq and Loop Genomics' LoopSeq combine 3[Formula: see text] barcoding with standard short-read sequencing to expand the range of linkage resolution from hundreds to tens of thousands of base-pairs. However, the lack of a 1:1 correspondence between a long fragment and a 3[Formula: see text] unique molecular identifier confounds the assignment of linkage between short reads. We introduce Ariadne, a novel assembly graph-based synthetic long read deconvolution algorithm, that can be used to extract single-species read-clouds from synthetic long read datasets to improve the taxonomic classification and de novo assembly of complex populations, such as metagenomes.
| S-EPMC10463629 | biostudies-literature