A molecular index for biological age identified from the metabolome and senescence-associated secretome in humans.

Unlike chronological age, biological age is a strong indicator of health of an individual. However, the molecular fingerprint associated with biological age is ill-defined. To define a high-resolution signature of biological age, we analyzed metabolome, circulating senescence-associated secretome (SASP)/inflammation markers and the interaction between them, from a cohort of healthy and rapid agers. The balance between two fatty acid oxidation mechanisms, ß-oxidation and ω-oxidation, associated with the extent of functional aging. Furthermore, a panel of 25 metabolites, Healthy Aging Metabolic (HAM) index, predicted healthy agers regardless of gender and race. HAM index was also validated in an independent cohort. Causal inference with machine learning implied three metabolites, ß-cryptoxanthin, prolylhydroxyproline, and eicosenoylcarnitine as putative drivers of biological aging. Multiple SASP markers were also elevated in rapid agers. Together, our findings reveal that a network of metabolic pathways underlie biological aging, and the HAM index could serve as a predictor of phenotypic aging in humans.

Microbiome depiction through user-adapted bioinformatic pipelines and parameters.

Introduction. The role of the microbiome in health and disease continues to be increasingly recognized. However, there is significant variability in the bioinformatic protocols for analysing genomic data. This, in part, has impeded the potential incorporation of microbiomics into the clinical setting and has challenged interstudy reproducibility. In microbial compositional analysis, there is a growing recognition for the need to move away from a one-size-fits-all approach to data processing.Gap Statement. Few evidence-based recommendations exist for setting parameters of programs that infer microbiota community profiles despite these parameters significantly impacting the accuracy of taxonomic inference.Aim. To compare three commonly used programs (DADA2, QIIME2, and mothur) and optimize them into four user-adapted pipelines for processing paired-end amplicon reads. We aim to increase the accuracy of compositional inference and help standardize microbiomic protocol.Methods. Two key parameters were isolated across four pipelines: filtering sequence reads based on a whole-number error threshold (maxEE) and truncating read ends based on a quality score threshold (QTrim). Closeness of sample inference was then evaluated using a mock community of known composition.Results. We observed that raw genomic data lost were proportionate to how stringently parameters were set. Exactly how much data were lost varied by pipeline. Accuracy of sample inference correlated with increased sequence read retention. Falsely detected taxa and unaccounted for microbial constituents were unique to pipeline and parameter. Implementation of optimized parameter values led to better approximation of the known mock community.Conclusions. Microbial compositions generated based on the 16S rRNA marker gene should be interpreted with caution. To improve microbial community profiling, bioinformatic protocols must be user-adapted. Analysis should be performed with consideration for the select target amplicon, pipelines and parameters used, and taxa of interest.

Detection of SARS-CoV-2 in high-efficiency particulate air (HEPA) filters of low-cost air purifiers in community-based organizations.

This study aims to investigate the presence of SARS-CoV-2 in public spaces and assess the utility of inexpensive air purifiers equipped with high-efficiency particulate air (HEPA) filters for viral detection. Samples were collected from six community-based organizations in underserved minority neighborhoods in Northwest Miami, Florida, from February to May 2022. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was used to detect SARS-CoV-2 in air purifier filters and surface swabs. Among 32 filters tested, three yielded positive results, while no positive surface swabs were found. Notably, positive samples were obtained exclusively from child daycare centers. These findings highlight the potential for airborne transmission of SARS-CoV-2 in indoor air, particularly in child daycare centers. Moreover, the study demonstrates the effectiveness of readily available HEPA filters in detecting the virus. Improving indoor ventilation and implementing air filtration systems are crucial in reducing COVID-19 transmission where people gather. Air filtration systems incorporating HEPA filters offer a valuable approach to virus detection and reducing transmission risks. Future research should explore the applicability of this technology for early identification and mitigation of viral outbreaks.

Attention Deficit Hyperactivity Disorder (ADHD) and the gut microbiome: An ecological perspective.

