The dynamic microbiome

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 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 were 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.

IMPORTANCE

We have developed a suite of tools and techniques capable of inferring interactions between microbiome entities. METALICA introduces novel techniques called unrolling and de-confounding that are employed to uncover multi-omic entities considered to be confounders for some of the relationships that may be inferred using standard causal inferencing tools. To evaluate our method, we conducted tests on the inflammatory bowel disease (IBD) dataset from the iHMP longitudinal study, which we pre-processed in accordance with our previous work. From this dataset, we generated various subsets, encompassing different combinations of metagenomics, metabolomics, and metatranscriptomics datasets. Using these multi-omics datasets, we demonstrate how the unrolling process aids in the identification of putative intermediaries (genes and/or metabolites) to explain the interactions between microbes. Additionally, the de-confounding process identifies potential common causes that may give rise to spurious relationships to be inferred. The most significant unrollings and de-confoundings were manually validated using the existing literature and databases.

Social and environmental predictors of gut microbiome age in wild baboons

Understanding why some individuals age faster than others is essential to evolutionary biology and geroscience, but measuring variation in biological age is difficult. One solution may lie in measuring gut microbiome composition because microbiota change with many age-related factors (e.g., immunity and behavior). Here we create a microbiome-based age predictor using 13,563 gut microbial profiles from 479 wild baboons collected over 14 years. The resulting "microbiome clock" predicts host chronological age. Deviations from the clock's predictions are linked to demographic and socio-environmental factors that predict baboon health and survival: animals who appear old-for-age tend to be male, sampled in the dry season (for females), and high social status (both sexes). However, an individual's "microbiome age" does not predict the attainment of developmental milestones or lifespan. Hence, the microbiome clock accurately reflects age and some social and environmental conditions, but not the pace of development or mortality risk.

The microbiome of the lichen Lobaria pulmonaria varies according to climate on a subcontinental scale

The Lobaria pulmonaria holobiont comprises algal, fungal, cyanobacterial and bacterial components. We investigated L. pulmonaria's bacterial microbiome in the adaptation of this ecologically sensitive lichen species to diverse climatic conditions. Our central hypothesis posited that microbiome composition and functionality aligns with subcontinental-scale (a stretch of ~1100 km) climatic parameters related to temperature and precipitation. We also tested the impact of short-term weather dynamics, sampling season and algal/fungal genotypes on microbiome variation. Metaproteomics provided insights into compositional and functional changes within the microbiome. Climatic variables explained 41.64% of microbiome variation, surpassing the combined influence of local weather and sampling season at 31.63%. Notably, annual mean temperature and temperature seasonality emerged as significant climatic drivers. Microbiome composition correlated with algal, not fungal genotype, suggesting similar environmental recruitment for the algal partner and microbiome. Differential abundance analyses revealed distinct protein compositions in Sub-Atlantic Lowland and Alpine regions, indicating differential microbiome responses to contrasting environmental/climatic conditions. Proteins involved in oxidative and cellular stress were notably different. Our findings highlight microbiome plasticity in adapting to stable climates, with limited responsiveness to short-term fluctuations, offering new insights into climate adaptation in lichen symbiosis.

mbtransfer: Microbiome intervention analysis using transfer functions and mirror statistics

Time series studies of microbiome interventions provide valuable data about microbial ecosystem structure. Unfortunately, existing models of microbial community dynamics have limited temporal memory and expressivity, relying on Markov or linearity assumptions. To address this, we introduce a new class of models based on transfer functions. These models learn impulse responses, capturing the potentially delayed effects of environmental changes on the microbial community. This allows us to simulate trajectories under hypothetical interventions and select significantly perturbed taxa with False Discovery Rate guarantees. Through simulations, we show that our approach effectively reduces forecasting errors compared to strong baselines and accurately pinpoints taxa of interest. Our case studies highlight the interpretability of the resulting differential response trajectories. An R package, mbtransfer, and notebooks to replicate the simulation and case studies are provided.

Forensic Microbiology: When, Where and How

Forensic microbiology is a relatively new discipline, born in part thanks to the development of advanced methodologies for the detection, identification and characterization of microorganisms, and also in relation to the growing impact of infectious diseases of iatrogenic origin. Indeed, the increased application of medical practices, such as transplants, which require immunosuppressive treatments, and the growing demand for prosthetic installations, associated with an increasing threat of antimicrobial resistance, have led to a rise in the number of infections of iatrogenic origin, which entails important medico-legal issues. On the other hand, the possibility of detecting minimal amounts of microorganisms, even in the form of residual traces (e.g., their nucleic acids), and of obtaining gene and genomic sequences at contained costs, has made it possible to ask new questions of whether cases of death or illness might have a microbiological origin, with the possibility of also tracing the origin of the microorganisms involved and reconstructing the chain of contagion. In addition to the more obvious applications, such as those mentioned above related to the origin of iatrogenic infections, or to possible cases of infections not properly diagnosed and treated, a less obvious application of forensic microbiology concerns its use in cases of violence or violent death, where the characterization of the microorganisms can contribute to the reconstruction of the case. Finally, paleomicrobiology, e.g., the reconstruction and characterization of microorganisms in historical or even archaeological remnants, can be considered as a sister discipline of forensic microbiology. In this article, we will review these different aspects and applications of forensic microbiology.

Human gut microbiome, diet, and mental disorders

Diet is one of the most important external factor shaping the composition and metabolic activities of the gut microbiome. The gut microbiome plays a crucial role in host health, including immune system development, nutrients metabolism, and the synthesis of bioactive molecules. In addition, the gut microbiome has been described as critical for the development of several mental disorders. Nutritional psychiatry is an emerging field of research that may provide a link between diet, microbial function, and brain health. In this study, we have reviewed the influence of different diet types, such as Western, Mediterranean, vegetarian, and ketogenic, on the gut microbiota composition and function, and their implication in various neuropsychiatric and psychological disorders.

Dynamic coexistence driven by physiological transitions in microbial communities

Microbial ecosystems are commonly modeled by fixed interactions between species in steady exponential growth states. However, microbes often modify their environments so strongly that they are forced out of the exponential state into stressed or non-growing states. Such dynamics are typical of ecological succession in nature and serial-dilution cycles in the laboratory. Here, we introduce a phenomenological model, the Community State model, to gain insight into the dynamic coexistence of microbes due to changes in their physiological states. Our model bypasses specific interactions (e.g., nutrient starvation, stress, aggregation) that lead to different combinations of physiological states, referred to collectively as "community states", and modeled by specifying the growth preference of each species along a global ecological coordinate, taken here to be the total community biomass density. We identify three key features of such dynamical communities that contrast starkly with steady-state communities: increased tolerance of community diversity to fast growth rates of species dominating different community states, enhanced community stability through staggered dominance of different species in different community states, and increased requirement on growth dominance for the inclusion of late-growing species. These features, derived explicitly for simplified models, are proposed here to be principles aiding the understanding of complex dynamical communities. Our model shifts the focus of ecosystem dynamics from bottom-up studies based on idealized inter-species interaction to top-down studies based on accessible macroscopic observables such as growth rates and total biomass density, enabling quantitative examination of community-wide characteristics.

