Meta-Apo improves accuracy of 16S-amplicon-based prediction of microbiome function

On the limits of 16S rRNA gene-based metagenome prediction and functional profiling

Molecular profiling techniques such as metagenomics, metatranscriptomics or metabolomics offer important insights into the functional diversity of the microbiome. In contrast, 16S rRNA gene sequencing, a widespread and cost-effective technique to measure microbial diversity, only allows for indirect estimation of microbial function. To mitigate this, tools such as PICRUSt2, Tax4Fun2, PanFP and MetGEM infer functional profiles from 16S rRNA gene sequencing data using different algorithms. Prior studies have cast doubts on the quality of these predictions, motivating us to systematically evaluate these tools using matched 16S rRNA gene sequencing, metagenomic datasets, and simulated data. Our contribution is threefold: (i) using simulated data, we investigate if technical biases could explain the discordance between inferred and expected results; (ii) considering human cohorts for type two diabetes, colorectal cancer and obesity, we test if health-related differential abundance measures of functional categories are concordant between 16S rRNA gene-inferred and metagenome-derived profiles and; (iii) since 16S rRNA gene copy number is an important confounder in functional profiles inference, we investigate if a customised copy number normalisation with the rrnDB database could improve the results. Our results show that 16S rRNA gene-based functional inference tools generally do not have the necessary sensitivity to delineate health-related functional changes in the microbiome and should thus be used with care. Furthermore, we outline important differences in the individual tools tested and offer recommendations for tool selection.

Community-scale models of microbiomes: Articulating metabolic modelling and metagenome sequencing

Building models is essential for understanding the functions and dynamics of microbial communities. Metabolic models built on genome-scale metabolic network reconstructions (GENREs) are especially relevant as a means to decipher the complex interactions occurring among species. Model reconstruction increasingly relies on metagenomics, which permits direct characterisation of naturally occurring communities that may contain organisms that cannot be isolated or cultured. In this review, we provide an overview of the field of metabolic modelling and its increasing reliance on and synergy with metagenomics and bioinformatics. We survey the means of assigning functions and reconstructing metabolic networks from (meta-)genomes, and present the variety and mathematical fundamentals of metabolic models that foster the understanding of microbial dynamics. We emphasise the characterisation of interactions and the scaling of model construction to large communities, two important bottlenecks in the applicability of these models. We give an overview of the current state of the art in metagenome sequencing and bioinformatics analysis, focusing on the reconstruction of genomes in microbial communities. Metagenomics benefits tremendously from third-generation sequencing, and we discuss the opportunities of long-read sequencing, strain-level characterisation and eukaryotic metagenomics. We aim at providing algorithmic and mathematical support, together with tool and application resources, that permit bridging the gap between metagenomics and metabolic modelling.

Microbial denitrification characteristics of typical decentralized wastewater treatment processes based on 16S rRNA sequencing

Despite the widespread application of decentralized wastewater treatment (WWT) facilities in China, relatively few research has used the multi-media biological filter (MMBF) facilities to investigate the microorganism characteristics. This study utilizes 16S rRNA high-throughput sequencing (HTS) technology to examine the microbial biodiversity of a representative wastewater treatment (WWT) system in an expressway service area. The pathways of nitrogen removal along the treatment route were analyzed in conjunction with water quality monitoring. The distribution and composition of microbial flora in the samples were examined, and the dominant flora were identified using LEfSe analysis. The FAPROTAX methodology was employed to investigate the relative abundance of genes associated with the nitrogen cycle and to discern the presence of functional genes involved in nitrogen metabolism. On average, the method has a high level of efficiency in removing COD, TN, NH3-N, and TP from the effluent. The analysis of the microbial community identified a total of 40 phyla, 111 classes, 143 orders, 263 families, and 419 genera. The phyla that were predominantly observed include Proteobacteria, Acidobacteria, Chloroflexi, Actinobacteria, Nitrospirae, Bacteroidetes. The results show that the system has achieved high performance in nitrogen removal, the abundance of nitrification genes is significantly higher than that of other nitrogen cycle genes such as denitrification, and there are six nitrogen metabolism pathways, primarily nitrification, among which Nitrospirae and Nitrospira are the core differentiated flora that can adapt to low temperature conditions and participate in nitrification, and are the dominant nitrogen removal flora in cold regions. This work aims to comprehensively investigate the diversity and functional properties of the bacterial community in decentralized WWT processes.

Comprehensive Assessment of 16S rRNA Gene Amplicon Sequencing for Microbiome Profiling across Multiple Habitats