Attention Deficit Hyperactivity Disorder (ADHD) is an increasingly prevalent neuropsychiatric disorder characterized by hyperactivity, inattention, and impulsivity. Symptoms emerge from underlying deficiencies in neurocircuitry, and recent research has suggested a role played by the gut microbiome. The gut microbiome is an ecosystem of interdependent taxa involved in an exponentially complex web of interactions, plus host gene and reaction pathways, some of which involve neurotransmitters with roles in ADHD neurocircuitry. Studies have analyzed the ADHD gut microbiome using macroscale metrics such as diversity and differential abundance, and have proposed several taxa as elevated or reduced in ADHD compared to Control. Few studies have delved into the complex underlying dynamics ultimately responsible for the emergence of such metrics, leaving a largely incomplete, sometimes contradictory, and ultimately inconclusive picture. We aim to help complete this picture by venturing beyond taxa abundances and into taxa relationships (i.e. cooperation and competition), using a publicly available gut microbiome dataset (targeted 16S, v3-4 region, qPCR) from an observational, case-control study of 30 Control (15 female, 15 male) and 28 ADHD (15 female, 13 male) undergraduate students. We first perform the same macroscale analyses prevalent in ADHD gut microbiome literature (diversity, differential abundance, and composition) to observe the degree of correspondence, or any new trends. We then estimate two-way ecological relationships by producing Control and ADHD Microbial Co-occurrence Networks (MCNs), using SparCC correlations (p ≤ 0.01). We perform community detection to find clusters of taxa estimated to mutually cooperate along with their centroids, and centrality calculations to estimate taxa most vital to overall gut ecology. We finally summarize our results, providing conjectures on how they can guide future experiments, some methods for improving our experiments, and general implications for the field.

Microbiome maps: Hilbert curve visualizations of metagenomic profiles.

Abundance profiles from metagenomic sequencing data synthesize information from billions of sequenced reads coming from thousands of microbial genomes. Analyzing and understanding these profiles can be a challenge since the data they represent are complex. Particularly challenging is their visualization, as existing techniques are inadequate when the taxa number is in the thousands. We present a technique, and accompanying software, for the visualization of metagenomic abundance profiles using a space-filling curve that transforms a profile into an interactive 2D image. We created Jasper, an easy to use tool for the visualization and exploration of metagenomic profiles from DNA sequencing data. It orders taxa using a space-filling Hilbert curve, and creates a "Microbiome Map", where each position in the image represents the abundance of a single taxon from a reference collection. Jasper can order taxa in multiple ways, and the resulting microbiome maps can highlight "hot spots" of microbes that are dominant in taxonomic clades or biological conditions. We use Jasper to visualize samples from a variety of microbiome studies, and discuss ways in which microbiome maps can be an invaluable tool to visualize spatial, temporal, disease, and differential profiles. Our approach can create detailed microbiome maps involving hundreds of thousands of microbial reference genomes with the potential to unravel latent relationships (taxonomic, spatio-temporal, functional, and other) that could remain hidden using traditional visualization techniques. The maps can also be converted into animated movies that bring to life the dynamicity of microbiomes.

Lower respiratory tract microbiome composition and community interactions in smokers.

