Functional sequencing read annotation for high precision microbiome analysis

RFW captures species-level metagenomic functions by integrating genome annotation information

Functional profiling of whole-metagenome shotgun sequencing (WMS) enables our understanding of microbe-host interactions. We demonstrate microbial functional information loss by current annotation methods at both the taxon and community levels, particularly at lower read depths. To address information loss, we develop a framework, RFW (reference-based functional profile inference on WMS), that utilizes information from genome functional annotations and taxonomic profiles to infer microbial function abundances from WMS. Furthermore, we provide an algorithm for absolute abundance change quantification between groups as part of the RFW framework. By applying RFW to several datasets related to autism spectrum disorder and colorectal cancer, we show that RFW augments downstream analyses, such as differential microbial function identification and association analysis between microbial function and host phenotype. RFW is open source and freely available at https://github.com/Xingyinliu-Lab/RFW.

CAMDA 2023: Finding patterns in urban microbiomes

The Critical Assessment of Massive Data Analysis (CAMDA) addresses the complexities of harnessing Big Data in life sciences by hosting annual competitions that inspire research groups to develop innovative solutions. In 2023, the Forensic Challenge focused on identifying the city of origin for 365 metagenomic samples collected from public transportation systems and identifying associations between bacterial distribution and other covariates. For microbiome classification, we incorporated both taxonomic and functional annotations as features. To identify the most informative Operational Taxonomic Units, we selected features by fitting negative binomial models. We then implemented supervised models conducting 5-fold cross-validation (CV) with a 4:1 training-to-validation ratio. After variable selection, which reduced the dataset to fewer than 300 OTUs, the Support Vector Classifier achieved the highest F1 score (0.96). When using functional features from MIFASER, the Neural Network model outperformed other models. When considering climatic and demographic variables of the cities, Dirichlet regression over Escherichia, Enterobacter, and Klebsiella bacteria abundances suggests that population increase is indeed associated with a rise in the mean of Escherichia while decreasing temperature is linked to higher proportions of Klebsiella. This study validates microbiome classification using taxonomic features and, to a lesser extent, functional features. It shows that demographic and climatic factors influence urban microbial distribution. A Docker container and a Conda environment are available at the repository: GitHub facilitating broader adoption and validation of these methods by the scientific community.

Alterations of oral microbiome and metabolic signatures and their interaction in oral lichen planus

Background

Oral lichen planus (OLP) is a chronic oral mucosal inflammatory disease with a risk of becoming malignant. Emerging evidence suggests that microbial imbalance plays an important role in the development of OLP. However, the association between the oral microbiota and the metabolic features in OLP is still unclear.

Methods

We conducted 16S rRNA sequencing and metabolomics profiling on 95 OLP patients and 105 healthy controls (HC).To study oral microbes and metabolic changes in OLP, we applied differential analysis, Spearman correlation analysis and four machine learning algoeithms.

Results

The alpha and beta diversity both differed between OLP and HC. After adjustment for gender and age, we found an increase in the relative abundance of Pseudomonas, Aggregatibacter, Campylobacter, and Lautropia in OLP, while 18 genera decreased in OLP. A total of 153 saliva metabolites distinguishing OLP from HC were identified. Notably, correlations were found between Oribacterium, specific lipid and amino acid metabolites, and OLP's clinical phenotype. Additionally, the combination of Pseudomonas, Rhodococcus and (±)10-HDoHE effectively distinguished OLP from HC.

Conclusions

Based on multi-omics data, this study provides comprehensive evidence of a novel interplay between oral microbiome and metabolome in OLP pathogenesis using the oral microbiota and metabolites of OLP patients.

Assembling bacterial puzzles: piecing together functions into microbial pathways