The 16S rRNA gene works as a rapid and effective marker for the identification of microorganisms in complex communities; hence, a huge number of microbiomes have been surveyed by 16S amplicon-based sequencing. The resolution of the 16S rRNA gene is always considered only at the genus level; however, it has not been verified on a wide range of microbes yet. To fully explore the ability and potential of the 16S rRNA gene in microbial profiling, here, we propose Qscore, a comprehensive method to evaluate the performance of amplicons by integrating the amplification rate, multitier taxonomic annotation, sequence type, and length. Our in silico assessment by a "global view" of 35,889 microbe species across multiple reference databases summarizes the optimal sequencing strategy for 16S short reads. On the other hand, since microbes are unevenly distributed according to their habitats, we also provide the recommended configuration for 16 typical ecosystems based on the Qscores of 157,390 microbiomes in the Microbiome Search Engine (MSE). Detailed data simulation further proves that the 16S amplicons produced with Qscore-suggested parameters exhibit high precision in microbiome profiling, which is close to that of shotgun metagenomes under CAMI metrics. Therefore, by reconsidering the precision of 16S-based microbiome profiling, our work not only enables the high-quality reusability of massive sequence legacy that has already been produced but is also significant for guiding microbiome studies in the future. We have implemented the Qscore as an online service at http://qscore.single-cell.cn to parse the recommended sequencing strategy for specific habitats or expected microbial structures. IMPORTANCE 16S rRNA has long been used as a biomarker to identify distinct microbes from complex communities. However, due to the influence of the amplification region, sequencing type, sequence processing, and reference database, the accuracy of 16S rRNA has not been fully verified on a global range. More importantly, the microbial composition of different habitats varies greatly, and it is necessary to adopt different strategies according to the corresponding target microbes to achieve optimal analytical performance. Here, we developed Qscore, which evaluates the comprehensive performance of 16S amplicons from multiple perspectives, thus providing the best sequencing strategies for common ecological environments by using big data.

Establishment and Validation of a New Analysis Strategy for the Study of Plant Endophytic Microorganisms

Amplicon sequencing of bacterial or fungal marker sequences is currently the main method for the study of endophytic microorganisms in plants. However, it cannot obtain all types of microorganisms, including bacteria, fungi, protozoa, etc., in samples, nor compare the relative content between endophytic microorganisms and plants and between different types of endophytes. Therefore, it is necessary to develop a better analysis strategy for endophytic microorganism investigation. In this study, a new analysis strategy was developed to obtain endophytic microbiome information from plant transcriptome data. Results showed that the new strategy can obtain the composition of microbial communities and the relative content between plants and endophytic microorganisms, and between different types of endophytic microorganisms from the plant transcriptome data. Compared with the amplicon sequencing method, more endophytic microorganisms and relative content information can be obtained with the new strategy, which can greatly broaden the research scope and save the experimental cost. Furthermore, the advantages and effectiveness of the new strategy were verified with different analysis of the microbial composition, correlation analysis, inoculant content test, and repeatability test.

Beyond Basic Diversity Estimates-Analytical Tools for Mechanistic Interpretations of Amplicon Sequencing Data

Understanding microbial ecology through amplifying short read regions, typically 16S rRNA for prokaryotic species or 18S rRNA for eukaryotic species, remains a popular, economical choice. These methods provide relative abundances of key microbial taxa, which, depending on the experimental design, can be used to infer mechanistic ecological underpinnings. In this review, we discuss recent advancements in in situ analytical tools that have the power to elucidate ecological phenomena, unveil the metabolic potential of microbial communities, identify complex multidimensional interactions between species, and compare stability and complexity under different conditions. Additionally, we highlight methods that incorporate various modalities and additional information, which in combination with abundance data, can help us understand how microbial communities respond to change in a typical ecosystem. Whilst the field of microbial informatics continues to progress substantially, our emphasis is on popular methods that are applicable to a broad range of study designs. The application of these methods can increase our mechanistic understanding of the ongoing dynamics of complex microbial communities.

The landscape in the gut microbiome of long-lived families reveals new insights on longevity and aging - relevant neural and immune function

Human longevity has a strong familial and genetic component. Dynamic characteristics of the gut microbiome during aging associated with longevity, neural, and immune function remained unknown. Here, we aim to reveal the synergistic changes in gut microbiome associated with decline in neural and immune system with aging and further obtain insights into the establishment of microbiome homeostasis that can benefit human longevity. Based on 16S rRNA and metagenomics sequencing data for 32 longevity families including three generations, centenarians, elderly, and young groups, we found centenarians showed increased diversity of gut microbiota, severely damaged connection among bacteria, depleted in microbial-associated essential amino acid function, and increased abundance of anti-inflammatory bacteria in comparison to young and elderly groups. Some potential probiotic species, such as Desulfovibrio piger, Gordonibacter pamelaeae, Odoribacter splanchnicus, and Ruminococcaceae bacterium D5 were enriched with aging, which might possibly support health maintenance. The level of Amyloid-β (Aβ) and brain-derived neurotrophic factor (BDNF) related to neural function showed increased and decreased with aging, respectively. The elevated level of inflammatory factors was observed in centenarians compared with young and elderly groups. The enriched Bacteroides fragilis in centenarians might promote longevity through up-regulating anti-inflammatory factor IL-10 expression to mediate the critical balance between health and disease. Impressively, the associated analysis for gut microbiota with the level of Aβ, BDNF, and inflammatory factors suggests Bifidobacterium pseudocatenulatum could be a particularly beneficial bacteria in the improvement of impaired neural and immune function. Our results provide a rationale for targeting the gut microbiome in future clinical applications of aging-related diseases and extending life span.Abbreviations: 16S rRNA: 16S ribosomal RNA; MAGs: Metagenome-assembled genomes; ASVs: Amplicon sequence variants; DNA: Deoxyribonucleic acid; FDR: False discovery rate: KEGG: Kyoto Encyclopedia of Genes and Genomes; PCoA: Principal coordinates analysis; PCR: Polymerase chain reaction; PICRUSt: Phylogenetic Investigation of Communities by Reconstruction of Unobserved States; Aβ: Amyloid-β (Aβ); BDNF: Brain-derived neurotrophic factor.