The lung microbiome impacts on lung function, making any smoking-induced changes in the lung microbiome potentially significant. The complex co-occurrence and co-avoidance patterns between the bacterial taxa in the lower respiratory tract (LRT) microbiome were explored for a cohort of active (AS), former (FS) and never (NS) smokers. Bronchoalveolar lavages (BALs) were collected from 55 volunteer subjects (9 NS, 24 FS and 22 AS). The LRT microbiome composition was assessed using 16S rRNA amplicon sequencing. Identification of differentially abundant taxa and co-occurrence patterns, discriminant analysis and biomarker inferences were performed. The data show that smoking results in a loss in the diversity of the LRT microbiome, change in the co-occurrence patterns and a weakening of the tight community structure present in healthy microbiomes. The increased abundance of the genus Ralstonia in the lung microbiomes of both former and active smokers is significant. Partial least square discriminant and DESeq2 analyses suggested a compositional difference between the cohorts in the LRT microbiome. The groups were sufficiently distinct from each other to suggest that cessation of smoking may not be sufficient for the lung microbiota to return to a similar composition to that of NS. The linear discriminant analysis effect size (LEfSe) analyses identified several bacterial taxa as potential biomarkers of smoking status. Network-based clustering analysis highlighted different co-occurring and co-avoiding microbial taxa in the three groups. The analysis found a cluster of bacterial taxa that co-occur in smokers and non-smokers alike. The clusters exhibited tighter and more significant associations in NS compared to FS and AS. Higher degree of rivalry between clusters was observed in the AS. The groups were sufficiently distinct from each other to suggest that cessation of smoking may not be sufficient for the lung microbiota to return to a similar composition to that of NS.

Diet Quality and Liver Health in People Living with HIV in the MASH Cohort: A Multi-Omic Analysis of the Fecal Microbiome and Metabolome.

The gut-liver axis has been recognized as a potential pathway in which dietary factors may contribute to liver disease in people living with HIV (PLWH). The objective of this study was to explore associations between dietary quality, the fecal microbiome, the metabolome, and liver health in PLWH from the Miami Adult Studies on HIV (MASH) cohort. We performed a cross-sectional analysis of 50 PLWH from the MASH cohort and utilized the USDA Healthy Eating Index (HEI)-2015 to measure diet quality. A Fibrosis-4 Index (FIB-4) score < 1.45 was used as a strong indication that advanced liver fibrosis was not present. Stool samples and fasting blood plasma samples were collected. Bacterial composition was characterized using 16S rRNA sequencing. Metabolomics in plasma were determined using gas and liquid chromatography/mass spectrometry. Statistical analyses included biomarker identification using linear discriminant analysis effect size. Compared to participants with FIB-4 ≥ 1.45, participants with FIB-4 < 1.45 had higher intake of dairy (p = 0.006). Fibrosis-4 Index score was inversely correlated with seafood and plant protein HEI component score (r = -0.320, p = 0.022). The relative abundances of butyrate-producing taxa Ruminococcaceae, Roseburia, and Lachnospiraceae were higher in participants with FIB-4 < 1.45. Participants with FIB-4 < 1.45 also had higher levels of caffeine (p = 0.045) and related metabolites such as trigonelline (p = 0.008) and 1-methylurate (p = 0.023). Dietary components appear to be associated with the fecal microbiome and metabolome, and liver health in PLWH. Future studies should investigate whether targeting specific dietary components may reduce liver-related morbidity and mortality in PLWH.

Microbiome Dysbiosis Shows Strong Association of Gut-Derived Altered Metabolomic Profile in Gulf War Chronic Multisymptom Illness Symptom Persistence Following Western Diet Feeding and Development of Obesity.

The pathophysiology of Gulf War Illness (GWI) remains elusive even after three decades. The persistence of multiple complex symptoms along with metabolic disorders such as obesity worsens the health of present Gulf War (GW) Veterans often by the interactions of the host gut microbiome and inflammatory mediators. In this study, we hypothesized that the administration of a Western diet might alter the host metabolomic profile, which is likely associated with the altered bacterial species. Using a five-month symptom persistence GWI model in mice and whole-genome sequencing, we characterized the species-level dysbiosis and global metabolomics, along with heterogenous co-occurrence network analysis, to study the bacteriome-metabolomic association. Microbial analysis at the species level showed a significant alteration of beneficial bacterial species. The beta diversity of the global metabolomic profile showed distinct clustering due to the Western diet, along with the alteration of metabolites associated with lipid, amino acid, nucleotide, vitamin, and xenobiotic metabolism pathways. Network analysis showed novel associations of gut bacterial species with metabolites and biochemical pathways that could be used as biomarkers or therapeutic targets to ameliorate symptom persistence in GW Veterans.

