MetaQUBIC: a computational pipeline for gene-level functional profiling of metagenome and metatranscriptome.
Ontology highlight
ABSTRACT:
SUBMITTER: Ma A
PROVIDER: S-EPMC6821265 | biostudies-literature | 2019 Nov
REPOSITORIES: biostudies-literature
Publications
MetaQUBIC: a computational pipeline for gene-level functional profiling of metagenome and metatranscriptome.
Ma Anjun A Sun Minxuan M McDermaid Adam A Liu Bingqiang B Ma Qin Q
Bioinformatics (Oxford, England) 20191101 21
<h4>Motivation</h4>Metagenomic and metatranscriptomic analyses can provide an abundance of information related to microbial communities. However, straightforward analysis of this data does not provide optimal results, with a required integration of data types being needed to thoroughly investigate these microbiomes and their environmental interactions.<h4>Results</h4>Here, we present MetaQUBIC, an integrated biclustering-based computational pipeline for gene module detection that integrates both ...[more]
Similar Datasets
Project description:BackgroundMetatranscriptomic landscapes can provide insights in functional relationships within natural microbial communities. Analysis of complex metatranscriptome datasets of these communities poses a considerable bioinformatic challenge since they are non-restricted with a varying number of participating strains and species. For RNA-Seq data a standard approach is to align the generated reads to a set of closely related reference genomes. This only works well for microbial communities for which a near complete catalogue of reference genomes is available at a small evolutionary distance. In this study, we focus on the design of a validated de novo metatranscriptome assembly pipeline for single-end Illumina RNA-Seq data to obtain functional and taxonomic profiles of murine microbial communities.ResultsThe here developed de novo assembly metatranscriptome pipeline combined rRNA removal, IDBA-UD assembler, functional annotation and taxonomic classification. Different assemblers were tested and validated using RNA-Seq data from an in silico generated mock community and in vivo RNA-Seq data from a restricted microbial community taken from a mouse model colonized with Altered Schaedler Flora (ASF). Precision and recall of resulting gene expression, functional and taxonomic profiles were compared to those obtained with a standard alignment method. The validated pipeline was subsequently used to generate expression profiles from non-restricted cecal communities of four C57BL/6J mice fed on a high-fat high-protein diet spiked with an RNA-Seq data set from a well-characterized human sample. The spike in control was used to estimate precision and recall at assembly, functional and taxonomic level of non-restricted communities.ConclusionsA generic de novo assembly pipeline for metatranscriptome data analysis was designed for microbial ecosystems, which can be applied for microbial metatranscriptome analysis in any chosen niche.
| S-EPMC4710500 | biostudies-literature
Project description:Many computational pipelines exist for the detection of differentially expressed genes. However, computational pipelines for functional gene detection rarely exist. We developed a new computational pipeline for functional gene identification from transcriptome profiling data. Key features of the pipeline include batch effect correction, clustering optimization by gap statistics, gene ontology analysis of clustered genes, and literature analysis for functional gene discovery. By leveraging this pipeline on RNA-seq datasets from two mouse retinal development studies, we identified 7 candidate genes involved in the formation of the photoreceptor outer segment. The expression of top three candidate genes (Pde8b, Laptm4b, and Nr1h4) in the outer segment of the developing mouse retina were experimentally validated by immunohistochemical analysis. This computational pipeline can accurately predict novel functional gene for a specific biological process, e.g., development of the outer segment and synapses of the photoreceptor cells in the mouse retina. This pipeline can also be useful to discover functional genes for other biological processes and in other organs and tissues.
| S-EPMC8651667 | biostudies-literature
Project description:Chronic obstructive pulmonary disease (COPD) is an inflammatory lung disorder characterized by the progressive obstruction of airflow and is currently the fourth leading cause of death in the world. The pathogenesis of COPD is thought to involve bacterial infections and inflammations. Owing to advancement in sequencing technology, evidence is emerging that supports an association between the lung microbiome and COPD. However, few studies have looked into the expression profile of the bacterial communities in the COPD lungs. In this study, we analyzed the sputum microbiome of four moderate and four severe COPD male patients both at the DNA and RNA level, using next generation sequencing technology. We found that bacterial composition determined by 16S rRNA gene sequencing may not directly translate to the set of actively expressing bacteria as defined by transcriptome sequencing. The two sequencing data agreed on Prevotella, Rothia, Neisseria, Porphyromonas, Veillonella, Fusobacterium and Streptococcus being among the most differentially abundant genera between the moderate and severe COPD samples, supporting their association with COPD severity. However, the two sequencing analyses disagreed on the relative abundance of these bacteria in the two COPD groups, implicating the importance of studying the actively expressing bacteria for enriching our understanding of COPD. Though we have described the metatranscriptome profiles of the lung microbiome in moderate and severe COPD, further investigations are required to determine the functional basis underlying the relationship between the microbial species in the lungs and pathogenesis of COPD.
| S-EPMC4948834 | biostudies-literature
Project description:BackgroundAlthough metatranscriptomics-the study of diverse microbial population activity based on RNA-seq data-is rapidly growing in popularity, there are limited options for biologists to analyze this type of data. Current approaches for processing metatranscriptomes rely on restricted databases and a dedicated computing cluster, or metagenome-based approaches that have not been fully evaluated for processing metatranscriptomic datasets. We created a new bioinformatics pipeline, designed specifically for metatranscriptome dataset analysis, which runs in conjunction with Metagenome-RAST (MG-RAST) servers. Designed for use by researchers with relatively little bioinformatics experience, SAMSA offers a breakdown of metatranscriptome transcription activity levels by organism or transcript function, and is fully open source. We used this new tool to evaluate best practices for sequencing stool metatranscriptomes.ResultsWorking with the MG-RAST annotation server, we constructed the Simple Annotation of Metatranscriptomes by Sequence Analysis (SAMSA) software package, a complete pipeline for the analysis of gut microbiome data. SAMSA can summarize and evaluate raw annotation results, identifying abundant species and significant functional differences between metatranscriptomes. Using pilot data and simulated subsets, we determined experimental requirements for fecal gut metatranscriptomes. Sequences need to be either long reads (longer than 100 bp) or joined paired-end reads. Each sample needs 40-50 million raw sequences, which can be expected to yield the 5-10 million annotated reads necessary for accurate abundance measures. We also demonstrated that ribosomal RNA depletion does not equally deplete ribosomes from all species within a sample, and remaining rRNA sequences should be discarded. Using publicly available metatranscriptome data in which rRNA was not depleted, we were able to demonstrate that overall organism transcriptional activity can be measured using mRNA counts. We were also able to detect significant differences between control and experimental groups in both organism transcriptional activity and specific cellular functions.ConclusionsBy making this new pipeline publicly available, we have created a powerful new tool for metatranscriptomics research, offering a new method for greater insight into the activity of diverse microbial communities. We further recommend that stool metatranscriptomes be ribodepleted and sequenced in a 100 bp paired end format with a minimum of 40 million reads per sample.
| S-EPMC5041328 | biostudies-literature
Project description:Analyses of metagenome data (MG) and metatranscriptome data (MT) are often challenged by a paucity of complete reference genome sequences and the uneven/low sequencing depth of the constituent organisms in the microbial community, which respectively limit the power of reference-based alignment and de novo sequence assembly. These limitations make accurate protein family classification and abundance estimation challenging, which in turn hamper downstream analyses such as abundance profiling of metabolic pathways, identification of differentially encoded/expressed genes, and de novo reconstruction of complete gene and protein sequences from the protein family of interest. The profile hidden Markov model (HMM) framework enables the construction of very useful probabilistic models for protein families that allow for accurate modeling of position specific matches, insertions, and deletions. We present a novel homology detection algorithm that integrates banded Viterbi algorithm for profile HMM parsing with an iterative simultaneous alignment and assembly computational framework. The algorithm searches a given profile HMM of a protein family against a database of fragmentary MG/MT sequencing data and simultaneously assembles complete or near-complete gene and protein sequences of the protein family. The resulting program, HMM-GRASPx, demonstrates superior performance in aligning and assembling homologs when benchmarked on both simulated marine MG and real human saliva MG datasets. On real supragingival plaque and stool MG datasets that were generated from healthy individuals, HMM-GRASPx accurately estimates the abundances of the antimicrobial resistance (AMR) gene families and enables accurate characterization of the resistome profiles of these microbial communities. For real human oral microbiome MT datasets, using the HMM-GRASPx estimated transcript abundances significantly improves detection of differentially expressed (DE) genes. Finally, HMM-GRASPx was used to reconstruct comprehensive sets of complete or near-complete protein and nucleotide sequences for the query protein families. HMM-GRASPx is freely available online from http://sourceforge.net/projects/hmm-graspx.
| S-EPMC4939949 | biostudies-literature
Project description:Although the composition of the human microbiome is now well-studied, the microbiota's >8 million genes and their regulation remain largely uncharacterized. This knowledge gap is in part because of the difficulty of acquiring large numbers of samples amenable to functional studies of the microbiota. We conducted what is, to our knowledge, one of the first human microbiome studies in a well-phenotyped prospective cohort incorporating taxonomic, metagenomic, and metatranscriptomic profiling at multiple body sites using self-collected samples. Stool and saliva were provided by eight healthy subjects, with the former preserved by three different methods (freezing, ethanol, and RNAlater) to validate self-collection. Within-subject microbial species, gene, and transcript abundances were highly concordant across sampling methods, with only a small fraction of transcripts (<5%) displaying between-method variation. Next, we investigated relationships between the oral and gut microbial communities, identifying a subset of abundant oral microbes that routinely survive transit to the gut, but with minimal transcriptional activity there. Finally, systematic comparison of the gut metagenome and metatranscriptome revealed that a substantial fraction (41%) of microbial transcripts were not differentially regulated relative to their genomic abundances. Of the remainder, consistently underexpressed pathways included sporulation and amino acid biosynthesis, whereas up-regulated pathways included ribosome biogenesis and methanogenesis. Across subjects, metatranscriptional profiles were significantly more individualized than DNA-level functional profiles, but less variable than microbial composition, indicative of subject-specific whole-community regulation. The results thus detail relationships between community genomic potential and gene expression in the gut, and establish the feasibility of metatranscriptomic investigations in subject-collected and shipped samples.
| S-EPMC4050606 | biostudies-literature
Project description:BACKGROUND:Recent evidence suggests that immunotherapy efficacy in melanoma is modulated by gut microbiota. Few studies have examined this phenomenon in humans, and none have incorporated metatranscriptomics, important for determining expression of metagenomic functions in the microbial community. METHODS:In melanoma patients undergoing immunotherapy, gut microbiome was characterized in pre-treatment stool using 16S rRNA gene and shotgun metagenome sequencing (n = 27). Transcriptional expression of metagenomic pathways was confirmed with metatranscriptome sequencing in a subset of 17. We examined associations of taxa and metagenomic pathways with progression-free survival (PFS) using 500 × 10-fold cross-validated elastic-net penalized Cox regression. RESULTS:Higher microbial community richness was associated with longer PFS in 16S and shotgun data (p < 0.05). Clustering based on overall microbiome composition divided patients into three groups with differing PFS; the low-risk group had 99% lower risk of progression than the high-risk group at any time during follow-up (p = 0.002). Among the species selected in regression, abundance of Bacteroides ovatus, Bacteroides dorei, Bacteroides massiliensis, Ruminococcus gnavus, and Blautia producta were related to shorter PFS, and Faecalibacterium prausnitzii, Coprococcus eutactus, Prevotella stercorea, Streptococcus sanguinis, Streptococcus anginosus, and Lachnospiraceae bacterium 3 1 46FAA to longer PFS. Metagenomic functions related to PFS that had correlated metatranscriptomic expression included risk-associated pathways of L-rhamnose degradation, guanosine nucleotide biosynthesis, and B vitamin biosynthesis. CONCLUSIONS:This work adds to the growing evidence that gut microbiota are related to immunotherapy outcomes, and identifies, for the first time, transcriptionally expressed metagenomic pathways related to PFS. Further research is warranted on microbial therapeutic targets to improve immunotherapy outcomes.
| S-EPMC6785875 | biostudies-other
Project description:The production of biogas by anaerobic digestion (AD) of agricultural residues, organic wastes, animal excrements, municipal sludge, and energy crops has a firm place in sustainable energy production and bio-economy strategies. Focusing on the microbial community involved in biomass conversion offers the opportunity to control and engineer the biogas process with the objective to optimize its efficiency. Taxonomic profiling of biogas producing communities by means of high-throughput 16S rRNA gene amplicon sequencing provided high-resolution insights into bacterial and archaeal structures of AD assemblages and their linkages to fed substrates and process parameters. Commonly, the bacterial phyla Firmicutes and Bacteroidetes appeared to dominate biogas communities in varying abundances depending on the apparent process conditions. Regarding the community of methanogenic Archaea, their diversity was mainly affected by the nature and composition of the substrates, availability of nutrients and ammonium/ammonia contents, but not by the temperature. It also appeared that a high proportion of 16S rRNA sequences can only be classified on higher taxonomic ranks indicating that many community members and their participation in AD within functional networks are still unknown. Although cultivation-based approaches to isolate microorganisms from biogas fermentation samples yielded hundreds of novel species and strains, this approach intrinsically is limited to the cultivable fraction of the community. To obtain genome sequence information of non-cultivable biogas community members, metagenome sequencing including assembly and binning strategies was highly valuable. Corresponding research has led to the compilation of hundreds of metagenome-assembled genomes (MAGs) frequently representing novel taxa whose metabolism and lifestyle could be reconstructed based on nucleotide sequence information. In contrast to metagenome analyses revealing the genetic potential of microbial communities, metatranscriptome sequencing provided insights into the metabolically active community. Taking advantage of genome sequence information, transcriptional activities were evaluated considering the microorganism's genetic background. Metaproteome studies uncovered enzyme profiles expressed by biogas community members. Enzymes involved in cellulose and hemicellulose decomposition and utilization of other complex biopolymers were identified. Future studies on biogas functional microbial networks will increasingly involve integrated multi-omics analyses evaluating metagenome, transcriptome, proteome, and metabolome datasets.
| S-EPMC5959977 | biostudies-literature
Project description:Metagenomics research has accelerated the studies of microbial organisms, providing insights into the composition and potential functionality of various microbial communities. Metatranscriptomics (studies of the transcripts from a mixture of microbial species) and other meta-omics approaches hold even greater promise for providing additional insights into functional and regulatory characteristics of the microbial communities. Current metatranscriptomics projects are often carried out without matched metagenomic datasets (of the same microbial communities). For the projects that produce both metatranscriptomic and metagenomic datasets, their analyses are often not integrated. Metagenome assemblies are far from perfect, partially explaining why metagenome assemblies are not used for the analysis of metatranscriptomic datasets.Here, we report a reads mapping algorithm for mapping of short reads onto a de Bruijn graph of assemblies. A hash table of junction k-mers (k-mers spanning branching structures in the de Bruijn graph) is used to facilitate fast mapping of reads to the graph. We developed an application of this mapping algorithm: a reference-based approach to metatranscriptome assembly using graphs of metagenome assembly as the reference. Our results show that this new approach (called TAG) helps to assemble substantially more transcripts that otherwise would have been missed or truncated because of the fragmented nature of the reference metagenome.TAG was implemented in C++ and has been tested extensively on the Linux platform. It is available for download as open source at http://omics.informatics.indiana.edu/TAG CONTACT: yye@indiana.edu.
| S-EPMC4896364 | biostudies-literature
Project description:Whole metagenome shotgun sequencing is a powerful approach for assaying the functional potential of microbial communities. We currently lack tools that efficiently and accurately align DNA reads against protein references, the technique necessary for constructing a functional profile. Here, we present PALADIN-a novel modification of the Burrows-Wheeler Aligner that provides accurate alignment, robust reporting capabilities and orders-of-magnitude improved efficiency by directly mapping in protein space.We compared the accuracy and efficiency of PALADIN against existing tools that employ nucleotide or protein alignment algorithms. Using simulated reads, PALADIN consistently outperformed the popular DNA read mappers BWA and NovoAlign in detected proteins, percentage of reads mapped and ontological similarity. We also compared PALADIN against four existing protein alignment tools: BLASTX, RAPSearch2, DIAMOND and Lambda, using empirically obtained reads. PALADIN yielded results seven times faster than the best performing alternative, DIAMOND and nearly 8000 times faster than BLASTX. PALADIN's accuracy was comparable to all tested solutions.PALADIN was implemented in C, and its source code and documentation are available at https://github.com/twestbrookunh/paladin.anthonyw@wildcats.unh.edu.Supplementary data are available at Bioinformatics online.
| S-EPMC5423455 | biostudies-literature