Functional metagenomics enables the study of unexplored bacterial diversity, gene families, and pathways essential to microbial communities. However, discovering biological insights with these data is impeded by the scarcity of quality annotations. Here, we use a co-occurrence-based analysis of predicted microbial protein functions to uncover pathways in genomic and metagenomic biological systems. Our approach, based on phylogenetic profiles, improves the identification of functional relationships, or participation in the same biochemical pathway, between enzymes over a comparable homology-based approach. We optimized the design of our profiles to identify potential pathways using minimal data, clustered functionally related enzyme pairs into multi-enzymatic pathways, and evaluated our predictions against reference pathways in the KEGG database. We then demonstrated a novel extension of this approach to predict inter-bacterial protein interactions amongst members of a marine microbiome. Most significantly, we show our method predicts emergent biochemical pathways between known and unknown functions. Thus, our work establishes a basis for identifying the potential functional capacities of the entire metagenome, capturing previously unknown and abstract functions into discrete putative pathways.

Learning from the unknown: exploring the range of bacterial functionality

Determining the repertoire of a microbe's molecular functions is a central question in microbial biology. Modern techniques achieve this goal by comparing microbial genetic material against reference databases of functionally annotated genes/proteins or known taxonomic markers such as 16S rRNA. Here, we describe a novel approach to exploring bacterial functional repertoires without reference databases. Our Fusion scheme establishes functional relationships between bacteria and assigns organisms to Fusion-taxa that differ from otherwise defined taxonomic clades. Three key findings of our work stand out. First, bacterial functional comparisons outperform marker genes in assigning taxonomic clades. Fusion profiles are also better for this task than other functional annotation schemes. Second, Fusion-taxa are robust to addition of novel organisms and are, arguably, able to capture the environment-driven bacterial diversity. Finally, our alignment-free nucleic acid-based Siamese Neural Network model, created using Fusion functions, enables finding shared functionality of very distant, possibly structurally different, microbial homologs. Our work can thus help annotate functional repertoires of bacterial organisms and further guide our understanding of microbial communities.

Deep learning of a bacterial and archaeal universal language of life enables transfer learning and illuminates microbial dark matter

The majority of microbial genomes have yet to be cultured, and most proteins identified in microbial genomes or environmental sequences cannot be functionally annotated. As a result, current computational approaches to describe microbial systems rely on incomplete reference databases that cannot adequately capture the functional diversity of the microbial tree of life, limiting our ability to model high-level features of biological sequences. Here we present LookingGlass, a deep learning model encoding contextually-aware, functionally and evolutionarily relevant representations of short DNA reads, that distinguishes reads of disparate function, homology, and environmental origin. We demonstrate the ability of LookingGlass to be fine-tuned via transfer learning to perform a range of diverse tasks: to identify novel oxidoreductases, to predict enzyme optimal temperature, and to recognize the reading frames of DNA sequence fragments. LookingGlass enables functionally relevant representations of otherwise unknown and unannotated sequences, shedding light on the microbial dark matter that dominates life on Earth.

Computational methods to analyse microbial systems rely on reference databases which do not capture their full functional diversity. Here the authors develop a deep learning model and apply it using transfer learning, creating biologically useful models for multiple different tasks.

DIAMOND+MEGAN: Fast and Easy Taxonomic and Functional Analysis of Short and Long Microbiome Sequences

One main approach to computational analysis of microbiome sequences is to first align against a reference database of annotated protein sequences (NCBI-nr) and then perform taxonomic and functional binning of the sequences based on the resulting alignments. For both short and long reads (or assembled contigs), alignment is performed using DIAMOND, whereas taxonomic and functional binning, followed by inter- active exploration and analysis, is performed using MEGAN. We provide two step-by-step descriptions of this approach: © 2021 The Authors. Basic Protocol 1: Taxonomic and functional analysis of short read microbiome sequences Support Protocol 1: Preprocessing Basic Protocol 2: taxonomic and functional analysis of assembled long read microbiome sequences Support Protocol 2: Taxonomic binning and CheckM.

Impact of vitamin A transport and storage on intestinal retinoid homeostasis and functions