Unfolding and De-confounding: Biologically meaningful causal inference from longitudinal multi-omic networks using METALICA.

A key challenge in the analysis of microbiome data is the integration of multi-omic datasets and the discovery of interactions between microbial taxa, their expressed genes, and the metabolites they consume and/or produce. In an effort to improve the state-of-the-art in inferring biologically meaningful multi-omic interactions, we sought to address some of the most fundamental issues in causal inference from longitudinal multi-omics microbiome data sets. We developed METALICA, a suite of tools and techniques that can infer interactions between microbiome entities. METALICA introduces novel unrolling and de-confounding techniques used to uncover multi-omic entities that are believed to act as confounders for some of the relationships that may be inferred using standard causal inferencing tools. The results lend support to predictions about biological models and processes by which microbial taxa interact with each other in a microbiome. The unrolling process helps to identify putative intermediaries (genes and/or metabolites) to explain the interactions between microbes; the de-confounding process identifies putative common causes that may lead to spurious relationships to be inferred. METALICA was applied to the networks inferred by existing causal discovery and network inference algorithms applied to a multi-omics data set resulting from a longitudinal study of IBD microbiomes. The most significant unrollings and de-confoundings were manually validated using the existing literature and databases.

Multiomic analysis reveals microbiome-related relationships between cocaine use and metabolites.

OBJECTIVE: Over 19 million individuals globally have a cocaine use disorder, a significant public health crisis. Cocaine has also been associated with a pro-inflammatory state and recently with imbalances in the intestinal microbiota as compared to nonuse. The objective of this pilot study was to characterize the gut microbiota and plasma metabolites in people with HIV (PWH) who use cocaine compared with those who do not. DESIGN: Cross-sectional study. METHODS: A pilot study in PWH was conducted on 25 cocaine users and 25 cocaine nonusers from the Miami Adult Studies on HIV cohort. Stool samples and blood plasma were collected. Bacterial composition was characterized using 16S rRNA sequencing. Metabolomics in plasma were determined using gas and liquid chromatography/mass spectrometry. RESULTS: The relative abundances of the Lachnopspira genus, Oscillospira genus, Bifidobacterium adolescentis species, and Euryarchaeota phylum were significantly higher in the cocaine- using PWH compared to cocaine-nonusing PWH. Cocaine-use was associated with higher levels of several metabolites: products of dopamine catabolism (3-methoxytyrosine and 3-methoxytyramine sulfate), phenylacetate, benzoate, butyrate, and butyrylglycine. CONCLUSIONS: Cocaine use was associated with higher abundances of taxa and metabolites known to be associated with pathogenic states that include gastrointestinal conditions. Understanding key intestinal bacterial functional pathways that are altered due to cocaine use in PWH will provide a better understanding of the relationships between the host intestinal microbiome and potentially provide novel treatments to improve health.

Prior exposure to microcystin alters host gut resistome and is associated with dysregulated immune homeostasis in translatable mouse models.

A strong association between exposure to the common harmful algal bloom toxin microcystin and the altered host gut microbiome has been shown. We tested the hypothesis that prior exposure to the cyanotoxin microcystin-LR may alter the host resistome. We show that the mice exposed to microcystin-LR had an altered microbiome signature that harbored antibiotic resistance genes. Host resistome genotypes such as mefA, msrD, mel, ant6, and tet40 increased in diversity and relative abundance following microcystin-LR exposure. Interestingly, the increased abundance of these genes was traced to resistance to common antibiotics such as tetracycline, macrolides, glycopeptide, and aminoglycosides, crucial for modern-day treatment of several diseases. Increased abundance of these genes was positively associated with increased expression of PD1, a T-cell homeostasis marker, and pleiotropic inflammatory cytokine IL-6 with a concomitant negative association with immunosurveillance markers IL-7 and TLR2. Microcystin-LR exposure also caused decreased TLR2, TLR4, and REG3G expressions, increased immunosenescence, and higher systemic levels of IL-6 in both wild-type and humanized mice. In conclusion, the results show a first-ever characterization of the host resistome following microcystin-LR exposure and its connection to host immune status and antimicrobial resistance that can be crucial to understand treatment options with antibiotics in microcystin-exposed subjects in clinical settings.