Parallel-Meta Suite: Interactive and rapid microbiome data analysis on multiple platforms

Massive microbiome sequencing data has been generated, which elucidates associations between microbes and their environmental phenotypes such as host health or ecosystem status. Outstanding bioinformatic tools are the basis to decipher the biological information hidden under microbiome data. However, most approaches placed difficulties on the accessibility to nonprofessional users. On the other side, the computing throughput has become a significant bottleneck of many analytical pipelines in processing large-scale datasets. In this study, we introduce Parallel-Meta Suite (PMS), an interactive software package for fast and comprehensive microbiome data analysis, visualization, and interpretation. It covers a wide array of functions for data preprocessing, statistics, visualization by state-of-the-art algorithms in a user-friendly graphical interface, which is accessible to diverse users. To meet the rapidly increasing computational demands, the entire procedure of PMS has been optimized by a parallel computing scheme, enabling the rapid processing of thousands of samples. PMS is compatible with multiple platforms, and an installer has been integrated for full-automatic installation.

Structure and Function of Oral Microbial Community in Periodontitis Based on Integrated Data

Objective

Microorganisms play a key role in the initiation and progression of periodontal disease. Research studies have focused on seeking specific microorganisms for diagnosing and monitoring the outcome of periodontitis treatment. Large samples may help to discover novel potential biomarkers and capture the common characteristics among different periodontitis patients. This study examines how to screen and merge high-quality periodontitis-related sequence datasets from several similar projects to analyze and mine the potential information comprehensively.

Methods

In all, 943 subgingival samples from nine publications were included based on predetermined screening criteria. A uniform pipeline (QIIME2) was applied to clean the raw sequence datasets and merge them together. Microbial structure, biomarkers, and correlation network were explored between periodontitis and healthy individuals. The microbiota patterns at different periodontal pocket depths were described. Additionally, potential microbial functions and metabolic pathways were predicted using PICRUSt to assess the differences between health and periodontitis.

Results

The subgingival microbial communities and functions in subjects with periodontitis were significantly different from those in healthy subjects. Treponema, TG5, Desulfobulbus, Catonella, Bacteroides, Aggregatibacter, Peptostreptococcus, and Eikenella were periodontitis biomarkers, while Veillonella, Corynebacterium, Neisseria, Rothia, Paludibacter, Capnocytophaga, and Kingella were signature of healthy periodontium. With the variation of pocket depth from shallow to deep pocket, the proportion of Spirochaetes, Bacteroidetes, TM7, and Fusobacteria increased, whereas that of Proteobacteria and Actinobacteria decreased. Synergistic relationships were observed among different pathobionts and negative relationships were noted between periodontal pathobionts and healthy microbiota.

Conclusion

This study shows significant differences in the oral microbial community and potential metabolic pathways between the periodontitis and healthy groups. Our integrated analysis provides potential biomarkers and directions for in-depth research. Moreover, a new method for integrating similar sequence data is shown here that can be applied to other microbial-related areas.

Hierarchical Meta-Storms enables comprehensive and rapid comparison of microbiome functional profiles on a large scale using hierarchical dissimilarity metrics and parallel computing

Functional beta-diversity analysis on numerous microbiomes interprets the linkages between metabolic functions and their meta-data. To evaluate the microbiome beta-diversity, widely used distance metrices only count overlapped gene families but omit their inherent relationships, resulting in erroneous distances due to the sparsity of high-dimensional function profiles. Here we propose Hierarchical Meta-Storms (HMS) to tackle such problem. HMS contains two core components: (i) a dissimilarity algorithm that comprehensively measures functional distances among microbiomes using multi-level metabolic hierarchy and (ii) a fast Principal Co-ordinates Analysis (PCoA) implementation that deduces the beta-diversity pattern optimized by parallel computing. Results showed HMS can detect the variations of microbial functions in upper-level metabolic pathways, however, always missed by other methods. In addition, HMS accomplished the pairwise distance matrix and PCoA for 20 000 microbiomes in 3.9 h on a single computing node, which was 23 times faster and 80% less RAM consumption compared to existing methods, enabling the in-depth data mining among microbiomes on a high resolution. HMS takes microbiome functional profiles as input, produces their pairwise distance matrix and PCoA coordinates.

Availability and implementation

It is coded in C/C++ with parallel computing and released in two alternative forms: a standalone software (https://github.com/qdu-bioinfo/hierarchical-meta-storms) and an equivalent R package (https://github.com/qdu-bioinfo/hrms).

Supplementary information

Supplementary data are available at Bioinformatics Advances online.