Lecithin:retinol acyltransferase and retinol-binding protein enable vitamin A (VA) storage and transport, respectively, maintaining tissue homeostasis of retinoids (VA derivatives). The precarious VA status of the lecithin:retinol acyltransferase-deficient (Lrat-/-) retinol-binding protein-deficient (Rbp-/-) mice rapidly deteriorates upon dietary VA restriction, leading to signs of severe vitamin A deficiency (VAD). As retinoids impact gut morphology and functions, VAD is often linked to intestinal pathological conditions and microbial dysbiosis. Thus, we investigated the contribution of VA storage and transport to intestinal retinoid homeostasis and functionalities. We showed the occurrence of intestinal VAD in Lrat-/-Rbp-/- mice, demonstrating the critical role of both pathways in preserving gut retinoid homeostasis. Moreover, in the mutant colon, VAD resulted in a compromised intestinal barrier as manifested by reduced mucins and antimicrobial defense, leaky gut, increased inflammation and oxidative stress, and altered mucosal immunocytokine profiles. These perturbations were accompanied by fecal dysbiosis, revealing that the VA status (sufficient vs. deficient), rather than the amount of dietary VA per se, is likely a major initial discriminant of the intestinal microbiome. Our data also pointed to a specific fecal taxonomic profile and distinct microbial functionalities associated with VAD. Overall, our findings revealed the suitability of the Lrat-/-Rbp-/- mice as a model to study intestinal dysfunctions and dysbiosis promoted by changes in tissue retinoid homeostasis induced by the host VA status and/or intake.

Guild-based analysis for understanding gut microbiome in human health and diseases

To demonstrate the causative role of gut microbiome in human health and diseases, we first need to identify, via next-generation sequencing, potentially important functional members associated with specific health outcomes and disease phenotypes. However, due to the strain-level genetic complexity of the gut microbiota, microbiome datasets are highly dimensional and highly sparse in nature, making it challenging to identify putative causative agents of a particular disease phenotype. Members of an ecosystem seldomly live independently from each other. Instead, they develop local interactions and form inter-member organizations to influence the ecosystem's higher-level patterns and functions. In the ecological study of macro-organisms, members are defined as belonging to the same "guild" if they exploit the same class of resources in a similar way or work together as a coherent functional group. Translating the concept of "guild" to the study of gut microbiota, we redefine guild as a group of bacteria that show consistent co-abundant behavior and likely to work together to contribute to the same ecological function. In this opinion article, we discuss how to use guilds as the aggregation unit to reduce dimensionality and sparsity in microbiome-wide association studies for identifying candidate gut bacteria that may causatively contribute to human health and diseases.

Metagenomics workflow for hybrid assembly, differential coverage binning, metatranscriptomics and pathway analysis (MUFFIN)

Metagenomics has redefined many areas of microbiology. However, metagenome-assembled genomes (MAGs) are often fragmented, primarily when sequencing was performed with short reads. Recent long-read sequencing technologies promise to improve genome reconstruction. However, the integration of two different sequencing modalities makes downstream analyses complex. We, therefore, developed MUFFIN, a complete metagenomic workflow that uses short and long reads to produce high-quality bins and their annotations. The workflow is written by using Nextflow, a workflow orchestration software, to achieve high reproducibility and fast and straightforward use. This workflow also produces the taxonomic classification and KEGG pathways of the bins and can be further used for quantification and annotation by providing RNA-Seq data (optionally). We tested the workflow using twenty biogas reactor samples and assessed the capacity of MUFFIN to process and output relevant files needed to analyze the microbial community and their function. MUFFIN produces functional pathway predictions and, if provided de novo metatranscript annotations across the metagenomic sample and for each bin. MUFFIN is available on github under GNUv3 licence: https://github.com/RVanDamme/MUFFIN.

A Cross-Sectional Study of Dairy Cattle Metagenomes Reveals Increased Antimicrobial Resistance in Animals Farmed in a Heavy Metal Contaminated Environment