Potential Autoimmunity Resulting from Molecular Mimicry between SARS-CoV-2 Spike and Human Proteins.

Molecular mimicry between viral antigens and host proteins can produce cross-reacting antibodies leading to autoimmunity. The coronavirus SARS-CoV-2 causes COVID-19, a disease curiously resulting in varied symptoms and outcomes, ranging from asymptomatic to fatal. Autoimmunity due to cross-reacting antibodies resulting from molecular mimicry between viral antigens and host proteins may provide an explanation. Thus, we computationally investigated molecular mimicry between SARS-CoV-2 Spike and known epitopes. We discovered molecular mimicry hotspots in Spike and highlight two examples with tentative high autoimmune potential and implications for understanding COVID-19 complications. We show that a TQLPP motif in Spike and thrombopoietin shares similar antibody binding properties. Antibodies cross-reacting with thrombopoietin may induce thrombocytopenia, a condition observed in COVID-19 patients. Another motif, ELDKY, is shared in multiple human proteins, such as PRKG1 involved in platelet activation and calcium regulation, and tropomyosin, which is linked to cardiac disease. Antibodies cross-reacting with PRKG1 and tropomyosin may cause known COVID-19 complications such as blood-clotting disorders and cardiac disease, respectively. Our findings illuminate COVID-19 pathogenesis and highlight the importance of considering autoimmune potential when developing therapeutic interventions to reduce adverse reactions.

A novel approach for combining the metagenome, metaresistome, metareplicome and causal inference to determine the microbes and their antibiotic resistance gene repertoire that contribute to dysbiosis.

The use of whole metagenomic data to infer the relative abundance of all its microbes is well established. The same data can be used to determine the replication rate of all eubacterial taxa with circular chromosomes. Despite their availability, the replication rate profiles (metareplicome) have not been fully exploited in microbiome analyses. Another relatively new approach is the application of causal inferencing to analyse microbiome data that goes beyond correlational studies. A novel scalable pipeline called MeRRCI (Metagenome, metaResistome, and metaReplicome for Causal Inferencing) was developed. MeRRCI combines efficient computation of the metagenome (bacterial relative abundance), metaresistome (antimicrobial gene abundance) and metareplicome (replication rates), and integrates environmental variables (metadata) for causality analysis using Bayesian networks. MeRRCI was applied to an infant gut microbiome data set to investigate the microbial community's response to antibiotics. Our analysis suggests that the current treatment stratagem contributes to preterm infant gut dysbiosis, allowing a proliferation of pathobionts. The study highlights the specific antibacterial resistance genes that may contribute to exponential cell division in the presence of antibiotics for various pathogens, namely Klebsiella pneumoniae, Citrobacter freundii, Staphylococcus epidermidis, Veilonella parvula and Clostridium perfringens. These organisms often contribute to the harmful long-term sequelae seen in these young infants.

Multi-Run Concrete Autoencoder to Identify Prognostic lncRNAs for 12 Cancers.