Methodological Considerations in Longitudinal Analyses of Microbiome Data: A Comprehensive Review

Biological processes underlying health and disease are inherently dynamic and are best understood when characterized in a time-informed manner. In this comprehensive review, we discuss challenges inherent in time-series microbiome data analyses and compare available approaches and methods to overcome them. Appropriate handling of longitudinal microbiome data can shed light on important roles, functions, patterns, and potential interactions between large numbers of microbial taxa or genes in the context of health, disease, or interventions. We present a comprehensive review and comparison of existing microbiome time-series analysis methods, for both preprocessing and downstream analyses, including differential analysis, clustering, network inference, and trait classification. We posit that the careful selection and appropriate utilization of computational tools for longitudinal microbiome analyses can help advance our understanding of the dynamic host-microbiome relationships that underlie health-maintaining homeostases, progressions to disease-promoting dysbioses, as well as phases of physiologic development like those encountered in childhood.

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.

[The Therapeutic Potential of Prebiotics and Probiotics in Child and Adolescent Psychiatric Disorders]

The Therapeutic Potential of Prebiotics and Probiotics in Child and Adolescent Psychiatric Disorders Abstract: This short review summarizes the literature available on therapeutic interventions with prebiotics and probiotics and their potential use in psychiatric disorders in childhood, adolescence, and adulthood. Most studies of children and adolescents are done on ADHD and autism spectrum disorders, whereas single reports exist largely on positive effects on cognitive symptoms and quality of life. Initial studies regarding anorexia nervosa point to a potential effect of weight gain and reduction of gastrointestinal symptoms. To date, the effects of prebiotics and probiotics in depression, bipolar disorder, anxiety disorders, and schizophrenia have been mainly investigated in adults. The best reported evidence exists for depression, whereas the effects on depressive symptomatology are small. Positive effects are seen on gastrointestinal symptoms in these disorders. Given these positive effects, the mixed literature reports may result from very heterogeneous study designs. Nevertheless, the high potential of prebiotics and probiotics may be seen for minors with mental health problems. Further studies that include child and adolescent psychiatric populations and reflect the complexity of the gut-brain axis are urgently needed.

The Vaginal Microbiota of Pregnant Women Varies with Gestational Age, Maternal Age, and Parity

The composition of the vaginal microbiota is heavily influenced by pregnancy and may factor into pregnancy complications, including spontaneous preterm birth. However, results among studies have been inconsistent due, in part, to variation in sample sizes and ethnicity. Thus, an association between the vaginal microbiota and preterm labor continues to be debated. Yet, before assessing associations between the composition of the vaginal microbiota and preterm labor, a robust and in-depth characterization of the vaginal microbiota throughout pregnancy in the specific study population under investigation is required. Here, we report a large longitudinal study (n = 474 women, 1,862 vaginal samples) of a predominantly African-American cohort-a population that experiences a relatively high rate of pregnancy complications-evaluating associations between individual identity, gestational age, and other maternal characteristics with the composition of the vaginal microbiota throughout gestation resulting in term delivery. The principal factors influencing the composition of the vaginal microbiota in pregnancy are individual identity and gestational age at sampling. Other factors are maternal age, parity, obesity, and self-reported Cannabis use. The general pattern across gestation is for the vaginal microbiota to remain or transition to a state of Lactobacillus dominance. This pattern can be modified by maternal parity and obesity. Regardless, network analyses reveal dynamic associations among specific bacterial taxa within the vaginal ecosystem, which shift throughout the course of pregnancy. This study provides a robust foundational understanding of the vaginal microbiota in pregnancy and sets the stage for further investigation of this microbiota in obstetrical disease. IMPORTANCE There is debate regarding links between the vaginal microbiota and pregnancy complications, especially spontaneous preterm birth. Inconsistencies in results among studies are likely due to differences in sample sizes and cohort ethnicity. Ethnicity is a complicating factor because, although all bacterial taxa commonly inhabiting the vagina are present among all ethnicities, the frequencies of these taxa vary among ethnicities. Therefore, an in-depth characterization of the vaginal microbiota throughout pregnancy in the specific study population under investigation is required prior to evaluating associations between the vaginal microbiota and obstetrical disease. This initial investigation is a large longitudinal study of the vaginal microbiota throughout gestation resulting in a term delivery in a predominantly African-American cohort, a population that experiences disproportionally negative maternal-fetal health outcomes. It establishes the magnitude of associations between maternal characteristics, such as age, parity, body mass index, and self-reported Cannabis use, on the vaginal microbiota in pregnancy.

Microeukaryotic predators shape the wastewater microbiome

The physicochemical parameters that shape the prokaryotic community composition in wastewater have been extensively studied. In contrast, it is poorly understood whether and how biotic interactions affect the prokaryotic community composition in wastewater. We used metatranscriptomics data from a bioreactor sampled weekly over 14 months to investigate the wastewater microbiome, including often neglected microeukaryotes. Our analysis revealed that while prokaryotes are unaffected by seasonal changes in water temperature, they are impacted by a seasonal, temperature-induced change in the microeukaryotic community. Our findings suggest that selective predation pressure exerted by microeukaryotes is a significant factor shaping the prokaryotic community in wastewater. This study underscores the importance of investigating the entire wastewater microbiome to develop a comprehensive understanding of wastewater treatment.

Caution regarding the specificities of pan-cancer microbial structure

Results published in an article by Poore et al. (Nature. 2020;579:567-574) suggested that machine learning models can almost perfectly distinguish between tumour types based on their microbial composition using machine learning models. Whilst we believe that there is the potential for microbial composition to be used in this manner, we have concerns with the paper that make us question the certainty of the conclusions drawn. We believe there are issues in the areas of the contribution of contamination, handling of batch effects, false positive classifications and limitations in the machine learning approaches used. This makes it difficult to identify whether the authors have identified true biological signal and how robust these models would be in use as clinical biomarkers. We commend Poore et al. on their approach to open data and reproducibility that has enabled this analysis. We hope that this discourse assists the future development of machine learning models and hypothesis generation in microbiome research.

The TaxUMAP atlas: Efficient display of large clinical microbiome data reveals ecological competition in protection against bacteremia