The use of heavy metals in economic and social development can create an accumulation of toxic waste in the environment. High concentrations of heavy metals can damage human and animal health, lead to the development of antibiotic resistance, and possibly change in bovine microbiota. It is important to investigate the influence of heavy metals in food systems to determine potential harmful effects environmental heavy metal contamination on human health. Because of a mining dam rupture, 43 million cubic meters of iron ore waste flowed into the Doce river basin surrounding Mariana City, Brazil, in 2015. Following this environmental disaster, we investigated the consequences of long-term exposure to contaminated drinking water on the microbiome and resistome of dairy cattle. We identified bacterial antimicrobial resistance (AMR) genes in the feces, rumen fluid, and nasopharynx of 16 dairy cattle 4 years after the environmental disaster. Cattle had been continuously exposed to heavy metal contaminated water until sample collection (A) and compared them to analogous samples from 16 dairy cattle in an unaffected farm, 356 km away (B). The microbiome and resistome of farm A and farm B differed in many aspects. The distribution of genes present in the cattle's nasopharynx, rumen, and feces conferring AMR was highly heterogeneous, and most genes were present in only a few samples. The relative abundance and prevalence (presence/absence) of AMR genes were higher in farm A than in farm B. Samples from farm A had a higher prevalence (presence) of genes conferring resistance to multiple drugs, metals, biocides, and multi-compound resistance. Fecal samples had a higher relative abundance of AMR genes, followed by rumen fluid samples, and the nasopharynx had the lowest relative abundance of AMR genes detected. Metagenome functional annotation suggested that selective pressures of heavy metal exposure potentially skewed pathway diversity toward fewer, more specialized functions. This is the first study that evaluates the consequences of a Brazilian environmental accident with mining ore dam failure in the microbiome of dairy cows. Our findings suggest that the long-term persistence of heavy metals in the environment may result in differences in the microbiota and enrichment of antimicrobial-resistant bacteria. Our results also suggest that AMR genes are most readily detected in fecal samples compared to rumen and nasopharyngeal samples which had relatively lower bacterial read counts. Since heavy metal contamination has an effect on the animal microbiome, environmental management is warranted to protect the food system from hazardous consequences.

Snow microbiome functional analyses reveal novel aspects of microbial metabolism of complex organic compounds

Microbes active in extreme cold are not as well explored as those of other extreme environments. Studies have revealed a substantial microbial diversity and identified cold-specific microbiome molecular functions. We analyzed the metagenomes and metatranscriptomes of 20 snow samples collected in early and late spring in Svalbard, Norway using mi-faser, our read-based computational microbiome function annotation tool. Our results reveal a more diverse microbiome functional capacity and activity in the early- vs. late-spring samples. We also find that functional dissimilarity between the same-sample metagenomes and metatranscriptomes is significantly higher in early than late spring samples. These findings suggest that early spring samples may contain a larger fraction of DNA of dormant (or dead) organisms, while late spring samples reflect a new, metabolically active community. We further show that the abundance of sequencing reads mapping to the fatty acid synthesis-related microbial pathways in late spring metagenomes and metatranscriptomes is significantly correlated with the organic acid levels measured in these samples. Similarly, the organic acid levels correlate with the pathway read abundances of geraniol degradation and inversely correlate with those of styrene degradation, suggesting a possible nutrient change. Our study thus highlights the activity of microbial degradation pathways of complex organic compounds previously unreported at low temperatures.

Computational Approaches for Unraveling the Effects of Variation in the Human Genome and Microbiome

The past two decades of analytical efforts have highlighted how much more remains to be learned about the human genome and, particularly, its complex involvement in promoting disease development and progression. While numerous computational tools exist for the assessment of the functional and pathogenic effects of genome variants, their precision is far from satisfactory, particularly for clinical use. Accumulating evidence also suggests that the human microbiome's interaction with the human genome plays a critical role in determining health and disease states. While numerous microbial taxonomic groups and molecular functions of the human microbiome have been associated with disease, the reproducibility of these findings is lacking. The human microbiome–genome interaction in healthy individuals is even less well understood. This review summarizes the available computational methods built to analyze the effect of variation in the human genome and microbiome. We address the applicability and precision of these methods across their possible uses. We also briefly discuss the exciting, necessary, and now possible integration of the two types of data to improve the understanding of pathogenicity mechanisms.

Hierarchically Labeled Database Indexing Allows Scalable Characterization of Microbiomes