BACKGROUND: Long non-coding RNA plays a vital role in changing the expression profiles of various target genes that lead to cancer development. Thus, identifying prognostic lncRNAs related to different cancers might help in developing cancer therapy. METHOD: To discover the critical lncRNAs that can identify the origin of different cancers, we propose the use of the state-of-the-art deep learning algorithm concrete autoencoder (CAE) in an unsupervised setting, which efficiently identifies a subset of the most informative features. However, CAE does not identify reproducible features in different runs due to its stochastic nature. We thus propose a multi-run CAE (mrCAE) to identify a stable set of features to address this issue. The assumption is that a feature appearing in multiple runs carries more meaningful information about the data under consideration. The genome-wide lncRNA expression profiles of 12 different types of cancers, with a total of 4768 samples available in The Cancer Genome Atlas (TCGA), were analyzed to discover the key lncRNAs. The lncRNAs identified by multiple runs of CAE were added to a final list of key lncRNAs that are capable of identifying 12 different cancers. RESULTS: Our results showed that mrCAE performs better in feature selection than single-run CAE, standard autoencoder (AE), and other state-of-the-art feature selection techniques. This study revealed a set of top-ranking 128 lncRNAs that could identify the origin of 12 different cancers with an accuracy of 95%. Survival analysis showed that 76 of 128 lncRNAs have the prognostic capability to differentiate high- and low-risk groups of patients with different cancers. CONCLUSION: The proposed mrCAE, which selects actual features, outperformed the AE even though it selects the latent or pseudo-features. By selecting actual features instead of pseudo-features, mrCAE can be valuable for precision medicine. The identified prognostic lncRNAs can be further studied to develop therapies for different cancers.

Mutation-induced changes in the receptor-binding interface of the SARS-CoV-2 Delta variant B.1.617.2 and implications for immune evasion.

Following the initial surges of the Alpha (B.1.1.7) and the Beta (B.1.351) variants, a more infectious Delta variant (B.1.617.2) is now surging, further deepening the health crises caused by the pandemic. The sharp rise in cases attributed to the Delta variant has made it especially disturbing and is a variant of concern. Fortunately, current vaccines offer protection against known variants of concern, including the Delta variant. However, the Delta variant has exhibited some ability to dodge the immune system as it is found that neutralizing antibodies from prior infections or vaccines are less receptive to binding with the Delta spike protein. Here, we investigated the structural changes caused by the mutations in the Delta variant's receptor-binding interface and explored the effects on binding with the ACE2 receptor as well as with neutralizing antibodies. We find that the receptor-binding ß-loop-ß motif adopts an altered but stable conformation causing separation in some of the antibody binding epitopes. Our study shows reduced binding of neutralizing antibodies and provides a possible mechanism for the immune evasion exhibited by the Delta variant.

Effect of metronidazole on vaginal microbiota associated with asymptomatic bacterial vaginosis.

Vaginal dysbiosis-induced by an overgrowth of anaerobic bacteria is referred to as bacterial vaginosis (BV). The dysbiosis is associated with an increased risk for acquisition of sexually transmitted infections. Women with symptomatic BV are treated with oral metronidazole (MET), but its effectiveness remains to be elucidated. This study used whole-genome sequencing (WGS) to determine the changes in the microbiota among women treated with MET. WGS was conducted on DNA obtained from 20 vaginal swabs collected at four time points over 12 months from five randomly selected African American (AA) women. The baseline visit included all women who were diagnosed with asymptomatic BV and were untreated. All subjects were tested subsequently once every 2 months and received a course of MET for each BV episode during the 12 months. The BV status was classified according to Nugent scores (NSs) of vaginal smears. The microbial and resistome profiles were analysed along with the sociodemographic metadata. Despite treatment, none of the five participants reverted to normal vaginal flora - two were consistently positive for BV, and the rest experienced episodic cases of BV. WGS analyses showed Gardnerella spp. as the most abundant organism. After treatment with MET, there was an observed decline of Lactobacillus and Prevotella species. One participant had a healthy vaginal microbiota based on NS at one follow-up time point. Resistance genes including tetM and lscA were detected. Though limited in subjects, this study shows specific microbiota changes with treatment, presence of many resistant genes in their microbiota, and recurrence and persistence of BV despite MET treatment. Thus, MET may not be an effective treatment option for asymptomatic BV, and whole metagenome sequence would better inform the choice of antibiotics.

Causal effects in microbiomes using interventional calculus.

Causal inference in biomedical research allows us to shift the paradigm from investigating associational relationships to causal ones. Inferring causal relationships can help in understanding the inner workings of biological processes. Association patterns can be coincidental and may lead to wrong conclusions about causality in complex systems. Microbiomes are highly complex, diverse, and dynamic environments. Microbes are key players in human health and disease. Hence knowledge of critical causal relationships among the entities in a microbiome, and the impact of internal and external factors on microbial abundance and their interactions are essential for understanding disease mechanisms and making appropriate treatment recommendations. In this paper, we employ causal inference techniques to understand causal relationships between various entities in a microbiome, and to use the resulting causal network to make useful computations. We introduce a novel pipeline for microbiome analysis, which includes adding an outcome or "disease" variable, and then computing the causal network, referred to as a "disease network", with the goal of identifying disease-relevant causal factors from the microbiome. Internventional techniques are then applied to the resulting network, allowing us to compute a measure called the causal effect of one or more microbial taxa on the outcome variable or the condition of interest. Finally, we propose a measure called causal influence that quantifies the total influence exerted by a microbial taxon on the rest of the microiome. Our pipeline is robust, sensitive, different from traditional approaches, and able to predict interventional effects without any controlled experiments. The pipeline can be used to identify potential eubiotic and dysbiotic microbial taxa in a microbiome. We validate our results using synthetic data sets and using results on real data sets that were previously published.

Dynamic Bayesian Networks for Integrating Multi-omics Time Series Microbiome Data.

A key challenge in the analysis of longitudinal microbiome data is the inference of temporal interactions between microbial taxa, their genes, the metabolites that they consume and produce, and host genes. To address these challenges, we developed a computational pipeline, a pipeline for the analysis of longitudinal multi-omics data (PALM), that first aligns multi-omics data and then uses dynamic Bayesian networks (DBNs) to reconstruct a unified model. Our approach overcomes differences in sampling and progression rates, utilizes a biologically inspired multi-omic framework, reduces the large number of entities and parameters in the DBNs, and validates the learned network. Applying PALM to data collected from inflammatory bowel disease patients, we show that it accurately identifies known and novel interactions. Targeted experimental validations further support a number of the predicted novel metabolite-taxon interactions.IMPORTANCE While a number of large consortia collect and profile several different types of microbiome and genomic time series data, very few methods exist for joint modeling of multi-omics data sets. We developed a new computational pipeline, PALM, which uses dynamic Bayesian networks (DBNs) and is designed to integrate multi-omics data from longitudinal microbiome studies. When used to integrate sequence, expression, and metabolomics data from microbiome samples along with host expression data, the resulting models identify interactions between taxa, their genes, and the metabolites that they produce and consume, as well as their impact on host expression. We tested the models both by using them to predict future changes in microbiome levels and by comparing the learned interactions to known interactions in the literature. Finally, we performed experimental validations for a few of the predicted interactions to demonstrate the ability of the method to identify novel relationships and their impact.

Inferring directional relationships in microbial communities using signed Bayesian networks.

BACKGROUND: Microbe-microbe and host-microbe interactions in a microbiome play a vital role in both health and disease. However, the structure of the microbial community and the colonization patterns are highly complex to infer even under controlled wet laboratory conditions. In this study, we investigate what information, if any, can be provided by a Bayesian Network (BN) about a microbial community. Unlike the previously proposed Co-occurrence Networks (CoNs), BNs are based on conditional dependencies and can help in revealing complex associations. RESULTS: In this paper, we propose a way of combining a BN and a CoN to construct a signed Bayesian Network (sBN). We report a surprising association between directed edges in signed BNs and known colonization orders. CONCLUSIONS: BNs are powerful tools for community analysis and extracting influences and colonization patterns, even though the analysis only uses an abundance matrix with no temporal information. We conclude that directed edges in sBNs when combined with negative correlations are consistent with and strongly suggestive of colonization order.