Longitudinal microbiome data provide valuable insight into disease states and clinical responses, but they are challenging to mine and view collectively. To address these limitations, we present TaxUMAP, a taxonomically informed visualization for displaying microbiome states in large clinical microbiome datasets. We used TaxUMAP to chart a microbiome atlas of 1,870 patients with cancer during therapy-induced perturbations. Bacterial density and diversity were positively associated, but the trend was reversed in liquid stool. Low-diversity states (dominations) remained stable after antibiotic treatment, and diverse communities had a broader range of antimicrobial resistance genes than dominations. When examining microbiome states associated with risk for bacteremia, TaxUMAP revealed that certain Klebsiella species were associated with lower risk for bacteremia localize in a region of the atlas that is depleted in high-risk enterobacteria. This indicated a competitive interaction that was validated experimentally. Thus, TaxUMAP can chart comprehensive longitudinal microbiome datasets, enabling insights into microbiome effects on human health.

Modeling interaction networks between host, diet, and bacteria predicts obesogenesis in a mouse model

Host-microbiome interactions are known to have substantial effects on human health, but the diversity of the human microbiome makes it difficult to definitively attribute specific microbiome features to a host phenotype. One approach to overcoming this challenge is to use animal models of host-microbiome interaction, but it must be determined that relevant aspects of host-microbiome interactions are reflected in the animal model. One such experimental validation is an experiment by Ridura et al. In that experiment, transplanting a microbiome from a human into a mouse also conferred the human donor's obesity phenotype. We have aggregated a collection of previously published host-microbiome mouse-model experiments and combined it with thousands of sequenced and annotated bacterial genomes and metametabolomic pathways. Three computational models were generated, each model reflecting an aspect of host-microbiome interactions: 1) Predict the change in microbiome community structure in response to host diet using a community interaction network, 2) Predict metagenomic data from microbiome community structure, and 3) Predict host obesogenesis from modeled microbiome metagenomic data. These computationally validated models were combined into an integrated model of host-microbiome-diet interactions and used to replicate the Ridura experiment in silico. The results of the computational models indicate that network-based models are significantly more predictive than similar but non-network-based models. Network-based models also provide additional insight into the molecular mechanisms of host-microbiome interaction by highlighting metabolites and metabolic pathways proposed to be associated with microbiome-based obesogenesis. While the models generated in this study are likely too specific to the animal models and experimental conditions used to train our models to be of general utility in a broader understanding of obesogenesis, the approach detailed here is expected to be a powerful tool of investigating multiple types of host-microbiome interactions.

MDITRE: Scalable and Interpretable Machine Learning for Predicting Host Status from Temporal Microbiome Dynamics

Longitudinal microbiome data sets are being generated with increasing regularity, and there is broad recognition that these studies are critical for unlocking the mechanisms through which the microbiome impacts human health and disease. However, there is a dearth of computational tools for analyzing microbiome time-series data. To address this gap, we developed an open-source software package, Microbiome Differentiable Interpretable Temporal Rule Engine (MDITRE), which implements a new highly efficient method leveraging deep-learning technologies to derive human-interpretable rules that predict host status from longitudinal microbiome data. Using semi-synthetic and a large compendium of publicly available 16S rRNA amplicon and metagenomics sequencing data sets, we demonstrate that in almost all cases, MDITRE performs on par with or better than popular uninterpretable machine learning methods, and orders-of-magnitude faster than the prior interpretable technique. MDITRE also provides a graphical user interface, which we show through case studies can be used to derive biologically meaningful interpretations linking patterns of microbiome changes over time with host phenotypes. IMPORTANCE The human microbiome, or collection of microbes living on and within us, changes over time. Linking these changes to the status of the human host is crucial to understanding how the microbiome influences a variety of human diseases. Due to the large scale and complexity of microbiome data, computational methods are essential. Existing computational methods for linking changes in the microbiome to the status of the human host are either unable to scale to large and complex microbiome data sets or cannot produce human-interpretable outputs. We present a new computational method and software package that overcomes the limitations of previous methods, allowing researchers to analyze larger and more complex data sets while producing easily interpretable outputs. Our method has the potential to enable new insights into how changes in the microbiome over time maintain health or lead to disease in humans and facilitate the development of diagnostic tests based on the microbiome.

Machine learning on the road to unlocking microbiota's potential for boosting immune checkpoint therapy

The intestinal microbiota is a complex and diverse ecological community that fulfills multiple functions and substantially impacts human health. Despite its plasticity, unfavorable conditions can cause perturbations leading to so-called dysbiosis, which have been connected to multiple diseases. Unfortunately, understanding the mechanisms underlying the crosstalk between those microorganisms and their host is proving to be difficult. Traditionally used bioinformatic tools have difficulties to fully exploit big data generated for this purpose by modern high throughput screens. Machine Learning (ML) may be a potential means of solving such problems, but it requires diligent application to allow for drawing valid conclusions. This is especially crucial as gaining insight into the mechanistic basis of microbial impact on human health is highly anticipated in numerous fields of study. This includes oncology, where growing amounts of studies implicate the gut ecosystems in both cancerogenesis and antineoplastic treatment outcomes. Based on these reports and first signs of clinical benefits related to microbiota modulation in human trials, hopes are rising for the development of microbiome-derived diagnostics and therapeutics. In this mini-review, we're inspecting analytical approaches used to uncover the role of gut microbiome in immune checkpoint therapy (ICT) with the use of shotgun metagenomic sequencing (SMS) data.

Small molecule modulation of microbiota: a systems pharmacology perspective

BACKGROUND

Microbes are associated with many human diseases and influence drug efficacy. Small-molecule drugs may revolutionize biomedicine by fine-tuning the microbiota on the basis of individual patient microbiome signatures. However, emerging endeavors in small-molecule microbiome drug discovery continue to follow a conventional "one-drug-one-target-one-disease" process. A systematic pharmacology approach that would suppress multiple interacting pathogenic species in the microbiome, could offer an attractive alternative solution.

RESULTS

We construct a disease-centric signed microbe-microbe interaction network using curated microbe metabolite information and their effects on host. We develop a Signed Random Walk with Restart algorithm for the accurate prediction of effect of microbes on human health and diseases. With a survey on the druggable and evolutionary space of microbe proteins, we find that 8-10% of them can be targeted by existing drugs or drug-like chemicals and that 25% of them have homologs to human proteins. We demonstrate that drugs for diabetes can be the lead compounds for development of microbiota-targeted therapeutics. We further show that the potential drug targets that specifically exist in pathogenic microbes are periplasmic and cellular outer membrane proteins.

CONCLUSION

The systematic studies of the polypharmacological landscape of the microbiome network may open a new avenue for the small-molecule drug discovery of the microbiome. We believe that the application of systematic method on the polypharmacological investigation could lead to the discovery of novel drug therapies.

Social complexity as a driving force of gut microbiota exchange among conspecific hosts in non-human primates

The emergent concept of the social microbiome implies a view of a highly connected biological world, in which microbial interchange across organisms may be influenced by social and ecological connections occurring at different levels of biological organization. We explore this idea reviewing evidence of whether increasing social complexity in primate societies is associated with both higher diversity and greater similarity in the composition of the gut microbiota. By proposing a series of predictions regarding such relationship, we evaluate the existence of a link between gut microbiota and primate social behavior. Overall, we find that enough empirical evidence already supports these predictions. Nonetheless, we conclude that studies with the necessary, sufficient, explicit, and available evidence are still scarce. Therefore, we reflect on the benefit of founding future analyses on the utility of social complexity as a theoretical framework.

Principal microbial groups: compositional alternative to phylogenetic grouping of microbiome data

Statistical and machine learning techniques based on relative abundances have been used to predict health conditions and to identify microbial biomarkers. However, high dimensionality, sparsity and the compositional nature of microbiome data represent statistical challenges. On the other hand, the taxon grouping allows summarizing microbiome abundance with a coarser resolution in a lower dimension, but it presents new challenges when correlating taxa with a disease. In this work, we present a novel approach that groups Operational Taxonomical Units (OTUs) based only on relative abundances as an alternative to taxon grouping. The proposed procedure acknowledges the compositional data making use of principal balances. The identified groups are called Principal Microbial Groups (PMGs). The procedure reduces the need for user-defined aggregation of $\textrm{OTU}$s and offers the possibility of working with coarse group of $\textrm{OTU}$s, which are not present in a phylogenetic tree. PMGs can be used for two different goals: (1) as a dimensionality reduction method for compositional data, (2) as an aggregation procedure that provides an alternative to taxon grouping for construction of microbial balances afterward used for disease prediction. We illustrate the procedure with a cirrhosis study data. PMGs provide a coherent data analysis for the search of biomarkers in human microbiota. The source code and demo data for PMGs are available at: https://github.com/asliboyraz/PMGs.

Compilation of longitudinal gut microbiome, serum metabolome, and clinical data in acute myeloid leukemia

Induction chemotherapy for patients with acute myeloid leukemia (AML) is a unique clinical scenario. These patients spend several weeks in the hospital, receiving multiple antibiotics, experiencing gastrointestinal mucosal damage, and suffering severe impairments in their immune system and nutrition. These factors cause major disruptions to the gut microbiota to a level rarely seen in other clinical conditions. Thus, the study of the gut microbiota in these patients can reveal novel aspects of microbiota-host relationships. When combined with the circulating metabolome, such studies could shed light on gut microbiota contribution to circulating metabolites. Collectively, gut microbiota and circulating metabolome are known to regulate host physiology. We have previously deposited amplicon sequences from 566 fecal samples from 68 AML patients. Here, we provide sample-level details and a link, using de-identified patient IDs, to additional data including serum metabolomics (260 samples from 36 patients) and clinical metadata. The detailed information provided enables comprehensive multi-omics analysis. We validate the technical quality of these data through 3 examples and demonstrate a method for integrated analysis.

Microbiota succession throughout life from the cradle to the grave

Associations between age and the human microbiota are robust and reproducible. The microbial composition at several body sites can predict human chronological age relatively accurately. Although it is largely unknown why specific microorganisms are more abundant at certain ages, human microbiota research has elucidated a series of microbial community transformations that occur between birth and death. In this Review, we explore microbial succession in the healthy human microbiota from the cradle to the grave. We discuss the stages from primary succession at birth, to disruptions by disease or antibiotic use, to microbial expansion at death. We address how these successions differ by body site and by domain (bacteria, fungi or viruses). We also review experimental tools that microbiota researchers use to conduct this work. Finally, we discuss future directions for studying the microbiota's relationship with age, including designing consistent, well-powered, longitudinal studies, performing robust statistical analyses and improving characterization of non-bacterial microorganisms.

Statistical challenges in longitudinal microbiome data analysis

The microbiome is a complex and dynamic community of microorganisms that co-exist interdependently within an ecosystem, and interact with its host or environment. Longitudinal studies can capture temporal variation within the microbiome to gain mechanistic insights into microbial systems; however, current statistical methods are limited due to the complex and inherent features of the data. We have identified three analytical objectives in longitudinal microbial studies: (1) differential abundance over time and between sample groups, demographic factors or clinical variables of interest; (2) clustering of microorganisms evolving concomitantly across time and (3) network modelling to identify temporal relationships between microorganisms. This review explores the strengths and limitations of current methods to fulfill these objectives, compares different methods in simulation and case studies for objectives (1) and (2), and highlights opportunities for further methodological developments. R tutorials are provided to reproduce the analyses conducted in this review.

Network-based approaches for the investigation of microbial community structure and function using metagenomics-based data

Network-based approaches offer a powerful framework to evaluate microbial community organization and function as it relates to a variety of environmental processes. Emerging studies are exploring network theory as a method for data integration that is likely to be critical for the integration of 'omics' data using systems biology approaches. Intricacies of network theory and methodological and computational complexities in network construction, however, impede the use of these tools for translational science. We provide a perspective on the methods of network construction, interpretation and emerging uses for these techniques in understanding host-microbiota interactions.

Calorie restriction improves metabolic state independently of gut microbiome composition: a randomized dietary intervention trial

Background

The gut microbiota has been suggested to play a significant role in the development of overweight and obesity. However, the effects of calorie restriction on gut microbiota of overweight and obese adults, especially over longer durations, are largely unexplored.

Methods

Here, we longitudinally analyzed the effects of intermittent calorie restriction (ICR) operationalized as the 5:2 diet versus continuous calorie restriction (CCR) on fecal microbiota of 147 overweight or obese adults in a 50-week parallel-arm randomized controlled trial, the HELENA Trial. The primary outcome of the trial was the differential effects of ICR versus CCR on gene expression in subcutaneous adipose tissue. Changes in the gut microbiome, which are the focus of this publication, were defined as exploratory endpoint of the trial. The trial comprised a 12-week intervention period, a 12-week maintenance period, and a final follow-up period of 26 weeks.

Results

Both diets resulted in ~5% weight loss. However, except for Lactobacillales being enriched after ICR, post-intervention microbiome composition did not significantly differ between groups. Overall weight loss was associated with significant metabolic improvements, but not with changes in the gut microbiome. Nonetheless, the abundance of the Dorea genus at baseline was moderately predictive of subsequent weight loss (AUROC of 0.74 for distinguishing the highest versus lowest weight loss quartiles). Despite the lack of consistent intervention effects on microbiome composition, significant study group-independent co-variation between gut bacterial families and metabolic biomarkers, anthropometric measures, and dietary composition was detectable. Our analysis in particular revealed associations between insulin sensitivity (HOMA-IR) and Akkermansiaceae, Christensenellaceae, and Tanerellaceae. It also suggests the possibility of a beneficial modulation of the latter two intestinal taxa by a diet high in vegetables and fiber, and low in processed meat.

Conclusions

Overall, our results suggest that the gut microbiome remains stable and highly individual-specific under dietary calorie restriction.

Trial registration

The trial, including the present microbiome component, was prospectively registered at ClinicalTrials.govNCT02449148 on May 20, 2015.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01030-0.

ARZIMM: A Novel Analytic Platform for the Inference of Microbial Interactions and Community Stability from Longitudinal Microbiome Study

Dynamic changes of microbiome communities may play important roles in human health and diseases. The recent rise in longitudinal microbiome studies calls for statistical methods that can model the temporal dynamic patterns and simultaneously quantify the microbial interactions and community stability. Here, we propose a novel autoregressive zero-inflated mixed-effects model (ARZIMM) to capture the sparse microbial interactions and estimate the community stability. ARZIMM employs a zero-inflated Poisson autoregressive model to model the excessive zero abundances and the non-zero abundances separately, a random effect to investigate the underlining dynamic pattern shared within the group, and a Lasso-type penalty to capture and estimate the sparse microbial interactions. Based on the estimated microbial interaction matrix, we further derive the estimate of community stability, and identify the core dynamic patterns through network inference. Through extensive simulation studies and real data analyses we evaluate ARZIMM in comparison with the other methods.

Conceptual strategies for characterizing interactions in microbial communities

Understanding the sets of inter- and intraspecies interactions in microbial communities is a fundamental goal of microbial ecology. However, the study and quantification of microbial interactions pose several challenges owing to their complexity, dynamic nature, and the sheer number of unique interactions within a typical community. To overcome such challenges, microbial ecologists must rely on various approaches to distill the system of study to a functional and conceptualizable level, allowing for a practical understanding of microbial interactions in both simplified and complex systems. This review broadly addresses the role of several conceptual approaches available for the microbial ecologist’s arsenal, examines specific tools used to accomplish such approaches, and describes how the assumptions, expectations, and philosophies underlying these tools change across scales of complexity.

Functional response regression model on correlated longitudinal microbiome sequencing data

Functional regression has been widely used on longitudinal data, but it is not clear how to apply functional regression to microbiome sequencing data. We propose a novel functional response regression model analyzing correlated longitudinal microbiome sequencing data, which extends the classic functional response regression model only working for independent functional responses. We derive the theory of generalized least squares estimators for predictors' effects when functional responses are correlated, and develop a data transformation technique to solve the computational challenge for analyzing correlated functional response data using existing functional regression method. We show by extensive simulations that our proposed method provides unbiased estimations for predictors' effect, and our model has accurate type I error and power performance for correlated functional response data, compared with classic functional response regression model. Finally we implement our method to a real infant gut microbiome study to evaluate the relationship of clinical factors to predominant taxa along time.

Microbial trend analysis for common dynamic trend, group comparison, and classification in longitudinal microbiome study

BACKGROUND

The human microbiome is inherently dynamic and its dynamic nature plays a critical role in maintaining health and driving disease. With an increasing number of longitudinal microbiome studies, scientists are eager to learn the comprehensive characterization of microbial dynamics and their implications to the health and disease-related phenotypes. However, due to the challenging structure of longitudinal microbiome data, few analytic methods are available to characterize the microbial dynamics over time.

RESULTS

We propose a microbial trend analysis (MTA) framework for the high-dimensional and phylogenetically-based longitudinal microbiome data. In particular, MTA can perform three tasks: 1) capture the common microbial dynamic trends for a group of subjects at the community level and identify the dominant taxa; 2) examine whether or not the microbial overall dynamic trends are significantly different between groups; 3) classify an individual subject based on its longitudinal microbial profiling. Our extensive simulations demonstrate that the proposed MTA framework is robust and powerful in hypothesis testing, taxon identification, and subject classification. Our real data analyses further illustrate the utility of MTA through a longitudinal study in mice.

CONCLUSIONS

The proposed MTA framework is an attractive and effective tool in investigating dynamic microbial pattern from longitudinal microbiome studies.

Towards a deeper understanding of microbial communities: integrating experimental data with dynamic models

Microbial communities and their functions are shaped by complex networks of interactions among microbes and with their environment. While the critical roles microbial communities play in numerous environments have become increasingly appreciated, we have a very limited understanding of their interactions and how these interactions combine to generate community-level behaviors. This knowledge gap hinders our ability to predict community responses to perturbations and to design interventions that manipulate these communities to our benefit. Dynamic models are promising tools to address these questions. We review existing modeling techniques to construct dynamic models of microbial communities at different scales and suggest ways to leverage multiple types of models and data to facilitate our understanding and engineering of microbial communities.

phyLoSTM: a novel deep learning model on disease prediction from longitudinal microbiome data

MOTIVATION

Research shows that human microbiome is highly dynamic on longitudinal timescales, changing dynamically with diet, or due to medical interventions. In this paper, we propose a novel deep learning framework "phyLoSTM", using a combination of Convolutional Neural Networks and Long Short Term Memory Networks (LSTM) for feature extraction and analysis of temporal dependency in longitudinal microbiome sequencing data along with host's environmental factors for disease prediction. Additional novelty in terms of handling variable timepoints in subjects through LSTMs, as well as, weight balancing between imbalanced cases and controls is proposed.

RESULTS

We simulated 100 datasets across multiple time points for model testing. To demonstrate the model's effectiveness, we also implemented this novel method into two real longitudinal human microbiome studies: (i) DIABIMMUNE three country cohort with food allergy outcomes (Milk, Egg, Peanut and Overall) (ii) DiGiulio study with preterm delivery as outcome. Extensive analysis and comparison of our approach yields encouraging performance with an AUC of 0.897 (increased by 5%) on simulated studies and AUCs of 0.762 (increased by 19%) and 0.713 (increased by 8%) on the two real longitudinal microbiome studies respectively, as compared to the next best performing method, Random Forest. The proposed methodology improves predictive accuracy on longitudinal human microbiome studies containing spatially correlated data, and evaluates the change of microbiome composition contributing to outcome prediction.

AVAILABILITY AND IMPLEMENTATION

https://github.com/divya031090/phyLoSTM.

Structural identifiability of the generalized Lotka-Volterra model for microbiome studies

Population dynamic models can be used in conjunction with time series of species abundances to infer interactions. Understanding microbial interactions is a prerequisite for numerous goals in microbiome research, including predicting how populations change over time, determining how manipulations of microbiomes affect dynamics and designing synthetic microbiomes to perform tasks. As such, there is great interest in adapting population dynamic theory for microbial systems. Despite the appeal, numerous hurdles exist. One hurdle is that the data commonly obtained from DNA sequencing yield estimates of relative abundances, while population dynamic models such as the generalized Lotka-Volterra model track absolute abundances or densities. It is not clear whether relative abundance data alone can be used to infer parameters of population dynamic models such as the Lotka-Volterra model. We used structural identifiability analyses to determine the extent to which a time series of relative abundances can be used to parametrize the generalized Lotka-Volterra model. We found that only with absolute abundance data to accompany relative abundance estimates from sequencing can all parameters be uniquely identified. However, relative abundance data alone do contain information on relative interaction strengths, which is sufficient for many studies where the goal is to estimate key interactions and their effects on dynamics. Using synthetic data of a simple community for which we know the underlying structure, local practical identifiability analysis showed that modest amounts of both process and measurement error do not fundamentally affect these identifiability properties.

Challenges, Strategies, and Perspectives for Reference-Independent Longitudinal Multi-Omic Microbiome Studies

In recent years, multi-omic studies have enabled resolving community structure and interrogating community function of microbial communities. Simultaneous generation of metagenomic, metatranscriptomic, metaproteomic, and (meta) metabolomic data is more feasible than ever before, thus enabling in-depth assessment of community structure, function, and phenotype, thus resulting in a multitude of multi-omic microbiome datasets and the development of innovative methods to integrate and interrogate those multi-omic datasets. Specifically, the application of reference-independent approaches provides opportunities in identifying novel organisms and functions. At present, most of these large-scale multi-omic datasets stem from spatial sampling (e.g., water/soil microbiomes at several depths, microbiomes in/on different parts of the human anatomy) or case-control studies (e.g., cohorts of human microbiomes). We believe that longitudinal multi-omic microbiome datasets are the logical next step in microbiome studies due to their characteristic advantages in providing a better understanding of community dynamics, including: observation of trends, inference of causality, and ultimately, prediction of community behavior. Furthermore, the acquisition of complementary host-derived omics, environmental measurements, and suitable metadata will further enhance the aforementioned advantages of longitudinal data, which will serve as the basis to resolve drivers of community structure and function to understand the biotic and abiotic factors governing communities and specific populations. Carefully setup future experiments hold great potential to further unveil ecological mechanisms to evolution, microbe-microbe interactions, or microbe-host interactions. In this article, we discuss the challenges, emerging strategies, and best-practices applicable to longitudinal microbiome studies ranging from sampling, biomolecular extraction, systematic multi-omic measurements, reference-independent data integration, modeling, and validation.

Exploring the Role and Potential of Probiotics in the Field of Mental Health: Major Depressive Disorder

The field of probiotic has been exponentially expanding over the recent decades with a more therapeutic-centered research. Probiotics mediated microbiota modulation within the microbiota–gut–brain axis (MGBA) have been proven to be beneficial in various health domains through pre-clinical and clinical studies. In the context of mental health, although probiotic research is still in its infancy stage, the promising role and potential of probiotics in various mental disorders demonstrated via in-vivo and in-vitro studies have laid a strong foundation for translating preclinical models to humans. The exploration of the therapeutic role and potential of probiotics in major depressive disorder (MDD) is an extremely noteworthy field of research. The possible etio-pathological mechanisms of depression involving inflammation, neurotransmitters, the hypothalamic–pituitary–adrenal (HPA) axis and epigenetic mechanisms potentially benefit from probiotic intervention. Probiotics, both as an adjunct to antidepressants or a stand-alone intervention, have a beneficial role and potential in mitigating anti-depressive effects, and confers some advantages compared to conventional treatments of depression using anti-depressants.

Dynamic Changes in Fecal Microbial Communities of Neonatal Dairy Calves by Aging and Diarrhea

Microbiota plays a critical role in the overall growth performance and health status of dairy cows, especially during their early life. Several studies have reported that fecal microbiome of neonatal calves is shifted by various factors such as diarrhea, antibiotic treatment, or environmental changes. Despite the importance of gut microbiome, a lack of knowledge regarding the composition and functions of microbiota impedes the development of new strategies for improving growth performance and disease resistance during the neonatal calf period. In this study, we utilized next-generation sequencing to monitor the time-dependent dynamics of the gut microbiota of dairy calves before weaning (1-8 weeks of age) and further investigated the microbiome changes caused by diarrhea. Metagenomic analysis revealed that continuous changes, including increasing gut microbiome diversity, occurred from 1 to 5 weeks of age. However, the composition and diversity of the fecal microbiome did not change after 6 weeks of age. The most prominent changes in the fecal microbiome composition caused by aging at family level were a decreased abundance of Bacteroidaceae and Enterobacteriaceae and an increased abundance of Prevotellaceae. Phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) analysis indicated that the abundance of microbial genes associated with various metabolic pathways changed with aging. All calves with diarrhea symptoms showed drastic microbiome changes and about a week later returned to the microbiome of pre-diarrheal stage regardless of age. At phylum level, abundance of Bacteroidetes was decreased (p = 0.09) and that of Proteobacteria increased (p = 0.07) during diarrhea. PICRUSt analysis indicated that microbial metabolism-related genes, such as starch and sucrose metabolism, sphingolipid metabolism, alanine aspartate, and glutamate metabolism were significantly altered in diarrheal calves. Together, these results highlight the important implications of gut microbiota in gut metabolism and health status of neonatal dairy calves.

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.

Compilation of longitudinal microbiota data and hospitalome from hematopoietic cell transplantation patients

The impact of the gut microbiota in human health is affected by several factors including its composition, drug administrations, therapeutic interventions and underlying diseases. Unfortunately, many human microbiota datasets available publicly were collected to study the impact of single variables, and typically consist of outpatients in cross-sectional studies, have small sample numbers and/or lack metadata to account for confounders. These limitations can complicate reusing the data for questions outside their original focus. Here, we provide comprehensive longitudinal patient dataset that overcomes those limitations: a collection of fecal microbiota compositions (>10,000 microbiota samples from >1,000 patients) and a rich description of the "hospitalome" experienced by the hosts, i.e., their drug exposures and other metadata from patients with cancer, hospitalized to receive allogeneic hematopoietic cell transplantation (allo-HCT) at a large cancer center in the United States. We present five examples of how to apply these data to address clinical and scientific questions on host-associated microbial communities.

Impact of Temporal pH Fluctuations on the Coexistence of Nasal Bacteria in an in silico Community

To manipulate nasal microbiota for respiratory health, we need to better understand how this microbial community is assembled and maintained. Previous work has demonstrated that the pH in the nasal passage experiences temporal fluctuations. Yet, the impact of such pH fluctuations on nasal microbiota is not fully understood. Here, we examine how temporal fluctuations in pH might affect the coexistence of nasal bacteria in in silico communities. We take advantage of the cultivability of nasal bacteria to experimentally assess their responses to pH and the presence of other species. Based on experimentally observed responses, we formulate a mathematical model to numerically investigate the impact of temporal pH fluctuations on species coexistence. We assemble in silico nasal communities using up to 20 strains that resemble the isolates that we have experimentally characterized. We then subject these in silico communities to pH fluctuations and assess how the community composition and coexistence is impacted. Using this model, we then simulate pH fluctuations-varying in amplitude or frequency-to identify conditions that best support species coexistence. We find that the composition of nasal communities is generally robust against pH fluctuations within the expected range of amplitudes and frequencies. Our results also show that cooperative communities and communities with lower niche overlap have significantly lower composition deviations when exposed to temporal pH fluctuations. Overall, our data suggest that nasal microbiota could be robust against environmental fluctuations.

Longitudinal assessment of the bovine ocular bacterial community dynamics in calves

BACKGROUND

Infectious Bovine Keratoconjunctivitis (IBK), commonly known as pinkeye, is one of the most significant diseases of beef cattle. As such, IBK costs the US beef industry at least 150 million annually. However, strategies to prevent IBK are limited, with most cases resulting in treatment with antibiotics once the disease has developed. Longitudinal studies evaluating establishment of the ocular microbiota may identify critical risk periods for IBK outbreaks or changes in the microbiota that may predispose animals to IBK.

RESULTS

In an attempt to characterize the establishment and colonization patterns of the bovine ocular microbiota, we conducted a longitudinal study consisting of 227 calves and evaluated the microbiota composition over time using amplicon sequence variants (ASVs) based on 16S rRNA sequencing data and culture-based approaches. Beef calves on trial consisted of both male (intact) and females. Breeds were composed of purebred Angus and composites with varying percentages of Simmental, Angus, and Red Angus breeds. Average age at the start of the trial was 65 days ±15.02 and all calves remained nursing on their dam until weaning (day 139 of the study). The trial consisted of 139 days with four sampling time points on day 0, 21, 41, and 139. The experimental population received three different vaccination treatments (autogenous, commercial (both inactivated bacteria), and adjuvant placebo), to assess the effectiveness of different vaccines for IBK prevention. A significant change in bacterial community composition was observed across time periods sampled compared to the baseline (p < 0.001). However, no treatment effect of vaccine was detected within the ocular bacterial community. The bacterial community composition with the greatest time span between sampling time periods (98d span) was most similar to the baseline sample collected, suggesting re-establishment of the ocular microbiota to baseline levels over time after perturbation. The effect of IgA levels on the microbial community was investigated in a subset of cattle within the study. However, no significant effect of IgA was observed. Significant changes in the ocular microbiota were identified when comparing communities pre- and post-clinical signs of IBK. Additionally, dynamic changes in opportunistic pathogens Moraxella spp. were observed and confirmed using culture based methods.

CONCLUSIONS

Our results indicate that the bovine ocular microbiota is well represented by opportunistic pathogens such as Moraxella and Mycoplasma. Furthermore, this study characterizes the diversity of the ocular microbiota in calves and demonstrates the plasticity of the ocular microbiota to change. Additionally, we demonstrate the ocular microbiome in calves is similar between the eyes and the perturbation of one eye results in similar changes in the other eye. We also demonstrate the bovine ocular microbiota is slow to recover post perturbation and as a result provide opportunistic pathogens a chance to establish within the eye leading to IBK and other diseases. Characterizing the dynamic nature of the ocular microbiota provides novel opportunities to develop potential probiotic intervention to reduce IBK outbreaks in cattle.

The gut microbiota is associated with immune cell dynamics in humans

The gut microbiota influences development1-3 and homeostasis4-7 of the mammalian immune system, and is associated with human inflammatory8 and immune diseases9,10 as well as responses to immunotherapy11-14. Nevertheless, our understanding of how gut bacteria modulate the immune system remains limited, particularly in humans, where the difficulty of direct experimentation makes inference challenging. Here we study hundreds of hospitalized-and closely monitored-patients with cancer receiving haematopoietic cell transplantation as they recover from chemotherapy and stem-cell engraftment. This aggressive treatment causes large shifts in both circulatory immune cell and microbiota populations, enabling the relationships between the two to be studied simultaneously. Analysis of observed daily changes in circulating neutrophil, lymphocyte and monocyte counts and more than 10,000 longitudinal microbiota samples revealed consistent associations between gut bacteria and immune cell dynamics. High-resolution clinical metadata and Bayesian inference allowed us to compare the effects of bacterial genera in relation to those of immunomodulatory medications, revealing a considerable influence of the gut microbiota-together and over time-on systemic immune cell dynamics. Our analysis establishes and quantifies the link between the gut microbiota and the human immune system, with implications for microbiota-driven modulation of immunity.

Emerging computational tools and models for studying gut microbiota composition and function

The gut microbiota and its metabolites play critical roles in human health and disease. Advances in high-throughput sequencing, mass spectrometry, and other omics assay platforms have improved our ability to generate large volumes of data exploring the temporal variations in the compositions and functions of microbial communities. To elucidate mechanisms, methods and tools are needed that can rigorously model the dependencies within time-series data. Longitudinal data are often sparse and unevenly sampled, and nontrivial challenges remain in determining statistical significance, normalization across different data types, and model validation. In this review, we highlight recent developments in models and software tools for the analysis of time series microbiome and metabolome data, as well as integration of these data.

Multidrug-resistant bacteria as intestinal colonizers and evolution of intestinal colonization in healthy university students in Portugal

Multidrug-resistant bacteria have been increasingly described in healthcare institutions, however community resistance also seems to be emerging. Escherichia coli an intestinal commensal bacteria, is also a pathogen and represents an important intestinal reservoir of resistance. Our aim was the study of the intestinal colonization and of the persistence of antibiotic resistant intestinal bacteria in healthy university students of Porto, in the north of Portugal. Samples from 30 university students were collected and analysed. Two E. coli isolates were randomly obtained from each student and Gram-negative bacilli resistant to antibiotics were studied. In addition, we evaluated changes in the Gram-negative intestinal colonization of ten university students in a short period of time. Molecular characterization showed a high presence of bla _TEM in commensal E. coli . Gram-negative bacteria with intrinsic and extrinsic resistance were isolated, namely Pseudomonas spp., Enterobacter spp. and Pantoea spp. We isolated three ESBL-producing E. coli from two students. These isolates showed bla _CTX-M group 1 (n=1), bla _CTX-M group 9 (n=2), bla _TEM (n=2), bla _SHV (n=1) and tetA (n=2) genes. Additionally, they showed specific virulence factors and conjugational transfer of antibiotic resistance and virulence genes. One Pseudomonas spp. isolate resistant to carbapenems was detected colonizing one student. Our results confirm that healthy young adults may be colonized with commensals showing clinically relevant antibiotic resistance mechanisms, creating a risk of silent spread of these bacteria in the community.

Zero-Inflated gaussian mixed models for analyzing longitudinal microbiome data

Motivation

The human microbiome is variable and dynamic in nature. Longitudinal studies could explain the mechanisms in maintaining the microbiome in health or causing dysbiosis in disease. However, it remains challenging to properly analyze the longitudinal microbiome data from either 16S rRNA or metagenome shotgun sequencing studies, output as proportions or counts. Most microbiome data are sparse, requiring statistical models to handle zero-inflation. Moreover, longitudinal design induces correlation among the samples and thus further complicates the analysis and interpretation of the microbiome data.

Results

In this article, we propose zero-inflated Gaussian mixed models (ZIGMMs) to analyze longitudinal microbiome data. ZIGMMs is a robust and flexible method which can be applicable for longitudinal microbiome proportion data or count data generated with either 16S rRNA or shotgun sequencing technologies. It can include various types of fixed effects and random effects and account for various within-subject correlation structures, and can effectively handle zero-inflation. We developed an efficient Expectation-Maximization (EM) algorithm to fit the ZIGMMs by taking advantage of the standard procedure for fitting linear mixed models. We demonstrate the computational efficiency of our EM algorithm by comparing with two other zero-inflated methods. We show that ZIGMMs outperform the previously used linear mixed models (LMMs), negative binomial mixed models (NBMMs) and zero-inflated Beta regression mixed model (ZIBR) in detecting associated effects in longitudinal microbiome data through extensive simulations. We also apply our method to two public longitudinal microbiome datasets and compare with LMMs and NBMMs in detecting dynamic effects of associated taxa.

Evolutionary ecology of Lyme Borrelia

The bacterial genus, Borrelia, is comprised of vector-borne spirochete species that infect and are transmitted from multiple host species. Some Borrelia species cause highly-prevalent diseases in humans and domestic animals. Evolutionary, ecological, and molecular research on many Borrelia species have resulted in tremendous progress toward understanding the biology and natural history of these species. Yet, many outstanding questions, such as how Borrelia populations will be impacted by climate and land-use change, will require an interdisciplinary approach. The evolutionary ecology research framework incorporates theory and data from evolutionary, ecological, and molecular studies while overcoming common assumptions within each field that can hinder integration across these disciplines. Evolutionary ecology offers a framework to evaluate the ecological consequences of evolved traits and to predict how present-day ecological processes may result in further evolutionary change. Studies of microbes with complex transmission cycles, like Borrelia, which interact with multiple vertebrate hosts and arthropod vectors, are poised to leverage the power of the evolutionary ecology framework to identify the molecular interactions involved in ecological processes that result in evolutionary change. Using existing data, we outline how evolutionary ecology theory can delineate how interactions with other species and the physical environment create selective forces or impact migration of Borrelia populations and result in micro-evolutionary changes. We further discuss the ecological and molecular consequences of those micro-evolutionary changes. While many of the currently outstanding questions will necessitate new experimental designs and additional empirical data, many others can be addressed immediately by integrating existing molecular and ecological data within an evolutionary ecology framework.

Skin microbiome of beluga whales: spatial, temporal, and health-related dynamics

BACKGROUND

Host-specific microbiomes play an important role in individual health and ecology; in marine mammals, epidermal microbiomes may be a protective barrier between the host and its aqueous environment. Understanding these epidermal-associated microbial communities, and their ecological- or health-driven variability, is the first step toward developing health indices for rapid assessment of individual or population health. In Cook Inlet, Alaska, an endangered population of beluga whales (Delphinapterus leucas) numbers fewer than 300 animals and continues to decline, despite more than a decade of conservation effort. Characterizing the epidermal microbiome of this species could provide insight into the ecology and health of this endangered population and allow the development of minimally invasive health indicators based on tissue samples.

RESULTS

We sequenced the hypervariable IV region of bacterial and archaeal SSU rRNA genes from epidermal tissue samples collected from endangered Cook Inlet beluga whales (n = 33) and the nearest neighboring population in Bristol Bay (n = 39) between 2012 and 2018. We examined the sequences using amplicon sequence variant (ASV)-based analyses, and no ASVs were associated with all individuals, indicating a greater degree of epidermal microbiome variability among beluga whales than in previously studied cetacean species and suggesting the absence of a species-specific core microbiome. Epidermal microbiome composition differed significantly between populations and across sampling years. Comparing the microbiomes of Bristol Bay individuals of known health status revealed 11 ASVs associated with potential pathogens that differed in abundance between healthy individuals and those with skin lesions or dermatitis. Molting and non-molting individuals also differed significantly in microbial diversity and the abundance of potential pathogen-associated ASVs, indicating the importance of molting in maintaining skin health.

CONCLUSIONS

We provide novel insights into the dynamics of Alaskan beluga whale epidermal microbial communities. A core epidermal microbiome was not identified across all animals. We characterize microbial dynamics related to population, sampling year and health state including level of skin molting. The results of this study provide a basis for future work to understand the role of the skin microbiome in beluga whale health and to develop health indices for management of the endangered Cook Inlet beluga whales, and cetaceans more broadly.

Fast zero-inflated negative binomial mixed modeling approach for analyzing longitudinal metagenomics data

MOTIVATION

Longitudinal metagenomics data, including both 16S rRNA and whole-metagenome shotgun sequencing data, enhanced our abilities to understand the dynamic associations between the human microbiome and various diseases. However, analytic tools have not been fully developed to simultaneously address the main challenges of longitudinal metagenomics data, i.e. high-dimensionality, dependence among samples and zero-inflation of observed counts.

RESULTS

We propose a fast zero-inflated negative binomial mixed modeling (FZINBMM) approach to analyze high-dimensional longitudinal metagenomic count data. The FZINBMM approach is based on zero-inflated negative binomial mixed models (ZINBMMs) for modeling longitudinal metagenomic count data and a fast EM-IWLS algorithm for fitting ZINBMMs. FZINBMM takes advantage of a commonly used procedure for fitting linear mixed models, which allows us to include various types of fixed and random effects and within-subject correlation structures and quickly analyze many taxa. We found that FZINBMM remarkably outperformed in computational efficiency and was statistically comparable with two R packages, GLMMadaptive and glmmTMB, that use numerical integration to fit ZINBMMs. Extensive simulations and real data applications showed that FZINBMM outperformed other previous methods, including linear mixed models, negative binomial mixed models and zero-inflated Gaussian mixed models.

AVAILABILITY AND IMPLEMENTATION

FZINBMM has been implemented in the R package NBZIMM, available in the public GitHub repository http://github.com//nyiuab//NBZIMM.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.