Increasingly available microbial reference data allow interpreting the composition and function of previously uncharacterized microbial communities in detail, via high-throughput sequencing analysis. However, efficient methods for read classification are required when the best database matches for short sequence reads are often shared among multiple reference sequences. Here, we take advantage of the fact that microbial sequences can be annotated relative to established tree structures, and we develop a highly scalable read classifier, PRROMenade, by enhancing the generalized Burrows-Wheeler transform with a labeling step to directly assign reads to the corresponding lowest taxonomic unit in an annotation tree. PRROMenade solves the multi-matching problem while allowing fast variable-size sequence classification for phylogenetic or functional annotation. Our simulations with 5% added differences from reference indicated only 1.5% error rate for PRROMenade functional classification. On metatranscriptomic data PRROMenade highlighted biologically relevant functional pathways related to diet-induced changes in the human gut microbiome.

Carnelian uncovers hidden functional patterns across diverse study populations from whole metagenome sequencing reads

Microbial populations exhibit functional changes in response to different ambient environments. Although whole metagenome sequencing promises enough raw data to study those changes, existing tools are limited in their ability to directly compare microbial metabolic function across samples and studies. We introduce Carnelian, an end-to-end pipeline for metabolic functional profiling uniquely suited to finding functional trends across diverse datasets. Carnelian is able to find shared metabolic pathways, concordant functional dysbioses, and distinguish Enzyme Commission (EC) terms missed by existing methodologies. We demonstrate Carnelian's effectiveness on type 2 diabetes, Crohn's disease, Parkinson's disease, and industrialized and non-industrialized gut microbiome cohorts.

Fingerprinting cities: differentiating subway microbiome functionality

BACKGROUND

Accumulating evidence suggests that the human microbiome impacts individual and public health. City subway systems are human-dense environments, where passengers often exchange microbes. The MetaSUB project participants collected samples from subway surfaces in different cities and performed metagenomic sequencing. Previous studies focused on taxonomic composition of these microbiomes and no explicit functional analysis had been done till now.

RESULTS

As a part of the 2018 CAMDA challenge, we functionally profiled the available ~ 400 subway metagenomes and built predictor for city origin. In cross-validation, our model reached 81% accuracy when only the top-ranked city assignment was considered and 95% accuracy if the second city was taken into account as well. Notably, this performance was only achievable if the similarity of distribution of cities in the training and testing sets was similar. To assure that our methods are applicable without such biased assumptions we balanced our training data to account for all represented cities equally well. After balancing, the performance of our method was slightly lower (76/94%, respectively, for one or two top ranked cities), but still consistently high. Here we attained an added benefit of independence of training set city representation. In testing, our unbalanced model thus reached (an over-estimated) performance of 90/97%, while our balanced model was at a more reliable 63/90% accuracy. While, by definition of our model, we were not able to predict the microbiome origins previously unseen, our balanced model correctly judged them to be NOT-from-training-cities over 80% of the time. Our function-based outlook on microbiomes also allowed us to note similarities between both regionally close and far-away cities. Curiously, we identified the depletion in mycobacterial functions as a signature of cities in New Zealand, while photosynthesis related functions fingerprinted New York, Porto and Tokyo.

CONCLUSIONS

We demonstrated the power of our high-speed function annotation method, mi-faser, by analysing ~ 400 shotgun metagenomes in 2 days, with the results recapitulating functional signals of different city subway microbiomes. We also showed the importance of balanced data in avoiding over-estimated performance. Our results revealed similarities between both geographically close (Ofa and Ilorin) and distant (Boston and Porto, Lisbon and New York) city subway microbiomes. The photosynthesis related functional signatures of NYC were previously unseen in taxonomy studies, highlighting the strength of functional analysis.

Supply of Methionine During Late-Pregnancy Alters Fecal Microbiota and Metabolome in Neonatal Dairy Calves Without Changes in Daily Feed Intake

To our knowledge, most studies demonstrating the role of manipulating maternal nutrition on hindgut (i.e., large intestine) microbiota in the offspring have been performed in non-ruminants. Whether this phenomenon exists in cattle is largely unknown. Therefore, the objectives of the current study were to evaluate the impact of maternal post-ruminal supply of methionine during late-pregnancy in dairy cows on fecal microbiota and metabolome in neonatal calves, and their association with body development and growth performance during the preweaning period. To achieve this, heifer calves, i.e., neonatal female offspring, born to Holstein cows receiving either a control (CON) diet (n = 13) or CON plus rumen-protected methionine (MET; Evonik Nutrition & Care GmbH) during the last 28 days of pregnancy were used. Fecal samples from heifers were collected from birth until 6 weeks of age, i.e., the preweaning period. Fecal microbiota was analyzed with QIIME 2 whereas fecal metabolites were measured using an untargeted LC-MS approach. At birth, MET heifers had greater (P ≤ 0.05) BW, HH, and WH. During the preweaning period, no differences between groups were detected for starter intake (P = 0.77). However, MET heifers maintained greater (P ≤ 0.05) BW, HH and tended (P = 0.06) to have greater WH and average daily gain (ADG) (P = 0.10). Fecal microbiota and metabolome profiles through 42 days of age in MET heifers indicated greater capacity for hindgut production of endogenous antibiotics and enhanced hindgut functionality and health. Enhancing maternal post-ruminal supply of methionine during late-gestation in dairy cows has a positive effect on hindgut functionality and health in their offspring through alterations in the fecal microbiota and metabolome without affecting feed intake. Those alterations could limit pathogen colonization of the hindgut while providing essential nutrients to the neonate. Together, such responses contribute to the ability of young calves to achieve better rates of nutrient utilization for growth.

Identification of Differentiating Metabolic Pathways between Infant Gut Microbiome Populations Reveals Depletion of Function-Level Adaptation to Human Milk in the Finnish Population

Knowing the limitations of taxonomy-based research, there is an emerging need for the development of higher-resolution techniques. The significance of this research is demonstrated by the novel method used for the analysis of function-level metagenomes. BiomeScout—the presented technology—utilizes proprietary algorithms for the detection of differences between functionalities present in metagenomic samples.

A variety of autoimmune and allergy events are becoming increasingly common, especially in Western countries. Some pieces of research link such conditions with the composition of microbiota during infancy. In this period, the predominant form of nutrition for gut microbiota is oligosaccharides from human milk (HMO). A number of gut-colonizing strains, such as Bifidobacterium and Bacteroides, are able to utilize HMO, but only some Bifidobacterium strains have evolved to digest the specific composition of human oligosaccharides. Differences in the proportions of the two genera that are able to utilize HMO have already been associated with the frequency of allergies and autoimmune diseases in the Finnish and the Russian populations. Our results show that differences in terms of the taxonomic annotation do not explain the reason for the differences in the Bifidobacterium/Bacteroides ratio between the Finnish and the Russian populations. In this paper, we present the results of function-level analysis. Unlike the typical workflow for gene abundance analysis, BiomeScout technology explains the differences in the Bifidobacterium/Bacteroides ratio. Our research shows the differences in the abundances of the two enzymes that are crucial for the utilization of short type 1 oligosaccharides.

IMPORTANCE Knowing the limitations of taxonomy-based research, there is an emerging need for the development of higher-resolution techniques. The significance of this research is demonstrated by the novel method used for the analysis of function-level metagenomes. BiomeScout—the presented technology—utilizes proprietary algorithms for the detection of differences between functionalities present in metagenomic samples.

clubber: removing the bioinformatics bottleneck in big data analyses

With the advent of modern day high-throughput technologies, the bottleneck in biological discovery has shifted from the cost of doing experiments to that of analyzing results. clubber is our automated cluster-load balancing system developed for optimizing these “big data” analyses. Its plug-and-play framework encourages re-use of existing solutions for bioinformatics problems. clubber’s goals are to reduce computation times and to facilitate use of cluster computing. The first goal is achieved by automating the balance of parallel submissions across available high performance computing (HPC) resources. Notably, the latter can be added on demand, including cloud-based resources, and/or featuring heterogeneous environments. The second goal of making HPCs user-friendly is facilitated by an interactive web interface and a RESTful API, allowing for job monitoring and result retrieval. We used clubber to speed up our pipeline for annotating molecular functionality of metagenomes. Here, we analyzed the Deepwater Horizon oil-spill study data to quantitatively show that the beach sands have not yet entirely recovered. Further, our analysis of the CAMI-challenge data revealed that microbiome taxonomic shifts do not necessarily correlate with functional shifts. These examples (21 metagenomes processed in 172 min) clearly illustrate the importance of clubber in the everyday computational biology environment.