Dataset Information

PhyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes.

ABSTRACT: The small-subunit rRNA (SSU rRNA) gene is the key marker in molecular ecology for all domains of life, but it is largely absent from metagenome-assembled genomes that often are the only resource available for environmental microbes. Here, we present phyloFlash, a pipeline to overcome this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based binning of full metagenomic assemblies. We show that a cleanup of artifacts is pivotal even with a curated reference database. With such a filtered database, the general-purpose mapper BBmap extracts SSU rRNA reads five times faster than the rRNA-specialized tool SortMeRNA with similar sensitivity and higher selectivity on simulated metagenomes. Reference-based targeted assemblers yielded either highly fragmented assemblies or high levels of chimerism, so we employ the general-purpose genomic assembler SPAdes. Our optimized implementation is independent of reference database composition and has satisfactory levels of chimera formation. phyloFlash quickly processes Illumina (meta)genomic data, is straightforward to use, even as part of high-throughput quality control, and has user-friendly output reports. The software is available at https://github.com/HRGV/phyloFlash (GPL3 license) and is documented with an online manual.IMPORTANCE To track organisms across all domains of life, the SSU rRNA gene is the gold standard. Many environmental microbes are known only from high-throughput sequence data, but the SSU rRNA gene, the key to visualization by molecular probes and link to existing literature, is often missing from metagenome-assembled genomes (MAGs). The easy-to-use phyloFlash software suite tackles this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based linking to MAGs. Starting from a cleaned reference database, phyloFlash profiles the taxonomic diversity and assembles the sorted SSU rRNA reads. The phyloFlash design is domain agnostic and covers eukaryotes, archaea, and bacteria alike. phyloFlash also provides utilities to visualize multisample comparisons and to integrate the recovered SSU rRNAs in a metagenomics workflow by linking them to MAGs using assembly graph parsing.

SUBMITTER: Gruber-Vodicka HR

PROVIDER: S-EPMC7593591 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

phyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes.

Gruber-Vodicka Harald R HR Seah Brandon K B BKB Pruesse Elmar E

mSystems 20201027 5

The small-subunit rRNA (SSU rRNA) gene is the key marker in molecular ecology for all domains of life, but it is largely absent from metagenome-assembled genomes that often are the only resource available for environmental microbes. Here, we present phyloFlash, a pipeline to overcome this gap with rapid, SSU rRNA-centered taxonomic classification, targeted assembly, and graph-based binning of full metagenomic assemblies. We show that a cleanup of artifacts is pivotal even with a curated referenc ...[more]

PMID: 33109753

Similar Datasets

Project description:Small subunit rRNA (SSU rRNA) amplicon sequencing can quantitatively and comprehensively profile natural microbiomes, representing a critically important tool for studying diverse global ecosystems. However, results will only be accurate if PCR primers perfectly match the rRNA of all organisms present. To evaluate how well marine microorganisms across all 3 domains are detected by this method, we compared commonly used primers with >300 million rRNA gene sequences retrieved from globally distributed marine metagenomes. The best-performing primers compared to 16S rRNA of bacteria and archaea were 515Y/926R and 515Y/806RB, which perfectly matched over 96% of all sequences. Considering cyanobacterial and chloroplast 16S rRNA, 515Y/926R had the highest coverage (99%), making this set ideal for quantifying marine primary producers. For eukaryotic 18S rRNA sequences, 515Y/926R also performed best (88%), followed by V4R/V4RB (18S rRNA specific; 82%)-demonstrating that the 515Y/926R combination performs best overall for all 3 domains. Using Atlantic and Pacific Ocean samples, we demonstrate high correspondence between 515Y/926R amplicon abundances (generated for this study) and metagenomic 16S rRNA (median R2 = 0.98, n = 272), indicating amplicons can produce equally accurate community composition data compared with shotgun metagenomics. Our analysis also revealed that expected performance of all primer sets could be improved with minor modifications, pointing toward a nearly completely universal primer set that could accurately quantify biogeochemically important taxa in ecosystems ranging from the deep sea to the surface. In addition, our reproducible bioinformatic workflow can guide microbiome researchers studying different ecosystems or human health to similarly improve existing primers and generate more accurate quantitative amplicon data. IMPORTANCE PCR amplification and sequencing of marker genes is a low-cost technique for monitoring prokaryotic and eukaryotic microbial communities across space and time but will work optimally only if environmental organisms match PCR primer sequences exactly. In this study, we evaluated how well primers match globally distributed short-read oceanic metagenomes. Our results demonstrate that primer sets vary widely in performance, and that at least for marine systems, rRNA amplicon data from some primers lack significant biases compared to metagenomes. We also show that it is theoretically possible to create a nearly universal primer set for diverse saline environments by defining a specific mixture of a few dozen oligonucleotides, and present a software pipeline that can guide rational design of primers for any environment with available meta'omic data.

Project description:Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth.

Project description:The structure of rhizobial communities nodulating Acacia in southeastern Australia from south Queensland to Tasmania was investigated by a molecular approach. A total of 118 isolates from nodule samples from 13 different Acacia species collected at 44 sites were characterized by small-subunit (SSU) ribosomal DNA (rDNA) PCR-restriction fragment length polymorphism analysis. Nine rhizobial genomospecies were identified, and these taxa corresponded to previously described genomospecies (B. Lafay and J. J. Burdon, Appl. Environ. Microbiol. 64:3989-3997, 1998). Eight of these genomospecies belonged to the Bradyrhizobium lineage and accounted for 96.6% of the isolates. The remaining genomospecies corresponded to Rhizobium tropici. For analysis of geographic patterns, results were grouped into five latitudinal regions regardless of host origin. In each region, as observed previously for rhizobial isolates taken from non-Acacia legumes (Lafay and Burdon, Appl. Environ. Microbiol. 64:3989-3997, 1998), rhizobial communities were dominated by one or two genomospecies, the identities of which varied from place to place. Despite this similarity in patterns, the most abundant genomospecies for Acacia isolates differed from the genomospecies found in the non-Acacia-derived rhizobial collection, suggesting that there is a difference in nodulation patterns of the Mimosoideae and the Papilionoideae. Only two genomospecies were both widespread and relatively abundant across the range of sites sampled. Genomospecies A was found in all regions except the most northern sites located in Queensland, whereas genomospecies B was not detected in Tasmania. This suggests that genomospecies A might be restricted to the more temperate regions of Australia, whereas in contrast, genomospecies B occurs in different climatic and edaphic conditions across the whole continent. The latter hypothesis is supported by the presence of genomospecies B in southwestern Australia, based on partial SSU rDNA sequence data (N. D. S. Marsudi, A. R. Glenn, and M. J. Dilworth, Soil Biol. Biochem. 31:1229-1238, 1998).

Dataset Information

PhyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes.

Publications

phyloFlash: Rapid Small-Subunit rRNA Profiling and Targeted Assembly from Metagenomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets